Free Access
Volume 628, August 2019
Article Number A97
Number of page(s) 5
Section Interstellar and circumstellar matter
Published online 13 August 2019

© ESO 2019

1 Introduction

Mid-infrared spectroscopy of young stellar objects provides important clues about their circumstellar envelopes and disks through the appearance of emission and absorption features that make it possible to study their mineralogy. The intensity and structure of those features reveal a wealth of information about aspects such as the composition of the dust, its growth, the location where its composing mineral species are formed, and their thermal history. The study of those features provides an extremely valuable link between the processes of dust formation in the winds of red giant stars, its presence and conditions in the surroundings of young stellar objects, the formation of planetary systems around other stars, and also about our own Solar System where comets are carriers of pristine material (see Henning 2010 for a review).

The source that powers the Herbig-Haro object HH 250, known as HH 250-IRS (=IRAS 19190+1048), is a binary system composed of two Class I young stellar objects separated by 053. The system lies in the outskirts of the Aquila Rift (Prato et al. 2008) and is embedded in the LDN 643 dark cloud. An investigation of the young stellar population in the Serpens-Aquila region lying in the same general direction of the sky using the Gaia Data Release 2 (DR2) catalog (Gaia Collaboration 2018) has revealed an aggregate of at least 15 young stars associated with LDN 643, including the previously identified binary T Tauri stars AS 353 and V536 Aql (Herczeg et al. 2019). Astrometry of the members of the aggregate yields a distance of 407 ± 16 pc, which is very similar to that of other aggregates in the Aquila Rift and the Serpens star forming clouds. Adopting this distance for HH 250-IRS translates into a projected separation of 216 AU1 between its components. Observations in the visible, infrared and submillimeter reported in Comerón et al. (2018, hereafter CRYC18) show that each of the binary components is surrounded by large amounts of circumstellar material and that the binary is itself surrounded by a marginally resolved circumbinary disk. HH 250 and HH 250-IRS are thus some of the rare examples of young stellar objects simultaneously displaying features of early stellar evolution such as an active outflow, a circumbinary disk, and two circumstellar disks that can be studied separately thanks to subarcsecond resolution instrumentation currently available at near-infrared and submillimeter wavelengths.

Despite similarities in the shape of the largely featureless near-infrared spectra of both sources presented by CRYC18, some clear differences appear at longer wavelengths. One of the components, HH 250-IRS SE, has overall bluer colors than the NW component, a marked drop in flux at 3.8 μm, and hints of an emission component in the N band possibly due to a silicate feature in emission. The NW component has a redder and smoother spectral energy distribution, in this case with hints of a silicate feature in absorption. Despite their similarities, small separation and coevality, the spectral energy distribution traced by the broadband photometry presented in CRYC18 suggests substantial differences between the circumstellar envelopes of both objects.

In this article we report the result of spectroscopic observations in the 10 μm atmosphericwindow exploring the silicate features of both sources. Besides confirming preliminary findings by CRYC18, we further discuss evidence for dust processing in the SE component, where evidence of crystallinity is found.

2 Observations

The observations reported here were obtained with VISIR (Lagage et al. 2004), the thermal infrared imager and spectrograph at ESO’s Very Large Telescope (VLT), in spectroscopic mode. A low-resolution grating was used covering the 8 to 13 μm interval at a resolution R = 270 with the 04 slit used. The slit was oriented so as to include both components of the binary in every exposure. A combination of chopping and nodding along the slit direction with an amplitude of 8 was used to cancel out the sky emission as well as the signature due to differences in the optical path when chopping.

The observations were obtained on the night of 22 to 23 June 2018. A total of 18 nodding cycles, each with 56 chopping cycles of the telescope secondary mirror at a frequency of 0.8 Hz, were performed. The total exposure time was 2140 s. However, approximately 20% of the set of exposures taken at given nod positions were rejected due to transient image quality degradations that caused overlap between the traces of the spectra of both components of the binary system. The K5III star HD 189695 was observed at a very similar airmass right after the exposures on our science targets to allow for the removal of telluric features.

Standard IRAF tasks were used for data reduction and calibration. Spectra at each of the two positions along the slit were extracted and stacked together for each target source and for the telluric correction star. The trace of the telluric correction spectrawas used to extract a spectrum of the sky at each of the two slit positions, which was used for wavelength calibration. The stacked, wavelength-calibrated spectra of the target object at each of the two slit positions were then combined, and the same was done for the telluric correction star. Telluric features removal was carried out by dividing the spectra of the target stars by that of the telluric correction star, and the result was multiplied by a template spectrum of the K3III star α Hya (Cohen et al. 1999) in order to get a relative flux calibration. The resulting spectra are shown in Fig. 1, where differences between both sources are visible.

thumbnail Fig. 1

Spectra of both components of HH 250-IRS binary in 8–13 μm region. The spectra are smoothed to the resolution λ∕Δλ ≃ 170 of the flux-calibrated spectrum of the telluric correction star (Cohen et al. 1999).

3 Results

The advent of mid-infrared spectrographs onboard space observatories such as ISO, Spitzer, and Herschel, and of sensitive mid-infrared spectrographs at large ground-based telescopes, has produced large collections of spectra with silicate features around young stellar objects of various masses, ranging from Herbig Ae/Be stars (van Boekel et al. 2005; Sicilia-Aguilar et al. 2011) down to lower mass young stellar objects (Kessler-Silacci et al. 2005), and all the way to the substellar limit (Merín et al. 2007), illustrating a rich diversity that can be interpreted in terms of mixtures of different species, grain sizes, locations, and processes of crystal formation. Such diversity is found even in our own planetary system, as illustrated by the silicate features observed in cometary comae (Wooden et al. 2007).

3.1 HH 250-IRS NW and the contribution of interstellar extinction

The 8 13 μm spectrum of HH 250-IRS NW (Fig. 1, bottom panel) shows a smooth, structureless feature in absorption with a minimum at 9.4 μm. The feature has the usual shape produced by small (~ 0.1 μm) grains of amorphous olivine and pyroxene silicate. The feature may be due either to absorption by unprocessed dust in the circumstellar envelope, or to interstellar grains in the intervening dust column. Images at visible wavelengths of the region where both sources are embedded (see Fig. 1 of CRYC18) clearly show that there is a considerable amount of extinction toward them, which should be accounted for at the time of discussing the intrinsic properties of the circumstellar environments of the HH 250-IRS binary components.

In Fig. 4 of CRYC18 we presented the L′-band spectrum toward both sources, which shows a prominent, broad absorption feature centered around 3.1 μm due to interstellar water ice grains. By a rough extrapolation of the continuum longward of the feature we estimate an optical depth at the core of the feature, which corresponds to a visible extinction AV ≃ 11.5 mag (Whittet et al. 2001). On the other hand, we estimate τSi ≃ 0.85 at the core of the silicate feature in the spectrum of HH 250-IRS NW by interpolating between the fluxes at 8 and 13 μm. This is most likely an underestimation, since the edges of the feature extend somewhat beyond both limits and the interpolated continuum should actually be higher than the level that we have assumed. The contribution of the foreground interstellar dust to the feature is estimated to be τSi, ISMAV ∕ 16.6 ≃ 0.7 using the AV -to-τSi, ISM proportionality of Rieke & Lebofsky (1985) and the value of AV derived from the depth of the 3.1 μm ice feature. Interstellar dust absorption thus seems to be insufficient to account for the full depth of the silicate feature in the spectrum of HH 250-IRS NW, leading us to conclude that the intrinsic spectrum of its envelope contains the signature of amorphous silicate absorption.

3.2 HH 250-IRS SE

No obvious feature of the same shape as the one dominating the 8 13 μm spectrum of HH 250-IRS NW appears in the spectrum of the SE component (Fig. 1, upper panel). However, being just 053 from the NW component and therefore equally embedded in the molecular cloud, the amount of foreground extinction is expected to be essentially identical. The most likely cause of its absence in the spectrum of HH 250-IRS SE is the existence of emission by the same species in its circumstellar envelope, with an intensity that approximately fills in the absorption caused by foreground dust. We will assume henceforth this to be the case.

The most obvious characteristic of the spectrum of HH 250-IRS SE is the flux increase that dominates the longer wavelength part of the interval shown in Fig. 1. Such a structured feature indicates the presence of emission due to crystalline dust species. The flux increase starts around 9.5 μm and reaches a plateau extending up to 11.2 μm. Figure 2 compares the structured emission feature with the mass-absorption coefficients of the two most common crystalline silicates in circumstellar disks, enstatite (MgSiO3) and forsterite (Mg2 SiO4), for two typical grain sizes of each species (van Boekel et al. 2005).

The main characteristic of the feature, namely the rise in flux between ~ 9.5 and ~ 10.2 μm, can be best reproduced with crystalline grains of forsterite as a main component. No emission features attributable to enstatite are obvious shortwards of 9.5 μm, whose most prominent feature appears near 9.4 μm (Koike et al. 2000; Molster et al. 2002; Olofsson et al. 2009), indicating that forsterite is the dominant form of crystalline silicates in the circumstellar enviroment of HH 250-IRS SE. As can be clearly seen in the comparison spectra of Fig. 2, the bands that are characteristic of the spectrum of forsterite become decreasingly prominent with increasing grain size, and dilute into a continuum above grain sizes ~ 1.6 μm (Lindsay et al. 2013). The fact that the signature of forsterite is still visible in our spectrum, although not in the form of well-defined peaks, indicates the dominance of rather large grains and that significant grain growth through coagulation has taken place, up to sizes near ~ 1 μm. On the other hand, Lindsay et al. (2013) have investigated the dependency of the features on the shapes of the forsterite crystals, which in turn is related to their temperature of formation (Kobatake et al. 2008). In the shape classification of Lindsay et al. (2013), both a-platelets and c-columns have the dominant 10 μm feature shifted to a peak at 10.2 μm, like what we see in the spectrum of HH 250-IRS SE. C-columns have the 11 μm feature shifted to wavelengths shorter than 11 μm, which we do not see in our spectra. We therefore interpret the crystalline features that we observe in HH 250-IRS SE as being predominantly due to micron-size, perhaps a-platelet-shaped grains of forsterite, although the latter conclusion is rather tentative and should be confirmed with higher signal-to-noise spectroscopy.

thumbnail Fig. 2

Silicate feature of HH 250-IRS SE (middle panel) compared to spectra of enstatite (top panel) and forsterite (bottom panel). The spectra of crystalline silicates are adapted from van Boekel et al. (2005). To facilitate the comparison a constant level has been subtracted from the spectrum of HH 250-IRS SE. The solid and dotted curves in the spectra of both silicates correspond to grains of 1 and 0.1 μm, in size, respectively.

4 Discussion

All the characteristics of HH 250-IRS NW noted in CRYC18 and here tend to confirm its classification as a relatively unevolved Class I object, surrounded by a thick envelope where dust is predominantly in the amorphous state. Since evidence for significant dust grain growth exists already during the relatively short-lived Class I phase (Sheehan & Eisner 2017), it is likely that it has already taken place in HH 250-IRS NW to some degree, but without producing yet detectable amounts of crystallinity, which is a slower process than coagulation (Suttner & Yorke 2001). An alternatively possibility might be that crystalline evolution has taken place in the inner disk of HH 250-IRS NW, but that the disk is seen close to edge-on from our vantage point, resulting in an optically thick outer layer of amorphous dust blocking our view of the central regions. This seems nevertheless unlikely, as computed model-spectral energy distributions predict a much deeper flux drop, by several orders of magnitude, in the 10 μm region (Robitaille et al. 2006) for objects in early evolutionary stages, even when the geometry departs from an edge-on orientation by a few tens of degrees.

Given the close proximity between the two components of HH 250-IRS, it is possible that the overall evolution of their circumstellar envelopes has been accelerated by the presence of the companion. Using submillimeter observations of a sample of Class 0 and Class I binaries in Ophiuchus with separations in the range 450–1100 AU, Patience et al. (2008) have shown that the frequency of circumstellar envelopes around the secondary component of the system drops abruptly from Class 0 to Class I. Their submillimeter observations are mainly sensitive to the evolution of the outer parts of envelopes, but demonstrate that the presence of a companion causes substantial evolution of the envelope even on timescales of ~ 105 yr or shorter. Such influence can be expected to be even more dramatic in the case of HH 250-IRS, in which the separation between components is between 2 and 5 times smaller than those probed by those authors. In agreement with the results of Patience et al. (2008), CRYC18 find only one compact source of dust continuum emission in the HH 250-IRS binary. The position of that source seems to be closer to HH 250-IRS SE, although this should be confirmed by more astrometrically accurate observations.

The presence of forsterite as the dominant crystalline species is a remarkable characteristic of the spectrum of HH 250-IRS SE, as the features of both species are most often observed simultaneously in the spectra of circumstellar environments. Forsterite formed by thermal annealing is the main form of silicate crystals in the innermost regions of circumstellar disks (van Boekel et al. 2005), reproducing model expectations (Gail 2004). Enstatite dominates the inner warm regions of the disk, but further outwards forsterite is found to be again the dominant form of crystalline silicates (Bouwman et al. 2008), in contrast with model predictions, perhaps indicating more complicated heating processes not considered in stationary disk models (Henning 2010), or as a consequence of the slow forsterite to enstatite conversion at low temperatures (Gail 2004).

We tentatively link the dearth of enstatite with the lack of dust in the inner disk that CRYC18 proposed as an explanation for the drop in L′-band flux of the spectral energy distribution of HH 250-IRS SE. The dust-devoid region may encompass the range of distances to the central star where enstatite would otherwise be the predominant crystalline silicate. The clearing of dust might be triggered bytidal interaction with HH 250-IRS NW. However, simulations of the interaction between circumstellar disks and binary companions at distances of several tens or hundreds of AU, which result in truncation of the outer disk and tidal tails, do not seem to be able to reproduce such a feature (Artymowicz & Lubow 1994). This has been observationally supported by Pascucci et al. (2008), who found no statistical difference in the degree of crystallinity and dust settling between a sample of single T Tauri stars, and another composed of binaries with separations between 10 and 450 AU, which suggests that the presence of companions in that range does not significantly affect the inner disk. An intriguing alternative possibility is that HH 250-IRS SE may itself be a tight binarysystem, in which the tidal interaction with the central binary is responsible for the inner truncation of the disk common to both components. Andrews et al. (2010) have studied two such cases in the TW Hya association, and Ruíz-Rodríguez et al. (2016) have shown that over one-third of a sample of young stellar objects in nearby star forming regions with spectral energy distributions typical of transitional disks, and lacking near-infrared excesses, are actually close binaries. If this were the case of HH 250-IRS SE as well, then HH 250-IRS would be a hierarchical triple system with two nested circumbinary disks, plus one circumstellar disk around the NW component. This is the expected outcome of the dynamical decay of a non-hierarchical triple system initially embedded in the dense interior of a cloud core (e.g. Reipurth et al. 2010; Reipurth & Mikkola 2015). Such an explanation is nevertheless speculative at this point in the absence of any other evidence of close binarity of HH 250-IRS SE.

The match between the structure of the crystalline silicate feature and the spectral characteristics of a-platelets also provides information on the temperature of formation. Kobatake et al. (2008) show that a-platelets are preferentially formed at temperatures between 970 and 1270 K, with more symmetric structures (equants in the terminology of Lindsay et al. 2013) formed at the highest temperatures and columnar structures formed below. As noted by Lindsay et al. (2013), the spectral signatures of platelet shapes are not detected in comets, implying that cometary dust in our Solar System condensed at higher temperature.

In view of the virtually certain coevality of the two components of the binary, the results obtained in CRYC18 and enlarged here by means of spectroscopy of the silicate feature caution against an interpretation of crystallinity as a straightforward indicator of the evolutionary stage of a circumstellar disk. SVS 20, another Class I binary discovered by Eiroa & Leinert (1987), has some features in common with HH 250-IRS, displaying similar spectral energy distributions with different mineralogy in the 10 μm band (Ciardi et al. 2005). However, the spectral energy distributions of both components of SVS 20 are much more similar to each other, and both components exhibit the crystalline emission features that are lacking in HH 250-IRS NW.

The strikingly different characteristics of the dust around the components of HH 250-IRS might be due to variations in the evolutionary timescales of circumstellar envelopes caused by the presence of companions, which in the case of HH 250-IRS-SE may either be the NW component or a hypothetical close companion responsible for producing a large gap in the inner disk. However, other causes are possible. Ábrahám et al. (2009) have found evidence for fast formation of crystalline features during the 2008 outburst of EX Lup, with strong indirect evidence that such features could disappear from the spectrum in a matter of decades, most likely due to radial transport of the freshly-formed crystals (Juhász et al. 2012). The existence of the Herbig-Haro object HH 250, with a kinematical age estimate of ~ 3500 yr, shows that at least one of the members of the HH 250-IRS system is at the jet-launching phase that usually coincides with the evolutionary stage at which FU Ori or EX Lup outbursts take place. Accretion variability, rather than steady disk evolution, could therefore be responsible for the derived crystallinity of HH 250-IRS-SE.

Since the depth of the 10 μm feature is a sensitive function of the disk orientation along the line of sight, it cannot be ruled out that non-coplanarity of the disks of both components also plays a role in the different appearance of their spectra. However, the moderate depth of the 10 μm feature in both objects indicates that neither of them is seen close to edge-on, as discussed above. We thus deem very unlikely that the possibly different orientations of their disks is a dominant factor causing the sharp differences between their spectra.

Our results on the silicate features study adds a new aspect to the wealth of structures and phenomena related to the earliest stages of stellar evolution that are found in the HH 250-IRS system, and presents it as a case study where the effects of multiplicity and accretion variability on the evolution of the mineralogy of the circumstellar environment may be investigated. Further observations, and particularly adaptive optics-assisted photometric and spectroscopic monitoring able to separately follow the evolution of each component, will be helpful in this regard.


We are grateful to Christian Hummel, from ESO’s User Support Department, for the careful verification of the Observation Blocks composing our program, and to the ESO staff astronomers in charge of Service Mode observations at the VLT when our program was executed. Comments by the anonymous referee greatly helped with the discussion of our results and are gratefully acknowledged.


  1. Ábrahám, P., Juhász, A., Dullemond, C. P., et al. 2009, Nature, 459, 224 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  2. Andrews, S. M., Czekala, I., Wilner, D. J., et al. 2010, ApJ, 710, 462 [NASA ADS] [CrossRef] [Google Scholar]
  3. Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bouwman, J., Henning, T., Hillenbrand, L. A., et al. 2008, ApJ, 683, 479 [NASA ADS] [CrossRef] [Google Scholar]
  5. Ciardi, D. R., Telesco, C. M., Packham, C., et al. 2005, ApJ, 629, 897 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cohen, M., Walker, R. G., Carter, B., et al. 1999, AJ, 117, 1864 [NASA ADS] [CrossRef] [Google Scholar]
  7. Comerón, F., Reipurth, B., Yen, H.-W., & Connelley, M. S. 2018, A&A, 612, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Eiroa, C., & Leinert, C. 1987, A&A, 188, 46 [NASA ADS] [Google Scholar]
  9. Gaia Collaboration (Brown, A. G. A., et al. ) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Gail, H.-P. 2004, A&A, 413, 571 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Henning, T. 2010, ARA&A, 48, 21 [NASA ADS] [CrossRef] [Google Scholar]
  12. Herczeg, G. J., Kuhn, M. A., Zhou, X., et al. 2019, ApJ, 878, 111 [NASA ADS] [CrossRef] [Google Scholar]
  13. Juhász, A., Dullemond, C. P., van Boekel, R., et al. 2012, ApJ, 744, 118 [NASA ADS] [CrossRef] [Google Scholar]
  14. Kessler-Silacci, J. E., Hillenbrand, L. A., Blake, G. A., & Meyer, M. R. 2005, ApJ, 622, 404 [Google Scholar]
  15. Kobatake, H., Tsukamoto, K., Nozawa, J., et al. 2008, Icarus, 198, 208 [NASA ADS] [CrossRef] [Google Scholar]
  16. Koike, C., Tsuchiyama, A., Shibai, H., et al. 2000, A&A, 363, 1115 [NASA ADS] [Google Scholar]
  17. Lagage, P. O., Pel, J. W., Authier, M., et al. 2004, The Messenger, 117, 12 [NASA ADS] [Google Scholar]
  18. Lindsay, S. S., Wooden, D. H., Harker, D. E., et al. 2013, ApJ, 766, 54 [NASA ADS] [CrossRef] [Google Scholar]
  19. Merín, B., Augereau, J.-C., van Dishoeck, E. F., et al. 2007, ApJ, 661, 361 [NASA ADS] [CrossRef] [Google Scholar]
  20. Molster, F. J., Waters, L. B. F. M., & Tielens, A. G. G. M. 2002, A&A, 382, 222 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Olofsson, J., Augereau, J.-C., van Dishoeck, E. F., et al. 2009, A&A, 507, 327 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Pascucci, I., Apai, D., Hardegree-Ullman, E. E., et al. 2008, ApJ, 673, 477 [NASA ADS] [CrossRef] [Google Scholar]
  23. Patience, J., Akeson, R. L., & Jensen, E. L. N. 2008, ApJ, 677, 616 [NASA ADS] [CrossRef] [Google Scholar]
  24. Prato, L., Rice, E. L., & Dame, T. M. 2008, Handbook of Star Forming Regions, ed. B. Reipurth (San Francisco: ASP), I, 18 [NASA ADS] [Google Scholar]
  25. Reipurth, B., & Mikkola, S. 2015, AJ, 149, 145 [NASA ADS] [CrossRef] [Google Scholar]
  26. Reipurth, B., Mikkola, S., Connelley, M., & Valtonen, M. 2010, ApJ, 725, L56 [NASA ADS] [CrossRef] [Google Scholar]
  27. Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618 [NASA ADS] [CrossRef] [Google Scholar]
  28. Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 256 [NASA ADS] [CrossRef] [Google Scholar]
  29. Ruíz-Rodríguez, D., Ireland, M., Cieza, L., & Kraus, A. 2016, MNRAS, 463, 3829 [NASA ADS] [CrossRef] [Google Scholar]
  30. Sheehan, P. D., & Eisner, J. A. 2017, ApJ, 851, 45 [NASA ADS] [CrossRef] [Google Scholar]
  31. Sicilia-Aguilar, A., Henning, T., Dullemond, C. P., et al. 2011, ApJ, 742, 39 [NASA ADS] [CrossRef] [Google Scholar]
  32. Suttner, G., & Yorke, H. W. 2001, ApJ, 551, 461 [CrossRef] [Google Scholar]
  33. van Boekel, R., Min, M., Waters, L. B. F. M., et al. 2005, A&A, 437, 189 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Whittet, D. C. B., Gerakines, P. A., Hough, J. H., & Shenoy, S. S. 2001, ApJ, 547, 872 [NASA ADS] [CrossRef] [Google Scholar]
  35. Wooden, D., Desch, S., Harker, D., Gail, H.-P., & Keller, L. 2007, Protostars and Planets V, eds. B. Reipurth, D. Jewitt, & K. Keil (Tucson, AZ: University of Arizona Press), 815 [Google Scholar]


The Gaia DR2 distance is almost twice the value adopted in Comerón et al. (2018), before the Gaia DR2 data became available. Some of the quantities discussed in that article need to be rescaled using the new distance, but the results obtained in the present work are independent of the distance.

All Figures

thumbnail Fig. 1

Spectra of both components of HH 250-IRS binary in 8–13 μm region. The spectra are smoothed to the resolution λ∕Δλ ≃ 170 of the flux-calibrated spectrum of the telluric correction star (Cohen et al. 1999).

In the text
thumbnail Fig. 2

Silicate feature of HH 250-IRS SE (middle panel) compared to spectra of enstatite (top panel) and forsterite (bottom panel). The spectra of crystalline silicates are adapted from van Boekel et al. (2005). To facilitate the comparison a constant level has been subtracted from the spectrum of HH 250-IRS SE. The solid and dotted curves in the spectra of both silicates correspond to grains of 1 and 0.1 μm, in size, respectively.

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.