Open Access
Issue
A&A
Volume 698, June 2025
Article Number A143
Number of page(s) 195
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202554411
Published online 19 June 2025

© The Authors 2025

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.

Open Access funding provided by Max Planck Society.

1 Introduction

Molecular clouds are the birthplaces of stars and their planetary systems. In these interstellar environments, star and planet formation results from an intricate interplay between physical and chemical processes. This work focuses on the chemistry that occurs in star-forming regions, in particular on the formation and evolution of complex organic molecules (COMs), which are defined as carbon bearing molecules that contain at least six atoms (Herbst & van Dishoeck 2009). Interstellar chemistry proceeds in both the gas phase and the solid phase, the latter corresponding to the surface and icy mantles of dust grains. While both phases are involved in the formation and evolution of COMs (e.g., Herbst & Leung 1986; Charnley et al. 1992; van Dishoeck & Blake 1998; Garrod et al. 2008), the major role of the solid phase in the production of COMs has become more obvious in the past decade, in particular since the detection of gas-phase saturated COMs at low temperatures in prestellar cores and protostellar envelopes (e.g., Öberg et al. 2010; Bacmann et al. 2012; Vastel et al. 2014) and, more recently, with the identification of several COMs in interstellar ices with the James Webb Space Telescope (JWST, e.g., Yang et al. 2022; Rocha et al. 2024; Nazari et al. 2024). While earlier kinetic models of interstellar chemistry relied on the increased mobility of radicals at the surface of dust grains during the warm-up phase of protostellar evolution to explain the production of many COMs (e.g., Garrod & Herbst 2006; Garrod et al. 2008), additional nondiffusive processes efficient at low temperatures on dust grains have been added to these models in order to account for the early formation of COMs in star-forming regions (e.g., Jin & Garrod 2020; Garrod et al. 2022).

One of our goals is to obtain observational constraints on the chemical and physical processes that play a role in the emergence of molecular complexity in star-forming regions. One way to achieve this is to test chemical models such as those mentioned above by comparing their predictions to the post-desorption chemical composition of the gas phase of hot molecular cores in star-forming regions.

Hot cores correspond to the dense and compact regions around nascent high-mass protostars where spectral emission from COMs is detected, which is associated with molecular desorption from dust grains upon heating by the central object (e.g., Walmsley et al. 1995). They have typical kinetic temperatures of 150−200 K. Their lifetime is on the order of 6 × 104 yr (e.g., Bonfand et al. 2017; Nony et al. 2024). In low mass star-forming regions, these objects are called hot corinos (Ceccarelli 2004). The chemical composition of hot cores can be determined by analyzing their spectrum in the millimeter/submillimeter wavelength range where molecules have most of their rotational transitions. With the advent of broadband instrumentation at single-dish telescopes and interferometers over the past two decades, spectral line surveys of multiple sources covering a broad frequency range have become more affordable (see, e.g., the review by Jørgensen et al. 2020). This means that the predictions of chemical models can be tested by determining the abundances of a large number of COMs for sample of sources that are located in different environments, have various masses and luminosities, or are in different evolutionary stages.

This work focuses on hot cores embedded in Sagittarius (Sgr) B2. This molecular cloud complex is located in the central molecular zone of our Galaxy, at a projected distance of about 100 pc from the galactic center, Sgr A*. Sgr B2 consists of several protoclusters, in particular the high mass star-forming region Sgr B2(N) that contains a population of hot cores and more evolved H II regions (e.g., Gaume et al. 1995; Schmiedeke et al. 2016; Bonfand et al. 2017; Sánchez-Monge et al. 2017). Since the 1970s, thanks to its high H2 column density, Sgr B2 has been a place of choice to search for interstellar molecules, in particular COMs (e.g., Ball et al. 1970; Rubin et al. 1971; Hollis et al. 2000; Belloche et al. 2008, 2022). It has been the target of multiple spectral line surveys, starting with single-dish telescopes in the 1980s and 1990s (e.g., Cummins et al. 1986; Turner 1989; Nummelin et al. 1998). Our own spectral line survey of Sgr B2(N) with the IRAM 30 m telescope at 3 mm, 2 mm, and 1.3 mm led to the detection of several new COMs (Belloche et al. 2008, 2009). Despite the limited angular resolution (halfpower beam width, HPBW, ∼25′′), we were able to identify the presence of two hot cores in Sgr B2(N) within the single-dish beam thanks to their different systemic velocities (Belloche et al. 2013). Our analysis of the IRAM 30 m spectrum showed that the spectral confusion, exacerbated by the presence of these two velocity components, was much higher at 1.3 mm than at 3 mm and that the opacity of the dust made it more difficult to derive reliable COM column densities at 1.3 mm in this source.

In order to make further progress with the determination of the chemical composition of Sgr B2(N)’s hot cores, we took advantage of the advent of the Atacama Large Millimeter/submillimeter Array (ALMA) in the early 2010s to perform an imaging spectral line survey of Sgr B 2(N) at 3 mm. Thanks to its angular resolution of ∼1.6′′, this survey, called Exploring Molecular Complexity with ALMA (EMoCA), resolved the two main hot cores into separate sources which we called Sgr B2(N1) and Sgr B2(N2) (Belloche et al. 2016). This gain in angular resolution reduced the spectral confusion considerably. The resolution was however still insufficient to resolve the internal structure of both hot cores, and the high dust optical depth of the main hot core, Sgr B2(N1), still made it difficult to explore its chemical composition even at 3 mm. This motivated a second high-sensitivity, high-angular-resolution (∼0.5′′) survey with ALMA, called Re-exploring Molecular Complexity with ALMA (ReMoCA) (Belloche et al. 2019). This survey allowed us to resolve the thermal structure of Sgr B2(N1), probe the transition between nonthermal and thermal desorption, and establish that a number of COMs co-desorb thermally with water in this source (Busch et al. 2022). Our series of spectral line surveys of Sgr B2(N) has demonstrated that a high angular resolution is key to reduce the spectral confusion of high mass star-forming regions by probing individual sources and pockets of gas with narrower line widths (Belloche et al. 2022).

In this article, we focus our attention on the secondary hot core Sgr B2(N2), which ReMoCA resolves into several sources, and we derive their individual chemical compositions. Möller et al. (2025) recently used ALMA at 1.2 mm to study the whole population of hot cores in Sgr B2(N) and Sgr B2(M). They derived the chemical composition of several dozen sources, albeit on the basis of a limited number of molecules, and focused their analysis on evaluating the evolutionary stages of the sources. Thanks to the lower degree of spectral confusion at 3 mm, we report here the detection of a much larger set of molecules and investigate the constraints that these observational results impose on chemical models. This article is structured as follows. We describe the observational setup and the method used to identify the molecules and measure their column densities in Sect. 2. Section 3 presents the results of this analysis and a comparison to other sources from the literature. In Sect. 4, we compare the observational results with the predictions of chemical models. We discuss the results obtained in this work in Sect. 5. Our conclusions are reported in Sect. 6.

2 Observations and radiative transfer modeling

2.1 ALMA observations

The imaging spectral line survey ReMoCA was performed with ALMA toward Sgr B2(N) at high angular resolution in the 3 mm atmospheric window between 84.1 and 114.4 GHz. The observations and data reduction were described in detail in Belloche et al. (2019). We summarize here the main features. The survey was performed in five parts that were tuned to different frequencies with a spectral resolution of 488.3 kHz(1.7−1.3 km s−1). We called these five spectral setups S1–S5. Each setup consists of four 1.88−GHz wide spectral windows labeled W0-W3 that together cover both sidebands. The five setups were observed on different days with various antenna configurations (see Table 1 of Belloche et al. 2019), which resulted in different angular resolutions. The spectral coverage, angular resolution, and noise level of each spectral window are listed in Table 1, which is slightly updated compared to the initial table published in Belloche et al. (2019). The survey has a median angular resolution of 0.6′′, which corresponds to 4900 au at the distance of Sgr B2 (8.2 kpc, Reid et al. 2019). The median rms noise level is 0.8 mJy beam−1, which corresponds to a brightness temperature noise level of 0.27 K for a half-power beam width (HPBW) of 0.6′′ at a frequency of 100 GHz.

Table 1

Beam sizes and noise levels.

Because of the high spectral line density, which translates into spectra close to the confusion limit, and the presence of several sources with different systemic velocities in the field of view, splitting the line and continuum emissions in the Fourier plane was not possible. This was done instead in the image plane as described in Belloche et al. (2019). This splitting was subsequently slightly improved, as we reported in Melosso et al. (2020a).

The primary beam of the ALMA 12 m antennas has a HPBW that varies between 69′′ at 84 GHz and 51′′ at 114 GHz (Remijan et al. 2015). The data cubes of the ReMoCA survey were not corrected for the primary beam response. However, this correction was applied to all spectra used to derive column densities and those displayed in the figures of this article. Figures showing maps are not corrected for the primary beam response. This has the advantage of keeping the noise level uniform. In these maps, the longest (shortest) angular distance to the phase center is 6.8′′(0.5′′), which corresponds to a correction factor of 1.05 (1.0003) at 84 GHz and 1.03 (1.0001) at 114 GHz. The distortion of the maps is thus marginal (smaller than 5%).

2.2 Radiative transfer modeling

We modeled the ReMoCA spectra under the assumption of local thermodynamic equilibrium (LTE), which is well justified given the high densities that characterize the hot cores in Sgr B2(N) (>1 × 107 cm−3, Bonfand et al. 2019). We used the astronomical software Weeds (Maret et al. 2011), which is part of the CLASS program of the GILDAS package1, to compute synthetic spectra. Weeds accounts for the finite angular resolution of the observations, the line optical depth, and the contribution of the background continuum emission to the equation of radiative transfer (see the Weeds documentation2 for more details). We emphasize that neglecting the strong continuum emission of hot cores in the equation of radiative transfer would lead to underestimating the molecular column densities.

We modeled the contribution of each molecule separately before adding them together. The synthetic spectrum of each species is defined by five parameters: size of the emitting region (θs), column density (N), temperature (Trot), line width (ΔV), and velocity offset (Voff) with respect to the assumed systemic velocity of the source. Many molecules contribute to the ReMoCA spectra with rotational transitions that originate not only from their vibrational ground state but also from various vibrationally excited states. In those cases, we modeled each vibrational state separately in order to be able to account for potential differences between vibrational and rotational temperatures.

When the vibrational and rotational temperatures of a given molecule are equal, the column density parameters N of all vibrational states of that molecule are identical in our Weeds models. This value corresponds to the total column density of the molecule, that is it includes the populations of the ground state and all vibrationally excited states. On the contrary, if the vibrational temperature of a given vibrational state V is higher than the rotational temperature because of, for example, radiative pumping, then, in order to fit the observed spectrum, the column density parameter NV of that state would have to be artificially increased to account for the overpopulation in that state. In such a case, the calculation of the total column density of the molecule would be tricky. In practice, we did not find any difference between the vibrational and rotational temperatures of the detected molecules.

We decided to keep the original angular resolutions of the 20 ReMoCA spectral windows despite their differences rather than degrade the data to the coarsest angular resolution (0.8′′). The reason for this is that the shortest separation between the sources studied in this article is 0.55′′. Such a small angle is resolved by setups S4 and S5 but not completely by the other setups. The LTE parameters derived for each source were optimized as much as possible on S4 and S5. Part of the discrepancies between the synthetic and observed spectra of setups S1–S3 result from their coarser angular resolution that leads to partial contamination between neighboring sources, an effect that cannot be accounted for by our simple modeling framework.

2.3 Spectroscopy

Table C.1 provides the list of the 448 spectroscopic entries that were used in this work to model the observed spectra or derive upper limits to the column densities of nondetected molecules. These spectroscopic entries were downloaded from the Cologne Database for Molecular Spectroscopy (CDMS3, Müller et al. 2001, 2005; Endres et al. 2016), the Lille Spectroscopic Database (LSD4), and the database for molecular spectroscopy of the Jet Propulsion Laboratory (JPL5, Pickett et al. 1998). They were all inserted into a local database connected to the radiative transfer software Weeds. In addition to these publicly available entries, our database also includes predictions that were provided to us by various collaborators over the past two decades and a few entries that we prepared ourselves on the basis of published studies.

Spectroscopic entries in all three public databases mentioned above have a specific five- or six-digit index called the tag. The third-to-last digit of the tag codes for the database (0 for JPL, 5 for CDMS, and 8 for LSD). The two or three digits that precede the database code represent the molecular weight of the molecule (e.g., 32 for CH3 OH or 103 for c-C6 H5CN). The last two digits distinguish molecules in the same database that have the same molecular weight. For a given molecule, we modeled the emission of vibrational excited states separately in order to account for potential deviations from LTE in the event that the vibrational temperature would be different from the rotational one (see Sect. 2.2). To do that, we had to split the original entries that contain several vibrational states in the public databases into multiple entries containing one state each, and we assigned new tags to the individual entries, with a database code different from 0 or 5. Unfortunately, we had used the database code 8 for some of our local entries before the LSD database was created. We have started to change these tags when new LSD entries are in conflict with our existing local entries. Some CDMS entries are provided with hyperfine structure in the documentation area of the CDMS website with no specific tag. Here again, we had to assign new tags to these CDMS entries with hyperfine structure.

We made a considerable effort to estimate a posteriori corrections to the column densities of the analyzed molecules when their partition functions did not include, or only partially, the contribution of vibrationally excited states or higher energy conformers. In most cases, the vibrational correction (Cvib) was computed in the harmonic approximation using the following equation: Cvib=i=1N11eEi/kT,Mathematical equation: $C_{\text{vib}}=\prod_{i=1}^{N} \frac{1}{1-e^{-E_{i} / k T}},$(1)

with k the Boltzmann constant, T the temperature, and Ei the energies of the N fundamental modes of vibration. In most cases, the conformational correction (Cconf) was computed assuming that a conformer j of energy Ej with respect to the energy of the lowest-energy conformer contributes to the total partition function of the molecule with a population equal to eEj/kT times the population of the lowest-energy conformer. We also accounted for the degeneracies of the conformers when they differed. For each spectroscopic entry, Table C.1 provides the chemical formula of the molecule, along with the vibrational state or name of the conformer when relevant, the tag, the database where it was downloaded from, and a list of relevant references (“Spectro” for the frequencies and dipole moment, “Cvib” for the vibrational energies, and “Cconf” for the energies of the conformers).

3 Results

3.1 Source selection and description

In this work, we focus our analysis on four positions located in Sgr B2(N2), the secondary hot core of Sgr B2(N). The first three positions correspond to the continuum sources AN06, AN03, and AN02 extracted by Sánchez-Monge et al. (2017) from their observations of Sgr B2(N) performed with ALMA in the 1.2 mm atmospheric window with an angular resolution of 0.4′′. These three sources represent all the continuum sources within Sgr B2(N2) in the list of Sánchez-Monge et al. (2017). Their spectral index between 211 and 275 GHz was found to be 2.4, 3.2, and 3.1, respectively, suggesting that their 1.2 mm ALMA continuum emission is dominated by dust emission. In addition to these three sources, we also analyzed the position called Sgr B2(N2b) where we recently detected normalpropanol and iso-propanol (Belloche et al. 2022). The spectral index of this position in the 1.2 mm spectral window is approximately 3 according to Fig. 7 of Sánchez-Monge et al. (2017). This suggests that the 1.2 mm ALMA continuum emission of N2b is dominated by dust as well. The coordinates and 1.2 mm spectral indices of all selected positions are listed in Table 2. We describe below the additional information that we derived from the ReMoCA survey on the physical structure or stage of evolution of these four sources.

Table 2

Properties of the positions analyzed within Sgr B2(N2).

3.1.1 AN03, a HCH II region

AN03 coincides with the hypercompact (HC) H II region K7 that was identified by De Pree et al. (2015) in Very Large Array (VLA) continuum observations at 7 mm. K7 is located at a J2000 equatorial position of 17h47m 19s.895, −2822 13′′.47, with uncertainties of 0.01 s and 0.1′′, and has a deconvolved full-width-at-half-maximum (FWHM) size of ∼0.08′′ at 7 mm (660 au at the distance of Sgr B2). Several recombination lines of hydrogen were covered by the ReMoCA survey and detected toward Sgr B2(N2). Figure 1 shows the ReMoCA spectrum at the frequencies of the H40α and H41α recombination lines toward the four positions of Sgr B2(N2) analyzed in this article. We selected these lines because they suffer little from contamination by rotational emission of molecules (see the molecular contribution predicted by our LTE model overlaid in blue in Fig. 1). In order to investigate the spatial morphology of the emission of these recombination lines in Sgr B2(N2), we integrated each line over a fixed velocity range, as highlighted in dark gray in Fig. 1. The resulting integrated intensity maps are shown in Fig. 2. The H40α line was independently covered by two setups with differing angular resolutions (S1W3 and S5W0). Figures 1 and 2 show the spectra and maps for both setups in order to illustrate the impact of the angular resolution.

The maps of Fig. 2 reveal a clear peak of H α emission at a J2000 equatorial position of 17h47m 199s.899, −2822 13′′.59 with an uncertainty of ∼0.1′′, as measured in the S5W0 map which has the highest angular resolution (HPBW 0.43′′ × 0.30′′). This H α emission peak coincides within the uncertainties with the positions of the 1.2 mm continuum source AN03 and the 7 mm HCH II region K7. Figure 2 also displays the positions of the water and methanol masers observed with the VLA by McGrath et al. (2004) and Lu et al. (2019), respectively. AN03 seems to be associated with the easternmost water maser which has a systemic velocity of 75.3 km s−1. Their angular separation is on the order of 0.15′′ (∼1200 au).

The recombination line maps of Fig. 2 also show an arc of extended emission that corresponds to K5, a large shell of ionized gas seen in the VLA continuum maps of Gaume et al. (1995) at 1.3 cm. The fact that the morphology of this extended structure is very similar in all three maps of Fig. 2 gives us confidence that the maps suffer little from contamination by molecular emission.

Thumbnail: Fig. 1 Refer to the following caption and surrounding text. Fig. 1

ALMA continuum-subtracted spectra toward the hot core positions AN06, AN03, AN02, and N2b (from left to right) at the frequencies of the hydrogen recombination lines H40α (top and middle rows, for two setups with different angular resolutions) and H41α (bottom row). In each panel, the vertical dashed line marks the systemic velocity of the source adopted for the LTE modeling of the molecular emission. The velocity axis refers to the rest frequency of the recombination line. The horizontal dashed line indicates the 3σ noise level. The blue spectrum represents the LTE model that includes the contribution of all molecules identified so far. The velocity range highlighted in dark gray is specific to each recombination line and represents the range of channels selected to compute the integrated intensity map of each recombination line shown in Fig. 2. The setup and spectral window are indicated in the left panel of each row along with the corresponding HPBW.

Thumbnail: Fig. 2 Refer to the following caption and surrounding text. Fig. 2

ALMA integrated intensity maps of the hydrogen recombination lines H41α (panel a) and H40α (panels b and c, for two setups with different angular resolutions). The intensity was integrated over the velocity range highlighted in dark gray in Fig. 1 in order to avoid contamination by molecular lines. The violet plus symbols mark the hot core positions AN06, AN03, AN02, and N2b. The green cross indicates the VLA position of the HCH II region K7 at 7 mm from De Pree et al. (2015). The blue triangles and violet square indicate the VLA positions of the water and Class II methanol masers reported by McGrath et al. (2004) and Lu et al. (2019), respectively. The setup and spectral window numbers are given in the bottom right corner of each panel along with the associated beam size (HPBW). The values of the noise level, σ, are 18, 21, and 15 mJy beam −1 km s−1, respectively. The contours start at 3σ (brown contour) and then increase by a factor of two at each step (black contours). The dotted blue contour, when present, shows the −3σ level.

Thumbnail: Fig. 3 Refer to the following caption and surrounding text. Fig. 3

ALMA continuum emission maps at 92.3, 98.9, and 99.6 GHz. The setup and spectral window numbers are given in the bottom right corner of each panel along with the associated beam size (HPBW). The symbols are the same as in Fig. 2. The values of the noise level, σ, are 0.73, 0.97, and 0.42 mJy beam −1, respectively. The contours start at 5σ (brown contour) and then increase by a factor of two at each step (black contours). The additional violet contour in panel c is at 60σ. It was added to emphasize the emission peak on AN02. Dotted blue contours, when present, indicate the −5σ and −10σ levels.

3.1.2 AN06, an isolated dust continuum source

Figure 3 shows the maps of continuum emission measured in the same ReMoCA spectral windows (S3W0, S1W3, and S5W0) as the recombination lines displayed in Fig. 2. These continuum emission maps trace a similar extended structure as the Hα maps of Fig. 2, indicating that the extended continuum emission detected at 3 mm with ReMoCA contains a significant contribution from thermal free-free emission of ionized gas. However, the contours of the highest angular resolution map (S5W0, Fig. 3c) clearly reveal 3 mm continuum peaks associated with the 1.2 mm continuum sources AN02 and AN06. This indicates that the compact 3 mm continuum emission around each source has a significant contribution from dust emission. AN06 is thus a bona fide dust continuum source that is not associated with any HCH II region because it has no compact counterpart in the ReMoCA maps of Hα emission and in the VLA maps of thermal free-free emission of De Pree et al. (2015).

The ReMoCA spectra of AN06 reveal for most molecules that we analyzed two velocity components that are clearly distinct. This is illustrated in Fig. 4 with two rotational transitions of ethanol that were carefully selected for their absence of contamination by other molecules. We divided the ethanol emission into three velocity ranges highlighted with the three darkest shades of gray in Fig. 4. The integrated intensity maps corresponding to these three velocity ranges are shown in Fig. 5. While the maps in the velocity ranges 71−77 km s−1 and 77−81 km s−1 do not show any particular structure toward AN06 (Figs. 5b, c, e, and f), this source is clearly associated with a compact emission peak in the ethanol maps integrated from 65 to 71 km s−1 (Figs. 5a and d). Therefore, we think that the dust continuum emission of AN06 is associated with the lower-velocity molecular component and we adopted a systemic velocity of 68.3 km s−1 for this source. The higher velocity component detected toward AN06, at a systemic velocity of about 73.5 km s−1 and with lower peak intensities than the lower velocity component, traces the edge of the Sgr B2(N2) dense core and is likely not related to AN06 itself (Figs. 5b and e). Both components have intrinsic line widths of about 3.5 km s−1. We conclude that AN06 is a bona fide continuum source that is likely not embedded in Sgr B2(N2) but rather lies in either the foreground or background of Sgr B2(N2).

3.1.3 AN02, a dust continuum source with two velocity components, possibly driving an outflow

AN02 is associated with a 3 mm continuum peak in Fig. 3c and has no compact counterpart in Hα emission at 3 mm (Fig. 2) or thermal free-free emission at 7 mm (De Pree et al. 2015). It thus represents a bona fide dust continuum source. Like for AN06, our LTE modeling of the molecular emission of AN02 required two velocity components to fit the asymmetric shape of the lines (see, e.g., Figs. 4c and g). However, with a velocity difference of only ∼4 km s−1 and intrinsic line widths of about 4 and 3 km s−1, these two components are not as well separated as in the case of AN06. They have similar peak temperatures and both have a local maximum at the position of AN02 in the integrated intensity maps shown in Fig. 5. Therefore, the data do not tell us whether the continuum source AN02 is associated with both components or only one of them. We assumed a systemic velocity of 72.9 km s−1 for AN02, which corresponds to the lower velocity component, but modeled both components together with velocity offsets of ∼0 and ∼4 km s−1.

Four water masers lie in the vicinity of AN02 toward the northwest, at angular distances of about 0.25′′ and 0.6–0.7′′. They are roughly aligned with AN02, suggesting that AN02 may drive an outflow that produces these masers. The easternmost maser, which we associated with AN03 in Sect. 3.1.1, is located at a position that is roughly symmetric to the position of the westernmost masers with respect to AN02. Therefore, the easternmost maser may equally well be associated with AN02 rather than AN03. The three westernmost masers have velocities of 75.3, 74.0, and 21.2 km s−1 (McGrath et al. 2004). The former velocities are both in between the velocities of the two molecular components associated with AN02 while the latter velocity is blueshifted by about −50 km s−1. The closest maser to AN02 has a velocity of 56.5 km s−1 that is blueshifted by about −16 km s−1.

A Class II methanol maser was detected by Lu et al. (2019) with the VLA in the vicinity of Sgr B2(N2) (see violet square in Fig. 3). It is located about 1′′ to the northwest of AN02. Mills et al. (2018) reported an astrometry offset affecting their VLA A-configuration data due to an error in the position of their phase calibrator in the VLA catalog. Such an offset may also affect the data set of Lu et al. (2019) which was calibrated with the same phase calibrator (X. Lu, priv. comm.), but the offset does not seem to be along the northwest direction (A. Ginsburg, priv. comm.). Therefore the source associated with the methanol maser is currently unknown.

We conclude that AN02 is a dust continuum source with two velocity components of molecular emission. This source is possibly driving an outflow along the northwest-southeast direction.

Thumbnail: Fig. 4 Refer to the following caption and surrounding text. Fig. 4

ALMA continuum-subtracted spectra toward the hot core positions AN06, AN03, AN02, and N2b (from left to right) at the frequencies of two contamination-free rotational transitions of ethanol indicated on the right along with the energy of the upper level in temperature unit. In each panel, the vertical dashed line marks the systemic velocity of the source adopted for the LTE modeling of the molecular emission. The velocity axis refers to the rest frequency of the ethanol transition. The horizontal dashed line indicates the 3σ noise level. The blue spectrum represents the LTE model that includes the contribution of all molecules identified so far. The red spectrum represents the LTE model of ethanol only. In the case of AN06, the synthetic ethanol spectra of the first and second velocity components are displayed in orange and red, respectively. The velocity ranges highlighted in the three darkest shades of gray are specific to each ethanol transition (but common to all sources, albeit with slight differences due to the finite spectral sampling) and represent the ranges of channels selected to compute the integrated intensity maps of each transition shown in Fig. 5. The setup and spectral window are indicated in the left panel of each row along with the corresponding HPBW.

Thumbnail: Fig. 5 Refer to the following caption and surrounding text. Fig. 5

ALMA integrated intensity maps of two contamination-free rotational transitions of ethanol with rest frequency and upper level energy indicated on the right of each row. The intensities were integrated over the velocity ranges highlighted in the three darkest shades of gray in Fig. 4. The integration range is indicated in the top left corner of each panel. The setup and spectral window numbers are given in the bottom right corner along with the associated beam size (HPBW). The symbols are the same as in Fig. 2. The noise level, σ, is indicated below the velocity range. The black contours start at 3σ and then increase by a factor of two at each step, except for the last contour of panel e which is at 320σ. The dotted blue contour, when present, indicates the −3σ level.

Thumbnail: Fig. 6 Refer to the following caption and surrounding text. Fig. 6

Maps of centroid velocity (left column) and line width (right column) of two contamination-free rotational transitions of ethanol with rest frequency and upper level energy indicated in the top left corner of each row. The kinematic information is plotted only for pixels with a signal-to-noise ratio in peak intensity higher than 10. The setup and window numbers are given in the bottom right corner of each panel along with the associated beam size (HPBW). The symbols are the same as in Fig. 2. These two transitions of ethanol have opacities on the order of 2 in N2b, which means that the intrinsic line widths around N2b are somewhat smaller than the plotted values.

3.1.4 N2b, a dust continuum source with narrow lines

Figure 6 shows maps of centroid velocity and line width obtained by fitting a Gaussian function to the spectra of both transitions of ethanol used in Figs. 4 and 5. Figures 6b and d reveal that the line width decreases toward the southwest of Sgr B2(N2), reaching values as low as ∼2.7 km s−1. Such narrow lines are advantageous for the search for new molecules but the narrowest line widths are reached at the edge of the detected molecular emission where the emission becomes faint (see Fig. 5). As a compromise between having narrow line widths and keeping high intensities, we selected position N2b for our previous work on propanol (Belloche et al. 2022). This position has intrinsic line widths on the order of 3.5 km s−1 and we adopted a systemic velocity of 74.2 km s−1. Figures 5b and e show that N2b is associated with a peak of integrated intensity of ethanol. It also corresponds to an extension of 3 mm continuum emission in Fig. 3c, which shows the map at highest angular resolution, while no such clear extension is visible in the H40α map of Fig. 2. This suggests that N2b may be associated with a compact dust core like the other three sources, even if it was not extracted as such by Sánchez-Monge et al. (2017). Their ALMA dust continuum map at 242 GHz also shows an extension at the position of N2b (see their Fig. 4).

Figure 6 shows other features that are consistent with the source descriptions provided in the previous sections. The blue area toward AN06 in panels a and c reveals the low velocity component that we exclusively assigned to AN06 and that dominates its emission. The overlap along the line of sight between this velocity component and the dense core of Sgr B2(N2) results in the sharp line width increase (∼8 km s−1) at the location of the red band that touches AN06 in panels b and d. The existence of two velocity components toward AN02 also results in a local increase in the line width (∼7 km s−1) and a slightly higher average centroid velocity (∼75 km s−1) than the rest of Sgr B2(N2), which has a systemic velocity of ∼74 km s−1. Finally, an unresolved pocket of gas with a systemic velocity of ∼75 km s−1 stands out at a bit less than 1′′ to the south of AN03. This region is not associated with any specific structure in the Hα and continuum maps shown in Figs. 2 and 3, therefore it is difficult to assess its nature, but it seems to be associated with a filamentary structure that is connected to AN02 in the ethanol maps integrated over the 77−81 km s−1 velocity range (Figs. 5c and f).

3.2 Chemical composition

The LTE parameters of all the molecules included in the model of at least one of the four positions presented in Sect. 3.1, obtained following the method described in Sect. 2.2, are listed in Tables E.1E.5. For a given catalog entry, each table indicates the number of spectral lines that we considered as detected, namely lines with a peak temperature above the 3σ noise level and for which the LTE model of this entry accounts for at least about two thirds of the detected signal. Exceptions to the latter criterion occur for the setups with the largest beams (e.g., S1) for which the LTE models, which were optimized for the setups with the highest angular resolution (S4 and S5), sometimes clearly underestimate the emission in the larger beams due to the nonuniformity of the molecular emission across the Sgr B2(N2) region. The excellent match between the LTE models and the ReMoCA spectra of N2b, AN02, AN03, and AN06 is illustrated in Fig. D.1 that also displays the ReMoCA spectrum of the ultracompact H II region K4 for comparison. K4 is located about 9′′ to the north of Sgr B2(N2).

Species that are considered as securely detected are labeled with a “d” in the third column of Tables E.1E.5. The number of lines required to qualify a catalog entry as detected depends on the signal-to-noise ratio of the lines and whether or not other related entries (e.g., other isotopologs or vibrational states) are detected with consistent parameters. Entries that are considered as tentatively detected with only a few lines are labeled with a “t.” Entries that are included in the LTE model because their predicted lines have peak temperatures close to or above the 3σ level but are heavily contaminated by emission or absorption of other species, and therefore cannot be considered as securely or tentatively identified, are labeled with a “c” to indicate that they merely contribute to the detected signal. They correspond typically to higher-energy vibrational states that are modeled with the same parameters as the lower-energy ones, or rarer isotopologs that are modeled with the parameters of the main species assuming a typical isotopic ratio. Species that are not detected are labeled with a “n” and the tables provide upper limits to their column densities.

Table 3 provides a concise overview of the detections obtained with the ReMoCA survey toward N2b, AN02, AN03, AN06, and AN06c2. This table contains 61 molecules. Thanks to the sensitivity of ALMA and the fact that the spectral confusion at 3 mm is not too severe, we were able to identify between 42 and 58 molecules (counting only the main isotopologs). Among these molecules, between 3 and 10 are detected only tentatively and between 11 and 14 are seen only in absorption. Furthermore, many COMs are detected: between 22 and 24 toward N2b, AN02, and AN03, and about 15 toward AN06 and AN06c2.

In addition to the molecules included in the LTE model of at least one of the four positions, we provide upper limits to the column density of 146 relevant molecules that were searched for in the ReMoCA survey but not detected (Tables E.6E.10). The upper limits were obtained assuming the same LTE parameters as related molecules that were detected or typical LTE parameters of the respective source, keeping only the column density as a free parameter. The molecules are grouped by classes (O-bearing, O+N-bearing, N-bearing, S-bearing, S+O-bearing, S+N-bearing, Cl-bearing, P-bearing, and hydrocarbons) and, within each class, by increasing entry index (tag) in our local Weeds database.

The column densities of 46 molecules extracted from Tables E.1E.5 and normalized to the column density of methanol are displayed in Fig. 7. We normalized to methanol because it is the most abundant detected COM and because deriving reliable H2 column densities is challenging. In order to compare the chemical compositions of the selected positions more easily, we show their abundances (relative to methanol) normalized to the abundances of N2b in Fig. F.1.

Table 3

List of molecules detected toward N2b, AN02, AN03, AN06, and AN06c2 with the ReMoCA survey.

Thumbnail: Fig. 7 Refer to the following caption and surrounding text. Fig. 7

Column densities derived toward N2b, AN02, AN03, and AN06 with our LTE modeling, normalized to the column density of methanol (see Tables E.1E.5). The column density of methanol is indicated in the top left corner of each panel. The chemical compositions of the two velocity components detected toward AN06 are displayed separately (AN06 in the top panel and the second velocity component in the panel labeled AN06c2). The panel of AN02 reports the sum of the column densities of its two velocity components. Hatched bars show tentative detections while empty bars with downward arrows indicate upper limits. Molecules with status “c” in Tables E.1E.5 are here represented as upper limits. The dashed lines indicate levels of 1% and 0.01% with respect to methanol.

3.3 Correlations of chemical composition between positions

Figure 8 shows correlation plots of the column densities normalized to methanol for various pairs of positions. The best correlation occurs between the chemical compositions of AN02 and AN03 (Fig. 8e). Most data points deviate by much less than a factor of 2 from the 1:1 relation. The only data points that deviate by slightly more than a factor of 2 are HC3N, HC5N (upper limit toward AN02), and SO2. C2H3CN, SO, and HNCS deviate by slightly less than a factor of 2. All six molecules are less abundant toward AN02 than toward AN03. The upper limits derived for HCCNC, HSCN, and E-HNCHCN toward AN02 and for n- and i−C3H7OH, and C2H5OCHO, toward both AN02 and AN03 are not stringent enough to tell whether or not these molecules deviate from the overall correlation.

The next pair of positions that show a high degree of correlation are the two velocity components of AN06 (Fig. 8h). In this case, however, the abundances (relative to methanol) of the second velocity component (AN06c2), which traces the gas at the edge of the dense core that contains AN03 and AN02, are systematically higher than the abundances of the main component (AN06), by roughly a factor of 2. There is one notable exception to the overall correlation: NH2CHO is remarkably underabundant toward AN06, with an abundance upper limit that is nearly two orders of magnitude below the abundance measured toward AN06c2. HNCO is offset on the same side of the overall correlation, but only by a factor of about 4, while HC3N, SO2, and CH3CCH are offset by a factor slightly larger than 2.

In contrast to the previous two pairs, the chemical composition of N2b correlates poorly with that of the other positions (Figs. F.1 and 8a−d). Most molecules have abundances relative to methanol distributed between the 1:1 and 10:1 lines in Figs. 8a and b, the abundances toward AN02 and AN03 being higher than toward N2b. The most extreme deviations from the 1:1 relation, by more than a factor of 20, occur for t-HCOOH (upper limit toward N2b) and NH2CN.

The correlation between N2b and AN06 is even poorer, with a spread of nearly two orders of magnitude (Figs. F.1a and 8d). Most species are less abundant (relative to methanol) in AN06 than in N2b. The most extreme cases are HNCO and NH2CHO which are almost two orders of magnitude (or more) less abundant, followed by the cyanides CH3CN, C2H5CN, C2H3CN, nand i-C3H7CN, and HC3N (about one order of magnitude). The only molecules that are more abundant in AN06 than in N2b by at least a factor of two are C2H5OCHO and CH3NH2.

AN06c2 correlates a bit better with N2b than AN06 (Figs. F.1b and 8c). The only species that deviate from the 1:1 relation by significantly more than a factor of 2 are CH3 NH2 (more abundant in AN06c2) as well as HNCO and C2 H5CN (more abundant in N2b).

Finally the chemical composition of both components of AN06 is poorly correlated with that of AN03 (Figs. 8f and g). All molecules are less abundant in AN06 than in AN03, by up to 2–3 orders of magnitude for HC3N, C2H3CN, HNCO, and NH2CHO. The situation is less extreme for AN06c2, with a handful of molecules close to the 1:1 relation, most molecules between the 1:1 and 1:10 relations, and only HC3N, C2H3CN, and HNCO being 15–30 times less abundant than in AN03.

In order to investigate if the lack of correlation between certain sources is related to systematic differences between classes of molecules, we show in Fig. F.2 the same correlation plots as in Fig. 8 but with colors coding for the atomic composition of the molecules: O-bearing in black, O+N-bearing in blue, N-bearing in red, pure hydrocarbon in green, and S-bearing in yellow. We notice that the correlations of the pairs N2b/AN02, N2b/AN06c2, N2b/AN06, and AN06/AN06c2 (Figs. F.2a, c, d, and h) are much tighter for the S-bearing molecules than for the full sample of molecules. This is not the case for the pairs of sources that involve AN03 (Figs. F.2b, e, f, and g). Furthermore, the class of O+N-bearing molecules stands out in Figs. F.2d and g: these molecules are located in the lower part of the correlation plot in both cases, revealing that they are underabundant in AN06 compared to N2b and AN03 (and AN02, given that it correlates well with AN03). In contrast to the S-bearing and O+N-bearing species, the O-bearing and N-bearing molecules do not show any obvious systematic differences between the five sources.

To summarize, the pairs of positions in Sgr B2(N2) that show a good correlation of their chemical composition (relative to methanol) are AN02/AN03 and AN06/AN06c2. The molecules that stand out in at least one panel of Fig. 8 are HC3N, HC5N, C2H3CN, C2H5CN, NH2CN, CH3NH2, NH2CHO, HNCO, t-HCOOH, C2 H5OCHO, and SO2. Finally, AN03 stands out for its S-bearing molecular content, which poorly correlates with the other sources, and AN06 for its underabundant O+N-bearing species compared to N2b, AN03, and AN02 (Fig. F.2).

Thumbnail: Fig. 8 Refer to the following caption and surrounding text. Fig. 8

Correlation plots of column densities normalized to methanol for various pairs of positions. The x- and y-axes of each panel correspond to the positions written in the bottom right and top left corners, respectively. The color coding of the molecules is the same as in Figs. 7 and F.1 and is indicated on the right. Latine and Greek letters, reported at the same abscissa as the corresponding molecule, were added to facilitate the identification of the data points. Filled data points represent firm and tentative detections while empty circles with arrows indicate upper limits toward at least one position. The dashed and dotted lines indicate deviations by a factor 10 and 2, respectively.

Thumbnail: Fig. 9 Refer to the following caption and surrounding text. Fig. 9

Abundance of protonated molecules (red) relative to their neutral form (black). Each panel corresponds to one position as labeled in the bottom right corner. Empty squares and downward arrows represent tentative detections and upper limits, respectively. Crosses of the same color represent predictions of the chemical model described in Sect. 4. The molecules are listed along the x axis below and above the bottom and top panels, respectively. The list of protonated molecules is indicated at the bottom right.

Thumbnail: Fig. 10 Refer to the following caption and surrounding text. Fig. 10

Same as Fig. 9, but for radicals (red) with respect to the hydrogenated form (black).

3.4 Correlations between families of molecules

Figures 915 explore the chemical composition of each selected position by plotting abundance ratios of specific pairs of species for different families of molecules.

Thumbnail: Fig. 11 Refer to the following caption and surrounding text. Fig. 11

Same as Fig. 9, but for molecules with a double (red) or triple (green) bond with respect to the saturated form with a single bond (black).

3.4.1 Protonated molecules

We compare in Fig. 9 six protonated molecules to their neutral form. None of these cations are detected toward any position. The most stringent constraint is obtained for formaldehyde toward N2b, with [H2COH+]/[H2CO] < 10−3.

3.4.2 Radicals

Radicals related to formaldehyde, methanol, acetaldehyde, methyl cyanide, cyanoacetylene, and formamide are not detected either (Fig. 10). The most stringent constraint is obtained for methyl cyanide with [CH2CN]/[CH3CN] < 10−3, again toward N2b.

3.4.3 Degree of bond saturation

Figure 11 shows the abundance ratios of saturated (single bond) and unsaturated (double or triple bond) molecules for nine families of molecules. In nearly all cases, the double-bond form is much less abundant than the single-bond one, by up to at least two orders of magnitude for the pairs C2 H5OH/C2H3 OH and n-C3H7CN/E-CH3CHCHCN in N2b. The only clear exception is the pair CH3NH2/CH2NH which has an abundance ratio of about 1 toward all positions, the unsaturated form even slightly dominating over the saturated one in AN02 and, marginally, N2b.

The three molecules C2 H5CN/C2H3CN/HC3N are detected toward all five positions, with the saturated form being the most abundant one by factors of 5 (in AN03) to 50 (in N2b). The relative ratios of the double-bond and triple-bond molecules depend on the position: they have a similar abundance in AN03 but the triple-bond one dominates over the double-bond one by a factor ∼2 in N2b, AN06, and AN06c2, while the opposite is true in AN02, with C2H3CN dominating over HC3N by a factor of 2. At the next stage in complexity in the alkyl cyanide family, the unsaturated forms of n- C3H7CN are less abundant than the saturated one by at least one order of magnitude in N2b, AN02, and AN03. The triple-bond species CH3 C3N is (tentatively) more abundant than the double-bond one, E-CH3CHCHCN, in AN02, which is opposite to what we found for the pair HC3N/C2H3CN in the same source.

Thumbnail: Fig. 12 Refer to the following caption and surrounding text. Fig. 12

Same as Fig. 9, but for molecules with O-bearing functional groups -CHO and -COOH with respect to molecules with a -CN functional group.

3.4.4 CN, CHO, and COOH functional groups

We examine in Fig. 12 how the abundances of five groups of molecules behave for molecules that share the same backbone but terminate with a nitrile functional group (-CN), an aldehyde group (-CHO), or a carboxylic group (-COOH). No obvious pattern is visible. Overall, the CN-bearing molecules are more abundant than their CHO- or COOH-bearing counterparts, with the notable exception of NH2CN which is about two orders of magnitude less abundant than NH2CHO in all three positions where it is detected (N2b, AN02, AN03). While acetaldehyde is at most a factor of 10 less abundant than methyl cyanide, we find that propanal is at least two orders of magnitude less abundant than ethyl cyanide in nearly all positions (the upper limit is inconclusive for AN06c2). In addition to the simple carboxylic acid HCOOH, only one complex carboxylic acid, CH3COOH, is (tentatively) detected toward AN02 and AN03, with an abundance about one order of magnitude lower than CH3CN. If this ratio also holds for the pair NH2 CH2CN/NH2CH2COOH, then glycine would have an abundance at least one order of magnitude lower than its current upper limit.

Thumbnail: Fig. 13 Refer to the following caption and surrounding text. Fig. 13

Same as Fig. 9, but for the reduced form (-CH3. functional group) with respect to the aldehyde form (-CHO functional group).

3.4.5 Oxydized and reduced forms

We show in Fig. 13 how the abundances behave when we go from an oxydized form (aldehyde functional group) to a reduced form (CH3 group). Here again, no obvious pattern is seen. While CH3OH and C2H5OH largely dominate over t-HCOOH and CH2(OH)CHO, respectively, in the sources where t-HCOOH or CH2(OH)CHO is detected, it is the opposite for CH3NH2 which is less abundant than NH2CHO in N2b, AN02, and AN03 (but maybe not in AN06c2). The pair CH3OCH3/CH3OCHO lies in between these two cases, with an abundance ratio close to 1 in all five sources. The upper limits obtained for CH3NHCH3 and C2H5OCH3 suggest that they behave like one of the latter two cases.

3.4.6 Size of carbon backbone

We investigate in Fig. 14 how the abundances vary with the size of the carbon backbone for 20 series of homologous molecules. H2CO, CH3OH, and NH2CHO show a similar pattern: an abundance drop of about one order of magnitude when a molecule has an additional CH2 group at one end of its heavy-atom backbone (CH3CHO, C2H5OH, and CH3NHCHO), and a further drop of about one order of magnitude for a second additional CH2 group (s-C2H5CHO, tentatively detected in ANO3, n-C3H7OH, detected in N2b, and C2H5NHCHO, with a stringent upper limit in AN02). This behavior is also seen for CH3OCHO for the case of one additional CH2 group, either with (tentative) detections or constraining upper limits of C2H5OCHO. We do not have results for the next degree of complexity in this case.

The same behavior is observed for the series of alkyl cyanides C2H5CN, n−C3H7CN (detection), and n-C4H9CN (upper limit), at least in N2b, but CH3CN does not fit into this pattern because it has a somewhat lower abundance than C2H5CN in nearly all five sources. Toward N2b, in addition to the previous cases, the derived upper limits shown in red also indicate an abundance drop of about one order of magnitude (at least) for one additional CH2 group for the following species: CH3OCH3, a-(CH2OH)2, CH2(OH)CHO, CH2CO, CH3NCO, CH3NC, C2H3CN, CH2NH, CH3SH, and CH3CCH. The upper limits are inconclusive for CH3C(O)NH2 and CH3NH2 in N2b, but they do imply an abundance drop by nearly an order of magnitude (at least) in AN02. The case of HC3N is more extreme, with an abundance drop of 2–3 orders of magnitude for CH3C3N in N2b, AN02, and AN03. For the two positions were HCCNC is tentatively detected, the upper limits derived for CH3CCNC are inconclusive.

The only clear exception to the general behavior noted in the previous paragraph is NH2CN, for which the next stage in backbone length, NH2CH2CN, is one order of magnitude more abundant in N2b. However, the latter molecule corresponds to an additional CH2 group between the nitrogen and carbon atoms of NH2CN. We do not have results for the molecule with an additional CH2 group at the end of the backbone (CH3NHCN).

Thumbnail: Fig. 14 Refer to the following caption and surrounding text. Fig. 14

Same as Fig. 9, but for series of molecules with an additional CH2 group in their backbone at each further step in complexity.

3.4.7 Structural isomers

Finally, we examine the relative abundances of structural isomers, that is molecules with the same elemental composition but a different arrangement of their constituent atoms. The results are displayed in Fig. 15 for 17 groups of species with a molecular size ranging from 4 to 12 atoms. Several pairs of isomers stand out for having similar abundances: C2H5OH/CH3OCH3, CH3C(O)NH2/CH3NHCHO, CH3COOH/CH2(OH)CHO (tentatively), n/i-C3H7OH, and n/i-C3H7CN. On the contrary, one isomer (given in parentheses) dominates by more than two orders of magnitude in the following groups: CHNO (HNCO), C2H3N(CH3CN), C3HN(HC3N), and C3H5N(C2H5CN). We note that methyl acetate (CH3C(O)OCH3), which was detected in Orion KL (Tercero et al. 2013, 2018), is missing in the C3H6O2 family shown in Fig. 15. We have a preliminary spectroscopic entry for this molecule but it suffers from several issues that currently prevent us from deriving reliable column densities or upper limits for the Sgr B2(N2) sources.

Focusing only on the detections and tentative detections displayed in Fig. 15, we do not find any significant difference in the isomer ratios between the five positions. These isomer ratios thus appear to be relatively robust with respect to the potential range of physical conditions or age spanned by this small sample of sources.

Thumbnail: Fig. 15 Refer to the following caption and surrounding text. Fig. 15

Same as Fig. 9, but for groups of structural isomers. The elemental composition of the groups is labeled along the x axis. The list of molecules belonging to each group is provided at the bottom of the figure with their respective colors as used in the plots.

3.5 Comparison of chemical composition of Sgr B2(N2) to other sources

Figures 1620 compare the column densities relative to methanol derived in Sect. 3.2 toward the sources of Sgr B2(N2) to the composition of sources that have been studied in detail in the literature: the hot core G31.41+0.31, hereafter G31.41 (GUAPOS interferometric survey, e.g., Mininni et al. 2020; López-Gallifa et al. 2024), the hot corino IRAS 16293–2422 B, hereafter IRAS 16293B (PILS interferometric survey, e.g., Jørgensen et al. 2016; Drozdovskaya et al. 2019), the (possibly shocked) molecular cloud region G+0.693–0.027, hereafter G+0.693 (Yebes and IRAM 30 m single-dish survey, e.g., Rivilla et al. 2022c; Jiménez-Serra et al. 2022)6, and the starless core TMC-1 in the Taurus molecular cloud (mainly the QUIJOTE single-dish survey, e.g., Cernicharo et al. 2021; Agúndez et al. 2025). The column densities of these sources compiled from numerous articles published by the teams of these surveys, as well as a few others in the case of TMC-1 (e.g., Gratier et al. 2016; Tennis et al. 2023), are listed in Table A.1. In the cases where this had not been done, we applied a posteriori vibrational or conformational corrections to the published column densities. We report in Table 4 the Pearson correlation coefficients of the abundance distributions shown inFigs. 1620, for all molecules and per class of molecules. The calculation does not take the upper limits into account and it was done only for samples with at least four items.

Jørgensen et al. (2020) reported a good correlation between the chemical compositions of IRAS 16293B (derived from the PILS survey) and Sgr B2(N2) (obtained from the EMoCA survey at an angular resolution of ∼1.6′′) after normalizing the column densities of the O-bearing and N-bearing molecules to CH3OH and HNCO (or CH3CN), respectively. Given that Sgr B2(N2) consists of several sources that were not resolved by the EMoCA survey, we revisit this correlation by taking advantage of the higher angular resolution of the ReMoCA survey. Figures 16b–18b show a clear differentiation of the four main classes of molecules (O-bearing, O+N-bearing, N-bearing, and S-bearing) between IRAS 16293B and N2b/AN02/AN03. For each class of molecules taken separately, there is a good correlation between IRAS 16293B and N2b/AN02/AN03, but the four groups are shifted with respect to each other: the O-bearing molecules occupy the top part of the distribution, the N-bearing species the bottom part, and the O+N-bearing and S-bearing molecules lie in between. The correlation coefficients listed in Table 4 confirm this visual impression: the coefficients of IRAS 16293B for the individual classes of molecules are systematically larger than for the whole sample of molecules. The tightest correlation is obtained for the S-bearing molecules, followed by the O-bearing ones. The N-bearing molecules are the ones that correlate the least, with correlation coefficients between 0.66 and 0.77. The O+N-bearing species have a coefficient that lies in between, except for N2b where it is slightly smaller than for the N-bearing molecules. With a coefficient larger than 0.8, the overall correlation seems to be stronger for AN06c2 and AN06 than for N2b/AN02/AN03 but this may be biased by the smaller sample of molecules (21 versus 25–28) and the upper limits of several O-bearing, N-bearing, and O+N-bearing species in Figs. 19b and 20b suggest significant deviations from the overall correlation.

Table 4

Correlations between the chemical composition of Sgr B2(N2) sources and sources from the literature.

Among the four sources from the literature displayed inFigs. 1618, the one that shows the best overall correlation with N2b, AN02, and AN03 is the hot core G31.41 (Table 4). One reason for this better overall correlation compared to IRAS 16293B is that the differentiation between the different classes of molecules is less pronounced. Still, like in the case of IRAS 16293B, the N-bearing species occupy the lower part of the distribution in Figs. 16a–18a: they are underabundant by about one order of magnitude in G31.41, especially compared to AN02 and AN03. However, in contrast to IRAS 16293B, the O-bearing molecules are well mixed with the O+N-bearing species and the S-bearing molecules in the overall distribution. Among the four classes of molecules, the tightest correlation occurs for the N-bearing molecules with a coefficient larger than 0.9. The O-bearing molecules are the ones that correlate the least. This is even more true for N2b when we consider the upper limits obtained for t-HCOOH and CH3COOH which are not taken into account in the calculation of the correlation coefficients (labels q and i in Fig. 16a). However, these two molecules do not stand out in the O-bearing class in the cases of ANO2 and AN03 where they are detected. There is also a strong overall correlation between G31.41 and AN06 (Fig. 20a), except for the class of O+N-bearing molecules, with HNCO, NH2CHO, and CH3C(O)NH2 being overabundant in G31.41 by about or more than one order of magnitude.

In contrast to G31.41 and IRAS 16293B, there is no overall correlation between G+0.693 and the Sgr B2(N2) sources, although they are all located in the same molecular cloud complex. The lack of correlation results mainly from the distribution of the N-bearing species which have correlation coefficients between 0.08 and 0.44 (Table 4). For each other class of molecules, we find a much higher degree of correlation, with all correlation coefficients being above 0.75. The S-bearing species are located above the O+N-bearing molecules and O-bearing species in Figs. 16c–20c, meaning that, relative to O+N-bearing molecules and O-bearing species, the S-bearing molecules are more prominent in G+0.693 than in the Sgr B2(N2) sources.

Finally, the situation is even more extreme for the starless core TMC-1. There is absolutely no overall correlation between this source and the Sgr B2(N2) sources, even for the O- or N-bearing molecules taken separately. The only class of molecules that shows some (low) degree of correlation is the S-bearing class but, as indicated by the smaller correlation coefficients (0.6–0.7), the dispersion of the data points in Figs. 16d–20d is larger than for the other three sources from the literature.

To summarize, the class of N-bearing species is the one that reveals the largest variance compared to the other classes of molecules: with respect to the Sgr B2(N2) sources, it shows the tightest correlation for G31.41, a poorer correlation for IRAS 16293B, and no correlation at all for G+0.693 and TMC-1. The N-bearing species are underabundant by 1−2 orders of magnitude in G31.41 and IRAS 16293B with respect to N2b, AN02, and AN03. This is also the case, albeit in a less pronounced manner, for the S-bearing and O+N-bearing molecules. In contrast to the N-bearing species, the class of S-bearing molecules has the smallest variance: it shows a high degree of correlation for G31.41, IRAS 16293B, and G+0.693 with respect to the Sgr B2(N2) sources (and a moderate degree of correlation for TMC-1). This is also the case, albeit to a lesser degree, for the class of O+N-bearing molecules. Finally, the O-bearing species stand out in G31.41, with a poorer correlation to the Sgr B2(N2) sources. Overall, the O-bearing species are the ones that are the closest to the 1:1 relation but this is biased by the fact that we analyzed the abundances relative to methanol.

Thumbnail: Fig. 16 Refer to the following caption and surrounding text. Fig. 16

Same as Fig. F.2 but comparing sources from the literature to N2b.

Thumbnail: Fig. 17 Refer to the following caption and surrounding text. Fig. 17

Same as Fig. F.2 but comparing sources from the literature to AN02.

Thumbnail: Fig. 18 Refer to the following caption and surrounding text. Fig. 18

Same as Fig. F.2 but comparing sources from the literature to AN03.

Thumbnail: Fig. 19 Refer to the following caption and surrounding text. Fig. 19

Same as Fig. F.2 but comparing sources from the literature to AN06c2.

Thumbnail: Fig. 20 Refer to the following caption and surrounding text. Fig. 20

Same as Fig. F.2 but comparing sources from the literature to AN06.

4 Comparison with chemical models

To aid in our interpretation of the varied molecular column densities determined toward each source, we made use of a preexisting grid of gas-grain chemical models. Shope et al. (2024) used the three-phase (gas/surface/bulk-ice) chemical kinetics model MAGICKAL to simulate hot-core chemistry under a range of physical parameters. Their chemical model and network were identical to the “final” model used by Garrod et al. (2022). Their physical treatment, consisting of a cold, isothermal collapse (stage 1) followed by a dynamically static warm-up phase (stage 2), was the same as that used by Garrod (2013), and similar to various related implementations. The Shope et al. grid varied several key parameters within their generic model: (i) the cosmic-ray ionization rate (CRIR) used throughout stages 1 and 2 of the model; (ii) the final gas density, nH, reached during stage 1, which carries through to stage 2; (iii) the warm-up timescale, twu, parameterized as the time taken to reach 200 K during stage 2; and (iv) the visual extinction, AV, used at the beginning of stage 1, which in turn scales the evolving AV value as collapse occurs. Parameters (i)–(iii) were varied logarithmically (nH = 2 × 106−2 × 1010 cm−3; ζ = 1.30 × 10−18−1.30 × 10−15 s −1; twu = 2 × 104−2 × 106 yr; see Table 1 of Shope et al. 2024). The initial visual extinction was tested with two values (2 mag and 3 mag). Due to the warm-up timescale only being relevant to stage 2, this resulted in a total of 84 stage-1 model runs, and 756 stage-2 runs, i.e. 9 stage-2 models for each preceding stage-1 run. Each stage-1/stage-2 combination is a single-point (0-D) representation of a hot core through its chemical evolution from cold, diffuse conditions to a dense, hot core that reaches a final temperature of 400 K.

Although the main focus of the Shope et al. (2024) models was a comparison with NGC 6334I, they also compared their model grid with the EMoCA-derived column densities for Sgr B2(N2) presented by Jørgensen et al. (2020) (see Belloche et al. 2016). This comparison involved finding the grid model with the best collective match to the observational column densities, each taken as a ratio with the column density of methanol, i.e. Robs,i = N(i)/N(CH3OH). Since the models produce only single-point fractional abundances, the peak value for each molecule (during stage 2) was taken in ratio with the peak methanol abundance, i.e. Rmod, i = Xpeak(i)/Xpeak(CH3OH), to represent the column density ratio (see Shope et al. 2024 and Belloche et al. 2019, for more technical discussions of the nuances of such comparisons). The logarithm of the quotient of the two ratios, i.e. log(Rmod,i/Robs,i) = mi, indicates the number of orders of magnitude that a modeled species diverges from its observed value. For species for which only an upper limit is determined, a modeled value that exceeds the observational value is treated in the same way, while a model value below the observed upper limit is treated as a perfect match, i.e. mi = 0. To determine the overall quality of match for a particular model, the root mean square of mi over all species i is calculated, with the best matching model having the lowest value. We note that this choice of matching parameter tends to favor models in which multiple species show a moderately good match, over models in which the match is very good for some species and very poor for others. In other words, species for which the model results are very divergent from the observations have a strong negative impact on the overall matching parameter.

Particularly noteworthy in the analysis of Shope et al. (2024) is their determination of the best-matching model for Sgr B2(N2) to have a warm-up timescale 2 × 106 yr and a cosmic-ray ionization rate of 1.3 × 10−17 s−1. The characteristic warm-up timescale should be considered simply a parameter, which is most important in determining how rapidly an already hot core becomes yet hotter; however, this timescale value is the longest of those tested in the grid. Meanwhile, the best-matching CRIR corresponds to the canonical interstellar value, ζ0, which is rather lower than might be expected for Sgr B2(N2), based on other recent studies (e.g., Bonfand et al. 2017). However, the determination of this match was based on the column densities of 11 chemical species, plus two upper limits, with several species omitted. Our present ReMoCA observations provide both a larger selection of molecules with which to compare, all of which are included in the model chemical network, as well as a sample of cores that the higher angular resolution of ReMoCA allowed us to distinguish. Collective fitting of this dataset may allow for a more accurate determination of the likely physical conditions in each core.

Because of these differences in the observational dataset, our method diverges slightly from that used by Shope et al. (2024). Crucially, we assumed that all five cores must experience the same CRIR. To impose this constraint, we determined firstly the best-matching model for each core within each fixed CRIR value; this produced a total of five match parameters for each CRIR value. We then took the root mean square of these five, so that each CRIR has an associated matching parameter that identifies the best collective value for the entire dataset (these values are shown in Fig. 21). We consider the best model for each source to be the one with the lowest matching parameter within the subset for which the CRIR has been collectively determined as the best match.

Some of the detected species are not included in our match calculations, as they are not present in the chemical network used in the model grid. However, following Shope et al. (2024), we also elected not to include certain species in the calculations even though they are present in the models. This includes species such as acetone (CH3C(O)CH3) and propanal (C2H5CHO), for which the network is sparse (although it has since been improved by Belloche et al. 2022). Others, like CH3NCO, have an uncertain chemistry that is not well reproduced by the models. Some omitted species, such as HC3N and HC5N, are not uniquely associated with the hot stage of the core, and therefore could be unduly influenced by chemistry in the lower density gas (either in the observations or during the early stages of the models). We also excluded the sulfur-bearing species OCS, H2CS, and CH3SH, as the dominant form of sulfur in dense clouds is uncertain and the chemistry is thus poorly defined. We retained SO and SO2 in the fit, as they are mainly produced in the gas phase in hot cores, following ice desorption, and they have been well reproduced by chemical models for a long time, given an appropriate initial sulfur abundance (here, Xinit(S) = 8 × 10−8). We note that the initial elemental abundances used in the models were not varied in the grid; Shope et al. (2024) used the values reported in Table 1 of Garrod (2013).

Table 5

Range of final (i.e. stage-2) gas densities from the model grid used in the comparison with each source, in the density-restricted setup.

There is a total of 37 species that could have in principle been included in the matching routine. Twelve were omitted (as described above), leaving 25, of which between 0 and 8 are observational upper limits. In contrast to Shope et al. (2024), for species bearing a nitrile (−CN) group the matching parameter is based on a ratio with CH3CN instead of CH3OH, in keeping with the comparison methods of Jørgensen et al. (2020). For methyl cyanide (CH3CN) itself, the ratio with methanol was retained.

Using this updated matching parameter method, we compared the observational data with the model grid in two ways. Firstly, we determined the best matching models based on an unrestricted comparison with all models; secondly, we restricted the comparison to a narrower range of gas densities that are reflective of observational values or estimates for each source (see Table 5). Values for AN02, AN03, and AN06 are based on determinations of n(H2) by Sánchez-Monge et al. (2017). The ReMoCA continuum emission is a factor of ∼3 fainter toward N2b than toward AN02 (see Fig. 3). Therefore we extended the restricted density range of this source to lower values. As we have seen in Sect. 3.1.2, two velocity components overlap along the line of sight to AN06. The secondary component, AN06c2, traces the edge of the dense core that contains AN02, AN03, and N2b, while the main component, AN06, is likely in the foreground or background. The continuum emission analyzed by Sánchez-Monge et al. (2017) contains the contribution of both components, yet we think that it is dominated by AN06. This motivated our decision to restrict the density range explored for AN06c2 to lower values.

Figure 21 indicates the collective matching parameter over all five sources, along with the individual subvalues for each source, at each CRIR value; panel (a) shows the results for the unrestricted setup, while panel (b) shows those of the density-restricted comparison. In both setups, CRIR values below the canonical interstellar value, ζ0, are disfavored. In the unrestricted setup, the best overall match occurs with ζ=10 ζ0. In the density-restricted setup, the best match is obtained at the maximum CRIR value tested, ζ=100 ζ0 = 1.3 × 10−15 s−1.

It is notable that the preference for the most extreme CRIR value is strongly associated with the variation in the ANO2 and AN03 matching parameters; indeed, the overall matching parameter is worse for all CRIR values in the restricted setup as a result of this dependence. In the unrestricted case, the best-matching models for both of these sources (with ζ = 10 ζ0) take stage-2 gas densities of nH = 2 × 106 cm−3, which are much lower than the observational estimates. Values of 2 × 108 cm−3 (with ζ = 100 ζ0) provide the best match in the restricted setup. Thus, at these higher densities for AN02 and AN03, a further elevated CRIR is required for the best match with observations, even though the other sources are provided a somewhat worse match under those conditions.

Our analysis continues with the exclusive use of the densityrestricted comparisons. Figure F.3 plots the matching parameter of each model for each observed source, ordered in bands corresponding to a fixed CRIR value. The bands are ordered from the best to the worst overall matching parameter. Within each band, and for each individual source, the matching parameter is plotted from best to worst. We note that there is no correspondence between model numbers among curves in the same band, as the results for each source-model combination are individually sorted. Furthermore, due to the density restrictions shown in Table 5, different sources have different numbers of grid models available to them. The plots indicate simply the degree of variation in matching parameter between different models for each source. Source AN06c2 universally shows the best matching parameter, comparing the curves at fixed sorted-model values. Matching the other sources is generally harder, i.e. there are fewer models that reproduce them well. This may be related to the fact that AN06c2 has the most upper limits, making a match easier. However, AN06 has the same number of upper limits, while being a worse match in all cases, when comparing best with best. The matching parameters for AN02 and AN03 generally appear to become the most rapidly divergent with increasing sorted-model value, i.e. they have a narrow range of physical conditions that produce a “good” match with observations.

Figure 22 shows the comparison between models and observations for each source; the best-matching model is shown in each case, for the best-matching CRIR of ζ = 1.3 × 10−15 s−1. Unfilled blue bars with an arrow indicate observational upper limits, while black bars indicate species that are omitted from the matching procedure, but which nevertheless exist in the models. With the CRIR value fixed, the key remaining parameters distinguishing the best-match models for each source are the gas density and the warm-up timescale. Particularly notable is the fact that, with the exception of AN06, the best match for each source is provided by the same model, taking values ζ = 100 ζ0, nH,final = 2 × 108 cm−3, twu = 2 × 104 yr, and AV,init = 3 mag. The best-matching model for AN06 has a factor 10 higher gas density, but is otherwise the same. We note that the best warm-up timescale for all sources is the shortest of those tested, and 2.5 × shorter than the “fast” warm-up timescale used in various past models by Garrod (2013) and others.

Many of the species omitted from the matching procedure are, unsurprisingly, not well reproduced. But in general, most common hot-core species included in the matching routine are reasonably well reproduced by the models, i.e. are within one order of magnitude (oom) of the observed values (as marked by the shaded region). For example, the structural isomers methyl formate and glycolaldehyde are both generally somewhat overproduced, although they are mostly within 1 oom of the observed values, and they also tend to scale with each other among the plots shown. However, due to the fact that four of the sources share a single best-matching model, any apparent variations in the quality of match between that (or any) pair of species among those four sources is solely due to variations in the observed values. In the two best-match models representative of either AN06 or collectively of AN02, AN03, AN06c2, and N2b, glycolaldehyde is overabundant compared with methyl formate.

H2CO, C2H5OH, CH3CHO, and CH2CO are consistently, if only modestly, overproduced for all five source-model matches. Species that appear consistently underproduced include CH3CN and C2H5CN, although this is not the case for the AN06 model. HNCO and NH2CHO are in some cases underproduced and in some cases overproduced, although they are mostly discrepant in the same sense; again, this apparent variation is essentially caused by matching a single model to a selection of sources that behave somewhat differently.

Figures. 915 show not only the observational ratios but the equivalent values from the best-matching chemical models for each source. Chemical model results are indicated by crosses. Although in our matching-parameter analysis we have only used molecules that are detected in at least one of our source datasets, it is valuable to compare the best-matching models with the observations for a much broader range of species, including those that have been omitted from the matching parameter.

The model abundances for protonated molecules with respect to their unprotonated forms (Fig. 9) are comfortably below all of the observational upper limits. Ion abundances in the models (during the hot stage, and under otherwise fixed CRIR conditions) tend to scale inversely with gas density, which influences collisional lifetimes. Within the range of densities tested in the restricted setup, the models would be unlikely to reach the observed upper limits. The model-source comparisons are therefore not a strong constraint on ion abundances.

Gas-phase destruction rates for radicals also tend to scale with gas density. The modeled values for radicals HCO, CH3O, and CH2OH with respect to H2CO and CH3OH (shown in Fig. 10) are substantially below the observational upper limits. The other four radicals plotted (CH2CHO, CH2CN, C3N, and NH2CO) are in some cases far above the upper limits. This may be caused by a lack of appropriate destruction mechanisms for these radicals; in particular, reactions with (abundant) atomic hydrogen could lead to their rapid destruction, as is the case for HCO, CH3O, and CH2OH (see, e.g., Tsang & Hampson 1986), leading to the production of H2 and another stable species. Other disproportionation and radical-pair production reactions are also possible, where energetically favorable. In the models, reactions with atomic H are frequently the most important destruction mechanisms for radicals that have them; more important even than sequential protonation and recombination, or CR-driven photodissociation. Atomic-H reactions for CH2CHO, CH2CN, NH2CO, and C3N are not currently present in the network (although see Bonfand et al., in prep.). In the latter case, however, such a pathway is unlikely, due to the absence of energetically favorable products. As such, the fact that the modeled C3N/HC3N ratio exceeds the observational upper limit for sources N2b, AN02, and AN03 indicates that a model with higher gas density could remedy this. However, the fact that HC3N was not included in the matching parameter could also affect the ratio.

Similarly, the ratios between low- and high-saturation species, C2H5CN/HC3CN and the equivalent larger nitriles (see Fig. 11) may also be affected by the difficulties of including those species in the matching parameters. The abundances of C2H3CN and C2H3OH are somewhat dependent on the branching ratio of the recombination of protonated ethyl cyanide and protonated ethanol, respectively, which could explain their overabundances for most model-source combinations. In Fig.12, -CHO and -COOH group-bearing species are somewhat overproduced relative to -CN, although relative to each other they are much closer to the observations. The NH2CHO/NH2CN ratio is higher than expected.

Modeled ratios of -CHO versus -CH3 group-bearing species shown in Fig. 13 are generally not far off the observations, although the CH2(OH)CHO/C2H5OH ratio is elevated in the case of N2b. As described by Shope et al. (2024), this may indicate that a higher gas density model is required.

As may be seen in Fig. 14, modeled ratios between the different nitrile homologs (CnH2n+1CN, n = 1−4) show a somewhat varied match with observations, with C2H5CN somewhat underproduced in all cases, while normal butyl cyanide is generally overproduced. The general underproduction of C2H5CN may be related to an overconversion to the larger cyanides in the ice mantles. For the CH3OCHO series, in cases where C2H5OCHO is detected, the ethyl to methyl formate abundance appears low in the models. As noted by Belloche et al. (2009), production of ethyl formate in the models is dominated by a methyl-radical addition to another radical whose precursor is methyl formate. In the present model grid, this occurs largely within the bulk ice, driven by CR-induced UV photoproduction of radicals. A number of different factors determining the availability of the CH3 radical (such as CH4 abundance and photodissociation efficiency) could be at play in determining the eventual modeled abundance of ethyl formate. Again, four of the sources are being matched with the same single model, so when they vary observationally, disagreements should become apparent. A higher resolution grid, with differentiated source-model matches for all sources, might lead to better individual matches for the molecules mentioned here.

Among the structural isomers shown in Fig. 15, CH3CN/CH3NC is not a good match for source AN06, although it is not bad for the other sources. In the present chemical network, which includes reactions used by Willis et al. (2020), the abundance of CH3NC is particularly sensitive to the reaction H+CH3NC → HCN+CH3, whose activation-energy barrier is highly uncertain. At the higher density experienced in the AN06 best-match model, this reaction would have a more pronounced effect. CH3NC was already omitted from the matching parameter due to the uncertainties in its chemistry. The poor match for AN06 could indicate that the barrier is too low in the network.

Figure 15 also shows the ratios of C2H4O2 species. CH3OCHO (methyl formate, MF) and CH2(OH)CHO (glycolaldehyde, GA) generally show roughly appropriate ratios with each other, although they are both very moderately overproduced in comparison with CH3COOH (acetic acid, AA) for sources where it is detected (AN02 and AN03). The production of acetic acid in the models occurs through somewhat different chemical pathways that involve acetaldehyde, CH3CHO (see Garrod et al. 2008, 2022). In the model representing N2b, AN02, AN03, and AN06c2, the MF/GA/AA ratio is 18/5.8/1. In the model representing AN06, a similar ratio 17/4.0/1 is achieved. These may be compared with the observed ratios of 29/1.6/<1 (N2b), 12/1.8/1 (AN02), and 14/1.5/1 (AN03). The modeled abundance of glycolaldehyde in particular, especially in comparison with source N2b, could be considered notably in excess of the observational value, when considering MF:GA ratios. Shope et al. (2024) considered the abundance of glycolaldehyde in particular and the possible explanations for extremely low ratios with methyl formate toward some sources. Those authors noted that a combination of high gas density and long timescale in the models was capable of reducing GA abundances during the warm-up period. Chemically, this occurred through the brief adsorption of gas-phase atomic H onto the grain surfaces, followed by reaction with a glycolaldehyde molecule that had recently become exposed on the surface but had not yet desorbed. This process was found to be effective under conditions where H adsorption would remain rapid for a long period of time.

The best-match models for the five sources studied in the present work achieve fairly high gas densities, but do not experience long enough periods for the effect to take hold strongly. In a less constrained physical model, or a more precise treatment of the warm-up process, GA destruction might be more efficient. As noted by Shope et al., the key period for GA destruction is not the entire warm-up period from ∼10−400 K, but the much shorter period between the start and end of major water-ice desorption. A somewhat longer period in this more limited temperature range would still be consistent with hot-core lifetimes on the order of a few 104 yr. More dynamically accurate simulations of hot-core chemical evolution are necessary to test these possibilities.

Thumbnail: Fig. 21 Refer to the following caption and surrounding text. Fig. 21

Overall matching parameter obtained for the comparison of the grid of chemical model results with observational abundances of all species across all sources, assuming a uniform CRIR. A lower matching parameter indicates a closer match between models and observations. Panel a: unrestricted density. Panel b: density restricted to the ranges indicated in Table 5.

Thumbnail: Fig. 22 Refer to the following caption and surrounding text. Fig. 22

Comparison of best-fit model abundances to observational results, using the density-restricted setup. Bars indicate the number of orders of magnitude by which each abundance ratio in the best-matching model exceeds (or otherwise) the observational value. Data correspond to the best-matching model for each individual source, within the subset of grid models with ζ = 100 ζ0, which is found to produce the best overall match across all sources. The shaded area represents values where the models and observations vary by 1 oom or less. Unfilled blue bars with an arrow indicate that the comparison is based on an observational upper limit. Black bars indicate observed species that were not included in the matching parameter analysis.

5 Discussion

5.1 Chemical evolution within Sgr B2(N2)

The ReMoCA survey has allowed us to derive the detailed chemical composition of four hot cores embedded in Sgr B2(N2) and a position, AN06c2, located at the edge of the dense core (Sect. 3.2). As demonstrated in Sect. 3.3, there is a tight correlation between the chemical compositions of AN02 and AN03, despite the fact that AN03 contains a HCH II region while AN02 does not (Sects. 3.1.1 and 3.1.3). While we could worry that the small angular separation between AN02 and AN03 (∼0.64′′) may have resulted in a mutual contamination of their observationally derived chemical compositions, we can exclude this interpretation because the correlation between N2b and AN02 is looser (Fig. 8a) despite their even smaller separation (∼0.55′′). The correlation between AN02 and AN03 is thus genuine. The presence of the HCH II region indicates that AN03 is more evolved than AN02 or more massive (hence more luminous). The similarity of the masses derived from dust emission by Sánchez-Monge et al. (2017) for AN02 and AN03 tends to discard the latter interpretation. According to the time sequence derived by Nony et al. (2024) for the high mass star-forming protocluster W49N, the lifetime of hot cores is 6 × 104 yr, including the phase when a hypercompact or ultracompact (UC) H II region coexists with the hot core emission. W49N contains a bit less than half as many hot cores with H/UCH II regions as hot cores without, which suggests a statistical lifetime of ∼2 × 104 yr for the former. Hot core lifetimes are similar in Sgr B2(N) and W49N (Bonfand et al. 2017) and so we conclude that the chemistry has not evolved significantly during a timescale of ∼2 × 104 yr around AN03 on a scale of 0.7′′(∼5700 au). With a size of ∼0.08′′, the HCH II region is one order of magnitude smaller than the molecular emission probed with ReMoCA toward AN03 (∼0.7′′−1′′). This is likely the reason why it has not yet had time to alter the chemistry of the hot core. As a matter of fact, there is no significant difference between AN02 and AN03 in terms of rotational temperature of their molecular emission (see Fig. F.4e).

The hot core N2b is most likely in an earlier evolutionary stage than AN02 and AN03 because, unlike these sources, it is neither associated with a HCH II region nor with masers. The rotational temperatures of its molecular emission are on average a bit lower than the temperatures of AN02 and AN03 by about 20 K (Fig. F.4a–b), which may support an earlier evolutionary stage. As reported in Sect. 3.1.4, N2b is associated with a compact dust core and its elevated rotational temperatures of ∼150 K suggest that it has already formed a protostar, unless its temperature structure is imposed by its neighbor AN02. If N2b contains a nascent protostar, then it is at most a few 104 yr younger than AN02. Therefore, it is remarkable that its chemical composition shows significant differences compared to that of AN02 and AN03. Overall, N2b’s molecular abundances relative to methanol are a factor ∼2−4 lower than those of AN02 and AN03, and their dispersion is larger than in the correlation plot of AN02 versus AN03 (compare panels a–b and e of Fig. 8). Two organic molecules stand out with abundances that are even lower by more than one order of magnitude in N2b : NH2CN and t-HCOOH. It is tempting to conclude that the chemistry evolves significantly during the initial phase (a few 104 yr) of thermal desorption at the hot core stage, which may suggest that gas-phase processes play an important role in reshaping the chemical composition after desorption of the molecules from the grains. In this context, it would be interesting to understand why NH2CN and t-HCOOH in particular appear to be so sensitive to this phase.

We proposed in Sect. 3.1.2 that AN06 is a hot core that is not embedded in the dense region that contains AN02, AN03, and N2b. Like N2b, it is not associated with a HCH if region or any masers. AN06 may thus be in an early evolutionary stage, maybe even earlier than N2b given its slightly lower rotational temperatures (see Fig. F.4d). Its chemical composition correlates with that of N2b, albeit with a much larger dispersion than AN02 versus AN03 (compare panels d and e of Fig. 8). If AN06 and N2b are at a similar stage of evolution, then the differences in their chemical compositions suggest that they reside in different environments. Alternatively, AN06 may be less massive than the other hot cores of Sgr B2(N2) and the chemical differences may result from a difference in evolutionary stage. Sánchez-Monge et al. (2017) derived a mass of AN06 that is 30% smaller than the mass of AN03. However, their dust continuum observations cannot disentangle the two velocity components that overlap along the line of sight. The mass calculated for AN06 may thus be largely overestimated. The most striking chemical difference between AN 06 and N2b is that HNCO and NH2CHO are nearly two orders of magnitude less abundant in the former. Understanding the cause of this huge deviation from the overall correlation may tell us what controls the larger dispersion of the correlation between AN06 and N2b. Finally, while the chemical compositions of AN06 and AN06c2 correlate well apart from an overall shift by a factor ∼2 (Fig. 8h), the abundance of formamide stands out by nearly two orders of magnitude. In this respect, the composition of AN06c2 is much more similar to that of the dense region that contains AN02, AN03, and N2b. This supports the idea that AN06c2 traces the edge of the region and AN06 overlaps along the line of sight but is physically unrelated to this region.

5.2 Tight correlation between Sgr B2(N2) and G31.41

Located in the Scutum spiral arm at a distance of 3.75 kpc (Immer et al. 2019), the hot core G31.41 is a small protocluster of luminosity ∼4.5 × 104 L that consists of four high-mass young stellar objects (YSOs) with gas masses of 15−26 M within a radius of 400-500 au (Beltrán et al. 2021). Each of them drives an outflow (Beltrán et al. 2022) but none of them harbors a HCH II region (Cesaroni et al. 2010; Beltrán et al. 2021). AN03 and AN02 have gas masses of about 500 M within a radius of about 2800 au (Sánchez-Monge et al. 2017). Assuming powerlaw density profiles with a power-law index in the range from −1.5 to −2, then the mass enclosed in a radius of 500 au is 40−90 M, implying that AN03 and AN02 may form highermass stars than G31.41’s YSOs by a factor 2–3, or even higher given that AN03 already contains a HCH II region. Among the four sources from the literature investigated in Sect. 3.5, G31.41 is the one with the tightest correlation of its overall chemical composition with that of N2b, AN02, and AN03. The overall chemical composition of hot cores thus seems to be relatively insensitive to the environmental features that distinguish the Galactic center region from the Galactic disk (e.g., level of turbulence, CRIR, gas temperature). The O-bearing species show a poorer correlation between G31.41 and N2b, AN02, and AN03 than the N-bearing ones. This suggests that the former may be more sensitive to the environmental conditions than the latter.

5.3 Chemical segregation between classes of molecules

Chemical segregation between O-bearing and N-bearing molecules has been reported in various star-forming regions since the 1990s (see, e.g., Qin et al. 2022 and short reviews on this topic in Sect. 4.3 of Jørgensen et al. 2020 and Sect. 4.7 of Busch et al. 2024). Among the four classes of molecules investigated in Sect. 3 (O-, N-, N+O-, and S-bearing), the N-bearing class stands out as the one with the largest variance across the sample of sources investigated in Sect. 3.5. For this class of molecules, the composition of Sgr B2(N2)’s hot cores correlates well with that of the other two hot cores/corinos G31.41 and IRAS 16293B but does not correlate at all with that of G+0.693 and TMC-1. This is in stark contrast to the class of S-bearing molecules which always shows a correlation with Sgr B2(N2)’s hot cores, strong in the case of G31.41, IRAS 16293B, and G+0.693, and moderate in the case of TMC-1. The O-bearing molecules lie in between, with a strong correlation in the case of G31.41, IRAS 16293B, and G+0.693, but a poor one in the case of TMC-1.

TMC-1 is a cold (∼10 K) starless core of moderate density (3−8 × 104 cm−3, see, e.g., Pratap et al. 1997; Lique et al. 2006) located in the low mass star-forming molecular cloud Taurus. Part of its molecular composition may be controlled by gasphase chemistry (e.g., Cernicharo et al. 2021), but pure gas-phase chemical models have difficulties in reproducing the abundances of certain molecules like acetaldehyde (e.g., Agúndez et al. 2025), suggesting that nonthermal desorption of molecules formed on the grains, triggered by, for instance, cosmic rays or the release of chemical energy, likely contributes to the composition of the gas phase as well (see, e.g., the case of CH3OH, C2H5OH, and C2H5CHO in Agúndez et al. 2023).

G+0.693 is also a starless region of moderate density (104− 105 cm−3, Zeng et al. 2020; Colzi et al. 2024) and low dust temperature (∼20 K, Etxaluze et al. 2013), but it is located in the Sgr B2 molecular cloud complex, about 54′′ away from Sgr B2(N2), which means that they both have been subject to similar levels of turbulence, cosmic rays, and radiation (UV photons, X-rays) during their evolution. The chemistry of G+0.693 is thought to be dominated by low-velocity shocks that eject molecules from the grain mantles through sputtering (e.g., Requena-Torres et al. 2006; Zeng et al. 2018). If we consider that most molecules are formed on the grains and that the molecular composition of Sgr B2(N2)’s hot cores was in large part inherited from the composition of the grain mantles formed during the prestellar phase, then differences in the compositions of G+0.693 and Sgr B2(N2)’s hot cores could result from the shocks that affect G+0.693. Alternatively, the composition of the hot cores may have been altered by processes on the grains during the warm-up phase. However, it is also likely that the large difference in gas density between G+0.693 and Sgr B2(N2)’s hot cores has a major effect on the post-desorption behavior of gas-phase molecules. Lower gas densities in G+0.693 should result in both a slower rate of destruction and a different balance of destructive ions, depending on how complete is the loss of ice due to sputtering during the passage of a shock (Willis et al., in prep.). The amount of ammonia in the gas can assist in the destruction of some molecules (especially amine group-bearing species) while enhancing the survival of others (Taquet et al. 2016; Garrod & Herbst 2023).

The CRIR in the Galactic center region is higher than in the Galactic disk by 1–2 orders of magnitude (e.g., Indriolo et al. 2015; Le Petit et al. 2016). The fact that N-bearing species of Sgr B2(N2)’s hot cores correlate well with hot cores/corinos in the Galactic disk (G31.41 and IRAS 16293B) but neither with G+0.693 nor with TMC-1 implies that the CRIR is not the dominant factor producing the large variance of the N-bearing class. We conclude that the class of N-bearing molecules reacts more sensitively to shocks (G+0.693), low-temperature gasphase chemistry subsequent to nonthermal desorption (TMC-1), or density (G+0.693 and TMC-1) than O-bearing and S-bearing species. This is in line with the findings of Busch et al. (2024) who reported an abundance enhancement of cyanopolyynes relative to methanol in post-shock gas in the outflow of Sgr B2(N)’s main hot core. In contrast, the strong correlation of the O-bearing content of Sgr B2(N2)’s hot cores with G+0.693 and the poor one with TMC-1 suggest that, in comparison to the processes that control the O-bearing content of hot cores, shocks do not greatly modify the relative abundances of the class of O-bearing species while low-temperature gas-phase chemistry does to some extent.

We showed in Sect. 3.5 that there is a tight correlation between the chemical compositions of the low-mass protostar IRAS 16293B and that of N2b, AN02, and AN03 after normalizing the column densities by class of molecules, but the classes are globally shifted in abundance with respect to each other: N-bearing molecules are underabundant by about 2 oom in IRAS 16293B, O+N-bearing and S-bearing species by about 1 oom. Similar shifts between classes of molecules also exist in the case of G31.41, but they are less pronounced: 1−1.5 oom for the N-bearing and <1 oom for the O+N-bearing and S-bearing species. G31.41 seems to be forming stars of lower mass than Sgr B2(N2)’s sources (see Sect. 5.2), but the difference may only be a factor of a few while there is a much larger gap in mass between G31.41 and IRAS 16293B (>1 oom). Therefore it seems unlikely that the abundance shifts between classes of molecules result primarily from physical differences between regions forming low-mass and high-mass stars, such as the temperature during the prestellar phase, the duration of the prestellar phase, or the level of UV radiation during the warm-up phase.

Shocks may play a role in hot cores and corinos, via, for example, accretion shocks of infalling material onto a circumstellar disk. This could in turn lead to somewhat different desorption behavior, including only partial desorption of ices, given a sufficiently weak shock. This could lead to only the outer layers being desorbed. These outer layers are often considered to be relatively water-poor, and have been associated with observed apolar spectral signatures of CO and CO2 (Garrod & Pauly 2011), while being richer in CO and its products such as methanol. It is likely that these layers would also be rich in O-bearing COMs (which are also ultimately related to CO), if they exist in the ice. On the other hand, the deeper ice layers are shown by models to have a much greater NH3 fraction. One might expect that amine-group bearing molecules would also be more abundant in those deeper ice layers, while their binding energies with water also tend to be greater than many of the CO-related species. Thus if only weak accretion shocks are active, only the weakly bound upper layers of O-bearing molecule-rich ice might be released into the gas. In hot cores, the accretion shocks may be stronger (leading to more complete ice desorption, either by sputtering or thermal desorption), and the heating of the grains by the protostar is likely to lead to a more global, sustained, and vigorous heating of the grains than in a hot corino. Hot cores may therefore undergo a more complete desorption of the ices than hot corinos, in general, leading to greater N-bearing species being released.

Other mechanisms may also lead to sporadic heating of the ices, such as accretion outbursts. Low-mass protostars may flashheat their surroundings multiple times (e.g. Vorobyov & Basu 2015). This could again lead to the desorption of only the outer and more weakly bound ice layers, which would be O-rich and N-poor. It is also suggested, based on simulations (Meyer et al. 2021), that outbursts are also common for high-mass protostars.

In principle, any of these processes that may lead to only partial ice desorption might well encourage the preferential release of O-rich material, while leaving deeper, N-rich material in place. Differences in molecular O/N ratios between sources would then be an indicator of how violently and completely the ice mantles are desorbed.

Unfortunately, while current chemical models can operationalize differences in binding energies between species, hot core models do not currently distinguish between ice layers beyond the typical bulk+surface-layer paradigm, and they are also not yet capable of taking local ice surface composition into account in the determination of binding energies. Testing the above idea will require substantial upgrades in model capabilities.

5.4 Chemical model behavior

The chemical model grid is capable of reproducing the observed abundances to a reasonable degree of accuracy (within 1 oom) in most cases. However, the rather coarse resolution of the grid means that we are unable to determine clearly the possible variations in physical conditions that would lead to the observed chemical distinctions between the five sources. The method imposes the requirement that a single value of CRIR be used in the matching routine, under the assumption that the same rate should be applicable to the entire region. A narrow range of gas density values is also imposed, based on observational constraints. Under these conditions, a single set of model parameters is found to be the best match to four of the sources. It is plausible that with a much higher resolution grid, better individual matches could be obtained for each source that would still be consistent with the imposed restrictions.

The matching method is nevertheless valuable in identifying two key parameter values: firstly, that the CRIR must necessarily be “high” to best match all five sources; and secondly, that the warm-up timescale for all best-matching models is 2 × 104 yr, the shortest value tested. Although the warm-up timescale in these models is a parameterized value that need not necessarily match the true dynamical timescale, the best-match value is both consistent across sources and in keeping with expected lifetimes for hot cores, on the order 6 × 104 yr (Nony et al. 2024; Bonfand et al. 2017). We note that our grid also includes values 3.56 and 6.32 × 104 yr, which did not produce the best match for any of the sources.

Again, with four out of five sources being reproduced by the same chemical model, the modeling comparison itself has limited value in determining the origins of the chemical differences between sources. However, the models can help to indicate what conditions could be adjusted to produce a closer match.

For example, the lower observed values of HNCO and NH2CHO toward AN06, as opposed to N2b, might indicate that different CRIR values are required between sources. The models indicate that the extreme ζ = 100 ζ0 value suppresses the abundances of HNCO and NH2CHO in the ices during the cold collapse stage, due to the enhanced UV field. Imposition of such a strict agreement in CRIR between Sgr B2(N2) sources might therefore be inappropriate; however, to account for a factor of nearly two orders of magnitude would require substantial variation, if the differences in these molecular abundances were caused by this alone.

Alternatively, other UV ionization sources could be active, requiring a slightly lower overall CRIR for all sources, and allowing more influence for variations in visual extinction between sources. Moreover, a more careful treatment of the 2-D/3-D structure during stages 1 and 2 would help to determine how important external UV may be to the eventual detected abundances of certain COMs, by allowing regions of the core with different initial visual extinctions to contribute to the molecular abundances that are ultimately observed. Using the simple physical model that we have employed here, it is not possible to determine which regions of a natal cold core may later contribute most strongly to the release of icy COMs into the gas phase. Also, as noted above, shocks or luminosity outbursts could be important to the release of the ice mantles, determining both the location of COM desorption and the degree of ice loss, which could in turn lead to variations in which molecules are released. Testing such ideas requires dedicated coupled chemical-dynamical modeling that cannot currently be carried out in high enough volumes to allow for a model-grid comparison similar to what is used here.

5.5 Comparison to ALMA 1.2 mm survey of Sgr B2(N)

Möller et al. (2025) derived the chemical composition of several dozen hot cores identified by Sánchez-Monge et al. (2017) on the basis of an imaging spectral line survey performed with ALMA between 211 and 275 GHz with a spectral resolution of 490 kHz (0.7−0.5 km s−1) and an angular resolution ranging from 0.4′′ to 0.7′′. They used super-resolution data cubes restored with a beam of 0.4′′. For each source, they modeled a spectrum averaged over a polygon (see the shaded blue polygons in their Fig. 2). Following Sánchez-Monge et al. (2017), they divided Sgr B2(N2) into three sources (AN02, AN03, and AN06) with polygon sizes on the order of 2′′ × 1′′. Their polygon around AN02 also includes the position that we named N2b.

Figure F.5 compares the column densities that we derived toward AN02, AN03, AN06, and AN06c2 to the column densities of the “Core Components” reported for AN02, AN03, and AN06 in Tables L.57, L.59, and L. 65 of Möller et al. (2025), respectively. In the case of AN02, most molecules were modeled by these authors with only one velocity component with a larger line width and thus we compared their 1.2 mm column densities to the sum of the column densities that we derived for the two velocity components that we identified in this source with ReMoCA. For both C2H5CN and CH3CN, we used only the main component reported by Möller et al. (2025) and ignored the other component that has a six orders of magnitude lower column density7. For CH3OH, we added the column densities of the two components modeled by Möller et al. (2025). In the case of AN03 (Table L.59 of Möller et al. 2025), we selected only the high-column-density components derived for CH3CN and CH3OH, and we ignored the HC3N component that has a line width smaller than 0.5 km s−1. Finally, Möller et al. (2025) modeled the spectrum of AN06 with two velocity components similar to ours only for C2H5CN and OCS. Most other molecules were modeled with a single velocity component with a larger line width. Therefore, apart from C2H5CN and OCS, Figs. F.5c and d display the same 1.2 mm column densities.

Given that Möller et al. (2025) used spectra averaged over large polygons while we analyzed spectra from a single pixel, we do not expect the ReMoCA column densities to match exactly the 1.2 mm column densities. While there is a decent overall correlation between the column densities of both studies, Fig. F.5a shows that the ReMoCA column densities of AN02 are generally higher (by a factor of ∼2) than the 1.2 mm ones. This must result in part from the fact that the 1.2 mm column densities were derived from spectra that were averaged over typical sizes of 2′′ × 1′′ while the ReMoCA ones are peak column densities. In addition to this, the 1.2 mm column densities of AN02 were derived assuming an emission size of 1.7′′ while, guided by their integrated intensity maps, we modeled most molecules with emission sizes (FWHM) of 0.7–1.0′′ Furthermore, we do not know if the column densities reported by Möller et al. (2025) account for vibrational and conformational corrections. These corrections can be as high as a factor 1.5–2 for some of the molecules displayed in Fig. F.5 (for instance C2H3CN, C2H5CN, HC3N, CH3OCHO). The situation is similar for AN 03, which was modeled with an emission size of 1.6′′ by Möller et al. (2025) versus 0.7−1.0′′ in our study, with a slightly better match between the two studies compared to AN02 (Fig. F.5b).

The case of AN06 looks different (Figs. F.5c and d). The 1.2 mm column densities are in many cases higher than the ReMoCA ones (by more than a factor of two). We think that this is due to the fact that the aperture used by Möller et al. (2025) includes more emission from the Sgr B2(N2) dense core, meaning more contamination from the AN03 region, than our single-pixel spectra. As a matter of fact, our Weeds models for AN06 were as far as possible optimized based on setups S4 and S5, which have the smallest beams in our survey (0.3–0.5′′). The line emission of AN06 in the ReMoCA setups with larger beams (S1–S3) is generally underestimated by our Weeds models, due to contamination by the AN03 environment, while it is well fitted for S4 and S5.

One molecule stands out in Fig. F.5. The column densities reported for NH2D by Möller et al. (2025) are much higher (up to two orders of magnitude for AN03) than the column densities derived from the ReMoCA survey. The ReMoCA column density toward AN03 was derived from four detected transitions while Möller et al. (2025) detected only one transition in some of their sources. In their Fig. F.20, they show the spectra toward AN05 and AN08 only, and so it is unclear to us if the NH2D column densities reported for AN02, AN03, and AN06 are reliable. Their Fig. F.20 indicates that there is significant contamination by other molecules at the frequency of the NH2D transition toward AN08. Therefore, it could well be that the emission detected toward AN02, AN03, and AN06 at the frequency of this NH2D transition is dominated by these contaminating molecules and does not trace NH2D.

The comparison presented in this section strengthens our confidence in the reliability of the column densities derived from the ReMoCA survey. They faithfully represent the genuine chemical composition of the four hot cores AN02, AN03, AN06, and N2b. Thanks to the lower degree of spectral confusion at 3 mm compared to the 1.2 mm range8, the ReMoCA survey has allowed us to detect many more molecules (up to 58 toward AN03), in particular many more COMs: while Möller et al. (2025) reported the detection of ten COMs toward AN02, AN03, and AN06, we managed to identify 22,24, and 17 COMs toward these sources, respectively (see Table 3).

5.6 Unidentified lines

There are still many unidentified lines in the ReMoCA spectra of the sources embedded in Sgr B2(N2) (see Fig. D.1). Table B.1 provides the list of unidentified lines that are brighter than 10 K (signal-to-noise ratio higher than 10–30 depending on the setup) in the spectrum of N2b. There are 255 such lines, which translates into about eight bright unidentified lines per GHz on average. As was already the case for our earlier singledish survey of Sgr B2(N) (Belloche et al. 2013), we think that most of these bright unidentified lines correspond to rotational transitions from within vibrationally excited states of molecules which have their vibrational ground state and in some cases a few vibrational states already included in our model but for which spectroscopic predictions of higher vibrational states are still missing in our database, for instance ethyl cyanide or ethanol.

6 Conclusions

We used the imaging spectral line survey ReMoCA performed with ALMA to probe the chemical composition of Sgr B2(N2), the secondary hot molecular core of the high mass star-forming protocluster Sgr B2(N). At the angular resolution of the survey, Sgr B2(N2) consists of four hot cores (N2b, AN02, AN03, and AN06). Two velocity components were detected toward AN06. The main component traces a hot core associated with a compact dust continuum source. It is likely in the foreground or background of the dense core Sgr B2(N2). The second velocity component, which we called AN06c2, traces the edge of Sgr B2(N2). AN02 is associated with a dust continuum source and possibly drives an outflow. It is likely younger than AN03 which is associated with a HCH II region. N2b is associated with a fainter dust continuum source, has narrow line widths, and is possibly younger than AN02 and AN03. All four sources have line-rich spectra. We identified up to 58 molecules, including up to 24 COMs. In addition to these molecules, many less abundant isotopologs were also detected. We derived the molecular composition of the four hot cores as well as AN06c2 under the assumption of LTE. The main results of our analysis are the following:

  1. The pairs of sources that show the best correlations of their chemical compositions are AN02/AN03 and AN06/AN06c2. AN03 stands out for its S-bearing molecular content that poorly correlates with that of the other sources, and AN06 for its underabundant O+N-bearing molecules.

  2. AN06 may be in a similar evolutionary stage as N2b. Their compositions correlate with each other, albeit with a larger dispersion than AN02 versus AN03. The most striking difference between AN06 and N2b is that HNCO and NH2CHO are almost two orders of magnitude less abundant in the former.

  3. Protonated molecules XH+and radicals X were not detected. The most stringent constraint was obtained in both cases for N2b, with abundance ratios XH+/X and X/XH lower than 10−3.

  4. Molecules with a double bond are in nearly all cases much less abundant than the corresponding single-bond species, by up to at least two orders of magnitude. The only clear exception is the pair CH2NH/CH3NH2 with a ratio close to unity.

  5. Only one complex carboxylic acid was (tentatively) detected. If the abundance ratio NH2CH2COOH/NH2CH2CN is similar to the ratio CH3COOH/CH3CN then the abundance of glycine must be at least one order of magnitude below its current ReMoCA upper limit.

  6. The abundances of series of homologous molecules drop by about one order of magnitude at each further step in complexity, except for NH2CN, which is less abundant than NH2CH2CN by one order of magnitude, and CH3CN, which has a similar abundance as C2H5CN.

    To gain more insights into the workings of interstellar chemistry, we compared the chemical composition of Sgr B2(N2)’s hot cores to that of other sources that have been studied in detail in the literature: the lower-density Sgr B2 source G+0.693, thought to be affected by shocks, the hot core G31.41 in the Galactic disk, the hot corino IRAS 16293B, and the cold starless core TMC-1, both located in the Solar neighborhood. We focused in particular on the behavior of four classes of molecules (O-bearing, N-bearing, O+N-bearing, and S-bearing). The main results of this comparison are the following:

  7. The source that shows the tightest correlation with N2b, AN02, and AN03 is G31.41. This implies that the overall chemical composition of hot cores is relatively insensitive to the environmental features that distinguish the Galactic center region from the Galactic disk.

  8. For each class of molecules taken separately, there is also a good correlation between IRAS 16293B and Sgr B2(N2)’s hot cores, but the four classes are shifted with respect to each other. Such a segregation also exists for G31.41 relative to Sgr B2(N2)’s sources but it is less pronounced.

  9. There is no overall correlation between G+0.693 and Sgr B2(N2)’s hot cores. The lack of correlation results mainly from the N-bearing species, while the other classes correlate better. The S-bearing species are more prominent in G+0.693.

  10. There is absolutely no correlation between the compositions of TMC-1 and Sgr B2(N2)’s hot cores.

  11. Among the four classes of molecules, the class of N-bearing species is the one that reveals the largest variance. Firstly, its abundance distribution with respect to Sgr B2(N2)’s hot cores goes from a tight correlation (G31.41) to no correlation (G+0.693 and TMC-1). Secondly, this class shows overall shifts with respect to the other classes that can be large, up to two orders of magnitude (IRAS 16293B).

  12. In contrast, the class of S-bearing molecules has the smallest variance, with a high degree of correlation for G31.41, IRAS 16293B, and G+0.693 with respect to Sgr B2(N2).

    We also compared the molecular composition of Sgr B2(N2)’s sources to a grid of generic hot-core models computed with the chemical kinetics code MAGICKAL. The main results of this comparison are the following:

  13. We confirm previous evidence for a strongly elevated CRIR for Sgr B2(N2). The models produce their best collective match with observational abundances for the five sources adopting a CRIR value ζ=1.3 × 10−15 s−1, which is the highest value tested. The shortest warm-up timescale (2 × 104 yr) is found to produce the optimum match with all sources.

  14. The chemical model grid is not able to distinguish adequately the chemical differences between sources. Four of the sources are best matched by a single model, under conditions where the gas density is constrained to an observational range of values. The models therefore provide little direct indication of the origins of chemical differences between the sources.

The model grid might provide more information given higher resolution in physical conditions. However, the match with observations may be much further improved by the use of an explicit treatment of gas dynamics in tandem with the chemistry.

We conclude from the observed behavior of the class of N-bearing molecules that this class reacts more sensitively to shocks (G+0.693), low-temperature gas phase chemistry after nonthermal desorption (TMC-1), or density (G+0.693 and TMC-1) than the classes of O-bearing and S-bearing species. A possible interpretation of the segregation between N-bearing and O-bearing molecules is that only partial ice desorption might encourage the preferential release of the outer ice layers that are rich in CO and related O-bearing species, while leaving deeper, N-rich material in place. The overall abundance shifts observed in the gas phase between the classes of N-bearing and O-bearing molecules may thus indicate how violently and completely the ice mantles are desorbed. Testing this idea will require substantial improvements of chemical models that do not currently distinguish between ice layers beyond the typical bulk+surface-layer paradigm.

Acknowledgements

AB thanks all the physicists who sent him spectroscopic predictions over the past two decades, in particular E. Alonso, J. L. Alonso, B. Arenas, L. Bizzocchi, L. Bonah, S. Brünken, C. Cabezas, Ning Chen, L. Coudert, C. Endres, Z. Fried, E. Gougoula, S. Gruet, B. Heyne, V. Ilyushin, Z. Kisiel, I. Kleiner, K. Kobayashi, L. Kolesniková, J. Koucký, L. Margulès, A. Maris, M.-A. Martin-Drumel, B. McGuire, C. Medcraft, I. Medvedev, M. Melosso, P. Misra, R. Motiyenko, M. Ordu, M. Sanz-Novo, M. Schnell, P. Stahl, O. Zingsheim, and Luyao Zou. We value the long-term efforts of the spectroscopic community that crucially feed the databases used by astrophysicists. RTG thanks the National Science Foundation for funding through the Astronomy & Astrophysics program (grant number 2206516). We thank B. Shope for advice on the chemical model grid. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.00074.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. The interferometric data are available in the ALMA archive at https://almascience.eso.org/aq/. Part of this work has been carried out within the Collaborative Research Center 956, subproject B3, funded by the Deutsche Forschungsgemeinschaft (DFG, the German research foundation) – project ID 184018867. HSPM acknowledges support by the DFG through the Collaborative Research Center 1601, sub-projects A4 and Inf – project ID 500700252.

Appendix A Column densities from the literature

Table A.1 reports column densities and rotational temperatures of G31.41, IRAS 16293B, G+0.693, and TMC-1 that we collected from the literature. These column density values were used inFigs. 1620.

Table A.1

Column densities of G31.41, IRAS 16293B, G+0.693, and TMC-1 collected from the literature.

Appendix B Bright unidentified lines

Table B.1 provides the list of unidentified lines that are brighter than 10 K in the ReMoCA spectrum of N2b.

Table B.1

List of unidentified lines brighter than 10 K in the ReMoCA spectrum of N2b.

Appendix C Spectroscopic references

Table C.1 provides the list of relevant references for the 448 spectroscopic entries that were used in this work to model the observed spectra or derive upper limits to the column densities of nondetected molecules.

Table C.1

Spectroscopic references of molecules analyzed in this work.

Appendix D Spectra

Figure D.1 shows the complete ReMoCA spectral survey toward five positions in Sgr B2(N): the four hot-core positions N2b, AN02, AN03, and AN06, and the ultracompact (UC) H II region K4. The continuum-subtracted spectra were shifted vertically by 0, 200,400,600, and 900 K, respectively. The continuum level is indicated in blue on the right below each spectrum. LTE synthetic spectra computed with Weeds and containing the contribution of all identified molecules are overlaid in red on the hot-core spectra. The parameters of these LTE models are listed in Tables E.1E.5. The molecular identification of a selection of spectral lines is indicated in various colors: magenta for rotational lines in the vibrational ground state, teal for rotational lines in vibrationally excited states (with the suffix v added to the label), and orange for lines seen in absorption. Green labels above the K4 spectrum mark the frequencies of hydrogen and helium recombination lines. To facilitate the visual matching between the spectral lines and their magenta or teal labels, short vertical bars are displayed in the same color below each hot-core spectrum. Each panel corresponds to half a spectral window (W) of a single spectral setup (S). The setup index (from 1 to 5) and the spectral window index (from 0 to 3) of each panel are indicated in the bottom left corner along with the HPBW. The spectra are displayed in brightness temperature scale and were corrected for primary beam attenuation. The panels are ordered in increasing frequency of the spectral windows. There is a substantial frequency overlap between the upper sideband of setups 1 and 2 and the lower sideband of setups 4 and 5. The angular resolution of these pairs of setups differs by a factor ∼1.5, which motivates our decision to display all spectral windows in Fig. D.1, despite their frequency overlap.

Thumbnail: Fig. D.1 Refer to the following caption and surrounding text. Fig. D.1

ALMA continuum-subtracted spectra of the hot core positions N2b, AN02, AN03, AN06 and the UCH II region K4 (black). LTE synthetic spectra are overlaid in red. See Appendix D for a detailed description of the figure.

Appendix E LTE model parameters and column density upper limits

Tables E.1E.5 list the parameters of our best-fit LTE models of N2b, AN02, AN03, the second velocity component toward AN06 (called AN06c2), and the first velocity component toward AN06, respectively. Tables E.6E.10 provide the column density upper limits derived for 146 molecules that were not detected in any of the investigated positions.

Table E.1

Parameters of our best-fit LTE model of N2b.

Table E.2

Parameters of our best-fit LTE model of AN02.

Table E.3

Parameters of our best-fit LTE model of AN03

Table E.4

Parameters of our best-fit LTE model of the second velocity component toward AN06 (AN06c2).

Table E.5

Parameters of our best-fit LTE model of AN06.

Table E.6

Column density upper limits toward N2b.

Table E.7

Column density upper limits toward N2b.

Table E.8

Column density upper limits toward AN03.

Table E.9

Column density upper limits toward the second velocity component toward AN06 (AN06c2).

Table E.10

Column density upper limits toward AN06.

Appendix F Additional figures

Figure F.1 shows the chemical composition (relative to methanol) of AN06, AN06c2, AN03, and AN02, normalized to the composition of N2b. Figure F.2 shows the same as Fig. 8 but with a different color scheme. Figure F.3 shows the matching parameters for all tested chemical models. Figure F.4 displays correlation plots of rotational temperatures for various pairs of positions in Sgr B2(N2). Figure F.5 compares the ReMoCA column densities to the column densities derived by Möller et al. (2025).

Thumbnail: Fig. F.1 Refer to the following caption and surrounding text. Fig. F.1

Same as for Fig. 7 but normalized to N2b. The dashed lines indicate values of 0.1, 1, and 10.

Thumbnail: Fig. F.2 Refer to the following caption and surrounding text. Fig. F.2

Same as Fig. 8 but with colors coding for the class of molecules (black for O-bearing, blue for O+N-bearing, red for N-bearing, green for pure hydrocarbon, and yellow for S-bearing).

Thumbnail: Fig. F.3 Refer to the following caption and surrounding text. Fig. F.3

Matching parameters for all tested models, using the density-restricted setup. Each curve shows the matching parameter for a single source as a function of all the models in the grid with which it is compared, using the same fixed CRIR. The results are arranged in bands, one for each CRIR. The bands are arranged left to right in order of best to worst overall matching parameter. Within each band, the models contributing to each curve are arranged in order of increasing (i.e. worse) matching parameter, as determined independently for each source comparison.

Thumbnail: Fig. F.4 Refer to the following caption and surrounding text. Fig. F.4

Correlation plots of rotational temperatures for various pairs of positions. The x- and y-axes of each panel correspond to the positions written in the bottom right and top left corners, respectively. The color coding of the molecules is the same as in Fig. 8, plus black for methanol, and is indicated on the right. The plain line indicates the one-to-one relation.

Thumbnail: Fig. F.5 Refer to the following caption and surrounding text. Fig. F.5

Comparison of the column densities obtained from the ReMoCA survey toward AN02, AN03, AN06, and AN06c2 to those derived from an ALMA survey of Sgr B2(N) at 1.2 mm by Möller et al. (2025). The ReMoCA column densities correspond to peak column densities while the 1.2 mm column densities were derived from spectra averaged over the shaded blue polygons shown in Fig. 2 of Möller et al. (2025). Panels c and d show the same 1.2 mm column densities except for C2H5CN and OCS which were modeled with two velocity components similar to the components that we extracted from the ReMoCA survey. The dashed, dotted, and plain lines and the color coding and labels of the molecules are the same as in Fig. 8.

References

  1. Agúndez, M., Loison, J. C., Hickson, K. M., et al. 2023, A&A, 673, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Agúndez, M., Molpeceres, G., Cabezas, C., et al. 2025, A&A, 693, L20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Ahrens, V., Lewen, F., Takano, S., et al. 2002, Z. Naturforsch. A, 57, 669 [NASA ADS] [CrossRef] [Google Scholar]
  4. Alekseev, E. A., Dyubko, S. F., Ilyushin, V. V., & Podnos, S. V., 1996, J. Mol. Spectr., 176, 316 [NASA ADS] [CrossRef] [Google Scholar]
  5. Allegrini, M., Johns, J. W. C., & McKellar, A. R. W., 1979, J. Chem. Phys., 70, 2829 [NASA ADS] [CrossRef] [Google Scholar]
  6. Alonso, E. R., Kolesniková, L., Peña, I., et al. 2015, J. Mol. Spectr., 316, 84 [Google Scholar]
  7. Alonso, E. R., Kolesniková, L., Tercero, B., et al. 2016, ApJ, 832, 42 [NASA ADS] [CrossRef] [Google Scholar]
  8. Alonso, E. R., Kolesniková, L., Białkowska-Jaworska, E., et al. 2018, ApJ, 861, 70 [Google Scholar]
  9. Alonso, E. R., McGuire, B. A., Kolesniková, L., et al. 2019, ApJ, 883, 18 [NASA ADS] [CrossRef] [Google Scholar]
  10. Alonso, E. R., Kolesniková, L., Belloche, A., et al. 2021, A&A, 647, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Amano, T., & Maeda, A., 2000, J. Mol. Spectr., 203, 140 [Google Scholar]
  12. Amano, T., Hirao, T., & Takano, J., 2005, J. Mol. Spectr., 234, 170 [Google Scholar]
  13. Amano, T., Hashimoto, K., & Hirao, T., 2006, J. Mol. Struct., 795, 190 [Google Scholar]
  14. Apponi, A. J., Sun, M., Halfen, D. T., Ziurys, L. M., & Müller, H. S. P., 2008, ApJ, 673, 1240 [NASA ADS] [CrossRef] [Google Scholar]
  15. Arenas, B. E., Gruet, S., Steber, A. L., Giuliano, B. M., & Schnell, M. 2017a, Phys. Chem. Chem. Phys., 19, 1751 [Google Scholar]
  16. Arenas, B. E., Gruet, S., Steber, A. L., & Schnell, M., 2017b, J. Mol. Spectr., 337, 9 [Google Scholar]
  17. Arenas, B. E., Fatima, M., Pérez, C., et al. 2021, ApJ, 912, 90 [Google Scholar]
  18. August, J., Kroto, H. W., McNaughton, D., Phillips, K., & Walton, D. R. M., 1988, J. Mol. Spectr., 130, 424 [Google Scholar]
  19. Austin, J. A., Levy, D. H., Gottlieb, C. A., & Radford, H. E., 1974, J. Chem. Phys., 60, 207 [NASA ADS] [CrossRef] [Google Scholar]
  20. Baba, T., Tanaka, T., Morino, I., Yamada, K. M. T., & Tanaka, K., 1999, J. Chem. Phys., 110, 4131 [NASA ADS] [CrossRef] [Google Scholar]
  21. Babinec, P., & Leszczynski, J., 2000, J. Mol. Struct.: THEOCHEM, 501, 277 [Google Scholar]
  22. Bacmann, A., Taquet, V., Faure, A., Kahane, C., & Ceccarelli, C., 2012, A&A, 541, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Badawi, H. M., 2005, J. Mol. Struct.: THEOCHEM, 726, 253 [Google Scholar]
  24. Bak, B., & Svanholt, H., 1980, Chem. Phys. Lett., 75, 528 [NASA ADS] [CrossRef] [Google Scholar]
  25. Bak, B., Nielsen, O. J., & Svanholt, H., 1978, Chem. Phys. Lett., 59, 330 [Google Scholar]
  26. Bakri, B., Demaison, J., Margulès, L., & Møllendal, H., 2001, J. Mol. Spectr., 208, 92 [NASA ADS] [CrossRef] [Google Scholar]
  27. Ball, J. A., Gottlieb, C. A., Lilley, A. E., & Radford, H. E., 1970, ApJ, 162, L203 [NASA ADS] [CrossRef] [Google Scholar]
  28. Bane, M. K., Thompson, C. D., Robertson, E. G., Appadoo, D. R. T., & McNaughton, D., 2011, Phys. Chem. Chem. Phys., 13, 6793 [Google Scholar]
  29. Bane, M. K., Thompson, C. D., Appadoo, D. R. T., & McNaughton, D., 2012, J. Chem. Phys., 137, 084306 [Google Scholar]
  30. Barlow, K. R., & Tschumper, G. S., 2023, Int. J. Quant. Chem., 123, e27075 [Google Scholar]
  31. Bastiansen, O., Bakken, P., Kloster-Jensen, E., Samdal, S., & Traetteberg, M., 1995, J. Mol. Struct., 352, 77 [Google Scholar]
  32. Bates, J. B., & Smith, W., 1970, Spectrochim. Acta Part A, 26, 455 [Google Scholar]
  33. Bauer, A., & Bogey, M., 1970, C. R. Acad. Sc. Paris B, 271, 892 [Google Scholar]
  34. Bauer, A., & Burie, J., 1969, C. R. Acad. Sc. Paris B, 268, 800 [Google Scholar]
  35. Bauer, A., Bogey, M., & Maes, S., 1971, J. Phys. France, 32, 763 [Google Scholar]
  36. Bauer, A., Tarrago, G., & Remy, A., 1975, J. Mol. Spectr., 58, 111 [Google Scholar]
  37. Baughcum, S. L., Duerst, R. W., Rowe, W. F., Smith, Z., & Wilson, E. B., 1981, J. Am. Chem. Soc., 103, 6296 [CrossRef] [Google Scholar]
  38. Behnke, M., Medvedev, I., Winnewisser, M., De Lucia, F. C., & Herbst, E., 2004, ApJS, 152, 97 [NASA ADS] [CrossRef] [Google Scholar]
  39. Belloche, A., Menten, K. M., Comito, C., et al. 2008, A&A, 482, 179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Belloche, A., Garrod, R. T., Müller, H. S. P., et al. 2009, A&A, 499, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Belloche, A., Müller, H. S. P., Menten, K. M., Schilke, P., & Comito, C., 2013, A&A, 559, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Belloche, A., Müller, H. S. P., Garrod, R. T., & Menten, K. M., 2016, A&A, 587, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Belloche, A., Meshcheryakov, A. A., Garrod, R. T., et al. 2017, A&A, 601, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Belloche, A., Garrod, R. T., Müller, H. S. P., et al. 2019, A&A, 628, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Belloche, A., Garrod, R. T., Zingsheim, O., Müller, H. S. P., & Menten, K. M., 2022, A&A, 662, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Belov, S. P., Tretyakov, M. Y., Kozin, I. N., et al. 1998, J. Mol. Spectr., 191, 17 [NASA ADS] [CrossRef] [Google Scholar]
  47. Beltrán, M. T., Rivilla, V. M., Cesaroni, R., et al. 2021, A&A, 648, A100 [Google Scholar]
  48. Beltrán, M. T., Rivilla, V. M., Cesaroni, R., et al. 2022, A&A, 659, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Benidar, A., Bégué, D., Richter, F., et al. 2015, ChemPhysChem, 16, 848 [Google Scholar]
  50. Bermudez, C., Bailleux, S., & Cernicharo, J., 2017, A&A, 598, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Bester, M., 1984, PhD thesis, University of Cologne, Germany [Google Scholar]
  52. Bester, M., Tanimoto, M., Vowinkel, B., Winnewisser, G., & Yamada, K., 1983, Z. Naturforsch. A, 38, 64 [NASA ADS] [Google Scholar]
  53. Bester, M., Yamada, K., Winnewisser, G., et al. 1984, A&A, 137, L20 [NASA ADS] [Google Scholar]
  54. Bestmann, G., & Dreizler, H., 1982, Z. Naturforsch. A, 37, 58 [NASA ADS] [Google Scholar]
  55. Bestmann, G., Lalowski, W., & Dreizler, H., 1985a, Z. Naturforsch. A, 40, 271 [Google Scholar]
  56. Bestmann, G., Dreizler, H., Vacherand, J. M., et al. 1985b, Z. Naturforsch. A, 40, 508 [Google Scholar]
  57. Bevan, J., Legon, A., & Rego, C., 1992, J. Chem. Soc., Faraday Trans., 88, 3119 [Google Scholar]
  58. Binev, Y. I., Tsenov, J. A., Juchnovski, I. N., & Binev, I. G., 2003, J. Mol. Struct.: THEOCHEM, 625, 207 [Google Scholar]
  59. Birk, M., & Winnewisser, M., 1986, Chem. Phys. Lett., 123, 382 [Google Scholar]
  60. Bist, H. D., Brand, J. C. D., & Williams, D. R., 1967, J. Mol. Spectr., 24, 402 [Google Scholar]
  61. Bizzocchi, L., Cludi, L., & Degli Esposti, C., 2003, J. Mol. Spectr., 218, 53 [Google Scholar]
  62. Bizzocchi, L., Lattanzi, V., Laas, J., et al. 2017a, A&A, 602, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Bizzocchi, L., Tamassia, F., Laas, J., et al. 2017b, ApJS, 233, 11 [Google Scholar]
  64. Bizzocchi, L., Prudenzano, D., Rivilla, V. M., et al. 2020, A&A, 640, A98 [EDP Sciences] [Google Scholar]
  65. Bizzocchi, L., Tonolo, F., Giuliano, B. M., et al. 2024, ApJ, 970, 26 [Google Scholar]
  66. Blackburn, P. B., Brown, R. D., Burden, F. R., Crofts, J. G., & Gillard, I. R., 1970, Chem. Phys. Lett., 7, 102 [Google Scholar]
  67. Blackman, G. L., Brown, R. D., Godfrey, P. D., & Gunn, H. I., 1976, Nature, 261, 395 [NASA ADS] [CrossRef] [Google Scholar]
  68. Blake, G. A., Sastry, K. V. L. N., & de Lucia, F. C., 1984, J. Chem. Phys., 80, 95 [NASA ADS] [CrossRef] [Google Scholar]
  69. Blanco, S., Lesarri, A., López, J. C., & Alonso, J. L., 2004, J. Am. Chem. Soc., 126, 11675 [Google Scholar]
  70. Blanco, S., López, J. C., Lesarri, A., & Alonso, J. L., 2006, J. Am. Chem. Soc., 128, 12111 [Google Scholar]
  71. Blom, C., & Bauder, A., 1981, Chem. Phys. Lett., 82, 492 [NASA ADS] [CrossRef] [Google Scholar]
  72. Blom, C. E., & Günthard, H. H., 1981, Chem. Phys. Lett., 84, 267 [NASA ADS] [CrossRef] [Google Scholar]
  73. Blom, C., Grassi, G., & Bauder, A., 1984, J. Am. Chem. Soc., 106, 7427 [Google Scholar]
  74. Blukis, U., Kasai, P. H., & Myers, R. J., 1963, J. Chem. Phys., 38, 2753 [Google Scholar]
  75. Bocquet, R., Demaison, J., Poteau, L., et al. 1996, J. Mol. Spectr., 177, 154 [Google Scholar]
  76. Bogey, M., Demuynck, C., & Destombes, J. L., 1982, Chem. Phys., 66, 99 [Google Scholar]
  77. Bogey, M., Demuynck, C., & Destombes, J. L., 1984, Can. J. Phys., 62, 1248 [NASA ADS] [CrossRef] [Google Scholar]
  78. Bogey, M., Demuynck, C., & Destombes, J. L. 1986a, Chem. Phys., 102, 141 [Google Scholar]
  79. Bogey, M., Demuynck, C., & Destombes, J. L. 1986b, Chem. Phys. Lett., 125, 383 [NASA ADS] [CrossRef] [Google Scholar]
  80. Bogey, M., Demuynck, C., Destombes, J. L., & Dubus, H., 1987, J. Mol. Spectr., 122, 313 [NASA ADS] [CrossRef] [Google Scholar]
  81. Bogey, M., Cordonnier, M., Demuynck, C., & Destombes, J. L., 1992, ApJ, 399, L103 [NASA ADS] [CrossRef] [Google Scholar]
  82. Bogey, M., Demuynck, C., Destombes, J. L., & Vallee, Y., 1995, J. Mol. Spectr., 172, 344 [Google Scholar]
  83. Bogey, M., Civiš, S., Delcroix, B., et al. 1997, J. Mol. Spectr., 182, 85 [NASA ADS] [CrossRef] [Google Scholar]
  84. Bohn, R. K., Pardus, J. L., August, J., Brupbacher, T., & Jäger, W., 1997, J. Mol. Struct., 413, 293 [Google Scholar]
  85. Bolton, K., Owen, N., & Sheridan, J., 1970, Spectrochim. Acta Part A, 26, 909 [Google Scholar]
  86. Bolton, K., Lister, D. G., & Sheridan, J., 1974, J. Chem. Soc., Faraday Trans., 70, 113 [Google Scholar]
  87. Bonah, L., Zingsheim, O., Müller, H. S. P., et al. 2022, J. Mol. Spectr., 388, 111674 [Google Scholar]
  88. Bonah, L., Schlemmer, S., Guillemin, J.-C., Harding, M. E., & Thorwirth, S., 2024, J. Phys. Chem. A, 128, 4859 [Google Scholar]
  89. Bonfand, M., Belloche, A., Menten, K. M., Garrod, R. T., & Müller, H. S. P., 2017, A&A, 604, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Bonfand, M., Belloche, A., Garrod, R. T., et al. 2019, A&A, 628, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Bossa, J. B., Ordu, M. H., Müller, H. S. P., Lewen, F., & Schlemmer, S., 2014, A&A, 570, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Botschwina, P., 2000, J. Mol. Spectr., 203, 203 [Google Scholar]
  93. Botschwina, P., & Heyl, A., 1999, Mol. Phys., 97, 209 [Google Scholar]
  94. Botschwina, P., Horn, M., Flügge, J., & Seeger, S., 1993, J. Chem. Soc., Faraday Trans., 89, 2219 [NASA ADS] [CrossRef] [Google Scholar]
  95. Bouchez, A., Walters, A., Müller, H. S. P., et al. 2012, J. Quant. Spec. Radiat. Transf., 113, 1148 [NASA ADS] [CrossRef] [Google Scholar]
  96. Bouchy, A., 1972, PhD thesis, Université de Nancy, France [Google Scholar]
  97. Bouchy, A., Demaison, J., Roussy, G., & Barriol, J., 1973, J. Mol. Struct., 18, 211 [Google Scholar]
  98. Braakman, R., & Blake, G. A., 2010a, J. Mol. Spectr., 262, 100 [Google Scholar]
  99. Braakman, R., & Blake, G. A., 2010b, J. Mol. Spectr., 262, 93 [Google Scholar]
  100. Braakman, R., Drouin, B. J., Widicus Weaver, S. L., & Blake, G. A., 2010, J. Mol. Spectr., 264, 43 [Google Scholar]
  101. Brock, L. R., & Rohlfing, E. A., 1997, J. Chem. Phys., 106, 10048 [Google Scholar]
  102. Brown, R. D., Godfrey, P. D., & Storey, J., 1975, J. Mol. Spectr., 58, 445 [Google Scholar]
  103. Brown, R. D., Godfrey, P. D., & Kleibömer, B., 1985, J. Mol. Spectr., 114, 257 [Google Scholar]
  104. Brown, R. D., Godfrey, P. D., McNaughton, D., Pierlot, A. P., & Taylor, W. H., 1990a, J. Mol. Spectr., 140, 340 [Google Scholar]
  105. Brown, R. D., Godfrey, P. D., Kleibomer, B., Pierlot, A. P., & McNaughton, D., 1990b, J. Mol. Spectr., 142, 195 [NASA ADS] [CrossRef] [Google Scholar]
  106. Brünken, S., Müller, H. S. P., Lewen, F., & Winnewisser, G., 2003, Phys. Chem. Chem. Phys., 5, 1515 [Google Scholar]
  107. Brünken, S., Gottlieb, C. A., McCarthy, M. C., & Thaddeus, P. 2009a, ApJ, 697, 880 [CrossRef] [Google Scholar]
  108. Brünken, S., Yu, Z., Gottlieb, C. A., McCarthy, M. C., & Thaddeus, P. 2009b, ApJ, 706, 1588 [CrossRef] [Google Scholar]
  109. Bteich, S. B., Tercero, B., Cernicharo, J., et al. 2016, A&A, 592, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  110. Bteich, S., Goubet, M., Motiyenko, R. A., Margulès, L., & Huet, T. R., 2018, J. Mol. Spectr., 348, 124 [Google Scholar]
  111. Bunn, H., Soliday, R. M., Sumner, I., & Raston, P. L., 2017, ApJ, 847, 67 [Google Scholar]
  112. Burenin, A. V., Karyakin, E. N., Krupnov, A. F., Shapin, S. M., & Val’dov, A. N., 1981a, J. Mol. Spectr., 85, 1 [Google Scholar]
  113. Burenin, A. V., Val’dov, A. N., Karyakin, E. N., Krupnov, A. F., & Shapin, S. M., 1981b, J. Mol. Spectr., 87, 312 [Google Scholar]
  114. Burrell, P. M., Bjarnov, E., & Schwendeman, R. H., 1980, J. Mol. Spectr., 82, 193 [Google Scholar]
  115. Busch, L. A., Belloche, A., Garrod, R. T., Müller, H. S. P., & Menten, K. M., 2022, A&A, 665, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  116. Busch, L. A., Belloche, A., Garrod, R. T., Müller, H. S. P., & Menten, K. M., 2024, A&A, 681, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  117. Butcher, S. S., & Wilson, E. Bright, J., 1964, J. Chem. Phys., 40, 1671 [NASA ADS] [CrossRef] [Google Scholar]
  118. Butler, R. A. H., De Lucia, F. C., Petkie, D. T., et al. 2001, ApJS, 134, 319 [CrossRef] [Google Scholar]
  119. Cabezas, C., Kawashima, Y., Bermúdez, C., et al. 2023, Spectrochim. Acta Part A, 291, 122353 [Google Scholar]
  120. Calabrese, C., Vigorito, A., Feng, G., et al. 2014, J. Mol. Spectr., 295, 37 [Google Scholar]
  121. Calcutt, H., Jørgensen, J. K., Müller, H. S. P., et al. 2018a, A&A, 616, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  122. Calcutt, H., Fiechter, M. R., Willis, E. R., et al. 2018b, A&A, 617, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  123. Calcutt, H., Willis, E. R., Jørgensen, J. K., et al. 2019, A&A, 631, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Caminati, W., 1981, J. Mol. Spectr., 86, 193 [Google Scholar]
  125. Carles, S., Møllendal, H., & Guillemin, J. C., 2013, A&A, 558, A6 [Google Scholar]
  126. Carles, S., Møllendal, H., Trolez, Y., & Guillemin, J. C., 2014, A&A, 564, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  127. Carroll, P. B., Drouin, B. J., & Widicus Weaver, S. L., 2010, ApJ, 723, 845 [Google Scholar]
  128. Carroll, P. B., McGuire, B. A., Zaleski, D. P., et al. 2013, J. Mol. Spectr., 284, 21 [Google Scholar]
  129. Carvajal, M., Kleiner, I., & Demaison, J., 2010, ApJS, 190, 315 [NASA ADS] [CrossRef] [Google Scholar]
  130. Casado, J., Nygaard, L., & Sørensen, G., 1971, J. Mol. Struct., 8, 211 [Google Scholar]
  131. Caselli, P., Myers, P. C., & Thaddeus, P., 1995, ApJ, 455, L77 [Google Scholar]
  132. Castro, M. E., Niño, A., & Muñoz-Caro, C., 2008, Theor. Chem. Acc., 119, 343 [Google Scholar]
  133. Cazzoli, G., & Puzzarini, C., 2008, A&A, 487, 1197 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  134. Cazzoli, G., Cludi, L., Cotti, G., Esposti, C. D., & Dore, L., 1994, J. Mol. Spectr., 167, 468 [Google Scholar]
  135. Cazzoli, G., Puzzarini, C., & Lapinov, A. V., 2004, ApJ, 611, 615 [NASA ADS] [CrossRef] [Google Scholar]
  136. Cazzoli, G., Puzzarini, C., & Gauss, J., 2005, ApJS, 159, 181 [NASA ADS] [CrossRef] [Google Scholar]
  137. Cazzoli, G., Cludi, L., & Puzzarini, C., 2006, J. Mol. Struct., 780, 260 [CrossRef] [Google Scholar]
  138. Cazzoli, G., Cludi, L., Buffa, G., & Puzzarini, C., 2012, ApJS, 203, 11 [NASA ADS] [CrossRef] [Google Scholar]
  139. Ceccarelli, C., 2004, in Star Formation in the Interstellar Medium: In Honor of David Hollenbach, eds. D. Johnstone, F. C. Adams, D. N. C. Lin, D. A. Neufeeld, & E. C. Ostriker, Astronomical Society of the Pacific Conference Series, 323, 195 [Google Scholar]
  140. Cesaroni, R., Hofner, P., Araya, E., & Kurtz, S., 2010, A&A, 509, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  141. Cernicharo, J., Kisiel, Z., Tercero, B., et al. 2016, A&A, 587, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  142. Cernicharo, J., Lefloch, B., Agúndez, M., et al. 2018, ApJ, 853, L22 [CrossRef] [Google Scholar]
  143. Cernicharo, J., Marcelino, N., Agúndez, M., et al. 2020a, A&A, 642, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  144. Cernicharo, J., Marcelino, N., Agúndez, M., et al. 2020b, A&A, 642, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  145. Cernicharo, J., Agúndez, M., Cabezas, C., et al. 2021, A&A, 647, L2 [EDP Sciences] [Google Scholar]
  146. Cernicharo, J., Tercero, B., Agúndez, M., et al. 2024a, A&A, 686, A139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  147. Cernicharo, J., Cabezas, C., Agúndez, M., et al. 2024b, A&A, 686, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  148. Chahbazian, R., Martin-Drumel, M.-A., & Pirali, O., 2024, J. Phys. Chem. A, 128, 370 [NASA ADS] [CrossRef] [Google Scholar]
  149. Charles, S., Jones, G., Owen, N., Cyvin, S., & Cyvin, B., 1973, J. Mol. Struct., 16, 225 [Google Scholar]
  150. Charles, S. W., Cullen, F. C., Owen, N. L., & Williams, G. A., 1987, J. Mol. Struct., 157, 17 [Google Scholar]
  151. Charnley, S. B., Tielens, A. G. G. M., & Millar, T. J., 1992, ApJ, 399, L71 [Google Scholar]
  152. Chen, W., Bocquet, R., Wlodarczak, G., & Boucher, D., 1991, Int. J. Infrared Millim. Waves, 12, 987 [NASA ADS] [CrossRef] [Google Scholar]
  153. Cherniak, E. A., & Costain, C. C., 1966, J. Chem. Phys., 45, 104 [Google Scholar]
  154. Chitarra, O., Martin-Drumel, M.-A., Gans, B., et al. 2020, A&A, 644, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  155. Chomiak, D., Taleb-Bendiab, A., Civis, S., & Amano, T., 1994, Can. J. Phys., 72, 1078 [Google Scholar]
  156. Christiansen, J. J., 2005, J. Mol. Spectr., 231, 131 [NASA ADS] [CrossRef] [Google Scholar]
  157. Christen, D., & Müller, H. S. P., 2003, Phys. Chem. Chem. Phys., 5, 3600 [NASA ADS] [CrossRef] [Google Scholar]
  158. Christen, D., Coudert, L. H., Suenram, R. D., & Lovas, F. J., 1995, J. Mol. Spectr., 172, 57 [Google Scholar]
  159. Christen, D., Coudert, L. H., Larsson, J. A., & Cremer, D., 2001, J. Mol. Spectr., 205, 185 [Google Scholar]
  160. Clark, W. W., & De Lucia, F. C., 1976, J. Mol. Spectr., 60, 332 [Google Scholar]
  161. Cohen, E. A., & Pickett, H. M., 1982, J. Mol. Spectr., 93, 83 [Google Scholar]
  162. Collier, B., Krueger, K., Miller, I., et al. 2021, ApJS, 253, 40 [Google Scholar]
  163. Colzi, L., Rivilla, V. M., Beltrán, M. T., et al. 2021, A&A, 653, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  164. Colzi, L., Martín-Pintado, J., Zeng, S., et al. 2024, A&A, 690, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  165. Cordonnier, M., & Coudert, L. H., 1996, J. Mol. Spectr., 178, 59 [Google Scholar]
  166. Cornet, R. A., Landsberg, B. M., & Winnewisser, G., 1980, J. Mol. Spectr., 82, 253 [Google Scholar]
  167. Coudert, L. H., Çarçabal, P., Chevalier, M., et al. 2002, J. Mol. Spectr., 212, 203 [NASA ADS] [CrossRef] [Google Scholar]
  168. Coutens, A., Jørgensen, J. K., van der Wiel, M. H. D., et al. 2016, A&A, 590, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  169. Coutens, A., Zakharenko, O., Lewen, F., et al. 2019, A&A, 623, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  170. Coxon, J. A.,& Hajigeorgiou, P. G., 2004, J. Chem. Phys., 121, 2992 [NASA ADS] [CrossRef] [Google Scholar]
  171. Craig, N. C., Groner, P., Conrad, A. R., Gurusinghe, R., & Tubergen, M. J., 2016, J. Mol. Spectr., 328, 1 [Google Scholar]
  172. Creswell, R. A., Pearson, E. F., Winnewisser, M., & Winnewisser, G., 1976, Z. Naturforsch. A, 31, 221 [Google Scholar]
  173. Creswell, R. A., & Schwendeman, R. H., 1977a, J. Mol. Spectr., 64, 295 [Google Scholar]
  174. Creswell, R. A., Winnewisser, G., & Gerry, M. C. L., 1977b, J. Mol. Spectr., 65, 420 [NASA ADS] [CrossRef] [Google Scholar]
  175. Crowder, G., 1989, J. Mol. Struct.: THEOCHEM, 200, 235 [Google Scholar]
  176. Crowder, G. A., & Carlisle, G. O., 1991, J. Comput. Chem., 12, 880 [Google Scholar]
  177. Császár, A. G., 1989, Chem. Phys. Lett., 162, 361 [Google Scholar]
  178. Császár, A. G., 1996, J. Phys. Chem., 100, 3541 [Google Scholar]
  179. Cummins, S. E., Linke, R. A., & Thaddeus, P., 1986, ApJS, 60, 819 [NASA ADS] [CrossRef] [Google Scholar]
  180. Cunningham, George L., J., Boyd, A. W., Myers, R. J., Gwinn, W. D., & Le van, W. I. 1951, J. Chem. Phys., 19, 676 [Google Scholar]
  181. Curl, R. F., J. 1959, J. Chem. Phys., 30, 1529 [Google Scholar]
  182. Custer, T., Szczepaniak, U., Gronowski, M., et al. 2016, J. Phys. Chem. A, 120, 5928 [Google Scholar]
  183. Dalbouha, S., Senent, M. L., Komiha, N., & Domínguez-Gómez, R., 2016, J. Chem. Phys., 145, 124309 [Google Scholar]
  184. Daly, A. M., Bermúdez, C., Kolesniková, L., & Alonso, J. L., 2015, ApJS, 218, 30 [Google Scholar]
  185. De Lucia, F. C., & Helminger, P., 1975, J. Mol. Spectr., 54, 200 [Google Scholar]
  186. De Pree, C. G., Peters, T., Mac Low, M. M., et al. 2015, ApJ, 815, 123 [Google Scholar]
  187. De Vos, J., Schröder, B., & Rauhut, G., 2023, Spectrochim. Acta Part A, 302, 123083 [Google Scholar]
  188. de Vreede, J. P. M., Gillis, M. P. W., & Dijkerman, H. A., 1988, J. Mol. Spectr., 128, 509 [Google Scholar]
  189. de Zafra, R. L., 1971, ApJ, 170, 165 [NASA ADS] [CrossRef] [Google Scholar]
  190. Degli Esposti, C., Dore, L., & Bizzocchi, L., 2014, A&A, 565, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  191. Degli Esposti, C., Melosso, M., Bizzocchi, L., Tamassia, F., & Dore, L., 2020, J. Mol. Struct., 1203, 127429 [Google Scholar]
  192. DeLeon, R. L., & Muenter, J. S., 1985, J. Chem. Phys., 82, 1702 [NASA ADS] [CrossRef] [Google Scholar]
  193. Demaison, J., & Rudolph, H., 1975, J. Mol. Struct., 24, 325 [NASA ADS] [CrossRef] [Google Scholar]
  194. Demaison, J., Dubrulle, A., Boucher, D., Burie, J., & Typke, V., 1979, J. Mol. Spectr., 76, 1 [NASA ADS] [CrossRef] [Google Scholar]
  195. Demaison, J., Boucher, D., Burie, J., & Dubrulle, A., 1983, Z. Naturforsch. A, 38, 447 [Google Scholar]
  196. Demaison, J., Boucher, D., Burie, J., & Dubrulle, A., 1984, Z. Naturforsch. A, 39, 560 [Google Scholar]
  197. Demaison, J., Pohl, I., & Rudolph, H. D., 1985, J. Mol. Spectrosc., 114, 210 [Google Scholar]
  198. Demaison, J., Maes, H., Van Eijck, B. P., Wlodarczak, G., & Lasne, M. C., 1987, J. Mol. Spectr., 125, 214 [NASA ADS] [CrossRef] [Google Scholar]
  199. Demaison, J., Burie, J., Boucher, D., & Wlodarczak, G., 1991, J. Mol. Spectr., 146, 455 [Google Scholar]
  200. Dixon, T. A., & Woods, R. C., 1977, J. Chem. Phys., 67, 3956 [CrossRef] [Google Scholar]
  201. Dobrowolski, J. C., Ostrowski, S., Kołos, R., & Jamróz, M. H., 2008, Vib. Spectrosc., 48, 82 [CrossRef] [Google Scholar]
  202. Dore, L., Cazzoli, G., Civiš, S., & Scappini, F., 1995, Chem. Phys. Lett., 244, 145 [Google Scholar]
  203. Dore, L., Puzzarini, C., & Cazzoli, G., 2001, Can. J. Phys., 79, 359 [Google Scholar]
  204. Dore, L., Bizzocchi, L., Degli Esposti, C., & Tinti, F., 2009, A&A, 496, 275 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  205. Dore, L., Bizzocchi, L., Degli Esposti, C., & Gauss, J., 2010, J. Mol. Spectr., 263, 44 [NASA ADS] [CrossRef] [Google Scholar]
  206. Dore, L., Bizzocchi, L., & Degli Esposti, C., 2012, A&A, 544, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  207. Draper, G. R., & Werner, R. L., 1974, J. Mol. Spectr., 50, 369 [Google Scholar]
  208. Drouin, B. J., Pearson, J. C., Walters, A., & Lattanzi, V., 2006, J. Mol. Spectr., 240, 227 [Google Scholar]
  209. Drozdovskaya, M. N., van Dishoeck, E. F., Jørgensen, J. K., et al. 2018, MNRAS, 476, 4949 [Google Scholar]
  210. Drozdovskaya, M. N., van Dishoeck, E. F., Rubin, M., Jørgensen, J. K., & Altwegg, K., 2019, MNRAS, 490, 50 [Google Scholar]
  211. Dubrulle, A., Demaison, J., Burie, J., & Boucher, D., 1980, Z. Naturforsch. A, 35, 471 [NASA ADS] [Google Scholar]
  212. Duncan, J. L., Ferguson, A. M., Harper, J., Tonge, K. H., & Hegelund, F., 1987, J. Mol. Spectr., 122, 72 [Google Scholar]
  213. Durig, J. R., & Compton, D. A. C., 1978, J. Chem. Phys., 69, 4713 [Google Scholar]
  214. Durig, J. R., Tong, C. K., Hawley, C. W., & Bragin, J., 1971, J. Phys. Chem., 75, 44 [Google Scholar]
  215. Durig, J., Guirgis, G., & Bell, S., 1989, J. Phys. Chem., 93, 3487 [Google Scholar]
  216. Durig, J. R., Guirgis, G. A., & Drew, A. S., 1994, J. Raman Spectrosc., 25, 907 [Google Scholar]
  217. Durig, J. R., Jin, Y., Phan, H., Liu, J., & Durig, D. T., 2002, Struct. Chem., 13, 1 [Google Scholar]
  218. Durig, J. R., Zhou, X., Durig, N. E., Nguyen, D., & Durig, D. T., 2009, J. Mol. Struct., 917, 37 [Google Scholar]
  219. Durig, J. R., Zhou, S. X., Zheng, C., & Durig, D. T., 2010, J. Mol. Struct., 971, 23 [NASA ADS] [CrossRef] [Google Scholar]
  220. Durig, J. R., Deeb, H., Darkhalil, I. D., et al. 2011, J. Mol. Struct., 985, 202 [NASA ADS] [CrossRef] [Google Scholar]
  221. Durig, J. R., Klaassen, J. J., Darkhalil, I. D., et al. 2012a, J. Mol. Struct., 1009, 30 [Google Scholar]
  222. Durig, J. R., Darkhalil, I. D., Klaassen, J. J., et al. 2012b, J. Raman Spectrosc., 43, 1329 [Google Scholar]
  223. Ebenstein, W. L., & Muenter, J. S., 1984, J. Chem. Phys., 80, 3989 [Google Scholar]
  224. Ehrlichmann, H., Grabow, J. U., Dreizler, H., et al. 1989, Z. Naturforsch. A, 44, 751 [Google Scholar]
  225. Ellingsen, B. H., Marstokk, K. M., & Møllendal, H., 1978, J. Mol. Struct., 48, 9 [Google Scholar]
  226. Endo, Y., & Nakajima, M., 2014, J. Mol. Spectr., 301, 15 [Google Scholar]
  227. Endo, Y., Saito, S., & Hirota, E., 1984, J. Chem. Phys., 81, 122 [Google Scholar]
  228. Endo, Y., Saito, S., & Hirota, E., 1985, J. Chem. Phys., 83, 2026 [Google Scholar]
  229. Endres, C. P., Drouin, B. J., Pearson, J. C., et al. 2009, A&A, 504, 635 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  230. Endres, C. P., Schlemmer, S., Schilke, P., Stutzki, J., & Müller, H. S. P., 2016, J. Mol. Spectrosc., 327, 95 [NASA ADS] [CrossRef] [Google Scholar]
  231. Etxaluze, M., Goicoechea, J. R., Cernicharo, J., et al. 2013, A&A, 556, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  232. Fabricant, B., Krieger, D., & Muenter, J. S., 1977, J. Chem. Phys., 67, 1576 [NASA ADS] [CrossRef] [Google Scholar]
  233. Fatima, M., Müller, H. S. P., Zingsheim, O., et al. 2023, A&A, 680, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  234. Fayt, A., Vigouroux, C., Willaert, F., et al. 2004, J. Mol. Struct., 695, 295 [NASA ADS] [Google Scholar]
  235. Fernández, J. M., Tejeda, G., Carvajal, M., & Senent, M. L., 2019, ApJS, 241, 13 [CrossRef] [Google Scholar]
  236. Fischer, E., & Botskor, I., 1982, J. Mol. Spectr., 91, 116 [NASA ADS] [CrossRef] [Google Scholar]
  237. Fischer, E., & Botskor, I., 1984, J. Mol. Spectr., 104, 226 [Google Scholar]
  238. Fisher, J., Paciga, G., Xu, L.-H., et al. 2007, J. Mol. Spectr., 245, 7 [Google Scholar]
  239. Flaud, J. M., Lafferty, W. J., Kwabia Tchana, F., Perrin, A., & Landsheere, X., 2012, J. Mol. Spectr., 271, 38 [Google Scholar]
  240. Fliege, E., Bestmann, G., Schwarz, R., & Dreizler, H., 1981, Z. Naturforsch. A, 36, 1124 [Google Scholar]
  241. Flygare, W. H., & Lowe, J. T., 1965, J. Chem. Phys., 43, 3645 [Google Scholar]
  242. Fortenberry, R. C., Huang, X., Francisco, J. S., Crawford, T. D., & Lee, T. J., 2012, J. Chem. Phys., 136 [Google Scholar]
  243. Fried, Z. T. P., El-Abd, S. J., Hays, B. M., et al. 2024, ApJ, 965, L23 [Google Scholar]
  244. Fuchs, U., Brünken, S., Fuchs, G., et al. 2004, Z. Naturforsch. A, 59, 861 [NASA ADS] [CrossRef] [Google Scholar]
  245. Fuentetaja, R., Bermúdez, C., Cabezas, C., et al. 2023, A&A, 671, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  246. Fuentetaja, R., Agúndez, M., Cabezas, C., et al. 2024, A&A, 688, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  247. Fusina, L., Di Lonardo, G., Johns, J. W. C., & Halonen, L., 1988, J. Mol. Spectr., 127, 240 [Google Scholar]
  248. Gabrys, C., Uy, D., Jagod, M.-F., Oka, T., & Amano, T. 1995, J. Phys. Chem., 99, 15611 [Google Scholar]
  249. Gadhi, J., Lahrouni, A., Legrand, J., & Demaison, J., 1995, J. Chim. Phys., 92, 1984 [Google Scholar]
  250. García de la Concepción, J., Colzi, L., Jiménez-Serra, I., et al. 2022, A&A, 658, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  251. Gardner, F. F., Godfrey, P. D., & Williams, D. R., 1980, MNRAS, 193, 713 [Google Scholar]
  252. Garrod, R. T., 2013, ApJ, 765, 60 [Google Scholar]
  253. Garrod, R. T., & Herbst, E., 2006, A&A, 457, 927 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  254. Garrod, R. T., & Pauly, T., 2011, ApJ, 735, 15 [Google Scholar]
  255. Garrod, R. T., & Herbst, E., 2023, Faraday Discuss., 245, 541 [NASA ADS] [CrossRef] [Google Scholar]
  256. Garrod, R. T., Widicus Weaver, S. L., & Herbst, E., 2008, ApJ, 682, 283 [Google Scholar]
  257. Garrod, R. T., Jin, M., Matis, K. A., et al. 2022, ApJS, 259, 1 [NASA ADS] [CrossRef] [Google Scholar]
  258. Gaume, R. A., Claussen, M. J., de Pree, C. G., Goss, W. M., & Mehringer, D. M., 1995, ApJ, 449, 663 [Google Scholar]
  259. Giesen, T. F., Harding, M. E., Gauss, J., Grabow, J.-U., & Müller, H. S. P., 2020, J. Mol. Spectr., 371, 111303 [Google Scholar]
  260. Godfrey, P., Firth, S., Hatherley, L., Brown, R., & Pierlot, A., 1993, J. Am. Chem. Soc., 115, 9687 [Google Scholar]
  261. Godon, M., 1974, PhD thesis, Université de Lille, France [Google Scholar]
  262. Golby, J. A., & Butcher, R. J., 1984, J. Mol. Spectr., 107, 292 [Google Scholar]
  263. Golubiatnikov, G. Y., Lapinov, A. V., Guarnieri, A., & Knöchel, R., 2005, J. Mol. Spectr., 234, 190 [Google Scholar]
  264. Gong, J. Z., Matthews, D. A., Changala, P. B., & Stanton, J. F., 2018, J. Chem. Phys., 149, 114102 [Google Scholar]
  265. Goorvitch, D., 1994, ApJS, 95, 535 [NASA ADS] [CrossRef] [Google Scholar]
  266. Gottlieb, C. A., Gottlieb, E. W., Thaddeus, P., & Kawamura, H., 1983, ApJ, 275, 916 [CrossRef] [Google Scholar]
  267. Gottlieb, C. A., Apponi, A. J., McCarthy, M. C., Thaddeus, P., & Linnartz, H., 2000, J. Chem. Phys., 113, 1910 [NASA ADS] [CrossRef] [Google Scholar]
  268. Grabow, J. U., & Stahl, W., 1990, Z. Naturforsch. A, 45, 1043 [Google Scholar]
  269. Grabska, J., Beć, K. B., Ozaki, Y., & Huck, C. W., 2017, J. Phys. Chem. A, 121, 1950 [Google Scholar]
  270. Gratier, P., Majumdar, L., Ohishi, M., et al. 2016, ApJS, 225, 25 [Google Scholar]
  271. Gregersen, E. M., & Evans, Neal J. I. 2001, ApJ, 553, 1042 [NASA ADS] [CrossRef] [Google Scholar]
  272. Gripp, J., Guarnieri, A., Stahl, W., & Lentz, D., 2000, J. Mol. Struct., 526, 81 [NASA ADS] [CrossRef] [Google Scholar]
  273. Groner, P., Gunde, R., Ha, T. K., & Bauder, A., 1987, Chem. Phys. Lett., 139, 479 [NASA ADS] [CrossRef] [Google Scholar]
  274. Guarnieri, A., & Huckauf, A., 2003, Z. Naturforsch. A, 58, 275 [Google Scholar]
  275. Guarnieri, A., Hinze, R., Krüger, M., et al. 1992, J. Mol. Spectr., 156, 39 [NASA ADS] [CrossRef] [Google Scholar]
  276. Guélin, M., Friberg, P., & Mezaoui, A., 1982, A&A, 109, 23 [Google Scholar]
  277. Gutowsky, H., & Germann, T., 1991, J. Mol. Spectr., 147, 91 [Google Scholar]
  278. Haekel, J., & Mäder, H., 1988, Z. Naturforsch. A, 43, 1111 [Google Scholar]
  279. Haekel, J., & Mäder, H., 1989, J. Quant. Spec. Radiat. Transf., 41, 9 [Google Scholar]
  280. Hahn, M., Bodenseh, H. K., & Ferner, M., 2004, J. Mol. Spectr., 223, 138 [Google Scholar]
  281. Halasa, A., Lapinski, L., Reva, I., et al. 2014, J. Phys. Chem. A, 118, 5626 [Google Scholar]
  282. Hamada, Y., Hashiguchi, K., Hirakawa, A. Y., et al. 1983, J. Mol. Spectr., 102, 123 [Google Scholar]
  283. Hamada, Y., Hashiguchi, K., Tsuboi, M., Koga, Y., & Kondo, S., 1984, J. Mol. Spectr., 105, 70 [Google Scholar]
  284. Hardy, J. A., Cox, A. P., Fliege, E., & Dreizler, H., 1982, Z. Naturforsch. A, 37, 1035 [NASA ADS] [CrossRef] [Google Scholar]
  285. Havenith, M., Zwart, E., Leo Meerts, W., & Ter Meulen, J. J., 1990, J. Chem. Phys., 93, 8446 [NASA ADS] [CrossRef] [Google Scholar]
  286. Haykal, I., Motiyenko, R. A., Margulès, L., & Huet, T. R., 2013, A&A, 549, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  287. Haykal, I., Carvajal, M., Tercero, B., et al. 2014, A&A, 568, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  288. Hays, B. M., Wehres, N., DePrince, B. A., et al. 2015, Chem. Phys. Lett., 630, 18 [Google Scholar]
  289. Heineking, N., Grabow, J. U., & Stahl, W. 1994, Mol. Phys., 81, 1177 [Google Scholar]
  290. Heise, H. M., Mäder, H., & Dreizler, H., 1976, Z. Naturforsch. A, 31, 1228 [Google Scholar]
  291. Heise, H. M., Winther, F., & Lutz, H., 1981, J. Mol. Spectr., 90, 531 [NASA ADS] [CrossRef] [Google Scholar]
  292. Helminger, P. A., & De Lucia, F. C., 1985, J. Mol. Spectr., 111, 66 [NASA ADS] [CrossRef] [Google Scholar]
  293. Henfrey, N. F., & Thrush, B. A., 1985, J. Mol. Spectr., 113, 426 [Google Scholar]
  294. Herbst, E., & Leung, C. M., 1986, MNRAS, 222, 689 [Google Scholar]
  295. Herbst, E., & van Dishoeck, E. F., 2009, ARA&A, 47, 427 [NASA ADS] [CrossRef] [Google Scholar]
  296. Herschbach, D. R., & Swalen, J. D., 1958, J. Chem. Phys., 29, 761 [NASA ADS] [CrossRef] [Google Scholar]
  297. Hinze, R., Zerbe-Foese, H., Doose, J., & Guarnieri, A., 1996, J. Mol. Spectr., 176, 133 [Google Scholar]
  298. Hirahara, Y., Ohshima, Y., & Endo, Y., 1993, ApJ, 403, L83 [CrossRef] [Google Scholar]
  299. Hirao, T., Okabayashi, T., & Tanimoto, M., 2001, J. Mol. Spectr., 208, 148 [Google Scholar]
  300. Hirata, Y., Kubota, S., Watanabe, S., Momose, T., & Kawaguchi, K., 2008, J. Mol. Spectr., 251, 314 [Google Scholar]
  301. Hirono, H., Shibano, J., Mitani, A., et al. 1987, J. Mol. Struct., 162, 359 [Google Scholar]
  302. Hirota, E., 1966, J. Chem. Phys., 45, 1984 [Google Scholar]
  303. Hirota, E., 1979, J. Phys. Chem., 83, 1457 [Google Scholar]
  304. Hirota, E., Sugisaki, R., Nielsen, C. J., & Sørensen, G. O., 1974, J. Mol. Spectr., 49, 251 [NASA ADS] [CrossRef] [Google Scholar]
  305. Hocking, W. H., 1976, Z. Naturforsch. A, 31, 1113 [NASA ADS] [CrossRef] [Google Scholar]
  306. Hocking, W. H., & Winnewisser, G., 1976a, Z. Naturforsch. A, 31, 422 [NASA ADS] [CrossRef] [Google Scholar]
  307. Hocking, W. H., & Winnewisser, G., 1976b, Z. Naturforsch. A, 31, 995 [Google Scholar]
  308. Hocking, W. H., Gerry, M. C. L., & Winnewisser, G., 1975, Can. J. Phys., 53, 1869 [NASA ADS] [CrossRef] [Google Scholar]
  309. Hoeft, J., Tiemann, E., & Törring, T., 1972, Z. Naturforsch. A, 27, 703 [Google Scholar]
  310. Hollenstein, H., Schär, R., Schwizgebel, N., Grassi, G., & Günthard, H., 1983, Spectrochim. Acta Part A, 39, 193 [Google Scholar]
  311. Hollis, J. M., Lovas, F. J., & Jewell, P. R., 2000, ApJ, 540, L107 [Google Scholar]
  312. Hopkinson, M. J., Kroto, H. W., Nixon, J. F., & Simmons, N. P. C., 1976, Chem. Phys. Lett., 42, 460 [Google Scholar]
  313. Horn, M., Oswald, M., Oswald, R., & Botschwina, P., 1995, Ber. Bunsenges. Phys. Chem., 99, 323 [Google Scholar]
  314. Hsu, S. L., & Flygare, W. H., 1971, J. Mol. Spectr., 37, 92 [Google Scholar]
  315. Hu, C. H., Shen, M., & Schaefer III, H. F. 1993, J. Am. Chem. Soc., 115, 2923 [Google Scholar]
  316. Huang, E., Liu, S., Zhou, X., et al. 2025, J. Mol. Struct., 1321, 139871 [Google Scholar]
  317. Hutchinson, M., Kroto, H. W., & Walton, D. R. M., 1980, J. Mol. Spectr., 82, 394 [Google Scholar]
  318. Ilyushin, V. V., Alekseev, E. A., Dyubko, S. F., Kleiner, I., & Hougen, J. T., 2004, J. Mol. Spectrosc., 227, 115 [Google Scholar]
  319. Ilyushin, V. V., Alekseev, E. A., Dyubko, S. F., Motiyenko, R. A., & Lovas, F. J., 2005, J. Mol. Spectr., 231, 15 [Google Scholar]
  320. Ilyushin, V., Kryvda, A., & Alekseev, E., 2009, J. Mol. Spectr., 255, 32 [Google Scholar]
  321. Ilyushin, V. V., Endres, C. P., Lewen, F., Schlemmer, S., & Drouin, B. J., 2013, J. Mol. Spectr., 290, 31 [NASA ADS] [CrossRef] [Google Scholar]
  322. Ilyushin, V. V., Zakharenko, O., Lewen, F., et al. 2020, Can. J. Phys., 98, 530 [NASA ADS] [CrossRef] [Google Scholar]
  323. Ilyushin, V. V., Margulès, L., Tercero, B., et al. 2021, J. Mol. Spectr., 379, 111454 [Google Scholar]
  324. Ilyushin, V. V., Müller, H. S. P., Drozdovskaya, M. N., et al. 2024, A&A, 687, A220 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  325. Immer, K., Li, J., Quiroga-Nuñez, L. H., et al. 2019, A&A, 632, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  326. Indriolo, N., Neufeld, D. A., Gerin, M., et al. 2015, ApJ, 800, 40 [NASA ADS] [CrossRef] [Google Scholar]
  327. Irvine, W. M., & Schloerb, F. P., 1984, ApJ, 282, 516 [Google Scholar]
  328. Jabri, A., Van, V., Nguyen, H. V. L., et al. 2016, A&A, 589, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  329. Jackels, C. F., 1982, J. Chem. Phys., 76, 505 [Google Scholar]
  330. Jaman, A. I., Chakraborty, S., & Chakraborty, R., 2015, J. Mol. Struct., 1079, 402 [Google Scholar]
  331. Jiménez-Serra, I., Rodríguez-Almeida, L. F., Martín-Pintado, J., et al. 2022, A&A, 663, A181 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  332. Jiménez-Serra, I., Codella, C., & Belloche, A., 2025, arXiv e-prints [arXiv:2503.17104] [Google Scholar]
  333. Jin, M., & Garrod, R. T., 2020, ApJS, 249, 26 [Google Scholar]
  334. Johns, J. W. C., Stone, J. M. R., & Winnewisser, G., 1972, J. Mol. Spectr., 42, 523 [Google Scholar]
  335. Johns, J. W. C., Nemes, L., Yamada, K. M. T., et al. 1992, J. Mol. Spectr., 156, 501 [Google Scholar]
  336. Johnson, H. R., & Strandberg, M. W. P., 1952, J. Chem. Phys., 20, 687 [Google Scholar]
  337. Johnson, R. D., & Hudgens, J. W., 1996, J. Phys. Chem., 100, 19874 [Google Scholar]
  338. Johnson, M., Alexander, M., Hertel, I., & Lineberger, W., 1984, Chem. Phys. Lett., 105, 374 [Google Scholar]
  339. Johnson, T. J., Sams, R. L., Profeta, L. T. M., et al. 2013, J. Phys. Chem. A, 117, 4096 [Google Scholar]
  340. Jørgensen, J. K., van der Wiel, M. H. D., Coutens, A., et al. 2016, A&A, 595, A117 [Google Scholar]
  341. Jørgensen, J. K., Müller, H. S. P., Calcutt, H., et al. 2018, A&A, 620, A170 [Google Scholar]
  342. Jørgensen, J. K., Belloche, A., & Garrod, R. T., 2020, ARA&A, 58, 727 [Google Scholar]
  343. Kakar, R. K., & Quade, C. R., 1980, J. Chem. Phys., 72, 4300 [Google Scholar]
  344. Kasten, W., & Dreizler, H., 1986a, Z. Naturforsch. A, 41, 637 [Google Scholar]
  345. Kasten, W., & Dreizler, H., 1986b, Z. Naturforsch. A, 41, 1173 [NASA ADS] [Google Scholar]
  346. Kaushik, V. K., & Woods, R. C., 1982, Z. Phys. Chem., 132, 117 [CrossRef] [Google Scholar]
  347. Kawashima, Y., Usami, T., Suenram, R. D., Golubiatnikov, G. Y., & Hirota, E., 2010, J. Mol. Spectr., 263, 11 [Google Scholar]
  348. Kawiecki, R., Devlin, F., Stephens, P., & Amos, R., 1991, J. Phys. Chem., 95, 9817 [Google Scholar]
  349. Keceli, M., Shiozaki, T., Yagi, K., & Hirata, S., 2009, Mol. Phys., 107, 1283 [CrossRef] [Google Scholar]
  350. Kewley, R., 1974, Can. J. Chem., 52, 509 [Google Scholar]
  351. Khlifi, M., & Raulin, F., 1991, Spectrochim. Acta Part A, 47, 171 [Google Scholar]
  352. Khlifi, M., Nollet, M., Paillous, P., et al. 1999, J. Mol. Spectr., 194, 206 [NASA ADS] [CrossRef] [Google Scholar]
  353. Kisiel, Z., 1985, Chem. Phys. Lett., 118, 334 [Google Scholar]
  354. Kisiel, Z., Białkowska-Jaworska, E., Desyatnyk, O., Pietrewicz, B. A., & Pszczółkowski, L., 2001, J. Mol. Spectr., 208, 113 [Google Scholar]
  355. Kisiel, Z., Martin-Drumel, M.-A., & Pirali, O., 2015, J. Mol. Spectr., 315, 83 [NASA ADS] [CrossRef] [Google Scholar]
  356. Kisiel, Z., Pszczółkowski, L., Białkowska-Jaworska, E., & Charnley, S. B., 2016, J. Mol. Spectr., 321, 13 [Google Scholar]
  357. Kisiel, Z., Kolesniková, L., Belloche, A., et al. 2022, A&A, 657, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  358. Klapper, G., Lewen, F., Gendriesch, R., Belov, S. P., & Winnewisser, G., 2000, J. Mol. Spectr., 201, 124 [Google Scholar]
  359. Klapper, G., Lewen, F., Gendriesch, R., Belov, S. P., & Winnewisser, G., 2001, Z. Naturforsch. A, 56, 329 [Google Scholar]
  360. Klapper, G., Surin, L., Lewen, F., et al. 2003, ApJ, 582, 262 [NASA ADS] [CrossRef] [Google Scholar]
  361. Klaus, T., Saleck, A. H., Belov, S. P., et al. 1996, J. Mol. Spectr., 180, 197 [NASA ADS] [CrossRef] [Google Scholar]
  362. Kleimeier, N. F., Eckhardt, A. K., & Kaiser, R. I., 2021, J. Am. Chem. Soc., 143, 14009 [CrossRef] [Google Scholar]
  363. Klisch, E., Klaus, T., Belov, S. P., Winnewisser, G., & Herbst, E., 1995, A&A, 304, L5 [NASA ADS] [Google Scholar]
  364. Kobayashi, K., Takamura, K., Sakai, Y., et al. 2013, ApJS, 205, 9 [Google Scholar]
  365. Kobayashi, K., Sakai, Y., Fujitake, M., et al. 2020, Can. J. Phys., 98, 551 [NASA ADS] [CrossRef] [Google Scholar]
  366. Koerber, M., Bisschop, S. E., Endres, C. P., et al. 2013, A&A, 558, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  367. Kolesniková, L., Daly, A. M., Alonso, J. L., Tercero, B., & Cernicharo, J., 2013, J. Mol. Spectr., 289, 13 [Google Scholar]
  368. Kolesniková, L., Tercero, B., Cernicharo, J., et al. 2014, ApJ, 784, L7 [Google Scholar]
  369. Kolesniková, L., Alonso, J. L., Bermúdez, C., et al. 2016, A&A, 591, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  370. Kolesniková, L., Alonso, E. R., Mata, S., & Alonso, J. L. 2017a, ApJS, 229, 26 [Google Scholar]
  371. Kolesniková, L., Alonso, E. R., Mata, S., Cernicharo, J., & Alonso, J. L. 2017b, ApJS, 233, 24 [Google Scholar]
  372. Kolesniková, L., Alonso, E. R., Tercero, B., Cernicharo, J., & Alonso, J. L. 2018a, A&A, 616, A173 [Google Scholar]
  373. Kolesniková, L., Peña, I., Alonso, E. R., et al. 2018b, A&A, 619, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  374. Kolesniková, L., Belloche, A., Koucký, J., et al. 2022, A&A, 659, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  375. Kołos, R., & Sobolewski, A. L., 2001, Chem. Phys. Lett., 344, 625 [Google Scholar]
  376. Kondo, S., & Hirota, E., 1970, J. Mol. Spectr., 34, 97 [Google Scholar]
  377. Koput, J., 1986, J. Mol. Spectr., 115, 131 [Google Scholar]
  378. Köroğlu, B., Loparo, Z., Nath, J., Peale, R. E., & Vasu, S. S., 2015, J. Quant. Spec. Radiat. Transf., 152, 107 [Google Scholar]
  379. Korppi-Tommola, J., 1978, Spectrochim. Acta Part A, 34, 1077 [Google Scholar]
  380. Koucký, J., Kolesniková, L., Luková, K., et al. 2022, A&A, 666, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  381. Kowal, A. T., 2006, J. Chem. Phys., 124, 014304 [Google Scholar]
  382. Kraśnicki, A., Kisiel, Z., Jabs, W., Winnewisser, B. P., & Winnewisser, M., 2011, J. Mol. Spectr., 267, 144 [CrossRef] [Google Scholar]
  383. Krause, H., & Sutter, D. H., 1989, Z. Naturforsch. A, 44, 1063 [Google Scholar]
  384. Kretschmer, U., Consalvo, D., Knaack, A., et al. 1996, Mol. Phys., 87, 1159 [Google Scholar]
  385. Krisher, L. C., & Saegebarth, E., 1971, J. Chem. Phys., 54, 4553 [Google Scholar]
  386. Krüger, M., Dreizler, H., Preugschat, D., & Lentz, D., 1991, Angew. Chem. Int. Ed. Engl., 30, 1644 [CrossRef] [Google Scholar]
  387. Krüger, M., Stahl, W., & Dreizler, H., 1993, J. Mol. Spectr., 158, 298 [CrossRef] [Google Scholar]
  388. Kryvda, A. V., Gerasimov, V. G., Dyubko, S. F., Alekseev, E. A., & Motiyenko, R. A., 2009, J. Mol. Spectr., 254, 28 [NASA ADS] [CrossRef] [Google Scholar]
  389. Kuczkowski, R. L., Lovas, F. J., Suenram, R. D., et al. 1981, J. Mol. Struct., 72, 143 [NASA ADS] [CrossRef] [Google Scholar]
  390. Kukolich, S. G., 1982, J. Chem. Phys., 76, 97 [NASA ADS] [CrossRef] [Google Scholar]
  391. Kukolich, S. G., & Nelson, A. C., 1971, Chem. Phys. Lett., 11, 383 [NASA ADS] [CrossRef] [Google Scholar]
  392. Kukolich, S., Nelson, A., & Yamanashi, B., 1971, J. Am. Chem. Soc., 93, 6769 [Google Scholar]
  393. Kukolich, S. G., Oates, D. E., & Wang, J. H. S., 1974, J. Chem. Phys., 61, 4686 [Google Scholar]
  394. Kurland, R. J., & Wilson, E. Bright, J., 1957, J. Chem. Phys., 27, 585 [Google Scholar]
  395. Kutzer, P., Weismann, D., Waßmuth, B., et al. 2016, J. Mol. Spectr., 329, 28 [NASA ADS] [CrossRef] [Google Scholar]
  396. Kuze, H., Kuga, T., & Shimizu, T., 1982, J. Mol. Spectr., 93, 248 [Google Scholar]
  397. Kuzmin, V., Novikov, V., Sagitova, E., et al. 2021, J. Mol. Struct., 1243, 130847 [Google Scholar]
  398. Kwon, C. H., Kim, H. L., & Kim, M. S., 2003, J. Phys. Chem. A, 107, 10969 [Google Scholar]
  399. Lafferty, W. J., 1968, J. Mol. Spectr., 25, 359 [Google Scholar]
  400. Landsberg, B. M., & Suenram, R. D., 1983, J. Mol. Spectr., 98, 210 [Google Scholar]
  401. Lapinov, A. V., Golubiatnikov, G. Y., Markov, V. N., & Guarnieri, A., 2007, Astron. Lett., 33, 121 [Google Scholar]
  402. Larsen, N. W., 1979, J. Mol. Struct., 51, 175 [Google Scholar]
  403. Larsen, N. W., & Nicolaisen, F. M., 1974, J. Mol. Struct., 22, 29 [Google Scholar]
  404. Larsen, N. W., & Winnewisser, B. P., 1974, Z. Naturforsch. A, 29, 1213 [NASA ADS] [Google Scholar]
  405. Lattanzi, V., Walters, A., Drouin, B. J., & Pearson, J. C., 2007, ApJ, 662, 771 [NASA ADS] [CrossRef] [Google Scholar]
  406. Laurie, V. W., 1961, J. Chem. Phys., 34, 1516 [NASA ADS] [CrossRef] [Google Scholar]
  407. Le Guennec, M., Wlodarczak, G., Burie, J., & Demaison, J., 1992, J. Mol. Spectr., 154, 305 [Google Scholar]
  408. Le Petit, F., Ruaud, M., Bron, E., et al. 2016, A&A, 585, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  409. Lee, T. J., & Schaefer, Henry F., I. 1986, J. Chem. Phys., 85, 3437 [NASA ADS] [CrossRef] [Google Scholar]
  410. Lees, R. M., 1971, Can. J. Phys., 49, 367 [Google Scholar]
  411. Lees, R. M., Reid, E. M., Xu, L.-H., & Billinghurst, B. E. 2020, Can. J. Phys., 98, 519 [Google Scholar]
  412. Lesarri, A. G., Cosleou, J., Li, X., Wlodarczak, G., & Demaison, J., 1995, J. Mol. Spectr., 172, 520 [Google Scholar]
  413. Lewis-Bevan, W., Gaston, R. D., Tyrrell, J., Stork, W. D., & Salmon, G. L., 1992, J. Am. Chem. Soc., 114, 1933 [Google Scholar]
  414. Li, J., Wang, J., Lu, X., et al. 2021, ApJ, 919, 4 [NASA ADS] [CrossRef] [Google Scholar]
  415. Lide, D. R.Jr., 1960, J. Chem. Phys., 33, 1514 [Google Scholar]
  416. Lide, D. R.Jr., & Mann, D. E., 1957, J. Chem. Phys., 27, 868 [Google Scholar]
  417. Lignell, A., Osadchuk, I., Räsänen, M., & Lundell, J., 2021, ApJ, 917, 68 [NASA ADS] [CrossRef] [Google Scholar]
  418. Ligterink, N. F. W., Coutens, A., Kofman, V., et al. 2017, MNRAS, 469, 2219 [Google Scholar]
  419. Ligterink, N. F. W., Calcutt, H., Coutens, A., et al. 2018a, A&A, 619, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  420. Ligterink, N. F. W., Terwisscha van Scheltinga, J., Taquet, V., et al. 2018b, MNRAS, 480, 3628 [Google Scholar]
  421. Lindenmaier, R., Tipton, N., Sams, R. L., et al. 2016, J. Phys. Chem. A, 120, 5993 [Google Scholar]
  422. Lique, F., Cernicharo, J., & Cox, P., 2006, ApJ, 653, 1342 [NASA ADS] [CrossRef] [Google Scholar]
  423. Lockley, T. J. L., Hearn, J. P. I., King, A. K., & Howard, B. J., 2002, J. Mol. Struct., 612, 199 [Google Scholar]
  424. Lopes, S., Nikitin, T., & Fausto, R., 2023, J. Photochem. Photobiol. A Chem., 438, 114487 [Google Scholar]
  425. López-Gallifa, Á., Rivilla, V. M., Beltrán, M. T., et al. 2024, MNRAS, 529, 3244 [CrossRef] [Google Scholar]
  426. Lovas, F. J., 1985, J. Phys. Chem. Ref. Data, 14, 395 [NASA ADS] [CrossRef] [Google Scholar]
  427. Lovas, F. J., & Suenram, R. D., 1987, J. Chem. Phys., 87, 2010 [Google Scholar]
  428. Lovas, F. J., & Groner, P., 2006, J. Mol. Spectr., 236, 173 [Google Scholar]
  429. Lovas, F. J., Suenram, R. D., Ogata, T., & Yamamoto, S., 1992, ApJ, 399, 325 [CrossRef] [Google Scholar]
  430. Lovas, F. J., Kawashima, Y., Grabow, J. U., et al. 1995, ApJ, 455, L201 [Google Scholar]
  431. Lovas, F. J., Plusquellic, D. F., Pate, B. H., et al. 2009, J. Mol. Spectr., 257, 82 [Google Scholar]
  432. Lowenthal, M., Khanna, R., & Moore, M. H., 2002, Spectrochim. Acta Part A, 58, 73 [Google Scholar]
  433. Lu, X., Mills, E. A. C., Ginsburg, A., et al. 2019, ApJS, 244, 35 [Google Scholar]
  434. Luková, K., Kolesniková, L., Koucký, J., et al. 2022, A&A, 665, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  435. Lykke, J. M., Coutens, A., Jørgensen, J. K., et al. 2017, A&A, 597, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  436. Mäder, H., Heise, H. M., & Dreizler, H., 1974, Z. Naturforsch. A, 29, 164 [Google Scholar]
  437. Maeda, A., Medvedev, I. R., De Lucia, F. C., & Herbst, E., 2006, ApJS, 166, 650 [NASA ADS] [CrossRef] [Google Scholar]
  438. Maeda, A., Medvedev, I. R., Winnewisser, M., et al. 2008, ApJS, 176, 543 [NASA ADS] [CrossRef] [Google Scholar]
  439. Maiwald, F., Lewen, F., Ahrens, V., et al. 2000, J. Mol. Spectr., 202, 166 [NASA ADS] [CrossRef] [Google Scholar]
  440. Maki, A. G., & Johnson, D. R., 1973, J. Mol. Spectr., 47, 226 [Google Scholar]
  441. Mallinson, P. D., & de Zafra, R. L., 1978, Mol. Phys., 36, 827 [NASA ADS] [CrossRef] [Google Scholar]
  442. Manigand, S., Coutens, A., Loison, J. C., et al. 2021, A&A, 645, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  443. Marcelino, N., Agúndez, M., Tercero, B., et al. 2020, A&A, 643, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  444. Marcelino, N., Tercero, B., Agúndez, M., & Cernicharo, J., 2021, A&A, 646, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  445. Marenich, A. V., & Boggs, J. E., 2005, Chem. Phys. Lett., 404, 351 [Google Scholar]
  446. Maret, S., Hily-Blant, P., Pety, J., Bardeau, S., & Reynier, E., 2011, A&A, 526, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  447. Margulès, L., Motiyenko, R., Demyk, K., et al. 2009a, A&A, 493, 565 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  448. Margulès, L., Perrin, A., Janečkovà, R., et al. 2009b, Can. J. Phys., 87, 425 [Google Scholar]
  449. Margulès, L., Motiyenko, R. A., Ilyushin, V. V., & Guillemin, J. C., 2015, A&A, 579, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  450. Margulès, L., Belloche, A., Müller, H. S. P., et al. 2016, A&A, 590, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  451. Margulès, L., McGuire, B. A., Senent, M. L., et al. 2017, A&A, 601, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  452. Margulès, L., Tercero, B., Guillemin, J. C., Motiyenko, R. A., & Cernicharo, J., 2018, A&A, 610, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  453. Margulès, L., Ilyushin, V. V., McGuire, B. A., et al. 2020a, J. Mol. Spectr., 371, 111304 [CrossRef] [Google Scholar]
  454. Margulès, L., McGuire, B. A., Motiyenko, R. A., et al. 2020b, A&A, 638, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  455. Margulès, L., Remijan, A., Belloche, A., et al. 2022, A&A, 663, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  456. Marimuthu, A. N., Huis in’t Veld, F., Thorwirth, S., Redlich, B., & Brünken, S. 2021, J. Mol. Spectr., 379, 111477 [Google Scholar]
  457. Maris, A., 2004, Phys. Chem. Chem. Phys., 6, 2611 [CrossRef] [Google Scholar]
  458. Markus, C. R., Thorwirth, S., Asvany, O., & Schlemmer, S., 2019, Phys. Chem. Chem. Phys., 21, 26406 [Google Scholar]
  459. Marstokk, K. M., & Møllendal, H., 1970, J. Mol. Struct., 5, 205 [Google Scholar]
  460. Marstokk, K. M., & Møllendal, H., 1973a, J. Mol. Struct., 15, 137 [NASA ADS] [CrossRef] [Google Scholar]
  461. Marstokk, K. M., & Møllendal, H., 1973b, J. Mol. Struct., 16, 259 [Google Scholar]
  462. Marstokk, K.-M., & Mollendal, H., 1985, Acta Chem. Scand. A, 39, 15 [Google Scholar]
  463. Marstokk, K. M., Møllendal, H., & Samdal, S., 2000, J. Mol. Struct., 524, 69 [Google Scholar]
  464. Martin-Drumel, M. A., Lee, K. L. K., Belloche, A., et al. 2019, A&A, 623, A167 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  465. Martin-Drumel, M.-A., Coutens, A., Loison, J.-C., Jørgensen, J. K., & Pirali, O., 2024, A&A, 687, A233 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  466. Martin-Drumel, M.-A., Spaniol, J.-T., Chitarra, O., Pirali, O., & Müller, H. S. P., 2025, J. Phys. Chem. A, 129, 754 [Google Scholar]
  467. Martinez, O., Lattanzi, V., Thorwirth, S., & McCarthy, M. C., 2013, J. Chem. Phys., 138, 094316 [NASA ADS] [CrossRef] [Google Scholar]
  468. Mbosei, L., Fayt, A., Dréan, P., & Cosléou, J., 2000, J. Mol. Struct., 517, 271 [CrossRef] [Google Scholar]
  469. McCarthy, M. C., Gottlieb, C. A., & Thaddeus, P., 1995a, J. Mol. Spectr., 173, 303 [Google Scholar]
  470. McCarthy, M. C., Gottlieb, C. A., Thaddeus, P., Horn, M., & Botschwina, P., 1995b, J. Chem. Phys., 103, 7820 [NASA ADS] [CrossRef] [Google Scholar]
  471. McCarthy, M. C., Lee, K. L. K., Carroll, P. B., et al. 2020, J. Phys. Chem. A, 124, 5170 [Google Scholar]
  472. McGrath, E. J., Goss, W. M., & De Pree, C. G., 2004, ApJS, 155, 577 [NASA ADS] [CrossRef] [Google Scholar]
  473. McGuire, B. A., Martin-Drumel, M.-A., Thorwirth, S., et al. 2016a, Phys. Chem. Chem. Phys., 18, 22693 [Google Scholar]
  474. McGuire, B. A., Carroll, P. B., Loomis, R. A., et al. 2016b, Science, 352, 1449 [Google Scholar]
  475. McGuire, B. A., Martin-Drumel, M.-A., Lee, K. L. K., et al. 2018, Phys. Chem. Chem. Phys., 20, 13870 [Google Scholar]
  476. McGuire, B. A., Burkhardt, A. M., Loomis, R. A., et al. 2020, ApJ, 900, L10 [Google Scholar]
  477. McKellar, A. R. W., & Appadoo, D. R. T., 2008, J. Mol. Spectr., 250, 106 [Google Scholar]
  478. McNaughton, D., & Evans, C. J., 1999, J. Mol. Spectr., 196, 274 [Google Scholar]
  479. McNaughton, D., & Robertson, E. G., 1994, J. Mol. Spectr., 163, 80 [NASA ADS] [CrossRef] [Google Scholar]
  480. McNaughton, D., Osman, O., & Kroto, H., 1988, J. Mol. Struct., 190, 195 [Google Scholar]
  481. McNaughton, D., Evans, C. J., Lane, S., & Nielsen, C. J., 1999, J. Mol. Spectr., 193, 104 [NASA ADS] [CrossRef] [Google Scholar]
  482. Medvedev, I. R., De Lucia, F. C., & Herbst, E., 2009, ApJS, 181, 433 [Google Scholar]
  483. Mehrotra, S. C., & Maeder, H., 1984, Can. J. Phys., 62, 1280 [Google Scholar]
  484. Melli, A., Melosso, M., Tasinato, N., et al. 2018, ApJ, 855, 123 [NASA ADS] [CrossRef] [Google Scholar]
  485. Melli, A., Melosso, M., Bizzocchi, L., et al. 2022, J. Phys. Chem. A, 126, 6210 [Google Scholar]
  486. Melosso, M., Melli, A., Puzzarini, C., et al. 2018, A&A, 609, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  487. Melosso, M., McGuire, B. A., Tamassia, F., Degli Esposti, C., & Dore, L., 2019, ACS Earth Space Chem., 3, 1189 [Google Scholar]
  488. Melosso, M., Belloche, A., Martin-Drumel, M. A., et al. 2020a, A&A, 641, A160 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  489. Melosso, M., Melli, A., Spada, L., et al. 2020b, J. Phys. Chem. A, 124, 1372 [NASA ADS] [CrossRef] [Google Scholar]
  490. Melosso, M., Bizzocchi, L., Gazzeh, H., et al. 2022, Chem. Commun., 58, 2750 [NASA ADS] [CrossRef] [Google Scholar]
  491. Merke, I., & Dreizler, H., 1987, Z. Naturforsch. A, 42, 1043 [Google Scholar]
  492. Meyer, D. M. A., Vorobyov, E. I., Elbakyan, V. G., et al. 2021, MNRAS, 500, 4448 [Google Scholar]
  493. Mills, E. A. C., Ginsburg, A., Clements, A. R., et al. 2018, ApJ, 869, L14 [NASA ADS] [CrossRef] [Google Scholar]
  494. Mininni, C., Beltrán, M. T., Rivilla, V. M., et al. 2020, A&A, 644, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  495. Mininni, C., Beltrán, M. T., Colzi, L., et al. 2023, A&A, 677, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  496. Missaoui, D., Brahem, S., Najar, F., Yazidi, O., & Senent, M. L., 2024, ACS Earth Space Chem., 8, 1236 [Google Scholar]
  497. Mito, A., Sakai, J., & Katayama, M., 1984a, J. Mol. Spectr., 103, 26 [Google Scholar]
  498. Mito, A., Sakai, J., & Katayama, M., 1984b, J. Mol. Spectr., 105, 410 [Google Scholar]
  499. Mladenović, M., & Lewerenz, M., 2008, Chem. Phys., 343, 129 [Google Scholar]
  500. Moïses, A., Boucher, D., Burie, J., Demaison, J., & Dubrulle, A., 1982, J. Mol. Spectr., 92, 497 [NASA ADS] [CrossRef] [Google Scholar]
  501. Møllendal, H., Samdal, S., Matrane, A., & Guille min, J.-C. 2011, J. Phys. Chem. A, 115, 7978 [Google Scholar]
  502. Møllendal, H., Margulès, L., Belloche, A., et al. 2012, A&A, 538, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  503. Möller, T., Schilke, P., Sánchez-Monge, Á., & Schmiedeke, A. 2025, A&A, 693, A160 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  504. Morino, Y., & Matsumura, C., 1967, Bull. Chem. Soc. Jpn., 40, 1095 [Google Scholar]
  505. Morino, I., Yamada, K. M. T., & Maki, A. G., 2000, J. Mol. Spectr., 200, 145 [Google Scholar]
  506. Moskienko, E. M., & Dyubko, S. F., 1991, Radiophys. Quant. Electron., 34, 181 [NASA ADS] [CrossRef] [Google Scholar]
  507. Motiyenko, R. A., Margulès, L., & Guillemin, J. C. 2012a, A&A, 544, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  508. Motiyenko, R. A., Tercero, B., Cernicharo, J., & Margulès, L. 2012b, A&A, 548, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  509. Motiyenko, R. A., Margulès, L., & Guillemin, J. C., 2013, A&A, 559, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  510. Motiyenko, R. A., Ilyushin, V. V., Drouin, B. J., Yu, S., & Margulès, L., 2014, A&A, 563, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  511. Motiyenko, R. A., Margulès, L., Alekseev, E. A., & Guillemin, J.-C., 2015, J. Phys. Chem. A, 119, 1048 [Google Scholar]
  512. Motiyenko, R. A., Margulès, L., Despois, D., & Guillemin, J.-C., 2018, Phys. Chem. Chem. Phys., 20, 5509 [NASA ADS] [CrossRef] [Google Scholar]
  513. Motiyenko, R. A., Armieieva, I. A., Margulès, L., Alekseev, E. A., & Guillemin, J. C., 2019, A&A, 623, A162 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  514. Motiyenko, R. A., Belloche, A., Garrod, R. T., et al. 2020, A&A, 642, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  515. Motoki, Y., Isobe, F., Ozeki, H., & Kobayashi, K., 2014, A&A, 566, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  516. Muenter, J. S., & Laurie, V. W., 1966, J. Chem. Phys., 45, 855 [Google Scholar]
  517. Müller, F., & Bauder, A., 1996, J. Mol. Spectr., 179, 61 [Google Scholar]
  518. Müller, H. S. P., & Brünken, S., 2005, J. Mol. Spectr., 232, 213 [CrossRef] [Google Scholar]
  519. Müller, H. S. P., & Christen, D., 2004, J. Mol. Spectr., 228, 298 [Google Scholar]
  520. Müller, H. S. P., & Lewen, F., 2017, J. Mol. Spectr., 331, 28 [CrossRef] [Google Scholar]
  521. Müller, H. S. P., Gendriesch, R., Margulès, L., et al. 2000a, Phys. Chem. Chem. Phys., 2, 3401 [CrossRef] [Google Scholar]
  522. Müller, H. S. P., Klaus, T., & Winnewisser, G. 2000b, A&A, 357, L65 [NASA ADS] [Google Scholar]
  523. Müller, H. S. P., Thorwirth, S., Bizzocchi, L., & Winnewisser, G., 2000c, Z. Naturforsch. A, 55, 491 [Google Scholar]
  524. Müller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G., 2001, A&A, 370, L49 [Google Scholar]
  525. Müller, H. S. P., Pracna, P., & Horneman, V. M., 2002, J. Mol. Spectr., 216, 397 [Google Scholar]
  526. Müller, H. S. P., Schlöder, F., Stutzki, J., & Winnewisser, G., 2005, J. Mol. Struct., 742, 215 [Google Scholar]
  527. Müller, H. S. P., Belloche, A., Menten, K. M., Comito, C., & Schilke, P., 2008, J. Mol. Spectr., 251, 319 [Google Scholar]
  528. Müller, H. S. P., Drouin, B. J., & Pearson, J. C., 2009, A&A, 506, 1487 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  529. Müller, H. S. P., Coutens, A., Walters, A., Grabow, J.-U., & Schlemmer, S. 2011, J. Mol. Spectr., 267, 100 [Google Scholar]
  530. Müller, H. S. P., Spezzano, S., Bizzocchi, L., et al. 2013, J. Phys. Chem. A, 117, 13843 [Google Scholar]
  531. Müller, H. S. P., Brown, L. R., Drouin, B. J., et al. 2015, J. Mol. Spectr., 312, 22 [CrossRef] [Google Scholar]
  532. Müller, H. S. P., Drouin, B. J., Pearson, J. C., et al. 2016a, A&A, 586, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  533. Müller, H. S. P., Belloche, A., Xu, L.-H., et al. 2016b, A&A, 587, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  534. Müller, H. S. P., Walters, A., Wehres, N., et al. 2016c, A&A, 595, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  535. Müller, H. S. P., Zingsheim, O., Wehres, N., et al. 2017, J. Phys. Chem. A, 121, 7121 [Google Scholar]
  536. Müller, H. S. P., Maeda, A., Thorwirth, S., et al. 2019, A&A, 621, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  537. Müller, H. S. P., Thorwirth, S., & Lewen, F., 2020, J. Mol. Struct., 1207, 127769 [Google Scholar]
  538. Müller, H. S. P., Guillemin, J.-C., Lewen, F., & Schlemmer, S. 2022, J. Mol. Spectr., 384, 111584 [Google Scholar]
  539. Müller, H. S. P., Garrod, R. T., Belloche, A., et al. 2023, MNRAS, 523, 2887 [CrossRef] [Google Scholar]
  540. Müller, H. S. P., Maeda, A., Lewen, F., et al. 2024a, Mol. Phys., 122, e2262057 [Google Scholar]
  541. Müller, H. S. P., Coutens, A., Jørgensen, J. K., et al. 2024b, A&A, 687, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  542. Müller, H. S. P., Ilyushin, V. V., Belloche, A., Lewen, F., & Schlemmer, S. 2024c, A&A, 688, A201 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  543. Nazari, P., Rocha, W. R. M., Rubinstein, A. E., et al. 2024, A&A, 686, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  544. Nemes, L., Luckhaus, D., Quack, M., & Johns, J. W. C., 2000, J. Mol. Struct., 517, 217 [Google Scholar]
  545. Nguyen, L., Walters, A., Margulès, L., et al. 2013, A&A, 553, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  546. Niedenhoff, M., Winnewisser, G., Yamada, K. M. T., & Belov, S. P., 1995a, J. Mol. Spectr., 169, 224 [Google Scholar]
  547. Niedenhoff, M., Yamada, K. M. T., Belov, S. P., & Winnewisser, G., 1995b, J. Mol. Spectr., 174, 151 [Google Scholar]
  548. Niedenhoff, M., Yamada, K. M. T., & Winnewisser, G., 1996, J. Mol. Spectr., 176, 342 [Google Scholar]
  549. Niedenhoff, M., Yamada, K. M. T., & Winnewisser, G., 1997, J. Mol. Spectr., 183, 176 [Google Scholar]
  550. Niide, Y., & Hayashi, M., 2003, J. Mol. Spectr., 220, 65 [Google Scholar]
  551. Nony, T., Galván-Madrid, R., Brouillet, N., et al. 2024, A&A, 687, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  552. Nummelin, A., Bergman, P., Hjalmarson, Å., et al. 1998, ApJS, 117, 427 [NASA ADS] [CrossRef] [Google Scholar]
  553. Öberg, K. I., Bottinelli, S., Jørgensen, J. K., & van Dishoeck, E. F. 2010, ApJ, 716, 825 [Google Scholar]
  554. Ohno, K., Matsuura, H., McNaughton, D., & Kroto, H. W., 1985, J. Mol. Spectr., 111, 415 [Google Scholar]
  555. Okabayashi, T., & Tanimoto, M., 1993, J. Chem. Phys., 99, 3268 [Google Scholar]
  556. Ordu, M. H., Müller, H. S. P., Walters, A., et al. 2012, A&A, 541, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  557. Ordu, M. H., Zingsheim, O., Belloche, A., et al. 2019, A&A, 629, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  558. Ozeki, H., Hirao, T., Saito, S., & Yamamoto, S., 2004, ApJ, 617, 680 [Google Scholar]
  559. Padovani, M., Walmsley, C. M., Tafalla, M., Galli, D., & Müller, H. S. P., 2009, A&A, 505, 1199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  560. Pascoli, G., & Lavendy, H., 1998, Int. J. Mass Spectrom., 181, 11 [Google Scholar]
  561. Pasinszki, T., & Westwood, N. P. C., 2001, J. Phys. Chem. A, 105, 1244 [Google Scholar]
  562. Pasinszki, T., Krebsz, M., Bazsó, G., & Tarczay, G., 2009, Chem.: Eur. J., 15, 6100 [Google Scholar]
  563. Patel, D., Margolese, D., & Dyke, T. R., 1979, J. Chem. Phys., 70, 2740 [NASA ADS] [CrossRef] [Google Scholar]
  564. Pearson, J. C., & Mueller, H. S. P., 1996, ApJ, 471, 1067 [NASA ADS] [CrossRef] [Google Scholar]
  565. Pearson, E. F., Creswell, R. A., Winnewisser, M., & Winnewisser, G., 1976, Z. Naturforsch. A, 31, 1394 [Google Scholar]
  566. Pearson, J. C., Sastry, K. V. L. N., Herbst, E., & Delucia, F. C., 1994, J. Mol. Spectr., 166, 120 [Google Scholar]
  567. Pearson, J. C., Brauer, C. S., & Drouin, B. J., 2008, J. Mol. Spectr., 251, 394 [Google Scholar]
  568. Pearson, J. C., Drouin, B. J., Maestrini, A., et al. 2011, Rev. Sci. Instrum., 82, 093105 [Google Scholar]
  569. Pearson, J. C., Yu, S., & Drouin, B. J., 2012, J. Mol. Spectr., 280, 119 [Google Scholar]
  570. Penn, R. E., & Curl, R. F., J. 1971, J. Chem. Phys., 55, 651 [NASA ADS] [CrossRef] [Google Scholar]
  571. Persson, M. V., Jørgensen, J. K., Müller, H. S. P., et al. 2018, A&A, 610, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  572. Peter, R., & Dreizler, H., 1965, Z. Naturforsch. A, 20, 301 [NASA ADS] [Google Scholar]
  573. Peterson, K. A., & Woods, R. C., 1988, J. Chem. Phys., 89, 4929 [Google Scholar]
  574. Pickett, H. M., & Scroggin, D. G., 1974, J. Chem. Phys., 61, 3954 [NASA ADS] [CrossRef] [Google Scholar]
  575. Pickett, H. M., & Boyd, T. L., 1978, Chem. Phys. Lett., 58, 446 [NASA ADS] [CrossRef] [Google Scholar]
  576. Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, J. Quant. Spec. Radiat. Transf., 60, 883 [Google Scholar]
  577. Plath, K. L., Axson, J. L., Nelson, G. C., et al. 2009, React. Kinet. Catal. Lett., 96, 209 [Google Scholar]
  578. Pliva, J., Le, L. D., Johns, J. W. C., Lu, Z., & Bernheim, R. A., 1995, J. Mol. Spectr., 173, 423 [Google Scholar]
  579. Plusquellic, D. F., Lovas, F. J., Pate, B. H., et al. 2009, J. Phys. Chem. A, 113, 12911 [Google Scholar]
  580. Powell, F. X., & Lide, David R., J. 1964, J. Chem. Phys., 41, 1413 [Google Scholar]
  581. Pracna, P., Müller, H. S. P., Urban, Š., Horneman, V. M., & Klee, S. 2009, J. Mol. Spectr., 256, 152 [Google Scholar]
  582. Pracna, P., Urban, J., Votava, O., et al. 2011, J. Phys. Chem. A, 115, 1063 [Google Scholar]
  583. Pratap, P., Dickens, J. E., Snell, R. L., et al. 1997, ApJ, 486, 862 [Google Scholar]
  584. Profeta, L. T., Sams, R. L., Johnson, T. J., & Williams, S. D., 2011, J. Phys. Chem. A, 115, 9886 [Google Scholar]
  585. Prudenzano, D., Laas, J., Bizzocchi, L., et al. 2018, A&A, 612, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  586. Puzzarini, C., 2015, J. Phys. Chem. A, 119, 11614 [CrossRef] [Google Scholar]
  587. Puzzarini, C., Dore, L., & Cazzoli, G., 2003, J. Mol. Spectrosc., 217, 19 [Google Scholar]
  588. Puzzarini, C., Biczysko, M., Bloino, J., & Barone, V., 2014, ApJ, 785, 107 [Google Scholar]
  589. Pyckhout, W., Van Alsenoy, C., Geise, H. J., et al. 1986, J. Mol. Struct., 147, 85 [Google Scholar]
  590. Qin, S.-L., Liu, T., Liu, X., et al. 2022, MNRAS, 511, 3463 [CrossRef] [Google Scholar]
  591. Radziszewski, J. G., Downing, J. W., Gudipati, M. S., et al. 1996, J. Am. Chem. Soc., 118, 10275 [Google Scholar]
  592. Randeil, J., Cox, A. P., & Dreizier, H., 1987, Z. Naturforsch. A, 42, 957 [Google Scholar]
  593. Raymonda, J., & Klemperer, W., 1971, J. Chem. Phys., 55, 232 [Google Scholar]
  594. Raymonda, J. W., Muenter, J. S., & Klemperer, W. A., 1970, J. Chem. Phys., 52, 3458 [NASA ADS] [CrossRef] [Google Scholar]
  595. Read, W. G., Cohen, E. A., & Pickett, H. M., 1986, J. Mol. Spectr., 115, 316 [Google Scholar]
  596. Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2019, ApJ, 885, 131 [Google Scholar]
  597. Reinartz, J., Meertz, W., & Dymanus, A., 1972, Chem. Phys. Lett., 16, 576 [Google Scholar]
  598. Remijan, A. J., Snyder, L. E., McGuire, B. A., et al. 2014, ApJ, 783, 77 [Google Scholar]
  599. Remijan, A., Adams, M., & Warmels, R., 2015, ALMA Cycle 3 Technical Handbook Version 1.0, Tech. rep., ALMA [Google Scholar]
  600. Requena-Torres, M. A., Martín-Pintado, J., Rodríguez-Franco, A., et al. 2006, A&A, 455, 971 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  601. Richard, C., Belloche, A., Margulès, L., et al. 2018, J. Mol. Spectr., 345, 51 [Google Scholar]
  602. Riveros, J. M., & Bright Wilson, E., J. 1967, J. Chem. Phys., 46, 4605 [Google Scholar]
  603. Rivilla, V. M., Martín-Pintado, J., Jiménez-Serra, I., et al. 2019, MNRAS, 483, L114 [NASA ADS] [CrossRef] [Google Scholar]
  604. Rivilla, V. M., Colzi, L., Jiménez-Serra, I., et al. 2022a, ApJ, 929, L11 [NASA ADS] [CrossRef] [Google Scholar]
  605. Rivilla, V. M., García De La Concepción, J., Jiménez-Serra, I., et al. 2022b, Front. Astron. Space Sci., 9, 829288 [NASA ADS] [CrossRef] [Google Scholar]
  606. Rivilla, V. M., Jiménez-Serra, I., Martín-Pintado, J., et al. 2022c, Front. Astron. Space Sci., 9, 876870 [NASA ADS] [CrossRef] [Google Scholar]
  607. Rocha, W. R. M., van Dishoeck, E. F., Ressler, M. E., et al. 2024, A&A, 683, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  608. Rodler, M., 1983, PhD thesis, ETH Zürich, Switzerland [Google Scholar]
  609. Rodler, M., 1985, J. Mol. Spectr., 114, 23 [Google Scholar]
  610. Rodler, M., Brown, R. D., Godfrey, P. D., & Tack, L. M., 1984, Chem. Phys. Lett., 110, 447 [Google Scholar]
  611. Rodler, M., Brown, R. D., Godfrey, P. D., & Kleibömer, B., 1986, J. Mol. Spectr., 118, 267 [NASA ADS] [CrossRef] [Google Scholar]
  612. Rodler, M., Jans-Bürli, S., Oldani, M., & Bauder, A., 1987, Chem. Phys. Lett., 142, 10 [Google Scholar]
  613. Rodríguez-Almeida, L. F., Jiménez-Serra, I., Rivilla, V. M., et al. 2021, ApJ, 912, L11 [Google Scholar]
  614. Römheld, M., 1978, PhD thesis, Universität Ulm, Germany [Google Scholar]
  615. Rosado, M. T. S., Duarte, M. L. R. S., & Fausto, R., 1997, J. Mol. Struct., 410, 343 [Google Scholar]
  616. Rubin, R. H., Swenson, Jr., G. W., Benson, R. C., Tigelaar, H. L., & Flygare, W. H. 1971, ApJ, 169, L39 [NASA ADS] [CrossRef] [Google Scholar]
  617. Saito, S., 1972, ApJ, 178, L95 [NASA ADS] [CrossRef] [Google Scholar]
  618. Saito, S., 1976, Chem. Phys. Lett., 42, 399 [Google Scholar]
  619. Saito, S., & Yamamoto, S., 1997, J. Chem. Phys., 107, 1732 [Google Scholar]
  620. Saito, S., Kawaguchi, K., Yamamoto, S., et al. 1987, ApJ, 317, L115 [NASA ADS] [CrossRef] [Google Scholar]
  621. Sakai, N., Saruwatari, O., Sakai, T., Takano, S., & Yamamoto, S., 2010, A&A, 512, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  622. Sakaizumi, T., Mure, H., Ohashi, O., & Yamaguchi, I., 1990, J. Mol. Spectr., 140, 62 [Google Scholar]
  623. Saleck, A. H., Simon, R., & Winnewisser, G., 1994, ApJ, 436, 176 [NASA ADS] [CrossRef] [Google Scholar]
  624. Sánchez-Monge, Á., Schilke, P., Schmiedeke, A., et al. 2017, A&A, 604, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  625. Sanz, M. E., McCarthy, M. C., & Thaddeus, P., 2005, J. Chem. Phys., 122, 194319 [Google Scholar]
  626. Sanz-Novo, M., Belloche, A., Alonso, J. L., et al. 2020a, A&A, 639, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  627. Sanz-Novo, M., León, I., Alonso, J. L., Largo, A., & Barrientos, C. 2020b, A&A, 644, A3 [Google Scholar]
  628. Sanz-Novo, M., Belloche, A., Rivilla, V. M., et al. 2022, A&A, 666, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  629. Sanz-Novo, M., Rivilla, V. M., Jiménez-Serra, I., et al. 2023, ApJ, 954, 3 [NASA ADS] [CrossRef] [Google Scholar]
  630. Sanz-Novo, M., Rivilla, V. M., Jiménez-Serra, I., et al. 2024a, ApJ, 965, 149 [Google Scholar]
  631. Sanz-Novo, M., Rivilla, V. M., Müller, H. S. P., et al. 2024b, ApJ, 965, L26 [Google Scholar]
  632. Sanz-Novo, M., Molpeceres, G., Rivilla, V. M., & Jiménez-Serra, I., 2025, A&A, 698, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  633. Sastry, K. V. L. N., Rao, V. M., & Dass, S. C., 1968, Can. J. Phys., 46, 959 [NASA ADS] [CrossRef] [Google Scholar]
  634. Sastry, K. V. L. N., Helminger, P., Charo, A., Herbst, E., & De Lucia, F. C., 1981, ApJ, 251, L119 [NASA ADS] [CrossRef] [Google Scholar]
  635. Sastry, K. V. L. N., Mukhopadhyay, I., Gupta, P. K., & VanderLinde, J., 1996, J. Mol. Spectr., 176, 38 [Google Scholar]
  636. Saykally, R. J., Szanto, P. G., Anderson, T. G., & Woods, R. C., 1976, ApJ, 204, L143 [NASA ADS] [CrossRef] [Google Scholar]
  637. Schmid-Burgk, J., Muders, D., Müller, H. S. P., & Brupbacher-Gatehouse, B., 2004, A&A, 419, 949 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  638. Schmiedeke, A., Schilke, P., Möller, T., et al. 2016, A&A, 588, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  639. Schwahn, G., Schieder, R., Bester, M., & Winnewisser, G., 1986, J. Mol. Spectrosc., 116, 263 [Google Scholar]
  640. Schäfer, E., & Winnewisser, M., 1983, Ber. Bunsenges. Phys. Chem., 87, 327 [Google Scholar]
  641. Senent, M. L., 2001, Mol. Phys., 99, 1311 [Google Scholar]
  642. Senent, M. L., 2004, J. Phys. Chem. A, 108, 6286 [Google Scholar]
  643. Senent, M. L., Puzzarini, C., Domínguez-Gómez, R., Carvajal, M., & Hochlaf, M., 2014, J. Chem. Phys., 140, 124302 [Google Scholar]
  644. Sharma, A. R., Wu, J., Braams, B. J., et al. 2006, J. Chem. Phys., 125, 224306 [Google Scholar]
  645. Shimanouchi, T., 1972, Tables of Molecular Vibrational Frequencies Consolidated Vol. I [Google Scholar]
  646. Shope, B. M., El-Abd, S. J., Brogan, C. L., et al. 2024, ApJ, 972, 146 [Google Scholar]
  647. Singh, A., Shivani, Misra, A., & Tandon, P. 2014, Res. Astron. Astrophys., 14, 275 [Google Scholar]
  648. Skatrud, D. D., De Lucia, F. C., Blake, G. A., & Sastry, K. V. L. N., 1983, J. Mol. Spectr., 99, 35 [Google Scholar]
  649. Smirnov, I. A., Alekseev, E. A., Piddyachiy, V. I., Ilyushin, V. V., & Motiyenko, R. A., 2013, J. Mol. Spectr., 293, 33 [Google Scholar]
  650. Smirnov, I. A., Alekseev, E. A., Ilyushin, V. V., et al. 2014, J. Mol. Spectr., 295, 44 [NASA ADS] [CrossRef] [Google Scholar]
  651. Smith, Z., Wilson, E., & Duerst, R. W., 1983, Spectrochim. Acta Part A, 39, 1117 [Google Scholar]
  652. Snels, M., Hollenstein, H., & Quack, M., 2006, J. Chem. Phys., 125, 194319 [Google Scholar]
  653. Spahn, H., Müller, H. S. P., Giesen, T. F., et al. 2008, Chem. Phys., 346, 132 [Google Scholar]
  654. Spezzano, S., Tamassia, F., Thorwirth, S., et al. 2012, ApJS, 200, 1 [Google Scholar]
  655. Steber, A. L., Harris, B. J., Neill, J. L., & Pate, B. H., 2012, J. Mol. Spectr., 280, 3 [Google Scholar]
  656. Stiefvater, O. L., 1975, J. Chem. Phys., 62, 233 [Google Scholar]
  657. Stolze, W. H., Sutter, D. H., & Wentrup, C., 1989, Z. Naturforsch. A, 44, 291 [Google Scholar]
  658. Sugie, M., Takeo, H., & Matsumura, C., 1985, J. Mol. Spectr., 111, 83 [Google Scholar]
  659. Sugisaki, R., Tanaka, T., & Hirota, E., 1974, J. Mol. Spectr., 49, 241 [Google Scholar]
  660. Sutter, D. H., & Dreizler, H., 2000, Z. Naturforsch. A, 55, 695 [Google Scholar]
  661. Suzuki, M., & Kozima, K., 1970, J. Mol. Spectr., 33, 407 [Google Scholar]
  662. Swalen, J. D., & Herschbach, D. R., 1957, J. Chem. Phys., 27, 100 [Google Scholar]
  663. Szalanski, L. B., Gerry, M. C. L., Winnewisser, G., Yamada, K., & Winnewisser, M., 1978, Can. J. Phys., 56, 1297 [NASA ADS] [CrossRef] [Google Scholar]
  664. Tack, L. M., & Kukolich, S. G., 1983, J. Chem. Phys., 78, 6512 [Google Scholar]
  665. Takano, M., Sasada, Y., & Satoh, T., 1968, J. Mol. Spectr., 26, 157 [NASA ADS] [CrossRef] [Google Scholar]
  666. Takano, S., Sugie, M., Sugawara, K.-i., et al. 1990, J. Mol. Spectr., 141, 13 [NASA ADS] [CrossRef] [Google Scholar]
  667. Takashi, R., Tanaka, K., & Tanaka, T., 1989, J. Mol. Spectr., 138, 450 [Google Scholar]
  668. Tamassia, F., Bizzocchi, L., Melosso, M., et al. 2022, J. Quant. Spec. Radiat. Transf., 279, 108044 [Google Scholar]
  669. Tanaka, K., Tanaka, T., & Suzuki, I., 1985, J. Chem. Phys., 82, 2835 [Google Scholar]
  670. Tanimoto, M., & Macdonald, J. N., 1979, J. Mol. Spectr., 78, 106 [Google Scholar]
  671. Tanimoto, M., Almond, V., Charles, S. W., Macdonald, J. N., & Owen, N. L., 1979, J. Mol. Spectr., 78, 95 [Google Scholar]
  672. Taquet, V., Wirström, E. S., & Charnley, S. B., 2016, ApJ, 821, 46 [Google Scholar]
  673. Tayyari, S., & Milani-Nejad, F., 1998, Spectrochim. Acta Part A, 54, 255 [Google Scholar]
  674. Tennis, J. D., Xue, C., Talbi, D., et al. 2023, MNRAS, 525, 2154 [Google Scholar]
  675. Tercero, B., Kleiner, I., Cernicharo, J., et al. 2013, ApJ, 770, L13 [NASA ADS] [CrossRef] [Google Scholar]
  676. Tercero, B., Cuadrado, S., López, A., et al. 2018, A&A, 620, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  677. Tercero, B., Marcelino, N., Roueff, E., et al. 2024, A&A, 682, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  678. Thackston, R., & Fortenberry, R. C., 2018, Icarus, 299, 187 [Google Scholar]
  679. Thomson, R., & Dalby, F. W., 1968, Can. J. Phys., 46, 2815 [Google Scholar]
  680. Thomas, L., Sherrard, E., & Sheridan, J., 1955, Trans. Faraday Soc., 51, 619 [Google Scholar]
  681. Thorwirth, S., Müller, H. S. P., Lewen, F., Gendriesch, R., & Winnewisser, G. 2000a, A&A, 363, L37 [NASA ADS] [Google Scholar]
  682. Thorwirth, S., Müller, H. S. P., & Winnewisser, G., 2000b, J. Mol. Spectr., 204, 133 [NASA ADS] [CrossRef] [Google Scholar]
  683. Thorwirth, S., Mueller, H. S. P., & Winnewisser, G., 2001, Phys. Chem. Chem. Phys., 3, 1236 [CrossRef] [Google Scholar]
  684. Thorwirth, S., Müller, H. S. P., Lewen, F., et al. 2003, ApJ, 585, L163 [NASA ADS] [CrossRef] [Google Scholar]
  685. Thorwirth, S., McCarthy, M. C., Dudek, J. B., & Thaddeus, P., 2004, J. Mol. Spectr., 225, 93 [Google Scholar]
  686. Tiemann, E., 1974, J. Phys. Chem. Ref. Data, 3, 259 [CrossRef] [Google Scholar]
  687. Tiemann, E., 1982, J. Mol. Spectr., 91, 60 [Google Scholar]
  688. Tinti, F., Bizzocchi, L., Degli Esposti, C., & Dore, L., 2007, ApJ, 669, L113 [NASA ADS] [CrossRef] [Google Scholar]
  689. Tsang, W. & Hampson, R. F., 1986, J. Phys. Chem. Ref. Data, 15, 1087 [NASA ADS] [CrossRef] [Google Scholar]
  690. Tsunekawa, S., Taniguchi, I., Tambo, A., et al. 1989, J. Mol. Spectr., 134, 63 [Google Scholar]
  691. Tudorie, M., Ilyushin, V., Vander Auwera, J., et al. 2012, J. Chem. Phys., 137, 064304 [CrossRef] [Google Scholar]
  692. Turner, B. E., 1989, ApJS, 70, 539 [NASA ADS] [CrossRef] [Google Scholar]
  693. Tyree, D. J., Chapman, T. W., Medvedev, I. R., & Kisiel, Z., 2022, J. Mol. Spectr., 390, 111706 [Google Scholar]
  694. Utkin, Y. G., Han, J.-X., Sun, F., et al. 2003, J. Chem. Phys., 118, 10470 [Google Scholar]
  695. van der Tak, F. F. S., Müller, H. S. P., Harding, M. E., & Gauss, J., 2009, A&A, 507, 347 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  696. van Dishoeck, E. F., & Blake, G. A., 1998, ARA&A, 36, 317 [Google Scholar]
  697. van Duk, F. A., & Dymanus, A., 1974, Chem. Phys., 6, 474 [Google Scholar]
  698. Vanek, M. D., Jennings, D. A., Wells, J. S., & Maki, A. G., 1989, J. Mol. Spectr., 138, 79 [Google Scholar]
  699. Vastel, C., Ceccarelli, C., Lefloch, B., & Bachiller, R., 2014, ApJ, 795, L2 [NASA ADS] [CrossRef] [Google Scholar]
  700. Vastel, C., Kawaguchi, K., Quénard, D., et al. 2018, MNRAS, 474, L76 [NASA ADS] [CrossRef] [Google Scholar]
  701. Vávra, K., Kolesniková, L., Belloche, A., et al. 2022, A&A, 666, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  702. Vigouroux, C., Fayt, A., Guarnieri, A., et al. 2000, J. Mol. Spectr., 202, 1 [NASA ADS] [CrossRef] [Google Scholar]
  703. Vormann, K., Andresen, U., Heineking, N., & Dreizler, H., 1988, Z. Naturforsch. A, 43, 283 [Google Scholar]
  704. Vorob’eva, E. M., & Dyubko, S. F., 1994, Radiophys. Quantum Electron., 37, 155 [Google Scholar]
  705. Vorobyov, E. I., & Basu, S., 2015, ApJ, 805, 115 [NASA ADS] [CrossRef] [Google Scholar]
  706. Vrtilek, J. M., Gottlieb, C. A., & Thaddeus, P., 1987, ApJ, 314, 716 [CrossRef] [Google Scholar]
  707. Walmsley, C. M., Cesaroni, R., Olmi, L., Churchwell, E., & Hofner, P., 1995, Ap&SS, 224, 173 [Google Scholar]
  708. Wang, J. H. S., Oates, D. E., Ben-Reuven, A., & Kukolich, S. G., 1973, J. Chem. Phys., 59, 5268 [Google Scholar]
  709. Wang, K., Shan, X., & Chen, X., 2009, J. Mol. Struct.: THEOCHEM, 909, 91 [Google Scholar]
  710. Ware, J. M., & Roberts, J. A., 1984, J. Chem. Phys., 81, 1215 [Google Scholar]
  711. Wehres, N., Hermanns, M., Wilkins, O. H., et al. 2018, A&A, 615, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  712. Widicus, S. L., Drouin, B. J., Dyl, K. A., & Blake, G. A., 2003, J. Mol. Spectr., 217, 278 [NASA ADS] [CrossRef] [Google Scholar]
  713. Widicus Weaver, S. L., Butler, R. A. H., Drouin, B. J., et al. 2005, ApJS, 158, 188 [Google Scholar]
  714. Wierzejewska, M., & Moc, J., 2003, J. Phys. Chem. A, 107, 11209 [NASA ADS] [CrossRef] [Google Scholar]
  715. Willis, E. R., Garrod, R. T., Belloche, A., et al. 2020, A&A, 636, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  716. Winnewisser, G., & Cook, R. L., 1968, J. Mol. Spectr., 28, 266 [NASA ADS] [CrossRef] [Google Scholar]
  717. Winnewisser, M., & Vogt, J., 1978, Z. Naturforsch. A, 33, 1323 [Google Scholar]
  718. Winnewisser, M., & Winnewisser, B. P., 1971, Z. Naturforsch. A, 26, 128 [Google Scholar]
  719. Winnewisser, M., Winnewisser, G., Honda, T., & Hirota, E., 1975, Z. Naturforsch. A, 30, 1001 [NASA ADS] [Google Scholar]
  720. Winnewisser, M., Pearson, E. F., Galica, J., & Winnewisser, B. P., 1982, J. Mol. Spectr., 91, 255 [Google Scholar]
  721. Winnewisser, M., Winnewisser, B. P., & Wentrup, C., 1984, J. Mol. Spectr., 105, 193 [CrossRef] [Google Scholar]
  722. Winnewisser, M., Winnewisser, B. P., & Winnewisser, G., 1985, in NATO Advanced Study Institute (ASI) Series C, 157, eds. G. H. F. Diercksen, W. F. Huebner, & P. W. Langhoff, 375 [Google Scholar]
  723. Winnewisser, M., Winnewisser, B. P., Stein, M., et al. 2002, J. Mol. Spectr., 216, 259 [NASA ADS] [CrossRef] [Google Scholar]
  724. Winnewisser, M., Medvedev, I. R., De Lucia, F. C., et al. 2005, ApJS, 159, 189 [NASA ADS] [CrossRef] [Google Scholar]
  725. Winton, R. S., & Gordy, W., 1970, Phys. Lett. A, 32, 219 [Google Scholar]
  726. Wlodarczak, G., Boucher, D., Bocquet, R., & Demaison, J., 1986, J. Mol. Spectr., 116, 251 [Google Scholar]
  727. Wlodarczak, G., Bocquet, R., Bauer, A., & Demaison, J., 1988, J. Mol. Spectr., 129, 371 [NASA ADS] [CrossRef] [Google Scholar]
  728. Wlodarczak, G., Burie, J., Demaison, J., Vormann, K., & Császár, A. G., 1989, J. Mol. Spectr., 134, 297 [Google Scholar]
  729. Wlodarczak, G., Demaison, J., Heineking, N., & Csaszar, A. G., 1994, J. Mol. Spectr., 167, 239 [Google Scholar]
  730. Wohlfart, K., Schnell, M., Grabow, J.-U., & Küpper, J., 2008, J. Mol. Spectr., 247, 119 [NASA ADS] [CrossRef] [Google Scholar]
  731. Wollrab, J. E., & Laurie, V. W., 1968, J. Chem. Phys., 48, 5058 [Google Scholar]
  732. Woon, D. E., 1995, Chem. Phys. Lett., 244, 45 [NASA ADS] [CrossRef] [Google Scholar]
  733. Woon, D. E., & Herbst, E., 2009, ApJS, 185, 273 [Google Scholar]
  734. Xu, L.-H., & Lovas, F. J. 1997, J. Phys. Chem. Ref. Data, 26, 17 [Google Scholar]
  735. Xu, L.-H., Fisher, J., Lees, R. M., et al. 2008, J. Mol. Spectr., 251, 305 [NASA ADS] [CrossRef] [Google Scholar]
  736. Xu, L.-H., Lees, R. M., Crabbe, G. T., et al. 2012, J. Chem. Phys., 137, 104313 [NASA ADS] [CrossRef] [Google Scholar]
  737. Yamada, K., 1977, J. Mol. Spectr., 68, 423 [Google Scholar]
  738. Yamada, K., & Winnewisser, M., 1975, Z. Naturforsch. A, 30, 672 [Google Scholar]
  739. Yamada, K., & Winnewisser, M., 1977, J. Mol. Spectr., 68, 307 [Google Scholar]
  740. Yamada, K. M. T., & Creswell, R. A., 1986, J. Mol. Spectr., 116, 384 [NASA ADS] [CrossRef] [Google Scholar]
  741. Yamada, K., Winnewisser, M., Winnewisser, G., Szalanski, L. B., & Gerry, M. C. L., 1979, J. Mol. Spectr., 78, 189 [Google Scholar]
  742. Yamada, K., Winnewisser, M., Winnewisser, G., Szalanski, L. B., & Gerry, M. C. L., 1980, J. Mol. Spectr., 79, 295 [Google Scholar]
  743. Yamada, K. M. T., Moravec, A., & Winnewisser, G., 1995, Z. Naturforsch. A, 50, 1179 [NASA ADS] [CrossRef] [Google Scholar]
  744. Yamada, K. M. T., Degli Esposti, C., Botschwina, P., et al. 2004, A&A, 425, 767 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  745. Yamaguchi, Y., Richards, Claude A., J., & Schaefer, Henry F., I. 1994, J. Chem. Phys., 101, 8945 [Google Scholar]
  746. Yang, Y.-L., Green, J. D., Pontoppidan, K. M., et al. 2022, ApJ, 941, L13 [NASA ADS] [CrossRef] [Google Scholar]
  747. Zakharenko, O., Bossa, J. B., Lewen, F., Schlemmer, S., & Müller, H. S. P., 2017, J. Mol. Spectr., 333, 23 [Google Scholar]
  748. Zakharenko, O., Ilyushin, V. V., Lewen, F., et al. 2019b, A&A, 629, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  749. Zakharenko, O., Lewen, F., Ilyushin, V. V., et al. 2019a, A&A, 627, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  750. Zaleski, D. P., Seifert, N. A., Steber, A. L., et al. 2013, ApJ, 765, L10 [Google Scholar]
  751. Zeisberger, E., & Botskor, I., 1983, J. Mol. Struct., 97, 323 [Google Scholar]
  752. Zeng, S., Jiménez-Serra, I., Rivilla, V. M., et al. 2018, MNRAS, 478, 2962 [Google Scholar]
  753. Zeng, S., Zhang, Q., Jiménez-Serra, I., et al. 2020, MNRAS, 497, 4896 [NASA ADS] [CrossRef] [Google Scholar]
  754. Zeng, S., Rivilla, V. M., Jiménez-Serra, I., et al. 2023, MNRAS, 523, 1448 [NASA ADS] [CrossRef] [Google Scholar]
  755. Zingsheim, O., Maßen, J., Müller, H. S. P., et al. 2022, A&A, 662, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  756. Zingsheim, O., Müller, H. S. P., Lewen, F., Jørgensen, J. K., & Schlemmer, S., 2017, J. Mol. Spectr., 342, 125 [NASA ADS] [CrossRef] [Google Scholar]

6

A comparison between G+0.693 and Sgr B(N2) on the basis of our previous survey EMoCA, which had a lower angular resolution (1.6′′) and could not disentangle the four sources embedded in Sgr B2(N2), was presented by Jiménez-Serra et al. (2025).

7

On top of its extremely low column density, the weak component of C2H5CN has a puzzling velocity offset of −159 km s−1 in Table L.57 of Möller et al. (2025). We do not know what this component represents.

8

As explained in Sect. 1, our previous single-dish survey of Sgr B2(N) at 3, 2, and 1.3 mm already showed that spectral confusion and optical depth of the continuum emission are severe issues at 1.3 mm for this source (Belloche et al. 2013).

All Tables

Table 1

Beam sizes and noise levels.

Table 2

Properties of the positions analyzed within Sgr B2(N2).

Table 3

List of molecules detected toward N2b, AN02, AN03, AN06, and AN06c2 with the ReMoCA survey.

Table 4

Correlations between the chemical composition of Sgr B2(N2) sources and sources from the literature.

Table 5

Range of final (i.e. stage-2) gas densities from the model grid used in the comparison with each source, in the density-restricted setup.

Table A.1

Column densities of G31.41, IRAS 16293B, G+0.693, and TMC-1 collected from the literature.

Table B.1

List of unidentified lines brighter than 10 K in the ReMoCA spectrum of N2b.

Table C.1

Spectroscopic references of molecules analyzed in this work.

Table E.1

Parameters of our best-fit LTE model of N2b.

Table E.2

Parameters of our best-fit LTE model of AN02.

Table E.3

Parameters of our best-fit LTE model of AN03

Table E.4

Parameters of our best-fit LTE model of the second velocity component toward AN06 (AN06c2).

Table E.5

Parameters of our best-fit LTE model of AN06.

Table E.6

Column density upper limits toward N2b.

Table E.7

Column density upper limits toward N2b.

Table E.8

Column density upper limits toward AN03.

Table E.9

Column density upper limits toward the second velocity component toward AN06 (AN06c2).

Table E.10

Column density upper limits toward AN06.

All Figures

Thumbnail: Fig. 1 Refer to the following caption and surrounding text. Fig. 1

ALMA continuum-subtracted spectra toward the hot core positions AN06, AN03, AN02, and N2b (from left to right) at the frequencies of the hydrogen recombination lines H40α (top and middle rows, for two setups with different angular resolutions) and H41α (bottom row). In each panel, the vertical dashed line marks the systemic velocity of the source adopted for the LTE modeling of the molecular emission. The velocity axis refers to the rest frequency of the recombination line. The horizontal dashed line indicates the 3σ noise level. The blue spectrum represents the LTE model that includes the contribution of all molecules identified so far. The velocity range highlighted in dark gray is specific to each recombination line and represents the range of channels selected to compute the integrated intensity map of each recombination line shown in Fig. 2. The setup and spectral window are indicated in the left panel of each row along with the corresponding HPBW.

In the text
Thumbnail: Fig. 2 Refer to the following caption and surrounding text. Fig. 2

ALMA integrated intensity maps of the hydrogen recombination lines H41α (panel a) and H40α (panels b and c, for two setups with different angular resolutions). The intensity was integrated over the velocity range highlighted in dark gray in Fig. 1 in order to avoid contamination by molecular lines. The violet plus symbols mark the hot core positions AN06, AN03, AN02, and N2b. The green cross indicates the VLA position of the HCH II region K7 at 7 mm from De Pree et al. (2015). The blue triangles and violet square indicate the VLA positions of the water and Class II methanol masers reported by McGrath et al. (2004) and Lu et al. (2019), respectively. The setup and spectral window numbers are given in the bottom right corner of each panel along with the associated beam size (HPBW). The values of the noise level, σ, are 18, 21, and 15 mJy beam −1 km s−1, respectively. The contours start at 3σ (brown contour) and then increase by a factor of two at each step (black contours). The dotted blue contour, when present, shows the −3σ level.

In the text
Thumbnail: Fig. 3 Refer to the following caption and surrounding text. Fig. 3

ALMA continuum emission maps at 92.3, 98.9, and 99.6 GHz. The setup and spectral window numbers are given in the bottom right corner of each panel along with the associated beam size (HPBW). The symbols are the same as in Fig. 2. The values of the noise level, σ, are 0.73, 0.97, and 0.42 mJy beam −1, respectively. The contours start at 5σ (brown contour) and then increase by a factor of two at each step (black contours). The additional violet contour in panel c is at 60σ. It was added to emphasize the emission peak on AN02. Dotted blue contours, when present, indicate the −5σ and −10σ levels.

In the text
Thumbnail: Fig. 4 Refer to the following caption and surrounding text. Fig. 4

ALMA continuum-subtracted spectra toward the hot core positions AN06, AN03, AN02, and N2b (from left to right) at the frequencies of two contamination-free rotational transitions of ethanol indicated on the right along with the energy of the upper level in temperature unit. In each panel, the vertical dashed line marks the systemic velocity of the source adopted for the LTE modeling of the molecular emission. The velocity axis refers to the rest frequency of the ethanol transition. The horizontal dashed line indicates the 3σ noise level. The blue spectrum represents the LTE model that includes the contribution of all molecules identified so far. The red spectrum represents the LTE model of ethanol only. In the case of AN06, the synthetic ethanol spectra of the first and second velocity components are displayed in orange and red, respectively. The velocity ranges highlighted in the three darkest shades of gray are specific to each ethanol transition (but common to all sources, albeit with slight differences due to the finite spectral sampling) and represent the ranges of channels selected to compute the integrated intensity maps of each transition shown in Fig. 5. The setup and spectral window are indicated in the left panel of each row along with the corresponding HPBW.

In the text
Thumbnail: Fig. 5 Refer to the following caption and surrounding text. Fig. 5

ALMA integrated intensity maps of two contamination-free rotational transitions of ethanol with rest frequency and upper level energy indicated on the right of each row. The intensities were integrated over the velocity ranges highlighted in the three darkest shades of gray in Fig. 4. The integration range is indicated in the top left corner of each panel. The setup and spectral window numbers are given in the bottom right corner along with the associated beam size (HPBW). The symbols are the same as in Fig. 2. The noise level, σ, is indicated below the velocity range. The black contours start at 3σ and then increase by a factor of two at each step, except for the last contour of panel e which is at 320σ. The dotted blue contour, when present, indicates the −3σ level.

In the text
Thumbnail: Fig. 6 Refer to the following caption and surrounding text. Fig. 6

Maps of centroid velocity (left column) and line width (right column) of two contamination-free rotational transitions of ethanol with rest frequency and upper level energy indicated in the top left corner of each row. The kinematic information is plotted only for pixels with a signal-to-noise ratio in peak intensity higher than 10. The setup and window numbers are given in the bottom right corner of each panel along with the associated beam size (HPBW). The symbols are the same as in Fig. 2. These two transitions of ethanol have opacities on the order of 2 in N2b, which means that the intrinsic line widths around N2b are somewhat smaller than the plotted values.

In the text
Thumbnail: Fig. 7 Refer to the following caption and surrounding text. Fig. 7

Column densities derived toward N2b, AN02, AN03, and AN06 with our LTE modeling, normalized to the column density of methanol (see Tables E.1E.5). The column density of methanol is indicated in the top left corner of each panel. The chemical compositions of the two velocity components detected toward AN06 are displayed separately (AN06 in the top panel and the second velocity component in the panel labeled AN06c2). The panel of AN02 reports the sum of the column densities of its two velocity components. Hatched bars show tentative detections while empty bars with downward arrows indicate upper limits. Molecules with status “c” in Tables E.1E.5 are here represented as upper limits. The dashed lines indicate levels of 1% and 0.01% with respect to methanol.

In the text
Thumbnail: Fig. 8 Refer to the following caption and surrounding text. Fig. 8

Correlation plots of column densities normalized to methanol for various pairs of positions. The x- and y-axes of each panel correspond to the positions written in the bottom right and top left corners, respectively. The color coding of the molecules is the same as in Figs. 7 and F.1 and is indicated on the right. Latine and Greek letters, reported at the same abscissa as the corresponding molecule, were added to facilitate the identification of the data points. Filled data points represent firm and tentative detections while empty circles with arrows indicate upper limits toward at least one position. The dashed and dotted lines indicate deviations by a factor 10 and 2, respectively.

In the text
Thumbnail: Fig. 9 Refer to the following caption and surrounding text. Fig. 9

Abundance of protonated molecules (red) relative to their neutral form (black). Each panel corresponds to one position as labeled in the bottom right corner. Empty squares and downward arrows represent tentative detections and upper limits, respectively. Crosses of the same color represent predictions of the chemical model described in Sect. 4. The molecules are listed along the x axis below and above the bottom and top panels, respectively. The list of protonated molecules is indicated at the bottom right.

In the text
Thumbnail: Fig. 10 Refer to the following caption and surrounding text. Fig. 10

Same as Fig. 9, but for radicals (red) with respect to the hydrogenated form (black).

In the text
Thumbnail: Fig. 11 Refer to the following caption and surrounding text. Fig. 11

Same as Fig. 9, but for molecules with a double (red) or triple (green) bond with respect to the saturated form with a single bond (black).

In the text
Thumbnail: Fig. 12 Refer to the following caption and surrounding text. Fig. 12

Same as Fig. 9, but for molecules with O-bearing functional groups -CHO and -COOH with respect to molecules with a -CN functional group.

In the text
Thumbnail: Fig. 13 Refer to the following caption and surrounding text. Fig. 13

Same as Fig. 9, but for the reduced form (-CH3. functional group) with respect to the aldehyde form (-CHO functional group).

In the text
Thumbnail: Fig. 14 Refer to the following caption and surrounding text. Fig. 14

Same as Fig. 9, but for series of molecules with an additional CH2 group in their backbone at each further step in complexity.

In the text
Thumbnail: Fig. 15 Refer to the following caption and surrounding text. Fig. 15

Same as Fig. 9, but for groups of structural isomers. The elemental composition of the groups is labeled along the x axis. The list of molecules belonging to each group is provided at the bottom of the figure with their respective colors as used in the plots.

In the text
Thumbnail: Fig. 16 Refer to the following caption and surrounding text. Fig. 16

Same as Fig. F.2 but comparing sources from the literature to N2b.

In the text
Thumbnail: Fig. 17 Refer to the following caption and surrounding text. Fig. 17

Same as Fig. F.2 but comparing sources from the literature to AN02.

In the text
Thumbnail: Fig. 18 Refer to the following caption and surrounding text. Fig. 18

Same as Fig. F.2 but comparing sources from the literature to AN03.

In the text
Thumbnail: Fig. 19 Refer to the following caption and surrounding text. Fig. 19

Same as Fig. F.2 but comparing sources from the literature to AN06c2.

In the text
Thumbnail: Fig. 20 Refer to the following caption and surrounding text. Fig. 20

Same as Fig. F.2 but comparing sources from the literature to AN06.

In the text
Thumbnail: Fig. 21 Refer to the following caption and surrounding text. Fig. 21

Overall matching parameter obtained for the comparison of the grid of chemical model results with observational abundances of all species across all sources, assuming a uniform CRIR. A lower matching parameter indicates a closer match between models and observations. Panel a: unrestricted density. Panel b: density restricted to the ranges indicated in Table 5.

In the text
Thumbnail: Fig. 22 Refer to the following caption and surrounding text. Fig. 22

Comparison of best-fit model abundances to observational results, using the density-restricted setup. Bars indicate the number of orders of magnitude by which each abundance ratio in the best-matching model exceeds (or otherwise) the observational value. Data correspond to the best-matching model for each individual source, within the subset of grid models with ζ = 100 ζ0, which is found to produce the best overall match across all sources. The shaded area represents values where the models and observations vary by 1 oom or less. Unfilled blue bars with an arrow indicate that the comparison is based on an observational upper limit. Black bars indicate observed species that were not included in the matching parameter analysis.

In the text
Thumbnail: Fig. D.1 Refer to the following caption and surrounding text. Fig. D.1

ALMA continuum-subtracted spectra of the hot core positions N2b, AN02, AN03, AN06 and the UCH II region K4 (black). LTE synthetic spectra are overlaid in red. See Appendix D for a detailed description of the figure.

In the text
Thumbnail: Fig. F.1 Refer to the following caption and surrounding text. Fig. F.1

Same as for Fig. 7 but normalized to N2b. The dashed lines indicate values of 0.1, 1, and 10.

In the text
Thumbnail: Fig. F.2 Refer to the following caption and surrounding text. Fig. F.2

Same as Fig. 8 but with colors coding for the class of molecules (black for O-bearing, blue for O+N-bearing, red for N-bearing, green for pure hydrocarbon, and yellow for S-bearing).

In the text
Thumbnail: Fig. F.3 Refer to the following caption and surrounding text. Fig. F.3

Matching parameters for all tested models, using the density-restricted setup. Each curve shows the matching parameter for a single source as a function of all the models in the grid with which it is compared, using the same fixed CRIR. The results are arranged in bands, one for each CRIR. The bands are arranged left to right in order of best to worst overall matching parameter. Within each band, the models contributing to each curve are arranged in order of increasing (i.e. worse) matching parameter, as determined independently for each source comparison.

In the text
Thumbnail: Fig. F.4 Refer to the following caption and surrounding text. Fig. F.4

Correlation plots of rotational temperatures for various pairs of positions. The x- and y-axes of each panel correspond to the positions written in the bottom right and top left corners, respectively. The color coding of the molecules is the same as in Fig. 8, plus black for methanol, and is indicated on the right. The plain line indicates the one-to-one relation.

In the text
Thumbnail: Fig. F.5 Refer to the following caption and surrounding text. Fig. F.5

Comparison of the column densities obtained from the ReMoCA survey toward AN02, AN03, AN06, and AN06c2 to those derived from an ALMA survey of Sgr B2(N) at 1.2 mm by Möller et al. (2025). The ReMoCA column densities correspond to peak column densities while the 1.2 mm column densities were derived from spectra averaged over the shaded blue polygons shown in Fig. 2 of Möller et al. (2025). Panels c and d show the same 1.2 mm column densities except for C2H5CN and OCS which were modeled with two velocity components similar to the components that we extracted from the ReMoCA survey. The dashed, dotted, and plain lines and the color coding and labels of the molecules are the same as in Fig. 8.

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.