# Diffuse interstellar bands as dust indicators: The contribution from 3D maps

R. Lallement, J. L. Vergely, N. L. J. Cox

► To cite this version:

R. Lallement, J. L. Vergely, N. L. J. Cox. Diffuse interstellar bands as dust indicators: The contribution from 3D maps. Astronomy & Astrophysics - A&A, 2024, 691, pp.A41. 10.1051/0004-6361/202451657. hal-04760132

HAL Id: hal-04760132

https://hal.science/hal-04760132v1

Submitted on 30 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Diffuse interstellar bands as dust indicators:
The contribution from 3D maps

R. Lallement<fnref n="1" />,*, J. L. Vergely<fnref n="2" />, and N. L. J. Cox<fnref n="3" />

Received 25 July 2024 / Accepted 3 September 2024

## ABSTRACT

**Context.** Three-dimensional (3D) distributions of the 862 nm diffuse interstellar band (DIB) carrier have been computed based on Gaia parallaxes and DIB catalogues, in parallel with 3D maps of dust extinction density. Three-dimensional maps provide local diagnostics and information on the distribution of structures in addition to line-of-sight (LOS) integrated quantities, and allow us to focus on poorly studied low-extinction areas. They make cross-matching with other catalogues possible through estimates of the DIB and extinction along any given path.

**Aims.** We re-examined the relationships between the density of DIB carriers and the absorption and emission properties of spatially co-located dust. Along with laboratory identifications of carriers, these properties may shed light on the formation and evolution of this organic matter. They may also help to model dust emission and absorption properties in a more detailed way.

**Methods.** We used the 3D maps of 862 nm DIBs and of dust extinction, as well as available DIB equivalent width (EW) catalogues and published measurements of parameters characterizing the dust extinction law and the dust emission. We studied the relationships between the extinction-normalized 862 nm DIB EW and the extinction level, the total-to-selective extinction ratio $R_V$, and the dust far-IR emission spectral index $\beta$. We re-visited the link between several DIBs and the UV absorption bump at 220 nm.

**Results.** The ratio of the 862 nm DIB carrier density to the optical extinction density ($DIB_{norm}^{862}$) is increasing in low-density clouds, confirming with local values the trend seen in the LOS data. In both cases, the coefficients of a fitted power law fall within the range of those measured towards SDSS high-latitude targets for 20 different bands, ranking this DIB among those with a high increase, above that of the broad 4430 Å DIB. This is consistent with the recent measurement of a larger scale height above the Plane for the 862 nm DIB compared to that of the 4430 Å DIB. Using map-integrated 862 nm DIB EWs and extinctions along the paths to APOGEE targets with published proxies $R'_V$ for the total-to-selective extinction ratio, we found that, despite a large scatter, $DIB_{norm}^{862}$ is positively correlated with $R'_V$ for those stars with low to moderate extinctions ($A_V = 0.2$ to 2–3 mag). Independently, using stars from the 862 nm DIB catalogue located outside the disk and for the same regime of extinction, $DIB_{norm}^{862}$ is found to be globally anti-correlated with the Planck opacity spectral index $\beta$. This is consistent with the observed anti-correlation between $\beta$ and $R'_V$. In the light of a recent result on the variability of the carbon/silicate ratio in dust grains as a source of this anti-correlation, it suggests that $DIB_{norm}^{862}$ increases with the fraction of carbonaceous to silicate grains in the co-located dust, in agreement with the carbonaceous nature of DIB carriers and recent evidences for a spatial correlation between $DIB_{norm}^{862}$ and the fluxes of carbon-rich ejecta of asymptotic giant branch (AGB) stars. At higher extinction both trends disappear, and there is evidence for a trend reversal. Regarding the link between the height of the 220 nm UV absorption bump and extinction-normalized EWs of DIBs, we found that two factors explain the absence of previous clear results: the correlation disappears when we move from $\sigma$-type to $\zeta$-type DIBs and/or from single-cloud LOSs to paths crossing multiple clouds distant from each other; $\zeta$-type bands can be used to predict low and high values of the bump height, provided one adds a correcting factor linked to the ambient radiation field (e.g. the 5780/5797 DIB ratio). We show examples of simple models of the bump height based on the 5780 Å band, the 5850 Å band and the 5780/5797 DIB ratio. We also found an anti-correlation between $DIB_{norm}$ and the width of the bump, which similarly disappears from $\sigma$-type to $\zeta$-type DIBs. This suggests that a fraction of the bump is generated outside the dense molecular clouds.

**Conclusions.** There are complex relationships between DIBs and dust; however, massive measurements of DIBs and extinction and the derived 3D maps may provide some constraints on the density, the nature, and the contribution to extinction and emission of the co-located dust. This requires large stellar spectroscopic surveys and space-based measurements of UV extinction.

**Key words.** ISM: abundances – ISM: clouds – dust, extinction – ISM: general – ISM: lines and bands – ISM: molecules

## 1. Introduction

Diffuse interstellar bands (DIBs) have already fascinated (and sometimes discouraged) two generations of specialists of the interstellar medium and astrochemists. These weak irregular absorptions are imprinted in the spectra of all astrophysical

sources located behind interstellar clouds. About 600 of these features of various intensities and widths have been observed in the visible and near-infrared (NIR), and discoveries of new bands are continuously made, especially for faint or very broad bands (Ebenbichler et al. 2024; Maíz Apellániz et al. 2021), and in the NIR (Hamano et al. 2016). Their width, shape, and strength strongly vary from one band to another (for catalogues and recent reviews, see Jenniskens & Desert 1994; Hobbs et al. 2008;

<footnote n="1">¹ GEPI, Observatoire de Paris, PSL University, CNRS, 5 Place Jules Janssen, 92190 Meudon, France</footnote>

<footnote n="2">² ACRI-ST, 260 Route du Pin Montard, BP234, 06904 Sophia-Antipolis, France</footnote>

<footnote n="3">³ ACRI-ST, Centre d'Etudes et de Recherche de Grasse (CERGA), 10 Av. Nicolas Copernic, 06130 Grasse, France</footnote>

<footnote>* Corresponding author; rosine.lallement@obspm.fr</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1. 862 nm DIB carrier and extinction volume densities. All voxels of the Cox et al. (2024) 3D maps have been used. Left: DIB density averaged in voxels with equal extinction density. An adjusted power law limited to the low-extinction regime is superimposed (see text). Right: ratio of the DIB EW density to the extinction density for the same voxels, and superimposed power law.</figcaption></figure>

a tool for estimates of integrated quantities along the LOS to any target, including those which were not part of the inverted catalogue used to produce the map, and have been observed for various purposes. Finally, a fourth advantage is the possibility of a documented selection of regions and directions appropriate to a specific study (e.g. selections of single-cloud LOSs or targets beyond the layer of dust). For this work we made use of all these advantages.

The 862 nm DIB is a broad (≈0.4 nm wide) DIB. It is not particularly strong; however, being spectrally close to ionized calcium (CaII) lines commonly used to constrain the stellar parameters, its study benefits from stellar spectroscopic surveys. It was intensively studied based on the spectra of the RAVE survey and was found to be well correlated with the reddening (Munari et al. 2008; Kos et al. 2013, 2014). This was also confirmed by data from the Gaia-ESO and GIBS surveys (Puspitarini et al. 2015; Zhao et al. 2021), and the 862 nm DIB was found useful as a tracer of the arms and of their dynamical properties. It was shown to be well correlated with the NIR DIB at 1.5273 µm (Zhao et al. 2021). It did not show any particular characteristics that set it apart from other bands, and, as a result, it is probably representative of most other DIBs.

Section 2 is devoted to the relationship between the 862 nm DIB and the extinction, and to the comparison with DIB-extinction relationships found for other DIBs. In Sect. 3, we detail the search for a link between the total-to-extinction ratio $R_V$ and the DIB strength normalized to the extinction ($DIB_{norm}^{862}$). In Sect. 4 we similarly detail the search for a link between $DIB_{norm}^{862}$ (for the 862 nm DIB) and the dust opacity spectral index $\beta$ measured by Planck. In Sect. 5, we revisit the studies between DIBs and the 220 nm bump of the extinction law. Section 6 discusses the various types of relationships, their limitations, and their potential use. We also discuss some implications.

## 2. Extinction-normalized DIB and extinction

Significant attention has been given to DIB variability in star-forming regions and, in particular, to DIB carrier depletion in dense UV-protected parts of dense clouds (the skin effect). It is due, in large part, to the observation preference for highly extinguished target stars. Conversely, much less attention has been

given to the low-extinction part of the DIB dependence on dust opacity. The 3D maps of the 862 nm DIB (Cox et al. 2024) are not of high enough spatial resolution to address the skin effect for individual clouds. However, for the first time, it is possible to use the maps to investigate the relationship between the local value of the DIB carrier density $\rho$DIB, measured as the spatial gradient of the DIB EW, and the local value of the extinction density $\rho A_V$ in wide low-dust areas, or similarly the relationship between the extinction-normalized DIB density $\rho DIB_{norm}^{862} = \rho$DIB/$\rho A_V$ and the extinction density $\rho A_V$.

To do so, we selected all $10 \times 10 \times 10$ pc$^3$ voxels of both DIB and extinction 3D maps containing enough target stars to perform the local inversion of DIB EWs (see details in Cox et al. 2024). We distributed the DIB densities $\rho$DIB in these voxels in bins of similar local extinction density $\rho A_V$, and computed their average in each bin. We note that densities $\rho A_V$ are lower than actual densities in moderate- or high-density clouds because the spatial resolution of the 3D maps used is poor, including that of the dust map that was especially computed based on the same series of LOS as those of the DIB catalogue, and there is a dilution of clumpy structures in the map voxels. We did the same sampling for the ratios $\rho$DIB / $\rho A_V$. Figure 1 displays the two distributions of bin-averaged DIB density $\rho$DIB on the one hand, and bin-averaged ratio $\rho$DIB / $\rho A_V$ (i.e. the averaged extinction-normalized DIB density $\rho$DIB$_{norm}^{862}$) as a function of the local extinction density $\rho A_V$. Superimposed are fitted power laws:

$$ \rho \text{DIB} = K (\rho A_V)^\gamma, \\ \rho \text{DIB}_{\text{norm}} = \frac{\rho \text{DIB}}{\rho A_V} = K (\rho A_V)^{\gamma-1} = K (\rho A_V)^{\alpha}. \quad (1) $$

We restricted the adjustment to the extinction density interval $0.00018 < \rho A_V < 0.0044 \text{ mag pc}^{-1}$. The lower limit ensures the exclusion of uncertain ratios corresponding to low statistics, while the upper limit excludes voxels containing dense clouds. The power-law coefficients $\gamma$ are found to be $\gamma = +0.74$ for the DIB versus extinction dependence and, accordingly, $\alpha = \gamma - 1 = -0.26$ for the normalized DIB dependence.

In addition to local values, we also used directly measured LOS-integrated EW values from two catalogues: the catalogue of individual LOSs (Gaia Collaboration 2023b) and the catalogue of averaged EWs among stars belonging to the same
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** 862 nm DIB EW and LOS-integrated extinction. Top: DIB EW from the *Gaia* DR3 catalogue as a function of the extinction. Average integrated from the Sun to the target star. A power law is fitted to the data from $A_V=0$ to $A_V=1.5$ mag. Isocontours of target number densities have been drawn to illustrate the various contributions. Bottom: same as top, but for the second catalogue of averaged EWs in spatial bins.</figcaption></figure>

spatial bin (Gaia Collaboration 2023a). For both datasets, we used the EXPLORE 3D extinction density maps (Vergely et al. 2022) to compute the corresponding integrated extinctions, from the Sun to each target star in the first case, and from the Sun to each centre of the spatial bin in the second case. Figure 2 displays the EW measurements for individual targets of the first catalogue (left), and for spatial bins from the second catalogue (right), both as a function of the corresponding integrated extinction. An adjusted power law is superimposed for the two datasets. We restricted the fit to the low-moderate extinction interval $0.12 \le A_V \le 1.5$ mag. The $\gamma$ coefficients were found to be 0.57 and 0.62, respectively, for the two catalogues. The two coefficients are similar, and are slightly lower than that found from the 3D maps (0.74). Taken as a whole, these results provide an order of magnitude of the uncertainty on the coefficients for such a simple relationship. One can see that they depend on the DIB extraction method and that they also differ depending on the use of direct or inverted measurements. However, the existence of similar trends for integrated values and for local values confirms that the increase, for lower values of the dust extinction density, in the extinction-normalized DIB strength is well established. In particular, the result based on local values, taken from 3D distributions, confirms that it is not due to an existing bias in the measurements of weak extinctions, or low-DIB EWs.

An interesting aspect of this exercise is the comparison with the results obtained by Lan et al. (2015) for a series of different DIBs. The authors used Sloane Digital Sky Survey (SDSS) spectra of halo stars, galaxies, and quasars at high Galactic latitudes, and extracted 20 different DIBs from stacked spectra in a large fraction of the high-latitude sky. They adjusted power-law functions to the DIB EWs as a function of the reddening, for the same extinction interval as the one we used above in the case of the *Gaia* catalogue (namely $0.12 \le A_V \le 1.5$ mag). Their Table 3 lists power-law coefficients $\gamma$ for the DIB dependence on $E(B-V)$. The whole set of $\gamma$ coefficients is in the range from 0.51 to 1.13 and six of the DIBs have a low $\gamma$ between 0.51 and 0.74. This demonstrates that the 862 nm DIB behaves as several other DIBs do, but is part of those bands with the highest relative increase with respect to extinction in dust-poor regions (i.e. the lowest $\gamma$). By comparison, for the strong 5780 Å DIB Lan et al. (2015) derived a $\gamma$ coefficient close to 1, which explains the linear correlation between its EW and the extinction in the regimes of low to moderate extinction. A comparison can also be made with the 9577 and 9632 Å $C_{60}^+$ DIBs. Schlarmann et al. (2021) inferred from published data a strong decrease (more than a factor of three) in the DIB/reddening ratio when the $E(B-V)$ reddening increases from 0.1 mag to 0.6 mag (quantities derived from their Fig. 3c). This would correspond roughly to a power-law coefficient $\gamma \approx 0.4$, which is lower than the lowest from the Lan et al. (2015) series (0.51). The stability of the $C_{60}$ structure may explain why it can persist in low-density, generally UV-permeated clouds.

Since the density of dust decreases with the distance to the Plane, the $\gamma$ coefficient is expected to have some link with the DIB carrier scale height: the lower the $\gamma$, the larger the scale height. A recent study by Zhao et al. (2023) derives a significant difference in scale heights for two DIBs, namely the 4430 Å DIB and the 862 nm DIB; the 862 nm DIB scale height is larger than that of the 4430 Å DIB. The Lan et al. (2015) $\gamma$ coefficient for the 4430 Å DIB is 0.89, above our three estimates of $\gamma$ for the 862 nm DIB, all three below 0.74. This is consistent with the heights measured by Zhao et al. (2023). Further measurements of DIB scale heights should help to confirm (or not) their direct link with the $\gamma$ coefficient. We note that the hierarchy of DIBs with respect to these Lan et al. (2015) $\gamma$ coefficients does not seem to correspond to the hierarchy of $\zeta$ to $\sigma$ DIBs, related to the UV radiation field (see e.g. the work of Ensor et al. 2017). This means that the phenomenon responsible for the increase in carrier density in low-dust clouds is not directly or simply related to the $\zeta$ versus $\sigma$ classification. This deserves further study, being, instead, a likely indication of the carrier formation process and site.

### 3. Extinction-normalized DIB and total-to-selective extinction ratio $R_V$

Schlafly et al. (2016, 2017) published a large catalogue of measurements of a proxy $R'_V$ for the total-to-selective extinction ratio $R_V = A_V/E(B-V)$, estimated for the LOSs to SDSS/APOGEE target stars. To do so, they combined accurate stellar atmospheric parameters derived from APOGEE spectra with available photometric data. They found wide-area variations in $R'_V$ on large angular scales ($\ge 30^\circ$). The lack of a link of these variations with the interstellar (IS) matter column density led the authors to conclude that they are not governed by grain size changes in dense clouds, and are due to other effects acting on large scales. They also found that $R'_V$ is globally higher outside the solar circle (at galactocentric distances larger than that of the Sun).
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** $R'_V$ values for APOGEE targets from Schlafly et al. (2016) (top) and extinction-normalized DIB for the same LOSs (bottom) in Galactic coordinates. The sizes of the markers were reduced for LOS extinctions beyond 4 mag. The limited resolution of the 3D maps is responsible for the slightly smoother image for the normalized DIB. There is a some similarity between the two independent measurements, with notable exceptions, e.g. in the Aquila region ($l = 0°, b = +5°, +15°$), and Taurus ($l = 160°, b = -15°$ or $b = -22°$) characterized by a strong opacity.</figcaption></figure>

Interestingly, a global increase in the extinction-normalized 862 nm DIB with the galactocentric radius is also one of the results of the 3D mapping presented in Cox et al. (2024), a radial increase potentially attributed to the increased prevalence of ejecta from C-rich AGBs over ejecta from O-rich stars. If the two phenomena were linked, it would imply the existence of a positive correlation between the extinction-normalized 862 nm DIB EW $DIB_{norm}^{862}$ and $R'_V$, which motivated the following study.

We used the catalogue of Schlafly et al. (2016), and for all LOSs to the APOGEE targets, we integrated along these paths within the 3D distribution of DIB density from Cox et al. (2024) and within the 3D distributions of extinction density from Vergely et al. (2022) to get estimates of the 862 nm DIB cumulative EW on the one hand and of the extinction Av of the target on the other hand. Their ratio $DIB_{norm}^{862}$ was then computed. Since the spatial extent of the Vergely et al. (2022) 3D maps is resolution-dependent, we used the series of extinction maps in a hierarchical mode. We started the integration in the highest resolution map and, if the map boundary is reached before the target star location, we continued to integrate in the second-best resolution map, and so on. Targets outside the 3D DIB density map, of smaller extent than the widest extinction map (+/- 4 kpc instead of +/- 5 kpc), were eliminated. As said in Sect. 1, this is one of the potential uses of 3D maps, namely the computation of DIB and extinction estimates for any target star that is not part of the inverted catalogue. The APOGEE catalogue is appropriate for our study for the following reasons: it is massive; it is homogeneously obtained; and most of the targets are distant, but the extinction is limited (i.e. very dense cloud cores are generally avoided). The last aspect is important since the resolution of the 3D maps is too poor to reconstruct very small structures. Instead, they are appropriate for the study of large-scale characteristics.

We made a first visual comparison of distributions over the sky of $R'_V$ measurements and computed extinction-normalized

DIB EWS $DIB_{norm}^{862}$. Figure 3 shows the locations of the APOGEE LOSs in Galactic coordinates, colour-coded according to $R'_V$ from Schlafly et al. (2016) (top) and to the normalized 862 nm DIB (bottom). There are some large-scale similarities, for example in the third quadrant where both quantities are particularly high. There are also discrepancies, especially in the highly opaque directions above the Galactic centre, or in the Taurus area where $R'_V$ is high, but the DIB-to-extinction ratio is notably low. As a result, these sky distributions suggest a link between the two quantities, but certainly different regimes must exist according to the extinction. We note that the $R'_V$ values from Schlafly et al. (2016) vary between 3.0 and 3.7, while $DIB_{norm}^{862}$ changes by a much larger factor, of the order of three.

A more direct comparison is presented in Fig. 4, where the extinction-normalized DIB for each star located within the 3D DIB map is displayed as a function of the inverse of the total-to-selective extinction ratio proxy $R'_V^{-1}$. Since all extinctions that have been inverted to produce the 3D extinction maps were estimated under the assumption of a unique fixed $R_V=3.1$ parameter, to avoid a bias and take into account the actual variability of $R_V$ we multiplied each individual extinction by a factor of $R_V/3.1$ before computing the DIB-to-extinction ratio. Despite the very large scatter, there is a clear decrease in $DIB_{norm}^{862}$ versus $R'_V^{-1}$ (i.e. an increase with $R'_V$). This decrease occurs in the low- to moderate-extinction regime $A_V \lesssim 3$ mag (extinctions are colour-coded), which is the range for most APOGEE targets used by Schlafly et al. (2016). Apart from the largest group of data points, one can distinguish several higher extinction data points at the lower left of the graph ($R'_V^{-1}$ between 0.2 and 0.26) that strongly depart from the general trend. This suggests that the trend disappears for some directions characterized by a high extinction. The decrease versus $R'_V^{-1}$ observed at low to moderate extinction, is opposite to what was found by Ramírez-Tannus et al. (2018) for $R_V$ based on targets of the M17 open cluster.
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Relationship between the extinction-normalized 862 nm DIB and the Schlafly et al. (2016) total to selective extinction $R'_V$ for APOGEE LOSs (dots). For low to moderate extinctions, despite the large scatter, $DIB_{norm}^{862}$ decreases with $R'_V^{-1}$. To illustrate the trend, iso-contours of data point numbers per 0.007 × 0.04 rectangular bins in $R'_V^{-1} \times DIB_{norm}^{862}$ are drawn for 40, 80, and 120 points. The trend is reversed at very high extinctions, as shown by the superimposed results of Ramírez-Tannus et al. (2018) for M17 (thick crosses). For these data points we neglected the difference between $R'_V$ and $R_V$. All markers are colour-coded according to the LOS extinction.</figcaption></figure>

However, the regime of extinction for the M17 stars considered by the authors is clearly above that of the APOGEE targets ($A_V$ varies between 5 and 9 mag for the M17 targets). We have superimposed their results in Fig. 4, using the same colour-coding representing the extinction as used for the APOGEE data. The M17 data points depart from the main distribution of the APOGEE data points to reach lower $R_V^{-1}$ and lower $DIB_{norm}^{862}$ values. From the combination of the two results, it appears that the DIB-$R_V$ positive relationship is suppressed when the extinction reaches levels of the order of ≈3 mag, with the appearance of the negative trend noted by Ramírez-Tannus et al. (2018). High $R_V$ values are generally associated with large grains formed by accretion or coagulation in UV-shielded molecular clouds, with high extinctions and DIB depletion. This is the trend observed by Ramírez-Tannus et al. (2018) above $A_V = 3$ mag. On the other hand, the opposite trend appearing at lower extinction is a different phenomenon.

**4. Extinction-normalized DIB and dust emission spectral index from Planck**

In the frame of their study, Schlafly et al. (2016) discovered a tight inverse relationship between $R'_V$ and the dust opacity spectral index $\beta$, the parameter used in the modified black-body representation of the Galactic dust far-IR emission. Properties or mechanisms at the origin of this trend have been searched for and, more recently, Zelko & Finkbeiner (2020) found that increasing the carbon-to-silicate ratio in dust grains may explain the simultaneous increase in $R_V$ and decrease in $\beta$. This effect comes in addition to the known increase in $R_V$ in opaque clouds,

attributed to the increased size of the grains. Such a result on the influence of the dust composition has important consequences in the context of extinction laws and also cosmological observations, through dust effects on the CMB spectral distortion (Zelko & Finkbeiner 2021). If the above positive correlation between $DIB_{norm}^{862}$ and $R'_V$ for low- to moderate-extinction APOGEE LOSs is due to the existence of a physical link between the DIB carrier density and the grain composition, it should, in turn, be accompanied by a negative correlation between $DIB_{norm}^{862}$ and the far-IR dust emission spectral index $\beta$. This expectation motivated the following study.

As already mentioned, contrary to $R_V$, $\beta$ is not a distance-limited quantity and target stars or locations in the 3D maps must be carefully selected to ensure that most of the IS matter in their direction is located in front of them, and there is negligible IS matter beyond. Here we used the first catalogue of individual 862 nm DIBs of Gaia Collaboration (2023b) due to the more extended distribution of directions and selected DIB data for which the quality flag is 0. For each direction we estimated the spectral index $\beta$ by interpolation in the high-resolution Planck map (Planck Collaboration Int. XLVIII 2016). We also used, in parallel, the more recent map of Casandjian et al. (2022). For each target, we computed the extinction $A_V$ by integrating in the highest resolution dust extinction density maps from Vergely et al. (2022), available from the EXPLORE facility. Similarly to the comparison with $R_V$, we first show, for a preliminary visual comparison, the distributions over the sky of the DIB-to-extinction ratio $DIB_{norm}^{862}$ and the two interpolated $\beta$ values (Fig. 5). We masked the Galactic latitudes lower than 6° since in those directions $\beta$ is too strongly influenced by dust beyond the limits of the extinction maps. We tuned the colour-coding of
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Images of 862 nm DIB<sup>862</sup><sub>norm</sub> (middle) and the two estimates of the dust opacity spectral index β (top: Planck Collaboration Int. XLVIII 2016 and bottom: Casandjian et al. 2022). There are strong similarities, again with notable exceptions in high-extinction regions, mainly the Aquila rift area (black ellipse).</figcaption></figure>

the β maps in order to obtain about the same contrast in β and DIB<sup>862</sup><sub>norm</sub>. There are obvious similarities of the patterns at large angular scale. We note that the optimal similarity is somewhat higher in the case of the Casandjian et al. (2022) estimates of β. In both cases, there are also obvious discrepancies, for example at ≈5–10° Galactic longitude and ≈20°–30° Galactic latitude in the Aquila Rift region.

Figure 6 illustrates the complex dependence on β of the extinction-normalized DIB, and especially the variability of this dependence according to the extinction regime. Here we used data from the second catalogue of Gaia Collaboration (2023a) for its reduced uncertainties on the EWs and selected flag=0 data. The extinction was computed by integration in the 3D maps from the Sun to the centre of the spatial bins. The figure displays DIB<sup>862</sup><sub>norm</sub> as a function of the LOS opacity. For each LOS the dust emission spectral index β interpolated from the Planck maps for the corresponding direction of each DIB bin centre is colour-coded. For Av values below ~3 mag, a decrease in the normalized DIB (going down along vertical lines in the figure) corresponds to an increase in β (colour variation from red to green). Above Av at ~3–4 mag, this trend is no longer visible, and there seems to be a weak reversed trend, namely β is decreasing (colours from violet to blue); however, the number of LOSs becomes low, something expected from the rare cases of high extinction and high Galactic latitude. This bi-modal behaviour, namely an anti-correlation between DIB<sup>862</sup><sub>norm</sub> and β for low to moderate extinctions, and a null-positive correlation with β in the high-opacity regime, is reminiscent of what we found for the relationship between the normalized DIB and R<sub>V</sub>. Since the R<sub>V</sub> and β data sources are completely independent, and since, as noted above, it has been shown that R<sub>V</sub> and beta are anti-correlated (Schlafly et al. 2016), this reinforces the two detected trends.

Given the complexity of the dependence on β, we performed comparisons in selected limited regimes of Galactic latitude and opacity. As an example of the various comparisons, Fig. 7 shows the dependence on β of the DIB-to-extinction ratio for Galactic latitudes between 10° and 20° and again for selected data (flag zero in the catalogue). The negative trend is clearly seen. Figure 8 shows another example of the dependence on β of the DIB-to-extinction ratio for a restricted range of opacity values, namely extinctions A<sub>V</sub> comprised between 1.5 and 2.5 mag. There is a significant negative slope for the DIB<sup>862</sup><sub>norm</sub>–β relationship in this regime of extinction. An orthogonal distance regression linear fit is superimposed. It assumes errors of 0.01 Å mag<sup>-1</sup> on DIB<sup>862</sup><sub>norm</sub> and of 0.02 on β.

Finally, Fig. 9 summarizes a series of results of the type shown in Fig. 8, here for 0.3 mag extinction intervals distributed between A<sub>V</sub> = 0.5 and A<sub>V</sub> = 4.5 mag, for the two DIB catalogues, and for the two β maps. At high extinction, the slope becomes very uncertain due to scarce data. At low extinction, the difference between individual LOSs and spatial bins has a strong effect. Between 1 and 3 magnitudes of extinction, the different determinations are all compatible and show a negative slope of the DIB<sup>862</sup><sub>norm</sub> versus β relationship. In summary, our results suggest that, despite a wide dispersion around the general trend, the as-yet-little-studied regime of low to moderate extinction is characterized by opposite dependences between DIB<sup>862</sup><sub>norm</sub> and R<sub>V</sub> and between DIB<sup>862</sup><sub>norm</sub> and β. Due to the very strong variability of the DIB-to-extinction ratio associated with the extinction density (Sect. 2), these trends are hard to detect unless massive datasets are available.

## 5. DIBs and 220 nm UV absorption bump

As mentioned in Sect. 1, contradictory conclusions were drawn on the relationship between the DIBs and the UV bump. We used the compilation of Xiang et al. (2017) and the 3D extinction maps in an attempt to re-examine this point. The compilation was initiated by Xiang et al. (2011), who concluded on the absence of a correlation between reddening-normalized DIBs and the reddening-normalized bump strength, based on poor Pearson correlation coefficients Pr. The strength of the bump was defined as the area above the linear profile in the classical description of Fitzpatrick & Massa (1990), namely πc<sub>3</sub>γ<sup>-1</sup>, where c<sub>3</sub> and γ are strength and width coefficients of the Drude profile fitted to the bump. The authors used the E(λ-V)/E(B-V) formulation of the extinction law, in the same way as Desert et al. (1995). Later on, and based on an augmented compilation of DIBs and extinction law parameters, Xiang et al. (2017) used the A(λ)/A(V) form of the extinction law and defined a new decomposition of the law. Based on this new definition, the authors confirmed the absence of correlation between extinction-normalized DIBs and the newly defined bump area. Here we re-examined the dependences on the bump in two parallel ways, using reddening-normalized DIBs and the E(λ-V)/E(B-V) formulation of the extinction law on the one hand (hereafter FORM1), and extinction-normalized DIBs and the A(λ)/A(V) form of the extinction law on the other hand (hereafter FORM2). The bump height is proportional to c<sub>3</sub>γ<sup>-2</sup> in both cases, but c<sub>3</sub> is not defined in the same way.

As we already discussed, LOS-integrated data cannot give rise to tight correlations, due to cloud multiplicity. Our goal was to use the LOS structure as it shows up in the 3D dust map to infer the potential level of degradation of a DIB-bump dependence due to the existence of clouds with different properties.
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Illustration of the dependence on the extinction regime of the relationship between the normalized DIB and $\beta$. Plotted are the normalized ratios as a function of the LOS extinction, color-coded for the $\beta$ value for each LOS. Marker sizes increase with target location distance to the Plane. The two estimates for $\beta$ are shown (left and right). Below $Av \approx 3$ mag (left part of the graphs), decreasing $DIB_{norm}^{862}$ (vertical descending lines) corresponds to an increase of $\beta$ (color from red to green). Above $\approx 3$ mag (right part of the graphs), there is an opposite trend: decreasing $DIB_{norm}^{862}$ (again vertical descending lines) corresponds to a decrease of $\beta$ (color from violet to green or cyan).</figcaption></figure>

<table><caption>**Table 1.** Line-of-sight flags F adopted for stars from the Xiang et al. (2017) catalogue.</caption><thead><tr><th>Star</th><th>Flag</th><th>Star</th><th>Flag</th><th>Star</th><th>Flag</th><th>Star</th><th>Flag</th></tr></thead><tbody><tr><td>HD 2905</td><td>0</td><td>HD 42087</td><td>0</td><td>HD 143018</td><td>1</td><td>HD 152247</td><td>0</td></tr><tr><td>HD 15558</td><td>0</td><td>HD 46056</td><td>1</td><td>HD 143275</td><td>1</td><td>HD 152248</td><td>0</td></tr><tr><td>HD 15570</td><td>0</td><td>HD 46150</td><td>1</td><td>HD 143567</td><td>1</td><td>HD 152249</td><td>0</td></tr><tr><td>HD 15629</td><td>0</td><td>HD 46202</td><td>1</td><td>HD 144217</td><td>1</td><td>HD 154445</td><td>1</td></tr><tr><td>HD 16691</td><td>0</td><td>HD 46223</td><td>1</td><td>HD 144470</td><td>1</td><td>HD 162978</td><td>1</td></tr><tr><td>HD 21291</td><td>0</td><td>HD 47129</td><td>1</td><td>HD 145502</td><td>1</td><td>HD 164492</td><td>1</td></tr><tr><td>HD 21483</td><td>1</td><td>HD 47240</td><td>1</td><td>HD 145554</td><td>1</td><td>HD 164794</td><td>1</td></tr><tr><td>HD 23060</td><td>0</td><td>HD 48099</td><td>1</td><td>HD 146001</td><td>1</td><td>HD 165052</td><td>1</td></tr><tr><td>HD 27778</td><td>1</td><td>HD 48434</td><td>0</td><td>HD 146029</td><td>1</td><td>HD 167971</td><td>1</td></tr><tr><td>HD 30614</td><td>1</td><td>HD 53974</td><td>0</td><td>HD 146416</td><td>1</td><td>HD 168076</td><td>1</td></tr><tr><td>HD+31643</td><td>1</td><td>BD+60497</td><td>0</td><td>HD 147165</td><td>1</td><td>HD 168112</td><td>1</td></tr><tr><td>HD 34078</td><td>0</td><td>BD+60594</td><td>0</td><td>HD 147701</td><td>1</td><td>HD 183143</td><td>0</td></tr><tr><td>HD 36879</td><td>0</td><td>HD 99872</td><td>1</td><td>HD 147888</td><td>1</td><td>HD 185418</td><td>0</td></tr><tr><td>HD 37367</td><td>1</td><td>HD 122879</td><td>1</td><td>HD 147889</td><td>1</td><td>HD 190603</td><td>1</td></tr><tr><td>HD 37903</td><td>1</td><td>HD 123008</td><td>0</td><td>HD 147933</td><td>1</td><td>HD 192281</td><td>1</td></tr><tr><td>HD 38087</td><td>1</td><td>HD 142096</td><td>1</td><td>HD 149757</td><td>1</td><td>HD 193322</td><td>1</td></tr><tr><td>HD 38131</td><td>0</td><td>HD 142165</td><td>1</td><td>HD 152233</td><td>0</td><td>HD 193682</td><td>0</td></tr><tr><td>HD 40893</td><td>0</td><td>HD 142315</td><td>1</td><td>HD 152236</td><td>0</td><td>HD 197770</td><td>0</td></tr><tr><td>HD 41117</td><td>0</td><td>HD 142378</td><td>1</td><td>HD 152246</td><td>0</td><td>HD 198478</td><td>1</td></tr></tbody></table>

For all targets from the Xiang et al. (2017) compilation, we inspected visually the series of clouds along the LOS, using the image of the extinction in a vertical plane containing the Sun and the target, computed from the 3D extinction distribution of Vergely et al. (2022). Admittedly somewhat arbitrarily, we defined a flag F according to the cloud distribution between the Sun and the star. If most of the extinction is produced by a single compact cloud or by a cloud group with limited extension (e.g. smaller than ≈100 pc), then F = 1. If, on the contrary, the extinction is produced by series of clouds distributed over long distances or to two groups separated by a large distance (e.g.

more than ≈500 pc), F = 0. This produces two groups of LOSs of about the same size. The series of targets and their attributed flags is listed in Table 1. Two examples are shown in Fig. 10 for the target star HD 192216 ($l = 88.92°, b = 3.02°, d = 1100$ pc) as an example of a flag F = 0, and for HD 143567 ($l = 350.87°, b = 22.68°, d = 135$ pc) as an example of a flag F = 1.

### 5.1. Bump height

We first selected the strong 5780 Å DIB. In Fig. 11 the reddening-normalized (respectively extinction-normalized) EWs are
<figure><img src="image_8.png" /></figure>

**Fig. 7.** Extinction-normalized 862 nm DIB EWs as a function of the spectral index $\beta$ from Casandjian et al. (2022), for DIB with quality flag 0 data from the second catalogue of Gaia Collaboration (2023a), and Galactic latitudes between 10° and 20°. Superimposed is a linear fit for this catalogue selection.

<figure><img src="image_9.png" /></figure>

**Fig. 8.** Extinction-normalized 862 nm DIB EW for selected data from the second catalogue (Gaia Collaboration 2023a) as a function of the sky-interpolated Planck spectral index $\beta$ (see text for selection). The colour scale refers to the distance to the Plane of the spatial bin, and excludes distances below 200 pc. Extinction values are restricted to the 1.5 to 2.5 mag interval. An orthogonal distance regression linear fit is superimposed (see text).

displayed as a function of the bump height in the FORM1 (respectively FORM2) representations of the extinction law. Adjusted linear fits are shown for all stars (top), flag F = 1 stars (middle), and flag F=0 stars (bottom). Uncertainties on the DIB equivalent widths from Xiang et al. (2017) are shown at the top. The displayed fitted line does not take any error bars into account, errors on the reddening (or on the extinction) being too uncertain. Instead, the dispersion around the fitted line is used as an estimate of the mean prevailing uncertainty. For this particular DIB, the Pearson correlation coefficients are 0.54 ± 0.09, 0.67 ± 0.10, and 0.01 ± 0.17 respectively for 94, 59, and

<figure><img src="image_10.png" /></figure>

**Fig. 9.** Ranges of gradients for extinction-normalized DIB EWs as a function of $\beta$ for data from the two catalogues and in various regimes of LOS extinctions. Errors are $1\sigma$ uncertainties on the slopes. There are differences at low extinctions due to the very different types of data (individual LOSs vs spatial bins) and at high extinctions due to the small number of data points. Between 1 and 3 magnitudes of extinction, all determinations are compatible and show a negative slope.

<figure><img src="image_11.png" /></figure>

**Fig. 10.** Illustration of the attribution of a flag F to a given star from the Xiang et al. (2017) list. Shown are images of the dust distribution in vertical planes containing the Sun and the star. The dark areas correspond to high dust opacity. The LOS to the target (green marker) is shown as a green dotted line. HD 143567 (top) is attributed a flag F=1 since most of the extinction is generated in a single very dense nearby cloud. The more distant HD 192216 (bottom) is attributed a flag F=0 since at least two dense clouds separated by a large distance contribute to the extinction.
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Comparisons between the σ-type 5780 Å DIB EW and the UV bump height, under formalism 1 of the extinction law (a–c) and formalism 2 (d–f). Data are from the compilation of Xiang et al. (2017). (a) and (d) show the bump height vs the reddening-normalized EW for all targets and the adjusted linear fit. The colour-coding refers to the extinction (from yellow for the lowest extinctions to black for the highest ones). Errors bars on EWs are indicated. (b and e) Same as (a) and (d), for flag F=1 stars only. (c and f) Same as (a) and (d) for flag F=0 only. The F=0 LOSs show no dependence, and lower the correlation based on all data. F=1 targets show a clearer positive correlation.</figcaption></figure>

35 stars for the first formalism. The estimated errors on the correlation coefficient are based on Student's t-distributions (i.e. assuming normal distributions for input data). In the case of the second formalism, these coefficients become 0.51 ± 0.09, 0.69 ± 0.10, and -0.14 ± 0.17. The slopes and $1\sigma$ uncertainties are 2.0 ± 0.3 × 10⁻³, 2.7 ± 0.4 × 10⁻³, and 0.02 ± 0.4 × 10⁻³ mag Å⁻¹, respectively, for formalism 1, and 2.0 ± 0.3 × 10⁻³, 2.6 ± 0.4 × 10⁻³, and -0.5 ± 0.6 × 10⁻³ mag Å⁻¹, respectively, for formalism 2. These results demonstrate that, in both cases, the selection of LOSs closer to the ideal case of a mono-cloud situation has a strong effect. For a LOS characterized by cloud multiplicity (F = 0), the UV bump height is about constant, and there is a reduced dispersion around the linear fit, which we interpret as a consequence of the averaging of the cloud properties. This series of data produces the central plateau in the figure displaying all the assembled data. On the contrary, restricting the data to those LOSs closer to the mono-cloud situation, there is a pronounced positive slope and a large variability around the linear fit. This shows the actual existence of a positive correlation between the strength of this particular DIB and the bump height.

While re-examining in the same way the potential relationships with the bump height for other DIBs than the 5780 Å band, we were surprised to find a decreasing correlation strength from the σ-type DIBs (e.g. like the 5780 Å DIB and the 6284 Å DIB) to ζ-type DIBs (e.g. like the 5797 Å DIB or the 5850 Å DIB). The former (σ) bands are well known for being significantly

reduced in dense, molecular clouds, while the latter are not. The σ-to-ζ ratio 5780/5797 is often used as an indicator of this difference, and is shown to increase with the UV radiation field prevailing in diffuse clouds (e.g. Krelowski & Walker 1987; Cami et al. 1997).

We show in Fig. 12, as an example, the results for the ζ-type 5850 Å DIB. The relationship between the reddening-normalized DIB EW and the bump height is dramatically different from that of the 5780 Å DIB, since no trend is detectable, whatever the LOS selection, in formalism 1, and a very weak trend is detected for formalism 2 (not shown here). Tables 2 and 3 list the correlation coefficients of the adjusted linear fits for the dependence of the bump height on reddening-normalized and extinction-normalized DIB EW. We limited the choice of DIBs to those with more than 12 LOSs. Coefficients are also listed for flag F=1 only and flag F=0 only. The estimated standard error on the correlation coefficient based on a Student's t-distribution (i.e. assuming normal distributions) is also listed. The DIBs are listed by decreasing Pearson correlation coefficient computed for all flags F=1 and F=0 LOS and formalism 1 (i.e. from 0.54 down to 0.07). We also see that the improvement (respectively degradation or total disappearance) of the correlation when considering only flag F=+1 (respectively flag F=0) we found for the 5780 Å DIB is also found for all other DIBs, which reinforces the case for actual correlations. We note that all correlations are stronger for formalism 2. However, for both formalisms, the sequence of
<table><caption>**Table 2.** Reddening-normalized DIB EW and bump height (formalism 1).</caption><thead><tr><th>DIB (Å)</th><th>N (ALL)</th><th>Pr</th><th>Err(Pr)</th><th>N (F = 1)</th><th>Pr</th><th>Err(Pr)</th><th>N (F = 0)</th><th>Pr</th><th>Err(Pr)</th></tr></thead><tbody><tr><td>5780</td><td>94</td><td>0.54</td><td>0.09</td><td>59</td><td>0.67</td><td>0.10</td><td>35</td><td>0.01</td><td>0.17</td></tr><tr><td>6284</td><td>64</td><td>0.52</td><td>0.11</td><td>41</td><td>0.57</td><td>0.13</td><td>23</td><td>0.26</td><td>0.21</td></tr><tr><td>6196</td><td>76</td><td>0.50</td><td>0.10</td><td>51</td><td>0.56</td><td>0.12</td><td>25</td><td>0.26</td><td>0.20</td></tr><tr><td>6614</td><td>58</td><td>0.56</td><td>0.11</td><td>44</td><td>0.63</td><td>0.12</td><td>14</td><td>0.16</td><td>0.28</td></tr><tr><td>6269</td><td>59</td><td>0.32</td><td>0.12</td><td>36</td><td>0.37</td><td>0.16</td><td>23</td><td>0.17</td><td>0.22</td></tr><tr><td>5797</td><td>90</td><td>0.16</td><td>0.11</td><td>57</td><td>0.17</td><td>0.13</td><td>33</td><td>-0.02</td><td>0.18</td></tr><tr><td>6379</td><td>43</td><td>0.25</td><td>0.15</td><td>34</td><td>0.27</td><td>0.17</td><td>9</td><td>0.06</td><td>0.38</td></tr><tr><td>5850</td><td>56</td><td>0.07</td><td>0.14</td><td>35</td><td>0.08</td><td>0.17</td><td>21</td><td>-0.17</td><td>0.22</td></tr></tbody></table>

**Notes.** N is the number of LOSs, Pr is the Pearson correlation coefficient, Err(Pr) is the uncertainty on the Pearson coefficient coefficient (see text). All is for flags 0 and 1, F = 0 for flag zero only, F = 1 for flag 1 only.

<table><caption>**Table 3.** Extinction-normalized DIB EW and bump height (formalism 2).</caption><thead><tr><th>DIB (Å)</th><th>N (ALL)</th><th>Pr</th><th>Err(Pr)</th><th>N (F = 1)</th><th>Pr</th><th>Err(Pr)</th><th>N (F = 0)</th><th>Pr</th><th>Err(Pr)</th></tr></thead><tbody><tr><td>5780</td><td>94</td><td>0.51</td><td>0.09</td><td>59</td><td>0.69</td><td>0.10</td><td>35</td><td>-0.14</td><td>0.17</td></tr><tr><td>6284</td><td>64</td><td>0.51</td><td>0.11</td><td>41</td><td>0.56</td><td>0.13</td><td>23</td><td>0.27</td><td>0.21</td></tr><tr><td>6196</td><td>76</td><td>0.58</td><td>0.09</td><td>51</td><td>0.66</td><td>0.11</td><td>25</td><td>0.20</td><td>0.20</td></tr><tr><td>6614</td><td>58</td><td>0.66</td><td>0.10</td><td>44</td><td>0.73</td><td>0.10</td><td>14</td><td>0.10</td><td>0.29</td></tr><tr><td>6269</td><td>59</td><td>0.42</td><td>0.12</td><td>36</td><td>0.44</td><td>0.15</td><td>23</td><td>0.26</td><td>0.21</td></tr><tr><td>5797</td><td>90</td><td>0.34</td><td>0.10</td><td>57</td><td>0.41</td><td>0.12</td><td>33</td><td>-0.13</td><td>0.18</td></tr><tr><td>6379</td><td>43</td><td>0.41</td><td>0.14</td><td>34</td><td>0.42</td><td>0.16</td><td>9</td><td>0.29</td><td>0.36</td></tr><tr><td>5850</td><td>56</td><td>0.29</td><td>0.13</td><td>35</td><td>0.35</td><td>0.16</td><td>21</td><td>-0.13</td><td>0.23</td></tr></tbody></table>

**Notes.** Same as Table 2 for formalism 2.

DIBs for decreasing trends shares similarities with the hierarchy of DIBs associated with the $\sigma$ to $\zeta$ types, or to the PC2 components of Ensor et al. (2017), found by the authors to represent the UV-field and is taken from their Fig. 7.

## 5.2. Estimates of bump height based on DIBs

The above results suggest that the response of a particular DIB to environmental conditions (here the strength of the UV radiation field) modifies the link with the bump height. We introduced a UV factor defined by the classical EW(DIB 5780)/EW(DIB 5797) ratio as an additional factor of the relationship between the bump and the reddening-normalized DIB. We show in Fig. 13a the relationship between the bump height and this UV factor. For flag F=0 stars there is no trend. For flag F=1 stars there may be some link, visible for extreme low and high values of the height; however, obviously a link far from a linear relationship. We then fitted the UV bump height (i.e. the $c_3\gamma^{-2}$ coefficient taken from the Xiang et al. 2017, Table 1) to the product of two power laws, one on the reddening-normalized DIB and the other on the UV factor (i.e. the $\lambda5780/\lambda5797$ ratio):

$$ \text{Height(model)} = c_3/\gamma^2 = K (\text{DIB}_{\text{norm}})^{A1} \cdot (\text{EW5780}/\text{EW5797})^{A2}. \quad (2) $$

Figure 13 shows the fitted model bump height for the case of the 5850 Å DIB, based on the use of all 5850 Å data with flags F=0 and F=1 (part b), and based on flag F=1 stars only (part c). The K, A1, and A2 coefficients are 0.35, 0.16, and 1.3 and 0.45, 0.23, and 0.88, respectively. Although this is certainly a far from perfect representation, we can see that this simplistic model reproduces the high and low values of the bump.

We pursued this exercise for other DIBs, namely the 5780, 5797, 6196, and 6614 Å DIBs. The results were approximately the same. We also tested combinations of power laws for several DIBs. The advantage is to increase the number of targets for the fitting and the comparison since EWs are not available for all bands. Finally, what we discovered from these various attempts is that $\sigma$ DIBs, such as 5780 Å, can be used directly and are a kind of reference; instead $\zeta$ DIBs, while corrected for their specific UV ratio ($\sigma/\zeta$), may also provide some estimates of the bump height. As an example, Fig. 13d shows the adjusted model based on a combination of the 5780 Å and 5850 Å DIBs. Its formulation is

$$ \text{Height(model)} = 0.58 \cdot (\text{DIB}(5780)_{\text{norm}})^{0.27} \cdot (5780/5797)^{0.073} \\ \cdot (\text{DIB}(5850)_{\text{norm}})^{0.026}. \quad (3) $$

This simplistic model reproduces the high and low values of the bump, and may provide a prior coefficient for the bump height. The $1\sigma$ deviation from the model is 0.5, compared with height values ranging from 2 to 5. Further extinction law measurements would help to refine and to correct this approach purely based on DIBs.

## 5.3. Bump width and bump area

We made the same comparisons and linear fits, this time as a function of the bump width and the bump area. The correlations with the bump width are displayed in Fig. 14 for the strong 5780 Å DIB. The Pr factors vary from -0.46 ± 0.09 to -0.42 ± 0.12 and reduced $\chi^2$ values vary from 8.9 to 7.4, before and after
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Same as Fig. 11 (a–c) but for the ζ-type 5850 Å DIB and formalism 1. Top: all flags. Middle: flag F=1. Bottom: flag F=0. No trend is detected, neither for all stars, nor for F=1 stars only.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Illustration of a composite model of the bump height based on DIBs and the UV ratio EW(5780)/EW(5797). (a) Bump height and UV ratio 5780/5797. (b) Composite power-law model (Eq. (2)) based on the 5850 Å DIB and the UV ratio, for all stars with DIB 5850 Å data (55 targets). (c) Same as panel b, for a model based on flag F=1 stars only. (d) Composite power-law model (Eq. (2)) based on DIB 5780, DIB 5850, and the UV ratio, adjusted on flag F=1 stars (89 targets). Flag F=0 stars are in red, Flag F=1 stars in violet. In panels b, c, and d, the *y* = *x* curves are superimposed.</figcaption></figure>

the LOS selection for formalism 1. The Pr factors are almost constant at -0.42 ± 0.09 and -0.38 ± 0.12 and reduced χ² constant at 7.4 for formalism 2. The slopes and their 1σ uncertainties are respectively -3.2 ± 0.6 × 10⁻⁴, -3.1 ± 0.9 × 10⁻⁴, and -3.6 ± 0.8 × 10⁻⁴ µm⁻¹ mag Å⁻¹ for 94, 59, and 35 targets with flags F=0 and F=1, only F=1, and only F=0 for formalism 1. The numbers become -9 ± 2 × 10⁻⁴, -8.4 ± 3 × 10⁻⁴, and -12.2 ± 3 × 10⁻⁴ µm⁻¹ mag Å⁻¹ for formalism 2. The uncertainties on the slopes are larger than in the case of the bump height, but all results favour an inverse trend between the reddening normalized DIB and the width. Results for correlation factors and uncertainties, for the same DIBs and in the same order as in Table 2 for the height, are listed in Table 4. We listed the results for formalism 1 only. Results for formalism 2 are similar. There is again a net decrease in the correlation strengths (here anti-correlation strengths) from σ to ζ DIBs. However, at variance with the case of the height, the strength of the correlation for the width does not seem to be linked to the LOS type.

The results for the bump area are shown in the case of formalism 1 in Fig. 15. The slopes and their 1σ uncertainties are 7.6 ± 3.7 × 10⁻⁴, 14 ± 5 × 10⁻⁴, -13 ± 5 × 10⁻⁴ mag µm⁻¹ Å⁻¹ respectively for 94, 59, 35 targets with flags F=0 and 1, F=1, and F=0 for formalism 1. The Pearson correlation coefficients are 0.21 ± 0.10, +0.39 ± 0.12, and -0.44 ± 0.16 respectively for the 94, 59, and 35 targets. The uncertainties on the slopes of the linear fits are much larger than in the case of the bump height and are still significantly larger than for the bump width. At variance with the case of the bump height and the bump width, there is no

constancy of the slope sign since there are two opposite trends for flag F=1 and flag F=0, and the links are unclear. Results for correlation factors and uncertainties, for the same DIBs and in the same order as in Table 2 for the height are listed in Table 5.

All together, the results are in agreement with those from Desert et al. (1995) for reddening-normalized DIBs, namely a
<table><caption>**Table 4. Reddening-normalized DIB EW and bump width.**</caption><thead><tr><th>DIB (Å)</th><th>N (ALL)</th><th>Pr</th><th>Err(Pr)</th><th>N (F = 1)</th><th>Pr</th><th>Err(Pr)</th><th>N (F = 0)</th><th>Pr</th><th>Err(Pr)</th></tr></thead><tbody><tr><td>5780</td><td>94</td><td>-0.46</td><td>0.09</td><td>59</td><td>-0.42</td><td>0.12</td><td>35</td><td>-0.59</td><td>0.14</td></tr><tr><td>6284</td><td>64</td><td>-0.50</td><td>0.11</td><td>41</td><td>-0.53</td><td>0.14</td><td>23</td><td>-0.41</td><td>0.20</td></tr><tr><td>6196</td><td>76</td><td>-0.23</td><td>0.11</td><td>51</td><td>-0.18</td><td>0.14</td><td>25</td><td>-0.42</td><td>0.19</td></tr><tr><td>6614</td><td>58</td><td>-0.23</td><td>0.13</td><td>44</td><td>-0.20</td><td>0.15</td><td>14</td><td>-0.72</td><td>0.20</td></tr><tr><td>6269</td><td>59</td><td>-0.10</td><td>0.13</td><td>36</td><td>-0.09</td><td>0.17</td><td>23</td><td>-0.22</td><td>0.21</td></tr><tr><td>5797</td><td>90</td><td>-0.13</td><td>0.10</td><td>57</td><td>-0.14</td><td>0.13</td><td>33</td><td>-0.13</td><td>0.18</td></tr><tr><td>6379</td><td>43</td><td>-0.06</td><td>0.16</td><td>34</td><td>0.02</td><td>0.18</td><td>9</td><td>-0.54</td><td>0.32</td></tr><tr><td>5850</td><td>56</td><td>-0.07</td><td>0.14</td><td>35</td><td>-0.08</td><td>0.17</td><td>21</td><td>-0.03</td><td>0.23</td></tr></tbody></table>

**Notes.** Same as Table 2, but for the bump width.

<table><caption>**Table 5. Reddening-normalized DIB EW and bump area.**</caption><thead><tr><th>DIB (Å)</th><th>N (ALL)</th><th>Pr</th><th>Err(Pr)</th><th>N (F = 1)</th><th>Pr</th><th>Err(Pr)</th><th>N (F = 0)</th><th>Pr</th><th>Err(Pr)</th></tr></thead><tbody><tr><td>5780</td><td>94</td><td>0.21</td><td>0.10</td><td>59</td><td>0.39</td><td>0.12</td><td>35</td><td>-0.44</td><td>0.16</td></tr><tr><td>6284</td><td>64</td><td>0.23</td><td>0.12</td><td>41</td><td>0.32</td><td>0.15</td><td>23</td><td>-0.11</td><td>0.22</td></tr><tr><td>6196</td><td>76</td><td>0.37</td><td>0.11</td><td>51</td><td>0.50</td><td>0.12</td><td>25</td><td>-0.13</td><td>0.21</td></tr><tr><td>6614</td><td>58</td><td>0.43</td><td>0.12</td><td>44</td><td>0.55</td><td>0.13</td><td>14</td><td>-0.50</td><td>0.25</td></tr><tr><td>6269</td><td>59</td><td>0.30</td><td>0.13</td><td>36</td><td>0.39</td><td>0.16</td><td>23</td><td>-0.02</td><td>0.22</td></tr><tr><td>5797</td><td>90</td><td>0.10</td><td>0.11</td><td>57</td><td>0.11</td><td>0.14</td><td>33</td><td>-0.12</td><td>0.18</td></tr><tr><td>6379</td><td>43</td><td>0.25</td><td>0.15</td><td>34</td><td>0.35</td><td>0.16</td><td>9</td><td>-0.38</td><td>0.35</td></tr><tr><td>5850</td><td>56</td><td>0.08</td><td>0.14</td><td>35</td><td>0.09</td><td>0.17</td><td>21</td><td>-0.13</td><td>0.23</td></tr></tbody></table>

**Notes.** Same as Table 2, but for the bump area.

positive correlation with the bump height and a negative correlation with the bump width, provided one uses σ-type DIBs and mono-cloud-type LOSs. They also agree with the absence of a clear global trend for the bump area found by Xiang et al. (2017) for extinction-normalized DIBs, and by Xiang et al. (2011) for reddening-normalized DIBs. This is consistent with the fact that the area is close to the product of the width and the height. At variance with the results of Xiang et al. (2017), we obtain some positive correlations for the area, although weaker than for the height, when we restrict to flag F=1 and if we use σ-type DIBs. We note that the area and the width are more difficult to measure than the height; in particular, they depend on the choice of the underlying linear part of the extinction law, and this may partly explain the results.

## 6. Summary and discussion

The existence of links between DIBs and dust properties has been the subject of numerous works, mainly based on spectra of hot strongly reddened stars. The results have been limited, and sometimes apparently contradictory. Several new kinds of studies are now permitted by the availability of 3D maps and massive datasets from Gaia for the 862 nm DIB, a diffuse band likely representative of most other DIBs: (i) the use of local diagnostics instead of LOS-integrated diagnostics, (ii) the use of massive information on LOSs with low extinction and on dust-poor areas, and (iii) the cross-match of the information contained in 3D maps with independent information from different surveys by means of integration in the 3D maps along any LOS. Moreover, for any DIB, the existence of 3D dust extinction density maps permits (iv) the use of the LOS structure to weigh correlative studies. We used these four new possibilities and, over

the course of searches for correlations, we obtained several new results.

In Sect. 2, we used (i) and (ii) and examined the relative increase in the 862 nm DIB density with respect to the extinction density occurring in low-extinction clouds, a trend detected in the 3D maps (Cox et al. 2024). We extended the study to LOS-integrated data of the two 862 nm DIB catalogues. We fitted the ratio of local (respectively LOS-integrated) DIB and extinction to a power law of the local (respectively integrated extinction) in the low-extinction regime ($A_V$ between 0.12 and 1.5 mag). The fitted power-law coefficients, comprised between 0.57 and 0.74 for LOS data from the two catalogues and for local (3D map) data, were compared with those found by Lan et al. (2015) for a series of 20 DIBs and for the same range of extinctions, based on SDSS stacked spectra for high-latitude directions. We found them in agreement with the coefficients of several other DIBs. This implies that the 862 nm DIB behaves similarly to several other DIBs in the low-extinction regime, with a power-law coefficient at the lower end of the Lan et al. (2015) series (i.e. presenting one of the highest relative increases with respect to extinction in dust-poor clouds, the lowest values of $\gamma$). This study, based on local values and massive data, validates the behaviour found with SDSS. The variability of the power-law $\gamma$ coefficients found by Lan et al. (2015), to date poorly discussed, contains information on the presence of the carriers far from the dense clouds. We note that the hierarchy of DIBs associated with these coefficients is significantly different from that associated with component 2 of the principal component analysis (PCA) of Ensor et al. (2017), assumed to represent the influence of the radiation field. It may be the manifestation of one of the unidentified components in their PCA analysis. Interestingly, the ranking of the 862 nm DIB, based on our results, among the 20 DIBs from the Lan et al. (2015) study, includes a comparison with the
<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Same as Fig. 11, but for the dependence on the bump width, for formalism 1 and 2. Top: all targets. Middle: Flag F=1 targets. Bottom: flag F=0 targets.</figcaption></figure>

broad 4430 Å DIB. This comparison implies that the 4430 Å DIB density increases more moderately in dust-poor clouds compared to the 862 nm DIB. This is consistent with the recent results of Zhao et al. (2023), who show a smaller extension of the 4430 Å DIB carriers at a large distance from the Plane in comparison with the 862 nm DIB carriers. This opens the possibility of using the Lan et al. (2015) coefficients to rank the scale heights of the various DIB carriers.

<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Same as Fig. 11, but for the dependence on the bump area and formalism 1. Top: all flags. Middle: F=1. Bottom: F=0.</figcaption></figure>

In Sect. 3, we used (iii), and investigated the link between the 862 nm DIB strength normalized to the extinction on the one hand, and the proxy $R'_V$ for the total-to-selective extinction ratio $R_V = A_V/E(B-V)$ derived by Schlafly et al. (2016) for a large number of SDSS/APOGEE target stars. One of the results from this work is the existence of large-scale variations of $R_V$ over the sky. For all the targets from the Schlafly et al. (2016) catalogue, we integrated from the Sun to the target within the 3D map of 862 nm DIB carrier density on the one hand, and within the corresponding 3D map of extinction density on the other hand (maps from Cox et al. 2024). The results were combined to obtain estimates of the extinction-normalized DIB$_{norm}^{862}$. Based on this dataset, and despite a large scatter due to the poor resolution of the 3D DIB map, we detected a positive correlation between DIB$_{norm}^{862}$ and $R'_V$ for the low- to moderate-reddening regime $A_V \lesssim 2.5$ mag. Combining with measurements of Ramírez-Tannus et al. (2018) and their estimates of $R_V$ for stars in the M17 cloud with $A_V \ge \approx 4$ mag, it appears that at higher extinctions the trend is suppressed and reversed. This type of bi-modal regime may explain the absence of previous clear results on the DIB-$R_V$ relationship when data from the two regimes are mixed together.

In Sect. 4, we also used (iii), and investigated the link between DIB$_{norm}^{862}$ and the dust emission spectral index $\beta$. To do so, we used Planck $\beta$ sky maps and selected targets from the two 862 nm DIB catalogues located outside the dust layer (i.e. with negligible interstellar dust beyond the star). This is necessary since the spectral index $\beta$ refers to the totality of the dust in the chosen direction. Such targets are preferentially found at high latitude. We integrated in the best resolution 3D extinction maps
to estimate the extinction of the targets. Despite a large scatter, we found a negative correlation between $DIB_{norm}^{862}$ and $\beta$ for $A_V \le 2.5$ mag, and a disappearance of the trend above this threshold (i.e. for high extinctions). There may be a hint for reversal; however, the number of appropriate targets with high extinctions and located beyond the dust layer becomes very small. The trend at low-moderate $A_V$ is more apparent when using recent estimates of $\beta$ by Casandjian et al. (2022).

The two trends, for $R'_V$ and $\beta$, in the low- to moderate-extinction regime, are found to be fully consistent with the anti-correlation between $\beta$ and $R'_V$ found by Schlafly et al. (2016), which reinforces the likelihood of these two independent results. The opposite trends we have derived are particularly interesting in view of the models and results of Zelko & Finkbeiner (2020), namely that the $\beta$ and $R'_V$ anti-correlation is induced by variations of the relative fractions of carbonaceous and silicate dust grains, with a higher $R'_V$ (and lower $\beta$) being associated with a higher fraction of carbonaceous grains. In this case, because it is now generally admitted that DIB carriers are carbonaceous macro-molecules, their density relatively to the total dust content, measured by $DIB_{norm}^{862}$, is naturally expected to correlate with $R'_V$ and anti-correlate with $\beta$, which is what we derive. The trends are no longer valid for LOSs crossing opaque clouds, where accretion onto grains and grain coagulation dominate, and the skin effect for the DIB carriers starts to play a role. Regarding the carbon to silicon abundance, we note that Cox et al. (2024) found evidence from the 3D maps that higher $DIB_{norm}^{862}$ values are spatially associated with regions where C-rich ejecta from AGBs increase relative to the O-rich fluxes. From the point of view of extinction and emission laws, this shows that useful predictions on the dust composition may be obtained based on DIB measurements.

In Sect. 5, we used (iv) and re-examined the link between the 220 nm UV bump and the DIBs, using data from the intensive compilation of Xiang et al. (2017). The results, so far, have been contradictory regarding the existence of a link. One of the reasons for the difficulty of such a study is the cloud multiplicity and the resulting mixing of the cloud (and DIB carrier) properties along the LOS of the target star, especially in the case of distant, highly reddened hot stars generally used for UV extinction law measurements and DIB extraction. Here we added to the previous searches an additional criterion based on the 3D dust maps and reflecting the cloud singleness (item iv in the above list). If the LOS is characterized by a unique dense cloud or group of clouds at short distance from each other, the flag is set to F=1. If, on the contrary, there are two or more groups of clouds contributing to the extinction (and to the DIB EW) separated by distances of the order of several hundred parsecs or more, one may expect different properties of the IS matter and the flag is set to F=0. This criterion is defined visually, and it should be used with limited ambition. However, our goal was to detect (or not) an improvement in the link between DIBs and the bump when selecting flag F=1 LOS by comparison with the full dataset or with flag F=0 LOS. This approach allowed us to demonstrate the existence of a positive relationship between the reddening-normalized (or extinction-normalized) DIB EWs and the height of the bump for those DIBs that are of $\sigma$ type, known to be favoured in UV-irradiated regions and to decrease in cloud cores (the UV field here is estimated by the 5780/5797 ratio). It confirmed, based on the larger number of LOSs from the Xiang et al. (2017) compilation, the previous conclusions drawn by Desert et al. (1995). The relationship disappears for DIBs such as $\lambda5797$ or $\lambda5850$ (i.e. of $\zeta$ type). This may explain the negative conclusion of Xiang et al. (2011) and Xiang et al. (2017)

regarding a global DIB-bump relationship associating all DIBs. We tested a simplistic predictive model of the bump height based on a combination of reddening-normalized DIBs and the 5780/5797 ratio. It provided some estimates of the bump height, especially in the case of particularly low, or particularly high bumps. On the other hand, we found a negative relationship between the reddening-normalized 5780 Å DIB EW and the bump width, similarly to Desert et al. (1995), although this trend is characterized by a larger scatter and a higher uncertainty on the slope compared to the case of the bump height. Finally, no general trend emerges from the comparisons between DIBs and the bump area, in agreement with the results of Xiang et al. (2017), although a weak positive trend is found for $\sigma$-type DIBs.

These results on the UV bump are surprising, in particular the fact that direct positive correlations between the bump height and DIBs, and direct anti-correlations between bump widths and DIBs, are found only for the DIBs that are the most prone to the skin effect (or $\sigma$ DIBs), the decrease in the carriers relative to the extinction in the dense cloud cores. Correlations may also be found in the case of $\zeta$ DIBs, but only after application of a corrective factor whose effect is to increase their strengths in LOS characterized by a high 5780 Å-to-5797 Å DIB ratio (or $\sigma$-to-$\zeta$ ratio). There may be evidence here for the building of the bump outside of the densest cores, and this point deserves further observations and modelling. In this respect, it is interesting to recall the result obtained by Fitzpatrick & Massa (1986) on the bump width. The authors inferred that the 220 nm bump is narrower in the diffuse ISM, and broader in dense quiescent clouds. Given that for $\sigma$-type DIBs, such as the 5780 Å DIB, the ratio $DIB_{norm}^{862}$ is found to increase when the bump width decreases (Fig. 14), the Fitzpatrick & Massa (1986) trend predicts, in turn, that $DIB_{norm}^{862}$ should increase from dense clouds to regions of lower density. This is what is statistically found based on 3D maps, and mono-cloud type LOSs, as discussed in Sect. 2, and adds a consistent link between the trends found between $DIB_{norm}^{862}$ and the bump on the one hand, and between $DIB_{norm}^{862}$ and the ISM type on the other hand.

It is beyond the scope of this observational work to discuss the implications in terms of the DIB carriers and species responsible for the UV bump. Our unique remark is that, despite the limited strength of all correlations, as expected from the integration along the LOSs, the results on bump height and width may be hints of a partial contribution to the bump of some specific molecular species and their specific lines, in addition to the global contribution of a large variety of species (e.g. PAHs). In this respect, it is interesting to note that recent works point to high densities of two-ring polycyclic aromatic hydrocarbons, their formation around AGBs, and their potential contribution to the bump. Stockett et al. (2023) found a mechanism that explains why cyanonaphthalene ($C_{10}H_7N$) can be orders of magnitude higher than predicted by models, following the discovery of quantities of this species and of indene in the Taurus Molecular Cloud with quantities orders of magnitude higher than expected (McGuire et al. 2021; Burkhardt et al. 2021). According to Stockett et al. (2023), their results challenge the widely accepted picture of rapid destruction of small PAHs in space. As a result of their advanced models, Krasnoukhov et al. (2022) found that under the conditions prevailing in the circumstellar envelopes of the asymptotic giant branch stars (AGBs), the formation of two-ring polycyclic aromatic hydrocarbons is favoured. On the other hand, based on a state-of-the-art laboratory experimental set-up, Fréreux et al. (2023) found that substituted naphthalene units, such as substructures of hydrogenated amorphous carbons,
are favoured in combustion processes and produce spectral sig-
natures around 217.5 nm. Further measurements of the bump
central wavelength and its potential association with the LOS
type may bring support to the existence of such contributions.

All together, these results based on 3D maps and surveys
reveal complex, but real relationships between the DIB carrier
densities and some dust extinction and emission parameters.
They suggest that future precise, massive measurements of var-
ious DIBs and of extinction, as well as massive measurements
of the UV part of the extinction law, may provide constraints on
the interplay between the carbonaceous macro-molecules at the
origin of DIBs and the size distribution of grains, on the spatial
variability of the dust composition, and in turn on the variabil-
ity of extinction and emission law parameters. DIBs are easy
to measure and their EWs are calibration-free measurements,
which are strong advantages. They are no longer restricted to
early-type stars, but extended to the far more numerous cool tar-
gets (e.g. Kos et al. 2013; Puspitarini et al. 2015; Zasowski et al.
2015), thanks to the huge progresses made over recent years
on stellar models and on their adjustments to data from recent
surveys and Gaia. Being distance-limited data, the use of dis-
tributed targets provides local measurements, allowing us in turn
to understand where and how the extinction and emission are
built, and to make predictions.

Acknowledgements. We thank our reviewer for his/her careful reading and help-
ful suggestions. This project has received funding from the European Union's
Horizon 2020 research and innovation programme under grant agreement
No 101004214 (EXPLORE project – https://explore-platform.eu). This
work has made use of data from the European Space Agency (ESA) mission Gaia
(https://www.cosmos.esa.int/gaia), processed by the Gaia Data Process-
ing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/
gaia/dpac/consortium). Funding for the DPAC has been provided by national
institutions, in particular the institutions participating in the Gaia Multilat-
eral Agreement. Funding for the SDSS and SDSS-II has been provided by
the Alfred P. Sloan Foundation, the Participating Institutions, the National Sci-
ence Foundation, the U.S. Department of Energy, the National Aeronautics and
Space Administration, the Japanese Monbukagakusho, the Max Planck Soci-
ety, and the Higher Education Funding Council for England. The SDSS Web
Site is http://www.sdss.org/. This research has made use of the VizieR cat-
alogue access tool, CDS, Strasbourg, France (DOI: 10.26093/cds/vizier).
The original description of the VizieR service is published in 2000, A&AS 143,
23.

References

Bertaux, J.-L., & Lallement, R. 2017, MNRAS, 469, S646

Burkhardt, A. M., Long Kelvin Lee, K., Bryan Changala, P., et al. 2021, ApJ,
913, L18

Cami, J., Sonnentrucker, P., Ehrenfreund, P., & Foing, B. H. 1997, A&A, 326,
822

Cami, J., Bernard-Salas, J., Peeters, E., & Malek, S. E. 2010, Science, 329, 1180

Campbell, E. K., Holz, M., Gerlich, D., & Maier, J. P. 2015, Nature, 523, 322

Casandjian, J.-M., Ballet, J., Grenier, I., & Remy, Q. 2022, ApJ, 940, 116

Cordiner, M. A., Linnartz, H., Cox, N. L. J., et al. 2019, ApJ, 875, L28

Cox, N. L. J., Cami, J., Kaper, L., et al. 2014, A&A, 569, A117

Cox, N. L. J., Cami, J., Farhang, A., et al. 2017, A&A, 606, A76

Cox, N. L. J., Vergely, J. L., & Lallement, R. 2024, A&A, 689, A38

Desert, F. X., Jenniskens, P., & Dennefeld, M. 1995, A&A, 303, 223

Ebenbichler, A., Smoker, J. V., Lallement, R., et al. 2024, A&A, 686, A50

Elyajouri, M., Lallement, R., Monreal-Ibero, A., Capitanio, L., & Cox, N. L. J.
2017, arXiv e-prints [arXiv:1702.00223]

Ensor, T., Cami, J., Bhatt, N. H., & Soddu, A. 2017, ApJ, 836, 162

Fan, H., Hobbs, L. M., Dahlstrom, J. A., et al. 2019, ApJ, 878, 151

Fitzpatrick, E. L., & Massa, D. 1986, ApJ, 307, 286

Fitzpatrick, E. L., & Massa, D. 1990, ApJS, 72, 163

Foing, B. H., & Ehrenfreund, P. 1994, Nature, 369, 296

Fréreux, J. N., Godard, M., Dartois, E., & Pino, T. 2023, A&A, 677, A149

Friedman, S. D., York, D. G., McCall, B. J., et al. 2011, ApJ, 727, 33

Gaia Collaboration (Schultheis, M., et al.) 2023a, A&A, 680, A38

Gaia Collaboration (Schultheis, M., et al.) 2023b, A&A, 674, A40

Geballe, T. R., Najarro, F., Figer, D. F., Schlegelmilch, B. W., & de La Fuente,
D. 2011, Nature, 479, 200

Hamano, S., Kobayashi, N., Kondo, S., et al. 2016, ApJ, 821, 42

Herbig, G. H. 1995, ARA&A, 33, 19

Hobbs, L. M., York, D. G., Snow, T. P., et al. 2008, ApJ, 680, 1256

Jenniskens, P., & Desert, F.-X. 1994, A&As, 106

Jones, A. P. 2016, Roy. Soc. Open Sci., 3, 160223

Kos, J., Zwitter, T., Grebel, E. K., et al. 2013, ApJ, 778, 86

Kos, J., Zwitter, T., Wyse, R., et al. 2014, Science, 345, 791

Krasnoukhov, V. S., Pivovarov, P. S., Zagidullin, M. V., et al. 2022, Astron. Rep.,
66, 811

Krelowski, J., & Walker, G. A. H. 1987, ApJ, 312, 860

Krelowski, J., & Westerlund, B. E. 1988, A&A, 190, 339

Krelowski, J., Snow, T. P., Seab, C. G., & Papaj, J. 1992, MNRAS, 258, 693

Kwok, S. 2023, Galaxies, 11, 104

Lan, T.-W., Ménard, B., & Zhu, G. 2015, MNRAS, 452, 3629

Li, K., Li, A., & Xiang, F. Y. 2019, MNRAS, 489, 708

Maŕz Apellániz, J., Barbá, R. H., Caballero, J. A., Bohlin, R. C., & Fariña, C.
2021, MNRAS, 501, 2487

McGuire, B. A., Loomis, R. A., Burkhardt, A. M., et al. 2021, Science, 371,
1265

Munari, U., Tomasella, L., Fiorucci, M., et al. 2008, A&A, 488, 969

Planck Collaboration Int. XLVIII. 2016, A&A, 596, A109

Puspitarini, L., Lallement, R., Babusiaux, C., et al. 2015, A&A, 573, A35

Ramírez-Tannus, M. C., Cox, N. L. J., Kaper, L., & de Koter, A. 2018, A&A,
620, A52

Sarre, P. J., Miles, J. R., Kerr, T. H., et al. 1995, MNRAS, 277, L41

Schlafly, E. F., Meisner, A. M., Stutz, A. M., et al. 2016, ApJ, 821, 78

Schlafly, E. F., Peek, J. E. G., Finkbeiner, D. P., & Green, G. M. 2017, ApJ, 838,
36

Schlarmann, L., Foing, B., Cami, J., & Fan, H. 2021, A&A, 656, L17

Stockett, M. H., Bull, J. N., Cederquist, H., et al. 2023, Nat. Commun., 14, 395

Vergely, J. L., Lallement, R., & Cox, N. L. J. 2022, A&A, 664, A174

Vos, D. A. I., Cox, N. L. J., Kaper, L., Spaans, M., & Ehrenfreund, P. 2011, A&A,
533, A129

Xiang, F. Y., Li, A., & Zhong, J. X. 2011, ApJ, 733, 91

Xiang, F. Y., Li, A., & Zhong, J. X. 2017, ApJ, 835, 107

Zasowski, G., Ménard, B., Bizyaev, D., et al. 2015, ApJ, 798, 35

Zelko, I. A., & Finkbeiner, D. P. 2020, ApJ, 904, 38

Zelko, I. A., & Finkbeiner, D. P. 2021, ApJ, 914, 68

Zhao, H., Schultheis, M., Rojas-Arriagada, A., et al. 2021, A&A, 654, A116

Zhao, H., Schultheis, M., Arentsen, A., et al. 2023, MNRAS, 519, 754