Open Access
Issue
A&A
Volume 671, March 2023
Article Number A2
Number of page(s) 30
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202244702
Published online 27 February 2023

© The Authors 2023

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

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

1 Introduction

On the 6 December 2020, the Hayabusa2 mission of the Japan Aerospace Exploration Agency (JAXA) returned samples collected on the Cb-type dark asteroid Ryugu to Earth (Tsuda et al. 2020). For the first time, Hayabusa2 brought back samples from the surface of a known carbonaceous asteroid. These Solar System objects preserved ancient material that escaped planetary accretion and differentiation, providing pristine mineral and organic phases formed within the first million years after the proto-Sun gravitational collapse. The amount of material collected by Hayabusa2 (5.4 g) was larger than the 100 mg that had been expected (Tachibana et al. 2014; Sawada et al. 2017). The infrared (IR) analyses of these samples provide unprecedented information on the formation and early evolution of the Solar System and allow us to make the link between primitive interplanetary material and astronomical observations. A key goal of this study between the Organic Macromolecule (Yabuta et al. 2022) and Stone (Nakamura et al. 2023) sub-teams in the Hayabusa2 Initial Analysis is to elucidate the distributions and chemical characteristics of macro-molecular organic materials and minerals in a C-type asteroid. In this work, we provide hyper-spectral analyses showing the characteristic content and distribution of both organics and minerals of Hayabusa2-returned particles that will enable the comparison with other types of asteroidal and cometary samples available in the laboratory.

Table 1

Ryugu samples analysed.

2 Experiments

2.1 Sample preparation

The Ryugu samples originated from chamber A and chamber C, respectively collected for the former at the surface of the first touchdown site and potentially containing sub-surface material ejected from the Hayabusa2-made artificial crater for the later (Tachibana et al. 2022). The studied samples include fragments of individual particles from A0064, C0002, C0040, C0046, and C0057 as well as aggregate samples from A0106-13, A0108-15, A0108-19, C0109-4, and C0109-10 (see Table 1). In total six fragments from chamber A and six from chamber C were prepared for Fourier Transform IR micro-spectroscopy (μ-FTIR), in particular for pre-characterisation before IR nanospectroscopy (AFM-IR Mathurin et al. 2022a,b; Dazzi & Prater 2017). Each fragment was micro-manipulated and transferred into a diamond compression cell made of two facing diamond windows. The fragments were flattened to a thickness of ~1 μm in order to provide an optimal thickness for IR spectroscopic transmission measurements. After compression in the diamond cell, the flattened sample was distributed generally in each diamond window, providing two complementary flat samples, named 'top' and 'bottom', depending on the diamond window in which the sample finally resided. In some cases the majority of the sample remained in one diamond window only.

2.2 Synchrotron MIR and FIR measurements

The transmission IR measurements were performed on the SMIS beamline at the synchrotron SOLEIL during three runs in July and October 2021, and in April 2022. Table 1 presents the names of the samples analysed both in the mid-IR (MIR, 4000700 cm−1) and the far-IR (FIR, 700100 cm−1), their approximate size once flattened, and the number of pixels in the recorded hyper-spectral IR maps. The synchrotron beam was coupled to a Nicolet continuum IR microscope, flushed under dry nitrogen at room temperature. For the measurements presented in this analysis, the IR spot size was optimised close to the diffraction limit, with a 6 μm × 6 μm sampling window, and maps were recorded using orientations adapted to the geometry of each Ryugu fragment. The map sampling step was 3 μm (Nyquist sampling). In order to better constrain the mineral components of the Ryugu samples, additional FIR spectra were recorded on a Nicolet Nic-Plan IR microscope using a bolometer, with spot sizes usually covering the entire size of the considered particles. Additional measurements over the entire size of the considered particles were also conducted on another Nicolet continuum IR microscope at ISMO, with a conventional Globar IR source covering an extended MIR range from 4000 cm−1 down to about 450 cm−1, ensuring the good behaviour in the overlapping 700−450 cm−1 spectral region. The FIR spectra were stiched to the MIR spectrum when necessary by applying a gain correction factor of typically a few percent, to build the global individual grain spectra covering the 4000−100 cm−1 range. The synchrotron IR beam μ-spectrometry was mandatory to sample the small components (e.g. small carbonate components) present in the samples, and to evaluate their spatial distribution and individual sizes, in order to build analytical chemical maps. The synchrotron hyper-spectral maps acquired on each Ryugu sample were used with the following complementary approaches: (i) the spectra corresponding to each map pixel were decomposed into individual contributions corresponding to specific vibrational modes (as shown in Figs. 1 and 2), to produce false-colour images of the spatial distribution of the different components (as shown, for instance, in Figs. 3 and 4); (ii) an average global spectrum (Fig. 5) was calculated for each map of the Ryugu samples combining the best signal-to-noise MIR (either synchrotron average map or Globar source) and FIR spectra, to estimate the CH2/CH3 ratios (Figs. 6-8) and the distribution of the different components in each map (Fig. 9), as well as the proportion of carbonates to phyllosilicates (Fig. 10).

thumbnail Fig. 1

Representative Ryugu sample A0108-15 transmittance spectra (6 × 6 μm spots) showing a carbonate-rich region (A), a region with inorganic 'opaques' (probably magnetite and/or Fe sulphides) (B), and the map average spectrum (green) corresponding to map 1 reported in Fig. 3. Vibrational modes of interest (including the ones used to build the coloured maps such as in Fig. 3) are shown with arrows.

2.3 AFM-IR measurements

An IR map at the tens of nanometres resolution was obtained using AFM-IR (shown in the lower left panel of Fig. 3), a technique combining an Atomic Force Microscope (AFM) to a tunable IR laser. It is based on the detection, by the AFM probe, of the photo-thermal expansion that occurs in the sample after IR absorption of the material located just under the AFM tip (Dazzi et al. 2005). The AFM-IR system was a NanoIR2 from Anasys Instruments (Bruker) coupled to a multi-chip Quantum Cascade Laser (QCL, Daylight Solutions) covering the MIR range from 900 cm−1 up to 1960 cm−1. An Au-coated silicon AFM tip was used (Budget sensors, tap 300). The AFM-IR was used in tapping mode (Mathurin et al. 2020, 2019). The IR map (AFM image of 30 μm × 30 μm, 1000 × 1000 pixels) was acquired at fixed tuned laser wavenumbers: 1020, 1420, and 1900 cm−1, corresponding to the absorption of silicate and carbonate chemical groups, and of the continuum, respectively. The IR mapping acquisitions were made using a 0.1 Hz scan rate with a laser power of 56.26% of the QCL chip, a repetition rate of 1340 kHz, and a pulse width of 120 ns. These maps were combined in a composite three-colour RGB image (green at 1020 cm−1, blue at 1420 cm−1, and red at 1900 cm−1), shown in the lower left panel of Fig. 3. Before combining the maps, minor AFM stripping was corrected by interpolation, the maps were realigned to the 1020 cm−1 map using their topographic AFM maps, and they were slightly smoothed using a 5-pixel kernel. To maximise the contrast in the composite image for bands with different intrinsic intensities, the intensity of each image recorded at a given wavenumber was first normalised to its maximum intensity over the entire map before the images were combined.

3 Analyses and results

Each spectrum in the synchrotron hyper-spectral maps was baseline corrected. The integrated optical depths over several absorption bands of interest were calculated, and divided by their corresponding integrated absorption cross section to estimate the column density measured in each pixel of the maps. The resulting maps with significant signal-to-noise ratios correspond to the silicate's SiO stretching mode, CO3 carbonates, OH stretching mode, and aliphatic CH. In addition, we built a map on phases that do not exhibit specific absorption bands but a wider light scatter and/or absorption behaviour that can be referred to as 'opaques'. For these opaque phases, instead of measuring a particular band, we monitored the non-baseline-corrected continuum drop by comparing the optical depth difference between the continuum on the right side of the silicates (i.e. at about 800 cm−1) and in an IR absorption-free region at about 2500 cm−1. These opaques maps betray the presence of species such as magnetite and/or sulphides, as these phases do not display specific signatures in the MIR range, but strongly absorb and/or scatter the IR beam.

The adopted integrated absorption cross sections (also called band strengths) are given in Table A.1. We detail below how they were determined from literature optical constants and/or spectra.

thumbnail Fig. 2

Close-up and deconvolution of the individual bands for the global spectra (corresponding to the maps of Figs. 3 and 4) recorded for sample A0108-15 (03-top-a), upper spectrum, with a high physisorbed and/or interlayer water contribution, and C0040-FC025 (bot-b), second spectrum from top, dominated by the OH hydroxyl contribution. Individual bands are represented by hatching: OH stretch (blue in the 3000–3700 cm−1 region), CH stretch (magenta), C=0 carbonyl stretch (brown), OH bending (blue around 1640 cm−1), aromatic C=C bending (black), carbonates (red), and silicates (green). Each wavenumber range is normalised with regard to the high wavenumber range by a factor (indicated in parenthesis) for clarity of the deconvolution. The spectrum of a particularly C-aromatic-rich grain, C0109-04 (top-a), is also shown. For comparison, a spectrum taken from Orgueil, which was prepared and recorded in the same conditions as the Ryugu samples, is shown at the bottom. See text for details.

thumbnail Fig. 3

Chamber A Ryugu sample #A0108-15 (03) (top diamond window). Upper left: Visible image recording the entire flattened grain with an extent of about 150 μm. The positions A and B correspond to the individual spectra shown in Fig. 1, showing a carbonate-rich region (A) and a region with inorganic opaques (probably magnetite and/or Fe sulphides) (B). The regions of interest mapped with the synchrotron measurements are numbered from 1 to 3. Top right and middle panel sets: Coloured phase maps for regions 1 to 3 corresponding, for each map, to silicates (a, green), carbonates (b, red), opaques, (c, yellow-orange), OH stretch (d, blue), CH stretch (e, magenta), and respective visible image (f). It should be noted that a few outlier (white) points with apparent high optical depth may be observed in the OH and CH stretch maps. They coincide with the presence of opaques or are on maps' borders. They are due to low signal-to-noise ratios, which lead to a higher uncertainty on the baseline determination (and correction). Lower panel: Images representing the same region mapped by AFM-IR (left, composite map of 30 μm × 30 μm, green, silicates; blue, carbonates; red, opaques), and mapped with synchrotron at the highest achievable spatial resolution (right,~6 μm aperture with a step size of ~3 μm, i.e. close to Nyquist sampling).

thumbnail Fig. 4

Chamber C Ryugu sample C0040-FC025 (bottom diamond window). Phase maps of silicates (a, green), carbonates (b, red), opaques (c, yellow-orange), OH stretch (d, blue), CH stretch (e, magenta), and corresponding visible images (f). Bottom: images representing the edge detection of individual ROIs based on the intensity map of carbonates (left), and of opaques (middle), and resulting size distribution histograms (right). The x-axis represents the projected surface of each region in μm2 (multiples of the spatial resolution of 9 μm2). As the samples are flattened to ~1 μm, this projected surface is directly proportional to the same volume value.

3.1 Band strengths

An integrated absorption cross section can be evaluated from the imaginary part of the optical constants for a given solid by integrating over a vibrational band: 𝒜(cm/[ group ])=v¯minv¯maxτ(v¯)dv¯N=v¯minv¯max4πkv¯×thickness×dv¯N=v¯minv¯max4πkv¯dv¯×MNA×ρ×ngroup,$\matrix{ {𝒜 \left( {{{{\rm{cm}}} \mathord{\left/ {\vphantom {{{\rm{cm}}} {\left[ {{\rm{group}}} \right]}}} \right. \kern-\nulldelimiterspace} {\left[ {{\rm{group}}} \right]}}} \right)} \hfill & { = \int_{{{\bar v}_{\min }}}^{{{\bar v}_{max}}} {{{\tau \left( {\bar v} \right){\rm{d}}\bar v} \over N}} } \hfill & { = \int_{{{\bar v}_{\min }}}^{{{\bar v}_{\max }}} {{{4\pi k\bar v \times {\rm{thickness}} \times {\rm{d}}\bar v} \over N}} } \hfill \cr {} \hfill & {} \hfill & {\, = \int_{{{\bar v}_{\min }}}^{{{\bar v}_{\max }}} {{{4\pi k\bar v{\rm{d}}\bar v \times M} \over {{N_{{\rm{A}} \times \rho \times {n_{group}}}}}},} } \hfill \cr }$(1)

where τ(v¯)$\tau \left( {\bar v} \right)$ is the optical depth of the vibrational band related to the group, N the column density of the considered group, M [g mol−1] is the molar mass, NA the Avogadro number, ρ(g cm−3 the volumetric mass density, ngroup the number of groups considered in the calculated molar mass to which 𝒜 is normalised (e.g. for dolomite, if M is calculated with the formula CaMg(CO3)2, and A is expressed per CO3 group, then nCO3=2${n_{{\rm{C}}{{\rm{O}}_{\rm{3}}}}} = 2$).

thumbnail Fig. 5

Global MIR to FIR spectra of the grains analysed, including large regions of interest for some of the grains from chamber A (left) and chamber C (right). Reference spectra for major species identified in the Ryugu samples spectra such as carbonates: dolomite, Roush (2021); magnesite, Roush (2021); breunnerite Mg0.78Fe0.22CO3, (this work) and phyllosilicates: saponite, serpentine, Glotch et al. (2007), are displayed for comparison (mineral compositions are given in Table A.1). For some grains with a large carbonate contribution to the spectrum, we also display the spectrum of dolomite, the carbonate the closest to the measured band positions and intensities. In addition, in the lower right panel are displayed the spectra of two CI1 meteorites closely related to the composition of the Ryugu samples, Orgueil (this work, measured under the same conditions in a diamond cell) and Ivuna (measured in a KBr pellet at 300°C, from Beck, Pierre (2012): MIR absorbance spectra of bulk CI chondrites in KBr pellets at ambient temperature, 150°C and 300°C. SSHADE/GhoSST, OSUG Data Center, https://doi.org/10.26302/SSHADE/EXPERIMENT_LB_20170719_001).

OH band strength

In Ryugu phyllosilicates, in addition to the OH relatively sharp hydroxyl group at about 3700 cm−1, some of the recorded spectra show a broad band centred at about 3300 cm−1 for interlayer water, with some contribution arising from residual ambient water physisorption in some cases. For these last modes, in the stretching mode region, we adopted the band strength of liquid water (see Table A.1). The estimate of the OH stretching was used to estimate the contribution in the MIR region of the OH bending mode, overlapping with aromatic C=C modes. In liquid water, the OH bending to OH stretching mode band strength ratio is about 0.072. This ratio was estimated also using natural saponite spectra recently published in a review (Kloprogge & Ponce 2021), leading to 0.071 ± 0.0017.

thumbnail Fig. 6

Aliphatic stretches' spectral decomposition used to calculate the methylene-to-methyl ratios. The spectra are fitted with five components corresponding to the asymmetric CH3 stretching mode (~2960 cm−1), the asymmetric CH2 stretch and its red shoulder (~2920−2900 cm−1), and the symmetric stretches of the same groups ~2880-2850 cm−1). The approximate positions of these modes are labelled above the upper left spectrum. The CH3 contributions are shown in red, and the CH2 ones in blue. In a few spectra, a contribution from dolomite affects the profile. For these, the expected dolomite spectrum contribution is shown in brown and the dolomite-subtracted spectrum is shown above. Spectra with a high fraction of dolomite contribution were not deconvolved. In the upper right corner are shown reference spectra of linear and branched alkanes, an alkyl benzene and a cycloalkane, with varying methylene-to-methyl ratios: isocetane, squalane, n-tetradecane (purchased from Merck), and n-tetracosane (provided by Sylvie Heron) were measured as received with our FTIR-microscope; n-hexane and n-decane spectra are from Dartois et al. (2004); triethylbenzene and dimethylcyclohexane spectra are from the NIST database: https://webbook.nist.gov/chemistry/. Although they are specific species, this series shows the expected asymmetric CH2/CH3 stretch band ratio trend.

thumbnail Fig. 7

Methylene-to-methyl (CH2/CH3) ratios (chamber A: triangles; chamber C: circles) expressed as the asymmetric stretching mode peak maximum intensity (blue) and as column density ratio (red), derived from the deconvolution of the aliphatic bands, as a function of the CH-to-Si column density ratio. See text for details.

Aliphatic CH band strength

The individual contributions to the observed profile in the CH stretching mode region for the methyl and methylene groups' asymmetric and symmetric vibrations have been discussed in the context of interstellar organics, based on laboratory experiments (Dartois et al. 2004, 2007) and literature data referenced therein. In addition, as discussed hereafter, a global band strength per CH can be calculated and used when the observed band is not deconvolved into its individual specific vibrations. They are given in Table A.1

C=C band strength

The (poly)aromatic C=C band strength may vary depending on their chemical speciation. The closest material to use as a reference is probably the insoluble organic matter (IOM) from meteorites (e.g. Kebukawa et al. 2011; Orthous-Daunay 2013). Using IOM spectra from Fig. 2 of Orthous-Daunay (2013), for which the H/C ratio has been measured independently (Alexander et al. 2007), provides a way to determine such C=C band strength. The IOM carbon content is dominated by carbon involved in aliphatic, aromatic, and carbonyl functional groups. The measured column density of C is thus given by the sum of these individual contributions: NC=NCaliph.+NCaro.+NCcarbony1    =γNHaliph.+τC=Cdv¯𝒜C=CAro.+τC=Odv¯𝒜C=O,$\matrix{ {{N_{\rm{C}}} = N_{\rm{C}}^{{\rm{aliph}}{\rm{.}}}} \hfill & { + N_{\rm{C}}^{{\rm{aro}}{\rm{.}}}} \hfill & { + N_{\rm{C}}^{{\rm{carbony1}}}} \hfill \cr {} \hfill & {} \hfill & {\,\,\,\, = \gamma N_H^{aliph.} + {{\int {{\tau _{C = C}}d\bar v} } \over {{𝒜}_{C = C}^{Aro.}}} + {{\int {{\tau _{C = O}}d\bar v} } \over {{{𝒜}_{C = O}}}},} \hfill \cr }$(2)

where 𝒜C=CAro.${𝒜}_{{\rm{C = C}}}^{{\rm{Aro}}{\rm{.}}}$ is the integrated absorption cross section for aromatic carbons, per C; 𝒜C=O${{𝒜}_{{\rm{C = O}}}}$ for the carbonyl group, per C. γ is the factor linking the CH aliphatic hydrogen column density to the aliphatic carbon one. This factor is linked to the methylene-to-methyl ratio β=NCH2/NCH3$\beta = {{{N_{{\rm{C}}{{\rm{H}}_{\rm{2}}}}}} \mathord{\left/ {\vphantom {{{N_{{\rm{C}}{{\rm{H}}_{\rm{2}}}}}} {{N_{{\rm{C}}{{\rm{H}}_{\rm{3}}}}}}}} \right. \kern-\nulldelimiterspace} {{N_{{\rm{C}}{{\rm{H}}_{\rm{3}}}}}}}$, defining the mean number of hydrogen per C in the aliphatic stretching mode region. From β, one derives that γ = (β + 1)/(2β + 3). To first order, one can neglect the carbonyl C=O contribution to the carbon budget, as the measured C=O band integrated optical depth (of the same order as the C=C band in the IOM spectra) is only due to an intrinsically high 𝒜C=O compared to that of C=C, implying a much lower C contribution. The γ parameter typically varies only from 0.40 to 0.46 for NCH2/NCH3=1${{{N_{{\rm{C}}{{\rm{H}}_{\rm{2}}}}}} \mathord{\left/ {\vphantom {{{N_{{\rm{C}}{{\rm{H}}_{\rm{2}}}}}} {{N_{{\rm{C}}{{\rm{H}}_{\rm{3}}}}}}}} \right. \kern-\nulldelimiterspace} {{N_{{\rm{C}}{{\rm{H}}_{\rm{3}}}}}}} = 1$ to 5, and, in the absence of the exact knowledge of the β ratio, one can adopt γ = 0.43 ± 0.03 to cover the range of most observed IOM values.

The IOM spectra show that the organic hydrogen content is dominated by aliphatic CHs, with methyl (CH3) and methylene (CH2) bonds. The measured aliphatics hydrogen column density on a spectrum is then given by: NHaliph=3NCH3+2NCH2=3τCH3asymdv¯𝒜CH3asym+2τCH2asymdv¯𝒜CH2asym.$N_{\rm{H}}^{{\rm{aliph}}} = 3{N_{{\rm{C}}{{\rm{H}}_3}}} + 2{N_{{\rm{C}}{{\rm{H}}_2}}} = 3{{\int {\tau _{{\rm{C}}{{\rm{H}}_3}}^{{\rm{asym}}}{\rm{d}}\bar v} } \over {{𝒜}_{{\rm{C}}{{\rm{H}}_3}}^{{\rm{asym}}}}} + 2{{\int {\tau _{{\rm{C}}{{\rm{H}}_2}}^{{\rm{asym}}}{\rm{d}}\bar v} } \over {{𝒜}_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{asym}}}}}.$(3)

To a good approximation (Dartois et al. 2004, 2007), 𝒜CH3asym/3𝒜CH2asym/24×1018cm/H=𝒜CH${𝒜}_{{\rm{C}}{{\rm{H}}_3}}^{{\rm{asym}}}/3 \approx {𝒜}_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{asym}}}/2 \approx 4 \times {10^{ - 18}}{{{\rm{cm}}} \mathord{\left/ {\vphantom {{{\rm{cm}}} {{\rm{H = }}{_{{\rm{CH}}}}}}} \right. \kern-\nulldelimiterspace} {{\rm{H = }}{{𝒜}_{{\rm{CH}}}}}}$(4)

and 𝒜CH3sy​m/𝒜CH3asy​m𝒜CH2sy m/𝒜CH2asy m0.2,${{{𝒜}_{{\rm{C}}{{\rm{H}}_3}}^{{\rm{sy}}{\rm{m}}}} \mathord{\left/ {\vphantom {{{𝒜}_{{\rm{C}}{{\rm{H}}_3}}^{{\rm{sy}}{\rm{m}}}} {{𝒜}_{{\rm{C}}{{\rm{H}}_3}}^{{\rm{asy}}{\rm{m}}} \approx {{{𝒜}_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{sy}}\,{\rm{m}}}} \mathord{\left/ {\vphantom {{_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{sy}}\,{\rm{m}}}} {_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{asy}}\,{\rm{m}}} \approx 0.2}}} \right. \kern-\nulldelimiterspace} {_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{asy}}\,{\rm{m}}} \approx 0.2}},}}} \right. \kern-\nulldelimiterspace} {_{{\rm{C}}{{\rm{H}}_3}}^{{\rm{asy}}{\rm{m}}} \approx {{_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{sy}}\,{\rm{m}}}} \mathord{\left/ {\vphantom {{_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{sy}}\,{\rm{m}}}} {_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{asy}}\,{\rm{m}}} \approx 0.2}}} \right. \kern-\nulldelimiterspace} {_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{asy}}\,{\rm{m}}} \approx 0.2}},}}$(5)

so that NHaliph.v¯2800v¯3000τaliphdv¯1.2𝒜CH.${\rm{N}}_{\rm{H}}^{{\rm{aliph}}{\rm{.}}} \approx {{\int_{{{\bar v}_{2800}}}^{{{\bar v}_{3000}}} {{\tau _{{\rm{aliph}}}}{\rm{d}}\bar v} } \over {1.2{{𝒜}_{{\rm{CH}}}}}}.$(6)

Then, with C/HNC/NHaliph.${{\rm{C}} \mathord{\left/ {\vphantom {{\rm{C}} {\rm{H}}}} \right. \kern-\nulldelimiterspace} {\rm{H}}} \approx {{{N_{\rm{C}}}} \mathord{\left/ {\vphantom {{{N_{\rm{C}}}} {N_{\rm{H}}^{{\rm{aliph}}{\rm{.}}}}}} \right. \kern-\nulldelimiterspace} {N_{\rm{H}}^{{\rm{aliph}}{\rm{.}}}}}$, the integrated absorption strength of aromatic C can be derived as 𝒜C=CAro.τC=Cdv¯v¯2800v¯3000τaliphdv¯×1.2𝒜CH/(CHγ).${{{𝒜}_{{\rm{C = C}}}^{{\rm{Aro}}{\rm{.}}} \approx {{\int {{\tau _{{\rm{C = C}}}}{\rm{d}}\bar v} } \over {\int_{{{\bar v}_{2800}}}^{{{\bar v}_{3000}}} {{\tau _{{\rm{aliph}}}}{\rm{d}}\bar v} }} \times 1.2{{𝒜}_{{\rm{CH}}}}} \mathord{\left/ {\vphantom {{_{{\rm{C = C}}}^{{\rm{Aro}}{\rm{.}}} \approx {{\int {{\tau _{{\rm{C = C}}}}{\rm{d}}\bar v} } \over {\int_{{{\bar v}_{2800}}}^{{{\bar v}_{3000}}} {{\tau _{{\rm{aliph}}}}{\rm{d}}\bar v} }} \times 1.2{_{{\rm{CH}}}}} {\left( {{{\rm{C}} \over {\rm{H}}} - \gamma } \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{{\rm{C}} \over {\rm{H}}} - \gamma } \right)}}.$(7)

Reanalysing data from several IOM measurements from Fig. 8 in Orthous-Daunay (2013), which reports the correlation between the ratio of integrated absorption of the CH stretching modes over the aromatic C=C one versus the H/C ratio, v¯2800v¯3000τaliphdv¯/τC=Cdv¯H/C$\int_{{{\bar v}_{2800}}}^{{{\bar v}_{3000}}} {{{{\tau _{{\rm{aliph}}}}{\rm{d}}\bar v} \mathord{\left/ {\vphantom {{{\tau _{{\rm{aliph}}}}{\rm{d}}\bar v} {\int {{\tau _{{\rm{C = C}}}}{\rm{d}}\bar v \approx {{\rm{H}} \mathord{\left/ {\vphantom {{\rm{H}} {\rm{C}}}} \right. \kern-\nulldelimiterspace} {\rm{C}}}} }}} \right. \kern-\nulldelimiterspace} {\int {{\tau _{{\rm{C = C}}}}{\rm{d}}\bar v \approx {{\rm{H}} \mathord{\left/ {\vphantom {{\rm{H}} {\rm{C}}}} \right. \kern-\nulldelimiterspace} {\rm{C}}}} }}} $, over the H/C range of about 0.35 to 0.8. Putting this correlation into Eq. (7) and normalising per carbon atom (i.e. specifying that each C=C contributes for two carbon atoms), one deduces: 𝒜C=CAro.(IOM)3.2±0.4×1018cm/C.${{{𝒜}_{{\rm{C = C}}}^{{\rm{Aro}}{\rm{.}}}\left( {{\rm{IOM}}} \right) \approx 3.2 \pm 0.4 \times {{10}^{ - 18}}{\rm{cm}}} \mathord{\left/ {\vphantom {{_{{\rm{C = C}}}^{{\rm{Aro}}{\rm{.}}}\left( {{\rm{IOM}}} \right) \approx 3.2 \pm 0.4 \times {{10}^{ - 18}}{\rm{cm}}} {\rm{C}}}} \right. \kern-\nulldelimiterspace} {\rm{C}}}.$(8)

Carbonate CO3 band strength

The CO3 asymmetric stretches from dolomite, magnesite, and calcite carbonates were calculated based on the optical constants from Roush (2021). We assumed a generic band strength of 2×10−16cm/CO3.

Silicate SiO band strength

The silicates' SiO stretching mode band strengths for saponite and serpentine were obtained from the analyses of optical constants from Glotch et al. (2007); Roush et al. (1991); Mooney & Knacke (1985). We assumed a generic band strength of 1.4×10−16cm/Si.

thumbnail Fig. 8

Comparison of methylene-to-methyl (CH2/CH3) intensity peak ratios for Ryugu bulk sample measurements (this work) with measurements from the literature for diffuse ISM observations (Sandford et al. 1991; Pendleton et al. 1994; Dartois et al. 2007; Godard et al. 2012), Stardust samples (comet 81P/Wild2, Keller et al. 2006), IDPs (Flynn et al. 2003), and meteorites bulk and IOM samples (Alexander et al. 2014; Briani et al. 2013; Dionnet et al. 2018; Kebukawa et al. 2010, 2011, 2019; Kerraouch et al. 2022; Merouane et al. 2012; Orthous-Daunay 2011, 2013; Quirico et al. 2018; Yesiltas et al. 2014, 2015, 2022; Yesiltas & Kebukawa 2016).

thumbnail Fig. 9

Relative proportions of the different functional groups (aliphatics, hydroxyl + interlayer and/or physisorbed water, carbonyl, carbonates, and silicates) from the FTIR analysis of the different samples, shown per functional group CH, OH, C=O, CO3, and Si, respectively.

thumbnail Fig. 10

Relative volume ratio of carbonates to (carbonates + phyllosili-cates) in the Ryugu samples (in %), calculated using Eq. (13). Samples from chamber A are plotted as blue triangles and chamber C samples are represented as blue dots (see text for details). Two measurements on polished sections of chamber C samples are reported for comparison as filled red dots (Table S6 of Nakamura et al. 2023).

3.2 Phase maps

The phase maps computed from the integration of band vibrations in the synchrotron hyper-spectral maps for one sample from chamber A (A0108-15) and one from chamber C (C0040-FC025) are shown in Figs. 3 and 4, respectively. All the other measured chamber A and C samples' phase maps are shown in the appendix. These two maps exemplify the typical information obtained with the performed synchrotron measurements. In all recorded maps, both the phyllosilicates and aliphatics components of the organics are widespread. They are present everywhere at the scale probed with the synchrotron (6 μm), sampled at 3 μm steps, and intimately mixed at such a scale. The organic aromatic C=C component is ubiquitous in the Ryugu samples (as seen in Raman spectra, e.g. Bonal et al. 2022). However, the corresponding signal in the FTIR spectra overlaps with that of the OH bending mode from interlayer water, which is significant in the case of Ryugu spectra. Therefore, no maps for the aromatic C=C were extracted in this work, and this component is discussed in more detail in Sect. 3.5. Physisorbed and interlayer water, observed through the OH stretching modes, is also observed at varying degrees. Carbonates and opaques (magnetite and/or sulphides) show a different behaviour. Image processing analyses are shown in Fig. 4 for the opaques and carbonate maps. Individual regions of interest (ROIs) were automatically detected by thresholding the images and deriving individual non-connected projected surfaces. These 2D projected surfaces correspond to probing initial 3D volumes flattened to about one micron in thickness. A histogram of the ROI size distribution is shown in the right panel of Fig. 4. The opaques distribution for these grains show a size distribution following, to first order, the inverse of their volume. The carbonate distribution is different, with fewer small grains than for opaques, and individual large carbonate grains with projected surfaces as large as a few thousand μm2. These correspond to individual carbonate grains of tens of microns in diameter prior to flattening, as also observed by Yokoyama et al. (2023) and Nakamura et al. (2023).

The AFM-IR measurements extend the analysis to small scales around a few tens of nanometres, which are several orders of magnitude more resolving than the synchrotron spatial scale, which is limited by the classical diffraction limit of a few μm. This can be seen from the comparison of the same 30 μm×30 μm region mapped by AFM-IR and with synchrotron IR micro-spectroscopy at the highest achievable spatial resolution (Fig. 3, lower panels; Yabuta et al. 2022; Mathurin et al. 2022a). The AFM-IR images show that phyllosilicates are widespread down to the finest scale (i.e. a few tens of nanometres). By contrast, isolated carbonates are visible and mixed into the phyllosilicates. The AFM-IR technique reveals the sizes of these phases below the synchrotron IR micro-spectroscopy diffraction limit. As an example, the large carbonate (about 10 μm large) that one can guess in the top right of the synchrotron IR micro-spectroscopy image (Fig. 3, bottom right), is clearly resolved and visible in the AFM-IR image (Fig. 3, bottom left). By contrast, the opaque phase detected in the centre of the synchrotron IR micro-spectroscopy image is also visible in the AFM-IR image but with contours somehow blurred, indicating that this opaque phase is most probably buried below the surface of the particle.

3.3 Composition of the Ryugu samples from infrared spectra

To present a chemical composition of the Ryugu samples, we derived the global spectrum for each sample, combining the respective MIR and FIR spectra (Fig. 5). Each spectrum was shifted by an arbitrary amount for clarity. In the upper right panel of Fig. 5, the spectra of saponite and serpentine are plotted for comparison (see Table A.1 and Glotch et al. 2007). In the top right of the figure, the reference spectra of three carbonates are shown: dolomite and magnesite (CaMg(CO3)2) and (MgCO3), see Table A.1 and Roush 2021), and breunnerite (Mg0.78Fe0.22CO3, this work). Figure 11 displays a close-up of selected spectra in the 1900−100 cm−1 spectral region, with reference spectra of possible mineral components. When the fraction of carbonates is high enough, the set of bands observed points to the identification of dolomite, as shown by the vertical dashed blue lines, with some bands in the FIR, such as at 160 cm−1, which are very characteristic. Dolomite is the dominant carbonate found in this set of samples, at the scale of the synchrotron IR microscopy. On the phyllosilicate side, the spectra represent a main component of saponite mixed with serpentine that gives rise to a shoulder on the red wing of the SiO stretching mode absorption. In order to determine the respective contributions of the components present in each sample, each global spectrum was deconvolved into individual components, as shown in Fig. 2. The absorption bands of the so-called opaque minerals (magnetite and iron sulphides) are responsible for a significant part of the extinction of the IR continuum in some spectra, as shown in Fig. 1. Their contributions are hard to retrieve from the IR spectra in Fig. 11, because their vibrational mode absorption bands occurring in the FIR (e.g. around ~570 and ~360 cm−1 for magnetite, and around ~270 cm−1 for FeS, Glotch & Rossman 2009; Querry 1985; Mukai 1990; Henning & Mutschke 1997) are accompanied by a high refractive index, and the expected resulting low-contrast bands overlap with the intense bands of the other major mineral components present.

The OH-CH stretch region is decomposed in two global components. The CH stretching mode contribution was obtained by subtracting a local baseline continuum. The OH contribution was then measured by integrating over the CH subtracted spectrum. In the 1800–1280 cm−1 region, the spectrum was deconvolved using five Gaussians: one for the carbonyl (C=O) contribution peaking at about 1710 cm−1, one for an OH bending peaking around 1640−1620 cm−1, one for an aromatic C=C peaking around 1600 cm−1, and two for bands peaking at 1500−1420 cm−1, defining the carbonates' contribution. The silicates' contribution was obtained by integrating the spectrum in the 1200−800 cm−1 range, after subtracting the potential O-C-O out-of-plane mode of the carbonates. In the case where physisorbed and interlayer water represents an important contribution to the spectrum, the deconvolution of the C=C from the overlapping OH bending and aromatic C=C bands becomes challenging. In such a case, an estimate of the contribution of the OH bending mode can be evaluated from the measured OH stretching-mode integrated optical depth, applying a known intensity factor of 0.072, as discussed above. We thus estimated the C=C contribution either by direct deconvolution or by subtracting the expected OH bending contribution. If the OH integrated optical depth contribution represented a high fraction, more than half of the OH + C=C value, the evaluated aromatic C=C contribution was considered only as an upper limit. The column densities were then evaluated using the band strength previouslydiscussed, andarereportedinTableB.1 andpresented in Fig. 9.

thumbnail Fig. 11

Close-up of selected characteristic spectra of the Ryugu samples in the 1900−100 cm−1 range, with pure mineral spectra shown for reference (dolomite and magnesite, Roush 2021; breunnerite Mg0.78Fe0.22CO3, this work; FeS Henning & Mutschke 1997; magnetite, Glotch & Rossman 2009; saponite and serpentine, Glotch et al. 2007). In the upper part of the figure, vertical dashed dark-green lines have been added to highlight the positions of specific small features of serpentine. In the lower part of the figure, vertical dashed blue lines are guides to the eye to show the characteristic positions of IR features in dolomite, which is the main carbonate identified in this work (see text for details).

3.4 CH2/CH3 ratios

Close-ups of the absorption in the aliphatics region of the analysed samples are shown in Fig. 6. The methylene-to-methyl ratio of the aliphatic component is evaluated in two usual ways for comparison to other planetary samples. One commonly found indicator in the literature is the peak intensity maximum of the asymmetric stretching modes ICH2peak/ICH3peak${{I_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{peak}}}} \mathord{\left/ {\vphantom {{I_{{\rm{C}}{{\rm{H}}_2}}^{{\rm{peak}}}} {I_{{\rm{C}}{{\rm{H}}_3}}^{{\rm{peak}}}}}} \right. \kern-\nulldelimiterspace} {I_{{\rm{C}}{{\rm{H}}_3}}^{{\rm{peak}}}}}$. This indicator is straightforward and does not require a deconvolution of the individual contributing bands. To gain insight into the column density ratio of the methylene-to-methyl ratio, a deconvolution of the aliphatics band is more informative. Such a deconvolution is applied in the analysis of interstellar dust grains made of hydrogenated amorphous carbon (e.g. Dartois et al. 2004, Dartois et al. 2007, and references therein). A minimum of five bands contribute to the aliphatic signature: the asymmetric CH3 peaking around 2960 cm−1, the asymmetric CH2 peaking around 2920 cm−1, a contribution from resonances of this last mode with overtones of the deformation modes, giving rise to an additional contribution in the 2900 cm−1 region, called the shoulder or Fermi resonance (this shoulder or resonance can be seen in the reference aliphatic spectra in the upper left of Fig. 6), the symmetric CH3 peaking around 2880 cm−1, and the symmetric CH2 peaking around 2860 cm−1. The spectra were deconvolved using five Gaussians assuming a common full width at half maximum for each band, and constraining the asymmetric to symmetric methyl band area to be 0.2, as measured in reference spectra (Dartois et al. 2004). This parametrisation, based on physico-chemistry constraints, avoids having to deal with too many free parameters that otherwise lead to a poor convergence of the fitting. The resulting fits are shown in Fig. 6, in red for methyl bands, and in blue for methylene bands. The overall fitted band is overlaid in cyan on the samples spectra. Care must be taken, in some cases, with the contribution of carbonate combination and overtone modes to this spectral region, when carbonates are abundant in the measured sample region of interest. This is shown for one sample in Fig. 6, where the contribution from dolomite to the spectrum is shown in brown, and the resulting dolomite-subtracted spectrum is shown just above. Once this deconvolution wais performed, we integrated the asymmetric bands optical depth τCH2dv¯$\int {{\tau _{{\rm{C}}{{\rm{H}}_{\rm{2}}}}}} {\rm{d}}\bar v$ and τCH3dv¯$\int {{\tau _{{\rm{C}}{{\rm{H}}_{\rm{3}}}}}} {\rm{d}}\bar v$ and, using the integrated absorption coefficients of each mode, we derived the methylene-to-methyl column density ratios given in Table 2. As peak intensity ratios are commonly used in the literature, both peak ratios and column density ratios are displayed in Fig. 7 as a function of the CH to Si column density ratio. The methylene-to-methyl ratio expressed in column density spans the 3–5.5 range. The methylene-to-methyl column density ratios seems to show a slight correlation: the more abundant the CH, the higher the CH2/CH3 ratio seems to be. In the samples analysed, the ones with the lowest CH/Si abundance are from chamber C. These samples seem also to be related to the ones showing more aromatic C=C contributions. The interstellar medium (ISM) value for the -(CH2)/-(CH3) ratio of hydrogenated amorphous carbons observed in the diffuse interstellar medium (DISM) is lower, around about 2−2.5 (e.g. Dartois et al. 2007; Sandford et al. 1991), or around 1−1.2 in peak intensity ratios. Figure 8 shows the comparison of CH2/CH3 peak ratios with literature data. The CH2/CH3 peak ratios in Ryugu are higher than the one typically observed for petrologic type 1 chondrites, and are comparable to the range observed for type 2 and 3 chondrites (Alexander et al. 2014; Briani et al. 2013; Dionnet et al. 2018; Kebukawa et al. 2010, 2011, 2019; Kerraouch et al. 2022; Merouane et al. 2012; Orthous-Daunay 2011, 2013; Quirico et al. 2018; Yesiltas et al. 2014, 2015, 2022; Yesiltas & Kebukawa 2016). This does not necessarily imply a higher level of heating of Ryugu than in C1 type chondrites, and should be constrained by other petrologic estimates. The Alexander et al. (2014) analysis of four lithologies in Tagish Lake IOM that have undergone different degrees of aqueous alteration suggest a negative correlation between the ratio of CH2/CH3 stretching bands with the (CH2+CH3)/C=C ratio, the least altered phase having the lowest CH2/CH3 ratio associated with the highest H/C ratio among the lithologies studied. The parent body aqueous alteration would, in this case, increase the CH2/CH3 ratio, but at the expense of H consumption as the carbon yield for the different lithologies varies much less. On the other hand, a recent report on the experimental aqueous alteration of the IOM of Paris at 150°C, seems to show a decrease in the CH2/CH3 ratio upon aqueous alteration (Laurent et al. 2022). If the same view is applied to Ryugu, then the initial organic matter H/C should have been even higher before aqueous alteration than that observed today. Possible interactions between the organic phase during aqueous alteration of the Ryugu matrix should be investigated in more detail. The high CH2/CH3 ratio in Ryugu might also arise from intrinsic initial formation processes of the organic phase, which favoured a high CH2/CH3 ratio, meaning longer chains than usual when compared to unheated carbonaceous chondrites.

Table 2

Aliphatic CH2/CH3 asymmetric vibrations' peak intensity ratios, asymmetric modes' integrated band ratios, and column density ratios.

3.5 Abundances

A summary of the column densities of the identified IR functional groups is given in Table B.1. Relative proportions of the different groups (aliphatics, hydroxyl + interlayer and/or physisorbed water, carbonyl, carbonates, and silicates) from FTIR analysis of the different samples, shown per group (CH, OH, C=O, CO3, and Si) are graphically presented in Fig. 9. The C=C contribution is not represented in the figure, as in many cases only upper limits could be calculated from the deconvolution due to the overlap with the OH bending mode in the samples (see below). Independent measurements of the Ryugu samples showed the presence of polycyclic aromatic hydrocarbons (PAHs) at the level of a few tens of nmol g−1 (a few ppm) during the extraction of soluble organic matter (Naraoka et al. 2022; Aponte et al. 2022). Thus, a possible contribution of aromatic C-H in the FTIR spectra at 3050 cm−1 at that abundance level is expected to be low.

The IR hyper-spectral maps show that Ryugu is dominated by phyllosilicates, with no identification of anhydrous silicate in our samples (Figs. 5 and 11). CI chondrites contain up to 94% phyllosilicates (Alfing et al. 2019; King et al. 2015), which are Fe-bearing Mg-rich serpentine and saponite (Tomeoka & Buseck 1988; Zolensky et al. 1993; Kerridge 1976; Osawa et al. 2005), except for Flensburg, which does not contain saponite (Bischoff et al. 2021). In Orgueil, the (most abundant) finegrained phyllosilicates are usually mixed with fine ferrihydrite particles, whereas coarser phyllosilicates are devoid of ferrihydrite (Tomeoka & Buseck 1988). Ivuna, however, like the Ryugu samples, does not show the presence of ferrihydrite (Zolensky et al. 1993; Nakamura et al. 2023). Infrared spectra of CI chon-drites show a sharp absorption at 3685 cm−1, which is attributed to the presence of serpentine (Osawa et al. 2005; Takir et al. 2013). This feature is also observed in the Ryugu spectra, supporting the important presence of serpentine in the analysed samples, but saponite is also present in the Ryugu samples (see Fig. 11). Yokoyama et al. (2023) suggest a serpentine-to-saponite molar ratio of 3:2 in the Ryugu samples, which is compatible with our analyses (see Figs. 5 and 11). In CI chondrites, there are variations in the composition of phyllosilicates, with Mg# (defined as the atomic ratio 100 × Mg/(Mg+Fe)) varying from ~50 to ~85 (Kerridge 1976; Zolensky et al. 1993; Bischoff et al. 2021). Such heterogeneity in compositions could also explain the variations in the phyllosilicate features in the IR spectra of Ryugu. Carbonates in the Ryugu samples are dominated by dolomite, as is the case for CI chondrites (Johnson & Prinz 1993; Endress & Bischoff 1996). We see no evidence for other types of carbonates in our analyses, although rare occurrences of breunnerite and calcite have been reported in Ryugu by other studies (Yokoyama et al. 2023; Nakamura et al. 2023). In CI chondrites, the abundance of carbonates is around 0.53 vol%, the abundance of magnetite varies from 5–10 vol%, and pyrrhotite (the main Fe sulphide) is present at the level of a few percent (Alfing et al. 2019; King et al. 2015). In the Ryugu crushed samples, we can sometimes optically identify the fram-boidal magnetite, but from the Ryugu hyper-spectral maps it is difficult to quantify the abundance of magnetite or pyrrhotite, as they have a very high refractive index. Nevertheless, it can be observed that the total projected surface of the opaques' ROIs (as shown, for instance, in the bottom panels of Figs. 4 or C.1) is always a few times higher than that of the carbonates in the vast majority of the hyper-spectral maps, except for the samples in which a large single carbonate is measured. Using the IR column densities (N) reported in Table B.1 (and represented in Fig. 9), the volume ratio of carbonates to phyllosilicates can be evaluated. We assumed a density of ρdolomite ~ 2.8 g cm−3 for a molar mass of Mdolomite ~184 g mol−1 (it should be noted that nCO3=2${n_{{\rm{C}}{{\rm{O}}_3}}} = 2$ as there are two CO3 groups for one Ca atom in dolomite). For phyllosilicates we assumed the Glotch et al. (2007) formula for saponite and serpentine (Table A.1), with molar masses of MSaponite ~376 g mol−1 and Mserpentine ~281 g mol−1, and for Si atoms, we assumed nSi = 3.96 and nSi = 1.8 in each formula, respectively. We assumed densities of ρsaponite ~ 2.3 g cm−3 and ρserpentine ~2.6 g cm−3, respectively. Assuming an average 3:2 serpentine-to-saponite ratio (Yokoyama et al. 2023), we can calculate: f=volumecarbonatesvolumephy llosilicates,$f = {{{\rm{volum}}{{\rm{e}}_{{\rm{carbonates}}}}} \over {{\rm{volum}}{{\rm{e}}_{{\rm{phy}}\,{\rm{llosilicates}}}}}},$(9) f=NCO3NSiMdolomite/nCO3ρdolomite ρphyllosilicatesMphyllosilicates/nSi ,$f = {{{N_{{\rm{C}}{{\rm{O}}_3}}}} \over {{N_{{\rm{Si}}}}}}{{{M_{{{{\rm{dolomite}}} \mathord{\left/ {\vphantom {{{\rm{dolomite}}} {{n_{{\rm{C}}{{\rm{O}}_3}}}}}} \right. \kern-\nulldelimiterspace} {{n_{{\rm{C}}{{\rm{O}}_3}}}}}}}} \over {{\rho _{{\rm{dolomite}}}}}}\left\langle {{{{\rho _{{\rm{phyllosilicates}}}}} \over {{M_{{{{\rm{phyllosilicates}}} \mathord{\left/ {\vphantom {{{\rm{phyllosilicates}}} {{n_{{\rm{Si}}}}}}} \right. \kern-\nulldelimiterspace} {{n_{{\rm{Si}}}}}}}}}}} \right\rangle ,$(10)

where Mdolomite/nCO3=MdolomitenCO3=Mdolomite2,${M_{{{{\rm{dolomite}}} \mathord{\left/ {\vphantom {{{\rm{dolomite}}} {{n_{{\rm{C}}{{\rm{O}}_3}}}}}} \right. \kern-\nulldelimiterspace} {{n_{{\rm{C}}{{\rm{O}}_3}}}}}}} = {{{M_{{\rm{dolomite}}}}} \over {{n_{{\rm{C}}{{\rm{O}}_3}}}}} = {{{M_{{\rm{dolomite}}}}} \over 2},$(11)

and Mphyllosilicates/nSiρphyllosilicates =35Mserpentineρserpentine/1.8+25Msaponiteρsaponite/3.96,$\left\langle {{{{M_{{{{\rm{phyllosilicates}}} \mathord{\left/ {\vphantom {{{\rm{phyllosilicates}}} {{n_{{\rm{Si}}}}}}} \right. \kern-\nulldelimiterspace} {{n_{{\rm{Si}}}}}}}}} \over {{\rho _{{\rm{phyllosilicates}}}}}}} \right\rangle = {3 \over 5}{{{M_{{\rm{serpentine}}}}} \over {{\rho _{{\rm{serpentine}}}}}}/1.8 + {2 \over 5}{{{M_{{\rm{saponite}}}}} \over {{\rho _{{\rm{saponite}}}}}}/3.96,$(12)

and then Volume(carbonatesphy llosilicates+carbonates)=(1+1/f)1.${\rm{Volume}}\left( {{{{\rm{carbonates}}} \over {{\rm{phy}}\,{\rm{llosilicates + carbonates}}}}} \right) = {\left( {{{1 + 1} \mathord{\left/ {\vphantom {{1 + 1} {\rm{f}}}} \right. \kern-\nulldelimiterspace} {\rm{f}}}} \right)^{ - 1}}.$(13)

In the Ryugu samples, the average relative volume of carbonates to phyllosilicate and carbonates using Eq. (13) is shown in Fig. 10. Error bars were calculated assuming the extreme cases of pure serpentine or pure saponite composition for the phyllosilicates. The volume ratios for most Ryugu data are at the high end of the range of volume ratios in CI chondrites reported in Alfing et al. (2019) and King et al. (2015). In a few cases, the volume ratios are very high because these samples contain large, single carbonates. It is important to note that the IR spectra are sensitive to the whole volume probed (column density), and thus the numbers are slightly different from the volume ratios estimated from the polished sections. The volume fractions of dolomite were reported as 1.8% and 2.9% in plates 5 and 6, respectively, of the C0002 samples, by Nakamura et al. (2023). These values translate to 2.1% and 3.2% when recalculated using Eq. (13), and these are plotted as red dots in Fig. 10.

The mean C/Si abundance ratio by weight for CI chondrites is 3.48 wt%/ 10.7 wt% (e.g. Palme et al. 2014; Lodders et al. 2009; Anders & Grevesse 1989), corresponding to a C/Si atomic ratio of about 0.76. For Ryugu the mean value measured for C is 4.63 wt%, with 3.08 wt% attributed to organic carbon (Yurimoto et al. 2022; Yokoyama et al. 2023), and thus the Ryugu organic fraction is in line with CI values. The CH to Si abundance we measured varies from 0.5 to 1.5 among the different samples measured (Table B.1). With the CH2/CH3 ratio measured (Fig. 7), it corresponds to an aliphatic C/Si atomic ratio of about half of this number, that is to say an atomic Caliphatics/Si ratio in the 0.25−0.75 range in our samples.

The mean carbonyl-to-CH ratio retrieved from the decon-volution is in the 0.013−0.064 range, with a mean around 0.034 ± 0.032 (2σ). The lowest values can, however, be slightly underestimated in the deconvolution if the OH bending mode contribution in the wing of the band is significant. Given the measured CH2/CH3 ratio of about two (Fig. 7), the C=O/CHx ratio (i.e. per CHx group) would be about twice the C=O/CH ratio, in the range 0.029−0.14. Previous nuclear magnetic resonance studies (Cody & Alexander 2005; Derenne & Robert 2010) measured a C=O/CHx ratio for Orgueil that is higher than this range, suggesting that the Ryugu samples could contain more CH than Orgueil. To confirm this statement, it would be interesting to compare the C=O/CHx ratio of Orgueil and Ryugu IOM residues using the same deconvolution.

As stated before, the backbone on which these aliphatics are attached also comprises aromatic carbons, which are observed also via Raman spectroscopy D and G bands (e.g. Bonal et al. 2022), and are widespread, as are the CH aliphatics in the IR measurements. Abundance with Raman is not easily accessible as this is emission spectroscopy and also because these C=C modes can be resonant. In the IR measurements, this aromatic C=C band contribution is accessible, although it is still less straightforward to derive in the whole-rock spectra as it can be overlapped with strong inter-layer bending modes from physisorbed water and saponite bands, possibly contributing to the spectra. We thus decided to calculate this contribution only as an upper limit to the aromatic C=C contribution, which is close to the true C=C column density for the spectra where these OH bending and/or saponite contributions are the least intense. This aromatic C column density is of the same order as the aliphatic C one. Some Ryugu fragments measured depart from the mean, such as the C0109 sample (see, lower spectrum of Fig. 2), for which the aromatic carbon contribution after spectral decomposition, as observed around 1600 cm−1, is higher. For these samples, the aromatic carbon column density is larger than the aliphatic carbon one, and thus has more polyaromatic moieties. The carbonate contribution is highly variable at the scale probed, and in a few samples we observed very large individual carbonate grains. In the measured samples, except for the carbonate grain we focused specifically on, the carbonate CO3 group variations observed with respect to the aliphatic carbon give an atomic Ccarbonates/Caliphatics ratio in the 0.03-0.34 range, with a median of 0.14 ± 0.08 The size distribution histogram of carbonates extracted from the maps show no particular behaviour. The biggest carbonates have large sizes of several thousands of square micrometres in projected size, and are isolated from a continuum of much smaller carbonates. The biggest one we measured (A0108-19 sample) is largely dominated by a dolomite composition from the positions of the bands in the FIR.

4 Summary and outlooks

Several Ryugu grains from chambers A and C were analysed by synchrotron IR micro-spectroscopy at the diffraction limit (~6 × 6 μm), providing hyper-spectral maps that record the spatial distributions of the main bands of organics, carbonates, and phyllosilicates. Strong absorption of the IR continuum also revealed the presence of minerals with a high refractive index in the samples, such as magnetite and Fe sulphides. Additional measurements allowed us to extend the wavelength range to the FIR at the scales of the grains studied, providing better assignments for the carbonates' and phyllosilicates' composition. The Ryugu IR spectra are compatible with CI chondrites spectra, with the silicates dominated by phyllosilicates. We did not identify any anhydrous silicates in the samples analysed. The carbonates are dominated by dolomite, as also observed in CI chondrites. The hyper-spectral maps show that the aliphatics organics are distributed all over the samples, and intimately mixed with the phyllosilicates at the scale of a few microns probed by the synchrotron beam. The methylene-to-methyl (CH2/CH3) column density ratio of the analysed samples varies between about 3.5 and five (or between about 1.5–2.5 in peak intensity ratio), which is compatible with what is typically observed in IOM from carbonaceous chondrites (Quirico et al. 2018; Kebukawa et al. 2011), and higher than the value determined in the DISM (Dartois et al. 2007; Sandford et al. 1991). The possible correlation observed between the CH2/CH3 ratio and the CH/Si ratio (Fig. 7) suggests a possible mixing of several organic phases, with an organic phase more ramified or with shorter chains associated with the phyllosilicates, and slightly longer chains, and therefore a CH2-richer phase, when the organics' abundance is higher locally.

Acknowledgements

This work was carried out in the frame of the Hayabusa2 Initial Analysis Teams, in the IOM sub-team led by Prof. H. Yabuta and the Stone sub-team led by Prof. T. Nakamura. Spectral data will be available on the Data ARchives and Transmission System (DARTS) of the Japan Aerospace Exploration Agency (JAXA). Part of the equipment used in this work has been financed by the French INSU-CNRS program "Physique et Chimie du Milieu Interstellaire" (PCMI). This work was supported in France by the Centre National d'Etudes Spatiales (CNES) for the participation to analyses of samples returned by the Hayabusa2 mission. We acknowledge synchrotron SOLEIL for provision of synchrotron radiation facilities for access to the SMIS beamline. We thank Sylvie Heron for providing us the n-tetracosane reference sample. We are grateful to Timothy D. Glotch for his constructive comments and recommendations.

Appendix A Band strengths

Table A.1

Integrated absorption coefficients (band strengths) 𝒜 (cm/[group]) of selected IR vibrations identified in the Ryugu samples, with their respective parameters and corresponding literature references.

Appendix B Chemical composition

Table B.1

Chemical composition of the Ryugu samples from FTIR

Appendix C Synchrotron maps

Synchrotron maps of the various samples are displayed in each figure for: a (silicates, green), b (carbonates, red), c (opaques, yellow-orange), d (OH stretch, blue), e (CH stretch, magenta), and f (corresponding visible image).

C.1. Chamber A

thumbnail Fig. C.1

Chamber A: A0064-FO013-17-top, ROI size distributions determined in a similar way as explained in caption of Fig. 4.

thumbnail Fig. C.2

Chamber A: A0108-15-03-bot-a

thumbnail Fig. C.3

Chamber A: A0108-15-03-bot-b

thumbnail Fig. C.4

Chamber A: A0108-15-03 Bottom

thumbnail Fig. C.5

Chamber A: A0108-19-03-bot-a

thumbnail Fig. C.6

Chamber A: A0108-19-03-top-a

thumbnail Fig. C.7

Chamber A: A0108-19b-bot-a

thumbnail Fig. C.8

Chamber A: A0108-19b-bot-b

thumbnail Fig. C.9

Chamber A: A0108-19b-top-a

thumbnail Fig. C10

Chamber A: A0108-19b-top-b

thumbnail Fig. C.11

Chamber A: A0108-19b-top-c

thumbnail Fig. C.12

Chamber A: A0106-13-bot-a

thumbnail Fig. C.13

Chamber A: A0106-13-top-a

thumbnail Fig. C.14

Chamber A: A0106-13-top-b

thumbnail Fig. C.15

Chamber A: A0064-F0020-bot-a

thumbnail Fig. C.16

Chamber A: A0064-F0020-top-a

thumbnail Fig. C.17

Chamber A: A0064-F0021-bot-a

thumbnail Fig. C.18

Chamber A: A0064-F0021-top-a

thumbnail Fig. C.19

Chamber C: C0046-FO004-010-topl

thumbnail Fig. C.20

Chamber C: C0046-FO004-010-bot-a

thumbnail Fig. C.21

Chamber C: C0002-FC016-top-a

thumbnail Fig. C.22

Chamber C: C0002-FC016-bot-a

thumbnail Fig. C.23

Chamber C: C0040-FC025-bot-a

thumbnail Fig. C.24

Chamber C: C0109-4-bot-a

thumbnail Fig. C.25

Chamber C: C0109-04-top-a

thumbnail Fig. C.26

Chamber C: C0109-10-top-a

thumbnail Fig. C.27

Chamber C: C0109-10-bot-a

thumbnail Fig. C.28

Chamber C: C0109-10 Bottom

thumbnail Fig. C.29

Chamber C: C0109-10 Bottom

thumbnail Fig. C.30

Chamber C: C0109-10 Bottom

thumbnail Fig. C.31

Chamber C: C0057-5-bot-a

thumbnail Fig. C.32

Chamber C: C0057-5-bot-b

thumbnail Fig. C.33

Chamber C: C0057-5-bot-c

thumbnail Fig. C.34

Chamber C: C0057-5-bot-d

References

  1. Alexander, C. M. O., Fogel, M., Yabuta, H., & Cody, G. D. 2007, Geochim. Cosmochim. Acta, 71, 4380 [NASA ADS] [CrossRef] [Google Scholar]
  2. Alexander, C. O., Cody, G., Kebukawa, Y., et al. 2014, Meteor. Planet. Sci., 49, 503 [CrossRef] [Google Scholar]
  3. Alfing, J., Patzek, M., & Bischoff, A. 2019, Geochemistry, 79, 125532 [NASA ADS] [CrossRef] [Google Scholar]
  4. Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 [Google Scholar]
  5. Anthony, J. W., Bideaux, R. A., Bladh, K. W., & Nichols, M. C. 2001, Handbook of Mineralogy, 1 (Mineral Data Publ.) [Google Scholar]
  6. Aponte, J. C., Dworkin, J., Glavin, D., et al. 2022, Earth, Planets and Space, submitted [Google Scholar]
  7. Bertie, J. E., & Lan, Z. 1996, Appl. Spectrosc., 50, 1047 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bischoff, A., Alexander, C. M. O., Barrat, J.-A., et al. 2021, Geochim. Cosmochim. Acta, 293, 142 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bonal, L., Quirico, E., Montagnac, G., et al. 2022, in 53rd Lunar and Planetary Science Conference, 2678, LPI Contributions, 1331 [NASA ADS] [Google Scholar]
  10. Briani, G., Quirico, E., Gounelle, M., et al. 2013, Geoch. Cosmochim. Acta, 122, 267 [NASA ADS] [CrossRef] [Google Scholar]
  11. Cody, G. D., & Alexander, C. M. O. D. 2005, Geochim. Cosmochim. Acta, 69, 1085 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dartois, E., Muñoz Caro, G. M., Deboffle, D., & d'Hendecourt, L. 2004, A & A, 423, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Dartois, E., Geballe, T. R., Pino, T., et al. 2007, A & A, 463, 635 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Dartois, E., Engrand, C., Brunetto, R., et al. 2013, Icarus, 224, 243 [NASA ADS] [CrossRef] [Google Scholar]
  15. Dazzi, A., & Prater, C. B. 2017, Chem. Rev., 117, 5146 [CrossRef] [Google Scholar]
  16. Dazzi, A., Prazeres, R., Glotin, F., & Ortega, J. 2005, Opt. Lett., 30, 2388 [NASA ADS] [CrossRef] [Google Scholar]
  17. Derenne, S., & Robert, F. 2010, Meteor. Planet. Sci., 45, 1461 [NASA ADS] [CrossRef] [Google Scholar]
  18. Dionnet, Z., Aléon-Toppani, A., Borondics, F., et al. 2018, Microsc. Microanal., 24, 2100 [NASA ADS] [CrossRef] [Google Scholar]
  19. Endress, M., & Bischoff, A. 1996, GCA, 60, 489 [NASA ADS] [Google Scholar]
  20. Flynn, G. J., Keller, L. P., Feser, M., Wirick, S., & Jacobsen, C. 2003, Geochim. Cosmochim. Acta, 67, 4791 [CrossRef] [Google Scholar]
  21. Glotch, T. D., & Rossman, G. R. 2009, Icarus, 204, 663 [NASA ADS] [CrossRef] [Google Scholar]
  22. Glotch, T. D., Rossman, G. R., & Aharonson, O. 2007, Icarus, 192, 605 [NASA ADS] [CrossRef] [Google Scholar]
  23. Godard, M., Geballe, T. R., Dartois, E., & Muñoz Caro, G. M. 2012, A & A, 537, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Henning, T., & Mutschke, H. 1997, A & A, 327, 743 [NASA ADS] [Google Scholar]
  25. Jiang, C.-L., Zeng, W., Liu, F.-S., Tang, B., & Liu, Q.-J. 2019, J. Phys. Chem. Solids, 131, 1 [NASA ADS] [CrossRef] [Google Scholar]
  26. Johnson, C. A., & Prinz, M. 1993, Geochim. Cosmochim. Acta, 57, 2843 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kebukawa, Y., Alexander, C. M. O. D., & Cody, G. D. 2011, Geochim. Cosmochim. Acta, 75, 3530 [NASA ADS] [CrossRef] [Google Scholar]
  28. Kebukawa, Y., Nakashima, S., Ishikawa, M., et al. 2010, Meteor. Planet. Sci., 45, 394 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kebukawa, Y., Alexander, C. M. O., & Cody, G. D. 2019, Meteor. Planet. Sci., 54, 1632 [NASA ADS] [CrossRef] [Google Scholar]
  30. Keller, L. P., Bajt, S., Baratta, G. A., et al. 2006, Science, 314, 1728 [CrossRef] [Google Scholar]
  31. Kerraouch, I., Kebukawa, Y., Bischoff, A., et al. 2022, Geochim. Cosmochim. Acta, 334, 155 [NASA ADS] [CrossRef] [Google Scholar]
  32. Kerridge, J. F. 1976, Earth Planet. Sci. Lett., 29, 194 [CrossRef] [Google Scholar]
  33. King, A., Schofield, P., Howard, K., & Russell, S. 2015, Geochim. Cosmochim. Acta, 165, 148 [NASA ADS] [CrossRef] [Google Scholar]
  34. Kloprogge, J. T., & Ponce, C. P. 2021, Minerals, 11, 112 [NASA ADS] [CrossRef] [Google Scholar]
  35. Laurent, B., Holin, M., Bernard, S., et al. 2022, in 53rd Lunar and Planetary Science Conference, 2678, LPI Contributions, 1432 [Google Scholar]
  36. Lodders, K., Palme, H., & Gail, H. P. 2009, Landolt-Börnstein (Springer-Verlag Berlin Heidelberg), 4B, 712 [Google Scholar]
  37. Mathurin, J., Dartois, E., Pino, T., et al. 2019, A & A, 622, A160 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Mathurin, J., Deniset-Besseau, A., & Dazzi, A. 2020, Acta Phys. Polon. A, 137, 29 [NASA ADS] [CrossRef] [Google Scholar]
  39. Mathurin, J., Dartois, E., Dazzi, A., et al. 2022a, in 53rd Lunar and Planetary Science Conference, 2678, LPI Contributions, 2142 [Google Scholar]
  40. Mathurin, J., Deniset-Besseau, A., Bazin, D., et al. 2022b, J. Appl. Phys., 131, 010901 [NASA ADS] [CrossRef] [Google Scholar]
  41. Merouane, S., Djouadi, Z., d'Hendecourt, L. L. S., Zanda, B., & Borg, J. 2012, ApJ, 756, 154 [NASA ADS] [CrossRef] [Google Scholar]
  42. Mooney, T., & Knacke, R. F. 1985, Icarus, 64, 493 [NASA ADS] [CrossRef] [Google Scholar]
  43. Mukai, T. 1990, Evolution of Interstellar Dust and Related Topics, 397 [Google Scholar]
  44. Nakamura, T., Matsumoto, M., Amano, K., et al. 2023, Science in press, https://doi.org/10.1126/science.abn8671 [Google Scholar]
  45. Naraoka, H., Takano, Y., Dworkin, J. P., et al. 2022, Science, submitted [Google Scholar]
  46. Orthous-Daunay, F.-R. 2011, Ph.D. thesis, Université de Grenoble, France [Google Scholar]
  47. Orthous-Daunay, F. R., Quirico, E., Beck, P., et al. 2013, Icarus, 223, 534 [NASA ADS] [CrossRef] [Google Scholar]
  48. Osawa, T., Kagi, H., Nakamura, T., & Noguchi, T. 2005, Meteor. Planet. Sci., 40, 71 [NASA ADS] [CrossRef] [Google Scholar]
  49. Palme, H., Lodders, K., & Jones, A. 2014, in Planets, Asteriods, Comets and The Solar System, 2, ed. A. M. Davis, 15 [NASA ADS] [Google Scholar]
  50. Pendleton, Y., Sandford, S., Allamandola, L., Tielens, A., & Sellgren, K. 1994, ApJ, 437, 683 [NASA ADS] [CrossRef] [Google Scholar]
  51. Phan, V. T. H., Quirico, E., Beck, P., et al. 2021, Spectrochim. Acta A: Mol. Spectrosc., 259, 119853 [NASA ADS] [CrossRef] [Google Scholar]
  52. Querry, M. R. 1985, Optical constants, Tech. rep. (Missouri Univ-Kansas City) [Google Scholar]
  53. Quirico, E., Bonal, L., Beck, P., et al. 2018, Geochim. Cosmochim. Acta, 241, 17 [CrossRef] [Google Scholar]
  54. Roush, T. L. 2021, Icarus, 354, 114056 [NASA ADS] [CrossRef] [Google Scholar]
  55. Roush, T., Pollack, J., & Orenberg, J. 1991, Icarus, 94, 191 [NASA ADS] [CrossRef] [Google Scholar]
  56. Sandford, S., Allamandola, L., Tielens, A., et al. 1991, ApJ, 371, 607 [NASA ADS] [CrossRef] [Google Scholar]
  57. Sawada, H., Okazaki, R., Tachibana, S., et al. 2017, Space Sci. Rev., 208, 81 [CrossRef] [Google Scholar]
  58. Tachibana, S., Abe, M., Arakawa, M., et al. 2014, Geochem. J., 48, 571 [NASA ADS] [CrossRef] [Google Scholar]
  59. Tachibana, S., Sawada, H., Okazaki, R., et al. 2022, Science, 375, 1011 [NASA ADS] [CrossRef] [Google Scholar]
  60. Takir, D., Emery, J. P., Mcsween, H. Y.Jr., et al. 2013, Meteor. Planet. Sci., 48, 1618 [NASA ADS] [Google Scholar]
  61. Tomeoka, K., & Buseck, P. R. 1988, GCA, 52, 1627 [NASA ADS] [Google Scholar]
  62. Tsuda, Y., Saiki, T., Terui, F., et al. 2020, Acta Astronaut., 171, 42 [NASA ADS] [CrossRef] [Google Scholar]
  63. Yabuta, H., Cody, C. D., Engrand, C., et al. 2022, Science, submitted [Google Scholar]
  64. Yesiltas, M., & Kebukawa, Y. 2016, Meteor. Planet. Sci., 51, 584 [NASA ADS] [CrossRef] [Google Scholar]
  65. Yesiltas, M., Kebukawa, Y., Peale, R. E., et al. 2014, Meteor. Planet. Sci., 49, 2027 [NASA ADS] [CrossRef] [Google Scholar]
  66. Yesiltas, M., Peale, R. E., Unger, M., Sedlmair, J., & Hirschmugl, C. J. 2015, Meteor. Planet. Sci., 50, 1684 [NASA ADS] [CrossRef] [Google Scholar]
  67. Yesiltas, M., Kebukawa, Y., Glotch, T. D., et al. 2022, Meteor. Planet. Sci., 57, 1665 [NASA ADS] [CrossRef] [Google Scholar]
  68. Yokoyama, T., Nagashima, K., Nakai, I., et al. 2023, Science, in press, https://doi.org/10.1126/science.abn7850 [Google Scholar]
  69. Yurimoto, H., Hayabusa2-Initial-Analysis Chemistry Team, & Hayabusa2-Initial-Analysis Core 2022, in 53rd Lunar and Planetary Science Conference, 2678, LPI Contributions, 1377 [Google Scholar]
  70. Zolensky, M., Barrett, R., & Browning, L. 1993, Geochim. Cosmochim. Acta, 57, 3123 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Ryugu samples analysed.

Table 2

Aliphatic CH2/CH3 asymmetric vibrations' peak intensity ratios, asymmetric modes' integrated band ratios, and column density ratios.

Table A.1

Integrated absorption coefficients (band strengths) 𝒜 (cm/[group]) of selected IR vibrations identified in the Ryugu samples, with their respective parameters and corresponding literature references.

Table B.1

Chemical composition of the Ryugu samples from FTIR

All Figures

thumbnail Fig. 1

Representative Ryugu sample A0108-15 transmittance spectra (6 × 6 μm spots) showing a carbonate-rich region (A), a region with inorganic 'opaques' (probably magnetite and/or Fe sulphides) (B), and the map average spectrum (green) corresponding to map 1 reported in Fig. 3. Vibrational modes of interest (including the ones used to build the coloured maps such as in Fig. 3) are shown with arrows.

In the text
thumbnail Fig. 2

Close-up and deconvolution of the individual bands for the global spectra (corresponding to the maps of Figs. 3 and 4) recorded for sample A0108-15 (03-top-a), upper spectrum, with a high physisorbed and/or interlayer water contribution, and C0040-FC025 (bot-b), second spectrum from top, dominated by the OH hydroxyl contribution. Individual bands are represented by hatching: OH stretch (blue in the 3000–3700 cm−1 region), CH stretch (magenta), C=0 carbonyl stretch (brown), OH bending (blue around 1640 cm−1), aromatic C=C bending (black), carbonates (red), and silicates (green). Each wavenumber range is normalised with regard to the high wavenumber range by a factor (indicated in parenthesis) for clarity of the deconvolution. The spectrum of a particularly C-aromatic-rich grain, C0109-04 (top-a), is also shown. For comparison, a spectrum taken from Orgueil, which was prepared and recorded in the same conditions as the Ryugu samples, is shown at the bottom. See text for details.

In the text
thumbnail Fig. 3

Chamber A Ryugu sample #A0108-15 (03) (top diamond window). Upper left: Visible image recording the entire flattened grain with an extent of about 150 μm. The positions A and B correspond to the individual spectra shown in Fig. 1, showing a carbonate-rich region (A) and a region with inorganic opaques (probably magnetite and/or Fe sulphides) (B). The regions of interest mapped with the synchrotron measurements are numbered from 1 to 3. Top right and middle panel sets: Coloured phase maps for regions 1 to 3 corresponding, for each map, to silicates (a, green), carbonates (b, red), opaques, (c, yellow-orange), OH stretch (d, blue), CH stretch (e, magenta), and respective visible image (f). It should be noted that a few outlier (white) points with apparent high optical depth may be observed in the OH and CH stretch maps. They coincide with the presence of opaques or are on maps' borders. They are due to low signal-to-noise ratios, which lead to a higher uncertainty on the baseline determination (and correction). Lower panel: Images representing the same region mapped by AFM-IR (left, composite map of 30 μm × 30 μm, green, silicates; blue, carbonates; red, opaques), and mapped with synchrotron at the highest achievable spatial resolution (right,~6 μm aperture with a step size of ~3 μm, i.e. close to Nyquist sampling).

In the text
thumbnail Fig. 4

Chamber C Ryugu sample C0040-FC025 (bottom diamond window). Phase maps of silicates (a, green), carbonates (b, red), opaques (c, yellow-orange), OH stretch (d, blue), CH stretch (e, magenta), and corresponding visible images (f). Bottom: images representing the edge detection of individual ROIs based on the intensity map of carbonates (left), and of opaques (middle), and resulting size distribution histograms (right). The x-axis represents the projected surface of each region in μm2 (multiples of the spatial resolution of 9 μm2). As the samples are flattened to ~1 μm, this projected surface is directly proportional to the same volume value.

In the text
thumbnail Fig. 5

Global MIR to FIR spectra of the grains analysed, including large regions of interest for some of the grains from chamber A (left) and chamber C (right). Reference spectra for major species identified in the Ryugu samples spectra such as carbonates: dolomite, Roush (2021); magnesite, Roush (2021); breunnerite Mg0.78Fe0.22CO3, (this work) and phyllosilicates: saponite, serpentine, Glotch et al. (2007), are displayed for comparison (mineral compositions are given in Table A.1). For some grains with a large carbonate contribution to the spectrum, we also display the spectrum of dolomite, the carbonate the closest to the measured band positions and intensities. In addition, in the lower right panel are displayed the spectra of two CI1 meteorites closely related to the composition of the Ryugu samples, Orgueil (this work, measured under the same conditions in a diamond cell) and Ivuna (measured in a KBr pellet at 300°C, from Beck, Pierre (2012): MIR absorbance spectra of bulk CI chondrites in KBr pellets at ambient temperature, 150°C and 300°C. SSHADE/GhoSST, OSUG Data Center, https://doi.org/10.26302/SSHADE/EXPERIMENT_LB_20170719_001).

In the text
thumbnail Fig. 6

Aliphatic stretches' spectral decomposition used to calculate the methylene-to-methyl ratios. The spectra are fitted with five components corresponding to the asymmetric CH3 stretching mode (~2960 cm−1), the asymmetric CH2 stretch and its red shoulder (~2920−2900 cm−1), and the symmetric stretches of the same groups ~2880-2850 cm−1). The approximate positions of these modes are labelled above the upper left spectrum. The CH3 contributions are shown in red, and the CH2 ones in blue. In a few spectra, a contribution from dolomite affects the profile. For these, the expected dolomite spectrum contribution is shown in brown and the dolomite-subtracted spectrum is shown above. Spectra with a high fraction of dolomite contribution were not deconvolved. In the upper right corner are shown reference spectra of linear and branched alkanes, an alkyl benzene and a cycloalkane, with varying methylene-to-methyl ratios: isocetane, squalane, n-tetradecane (purchased from Merck), and n-tetracosane (provided by Sylvie Heron) were measured as received with our FTIR-microscope; n-hexane and n-decane spectra are from Dartois et al. (2004); triethylbenzene and dimethylcyclohexane spectra are from the NIST database: https://webbook.nist.gov/chemistry/. Although they are specific species, this series shows the expected asymmetric CH2/CH3 stretch band ratio trend.

In the text
thumbnail Fig. 7

Methylene-to-methyl (CH2/CH3) ratios (chamber A: triangles; chamber C: circles) expressed as the asymmetric stretching mode peak maximum intensity (blue) and as column density ratio (red), derived from the deconvolution of the aliphatic bands, as a function of the CH-to-Si column density ratio. See text for details.

In the text
thumbnail Fig. 8

Comparison of methylene-to-methyl (CH2/CH3) intensity peak ratios for Ryugu bulk sample measurements (this work) with measurements from the literature for diffuse ISM observations (Sandford et al. 1991; Pendleton et al. 1994; Dartois et al. 2007; Godard et al. 2012), Stardust samples (comet 81P/Wild2, Keller et al. 2006), IDPs (Flynn et al. 2003), and meteorites bulk and IOM samples (Alexander et al. 2014; Briani et al. 2013; Dionnet et al. 2018; Kebukawa et al. 2010, 2011, 2019; Kerraouch et al. 2022; Merouane et al. 2012; Orthous-Daunay 2011, 2013; Quirico et al. 2018; Yesiltas et al. 2014, 2015, 2022; Yesiltas & Kebukawa 2016).

In the text
thumbnail Fig. 9

Relative proportions of the different functional groups (aliphatics, hydroxyl + interlayer and/or physisorbed water, carbonyl, carbonates, and silicates) from the FTIR analysis of the different samples, shown per functional group CH, OH, C=O, CO3, and Si, respectively.

In the text
thumbnail Fig. 10

Relative volume ratio of carbonates to (carbonates + phyllosili-cates) in the Ryugu samples (in %), calculated using Eq. (13). Samples from chamber A are plotted as blue triangles and chamber C samples are represented as blue dots (see text for details). Two measurements on polished sections of chamber C samples are reported for comparison as filled red dots (Table S6 of Nakamura et al. 2023).

In the text
thumbnail Fig. 11

Close-up of selected characteristic spectra of the Ryugu samples in the 1900−100 cm−1 range, with pure mineral spectra shown for reference (dolomite and magnesite, Roush 2021; breunnerite Mg0.78Fe0.22CO3, this work; FeS Henning & Mutschke 1997; magnetite, Glotch & Rossman 2009; saponite and serpentine, Glotch et al. 2007). In the upper part of the figure, vertical dashed dark-green lines have been added to highlight the positions of specific small features of serpentine. In the lower part of the figure, vertical dashed blue lines are guides to the eye to show the characteristic positions of IR features in dolomite, which is the main carbonate identified in this work (see text for details).

In the text
thumbnail Fig. C.1

Chamber A: A0064-FO013-17-top, ROI size distributions determined in a similar way as explained in caption of Fig. 4.

In the text
thumbnail Fig. C.2

Chamber A: A0108-15-03-bot-a

In the text
thumbnail Fig. C.3

Chamber A: A0108-15-03-bot-b

In the text
thumbnail Fig. C.4

Chamber A: A0108-15-03 Bottom

In the text
thumbnail Fig. C.5

Chamber A: A0108-19-03-bot-a

In the text
thumbnail Fig. C.6

Chamber A: A0108-19-03-top-a

In the text
thumbnail Fig. C.7

Chamber A: A0108-19b-bot-a

In the text
thumbnail Fig. C.8

Chamber A: A0108-19b-bot-b

In the text
thumbnail Fig. C.9

Chamber A: A0108-19b-top-a

In the text
thumbnail Fig. C10

Chamber A: A0108-19b-top-b

In the text
thumbnail Fig. C.11

Chamber A: A0108-19b-top-c

In the text
thumbnail Fig. C.12

Chamber A: A0106-13-bot-a

In the text
thumbnail Fig. C.13

Chamber A: A0106-13-top-a

In the text
thumbnail Fig. C.14

Chamber A: A0106-13-top-b

In the text
thumbnail Fig. C.15

Chamber A: A0064-F0020-bot-a

In the text
thumbnail Fig. C.16

Chamber A: A0064-F0020-top-a

In the text
thumbnail Fig. C.17

Chamber A: A0064-F0021-bot-a

In the text
thumbnail Fig. C.18

Chamber A: A0064-F0021-top-a

In the text
thumbnail Fig. C.19

Chamber C: C0046-FO004-010-topl

In the text
thumbnail Fig. C.20

Chamber C: C0046-FO004-010-bot-a

In the text
thumbnail Fig. C.21

Chamber C: C0002-FC016-top-a

In the text
thumbnail Fig. C.22

Chamber C: C0002-FC016-bot-a

In the text
thumbnail Fig. C.23

Chamber C: C0040-FC025-bot-a

In the text
thumbnail Fig. C.24

Chamber C: C0109-4-bot-a

In the text
thumbnail Fig. C.25

Chamber C: C0109-04-top-a

In the text
thumbnail Fig. C.26

Chamber C: C0109-10-top-a

In the text
thumbnail Fig. C.27

Chamber C: C0109-10-bot-a

In the text
thumbnail Fig. C.28

Chamber C: C0109-10 Bottom

In the text
thumbnail Fig. C.29

Chamber C: C0109-10 Bottom

In the text
thumbnail Fig. C.30

Chamber C: C0109-10 Bottom

In the text
thumbnail Fig. C.31

Chamber C: C0057-5-bot-a

In the text
thumbnail Fig. C.32

Chamber C: C0057-5-bot-b

In the text
thumbnail Fig. C.33

Chamber C: C0057-5-bot-c

In the text
thumbnail Fig. C.34

Chamber C: C0057-5-bot-d

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.