Evolution of Chemistry in the envelope of HOt corinoS (ECHOS) I. Extremely young sulphur chemistry in the isolated Class 0 object B 335

Context. Within the project Evolution of Chemistry in the envelope of HOt corinoS (ECHOS), we present a study of sulphur chemistry in the envelope of the Class 0 source B335 through observations in the spectral range λ = 7, 3, and 2 mm. Aims. Our goal is to characterise the sulphur chemistry in this isolated protostellar source and compare it with other Class 0 objects to determine the environmental and evolutionary effects on the sulphur chemistry in these young sources. Methods. We have modelled observations and computed column densities assuming local thermodynamic equilibrium and large velocity gradient approximation. We have also used the code Nautilus to study the time evolution of sulphur species, as well as of several sulphur molecular ratios. Results. We have detected 20 sulphur species in B335 with a total gas-phase S abundance similar to that found in the envelopes of other Class 0 objects, but with significant differences in the abundances between sulphur carbon chains and sulphur molecules containing oxygen and nitrogen. Our results highlight the nature of B335 as a source especially rich in sulphur carbon chains unlike other Class 0 sources. The low presence or absence of some molecules, such as SO and SO + , suggests a chemistry not particularly influenced by shocks. We, however, detect a large presence of HCS + that, together with the low rotational temperatures obtained for all the S species ( < 15K), reveals the moderate or low density of the envelope of B335. Model results also show the large influence of the cosmic ray ionisation rate and density variations on the abundances of some S species (e.g. SO, SO 2 , CCS, and CCCS) with differences of up to ∼ 4 orders of magnitude. We also find that observations are better reproduced by models with a sulphur depletion factor of 10 with respect to the sulphur cosmic elemental abundance. Conclusions. The comparison between our model and observational results for B335 reveals an age of 10 4 < t < 10 5 yr, which highlights the particularly early evolutionary stage of this source. B335 presents a different chemistry compared to other young protostars that have formed in dense molecular clouds, which could be the result of accretion of surrounding material from the diffuse cloud onto the protostellar envelope of B335. In addition, the theoretical analysis and comparison with observations of the SO 2 /C 2 S, SO/CS, and HCS + /CS ratios within a sample of prestellar cores and Class 0 objects show that they could be used as good chemical evolutionary indicators of the prestellar to protostellar transition.


Introduction
Hot corinos are regions that arise during the formation of solar-type protostars.The progressive collapse of a prestellar core leads to the heating of the infalling envelope of dust and gas to temperatures of hundreds of Kelvin, while the densities increase up to ∼10 8 -10 9 cm −3 in the inner ∼100 au around the central protostar.These are also the scales at which protoplanetary disks are expected to arise owing to the conservation of angular momentum.As the temperatures increase, the water-rich ice mantles sublimate, injecting molecules into the gas phase.Hot corinos are therefore chemically rich regions in Solar-like young protostars (Ceccarelli 2004;Ceccarelli et al. 2007b).In particular, these regions are relatively rich in interstellar com-plex organic molecules (Herbst & van Dishoeck 2009).Studying, therefore, the inventory and abundances of molecules in the earliest stages of star formation, that is in hot corinos, is essential to probe the initial conditions of planet and star formation processes.In spite of the detection of various hot corinos in the last decades (e.g.Cazaux et al. 2003;Okoda et al. 2023), just a few of them (such as IRAS 16293-2422, NGC 1333, and IRAS 4A van Dishoeck et al. 1995;Ceccarelli et al. 2007a;Jaber et al. 2014;Caselli & Ceccarelli 2012;Jørgensen et al. 2016) have been investigated in detail until now.Many questions, therefore, still remain open about their nature and molecular composition: for instance, chemical differences between the earliest stages of star formation (Class 0 versus Class I objects), the role of the A&A proofs: manuscript no.paper_B335 chemical characteristics inherited from the parent cloud, as well as the environmental effect on the chemical composition of hot corinos.
One way to understand the dynamics and evolution of the earliest stages of star formation is through molecular depletion.At temperatures of roughly 10 K and densities above 10 4 cm −3 , several molecules, such as CO and CS, condense out onto dust grain surfaces (e.g.Caselli et al. 1999;Crapsi et al. 2005), and depletion increases with time (e.g.Bergin & Langer 1997;Aikawa et al. 2003).Therefore, depletion can be used as a clock, since evolved cores should be more depleted of certain species than younger cores.Regarding depletion, sulphur plays a fundamental role since, although sulphur is one of the most abundant species in the Universe and plays a crucial role in biological systems on Earth (e.g.Leustek 2002;Francioso et al. 2020), S-bearing molecules are not as abundant as expected in the interstellar medium (ISM).In particular, unlike the diffuse ISM and photon-dominated regions where the observed gaseous sulphur accounts (e.g.Goicoechea et al. 2006;Howk et al. 2006) for its total solar abundance (S/H∼1.5×10−5 ; Asplund et al. 2009), sulphur is strongly depleted in more dense molecular gas, since the sum of the detectable S-bearing molecules only accounts for a very small fraction of the elemental S abundance (Tieftrunk et al. 1994).In fact, to reproduce observations in hot corinos and hot cores, one needs to assume a significant sulphur depletion of at least one order of magnitude lower than the solar elemental sulphur abundance (e.g.Wakelam et al. 2004;Esplugues et al. 2014;Crockett et al. 2014;Vastel et al. 2018;Bulut et al. 2021;Navarro-Almaida et al. 2021;Esplugues et al. 2022;Hily-Blant et al. 2022;Fuente et al. 2023).Most of the sulphur seems to be locked on the icy grain mantles (e.g.Millar & Herbst 1990;Ruffle et al. 1999;Vidal et al. 2017;Laas & Caselli 2019) in the form of organo-sulphur molecules, also present in the comet 67P (Calmonte et al. 2016) and in meteoritic material (e.g.Naraoka et al. 2022;Ruf et al. 2021).However, OCS and SO 2 are the only S-bearing molecules unambiguously detected in ice mantles to date (Geballe et al. 1985;Palumbo et al. 1995;McClure et al. 2023), while H 2 S has not been detected (Yang et al. 2022;McClure et al. 2023).Theoretical models also suggest that sulphur could be locked into pure allotropes of sulphur, especially S 8 (Shingledecker et al. 2020).In any case, sulphur reservoirs in dense regions remain as an open question, as well as the time at which sulphur depletion occurs.
Given the low number of discovered hot corinos and, therefore, the limited knowledge about them, the project Evolution of Chemistry in the envelope of HOt corinoS (ECHOS) aims to provide a complete inventory of the chemistry of hot corinos (with a sample including IRAS16293-2422, NGC 1333 IRAS 4A, HH 212, andB 335).ECHOS is carried out through homogeneous and systematic spectral surveys between 30 and 200 GHz using the Yebes-40m and IRAM-30m telescopes.In this paper, we focus on the analysis of B 335 and present a study of sulphurbearing molecules in this region.Observations and a description of the source are presented in Sects. 2 and 3.In Sect. 4 we present the data and calculate column densities and abundances using different methods, depending on the number of detected transitions for each species, and also depending on the availability of collisional coefficients.A discussion of the results (including a comparison with other similar objects, as well as the use of a chemical code to study the fractional abundance evolution of sulphur species) is presented in Sect. 5. We finally summarise our conclusions in Sect.6.

Observations
The observations were carried out with the Yebes-40m and IRAM-30m telescopes pointing to the position with coordinates α J2000.0 = 19 h 37 m 0.93 s , δ J2000.0 = +07 o 34 ′ 09.9 ′′ .For both telescopes, data were reduced and processed by using the CLASS and GREG packages provided within the GILDAS software1 , developed by the IRAM institute.

Yebes-40m telescope
The observations from the Yebes-40m radiotelescope located in Yebes (Guadalajara, Spain) were obtained using a receiver that consists of two cold high electron mobility transistor amplifiers covering the 31.0-50.3GHz Q-band with horizontal and vertical polarisations.The backends are 2×8×2.5GHz fast Fourier transform spectrometers with a spectral resolution of 38.15 kHz, providing the whole coverage of the Q-band in both polarisations.The observations were performed in the positionswitching mode with (-400 ′′ , 0 ′′ ) as the reference position.The main beam efficiency varies from 0.6 at 32 GHz to 0.43 at 50 GHz.Pointing corrections were derived from nearby quasars and SiO masers, and errors remained within 2-3 ′′ .

IRAM-30m telescope
The observations from the IRAM-30m radio telescope located in Pico Veleta (Granada, Spain) cover the 73-118 GHz and 130-175 GHz spectral ranges.These observations were carried out in one session in July 2021 with the Wobbler-Switching mode (±120 ′′ ), using the broad-band EMIR (E090 and E150) receivers and the fast Fourier transform spectrometer in its 200 kHz of spectral resolution.The intensity scale in antenna temperature (T ⋆ A , which is corrected for atmospheric absorption and for antenna ohmic and spillover losses) was calibrated using two absorbers at different temperatures and the atmospheric transmission model (ATM; Cernicharo et al. 1985;Pardo et al. 2001).Calibration uncertainties were adopted to be 10%.
To convert to main beam brightness temperature (T MB ), we used the expression where F eff is the telescope forward efficiency and B eff is the main beam efficiency 2,3 .See Table 1 for a summary of the telescope beam sizes depending on frequency.

The source
B 335 is an isolated dense globule associated with the embedded far-infrared source IRAS 19347+0727 (Keene et al. 1980(Keene et al. , 1983)), located at a distance of ∼100 pc (Olofsson & Olofsson 2009).The central source is a Class 0 source with a luminosity of 0.72 L ⊙ (Evans et al. 2015), whose mass is estimated to be 0.08 M ⊙ (Yen et al. 2015).Complex organic molecules (COMs), such as acetaldehyde (CH 3 CHO), methyl formate (HCOOCH 3 ), and formamide (NH 2 CHO), have been detected (Imai et al. 2016;Okoda et al. 2022) in the vicinity of the protostar (r∼10 au).It reveals that B 335 harbors a hot corino (Cazaux et al. 2003;Herbst & van Dishoeck 2009), which warms up the surrounding regions leading to the sublimation of grain mantles and injecting complex organic molecules into the gas-phase (e.g.Ceccarelli 2004;Caselli & Ceccarelli 2012).Emission of CCH and c-C 3 H 2 has also been detected around the protostar, which indicates the presence of carbon-chain molecules in this source (Sakai et al. 2008;Sakai & Yamamoto 2013).This carbon-chain molecules could suggest the presence of a so-called warm carbon-chain chemistry (WCCC, Sakai et al. 2008), characterised by various carbon-chain molecules being produced by reactions of CH 4 (evaporated from grain mantles in a lukewarm region) and where the formation of HC n N is slower than that of C n H m molecules (Sakai et al. 2008(Sakai et al. , 2009;;Oya et al. 2017).
Particularly interesting are the single-dish observations of CS and H 2 CO in B 335 (Zhou et al. 1993;Choi et al. 1995), since their optically thick lines show asymmetric profiles.This implies the presence of infalling motions in the envelope around B335.A more detailed analysis with C 18 O, SO, HCN, and HCO + ALMA observations (Evans et al. 2015;Yen et al. 2015;Evans et al. 2023) confirms an infalling gas motion of the inner part of the protostellar core, but rotation motion has not been detected at a spatial resolution of 0.3 ′′ (∼30 au).A rotation structure is only found at a scale of about 10 au from CH 3 OH and HCOOH emission (Imai et al. 2019).Single-dish and interferometric 12 CO (1-0) and 13 CO (1-0) observations (e.g.Hirano et al. 1988Hirano et al. , 1992;;Stutz et al. 2008;Bjerkeli et al. 2019) show also the presence of an east-west molecular outflow on ∼0.2 pc (Cabrit et al. 1988;Moriarty-Schieven & Snell 1989) with an inclination between 3 o -10 o on the plane of the sky and an outflow opening angle of 45 o (Hirano et al. 1988(Hirano et al. , 1992;;Stutz et al. 2008;Yen et al. 2010).A dynamical analysis of the outflow molecular structure (Yıldız et al. 2015) implies that dynamical timescales for the CO emitting gas are of the order of 10 4 yr.This result, along with the absence of a rotationally supported disk on scales greater than ∼10 au in size (Evans et al. 2015;Yen et al. 2015), suggest that the B 335 system is a particularly young source compared to other Class 0 objects, such as L 1527 which contains an embedded IRAS source and has a dynamical age of 2×10 4 -10 5 yr (Tamura et al. 1996;André et al. 2000;Agúndez et al. 2019).The detection of abundant CCS (see Sect. 4.3) and other carbon-chain molecules in B 335 also confirms the early stage of this source since it is generally recognised that CCS and other carbon-chain molecules are abundant 2 https://rt40m.oan.es/rayo/index.php 3http://www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies in the early stage of chemical evolution and become deficient at the advanced stage (Sakai et al. 2008).
So far, B 335 is one of the few hot corinos identified in a Bok globule, isolated from a large molecular cloud complex.This is therefore an ideal source for the study of the earliest stages of star formation with a chemical composition that could be regarded as a standard template for isolated protostellar cores.

Data analysis and results
We have detected 81 emission lines from 21 different sulphur-bearing species (including isotopologs) in B 335.In particular, we have detected three transitions of CS, three of 13 CS, three of C 34 S, three of C 33 S, nineteen of CCS, three of CC 34 S, seven of CCCS, three of CCC 34 S, five of OCS, three of HCS + , two of HC 34 S + , seven of SO, one of SO 2 , one of 34 SO, one of H 2 S, one of H 2 34 S, ten of H 2 CS, one of HSCN, one of HNCS, and six of NS (see .These transitions span an energy range of E up =2.0-65.3K.

Line profiles
All the lines are observed in emission.We first fitted the observed lines with Gaussian profiles using the CLASS software to derive the radial velocity (v LSR ), the line width, and the intensity for each line.Results are given in Table A.1.We observed that most of the detected molecules show narrow line profiles (∆v 1.5 km s −1 ) that can be adjusted by a single velocity component.
It is well known that the gravitational collapse of a molecular cloud when a new star is formed is accompanied by the development of highly supersonic outflows (e.g.Bachiller & Pérez Gutiérrez 1997;Holdship et al. 2019).The material ejected through these jets collides with the surrounding cloud, compressing and heating the gas, which leads to a drastic alteration of the chemistry as endothermic reactions become efficient and dust grains are partially destroyed.Many observations have shown that several sulphur-bearing species are enhanced in the shocked regions with respect to other areas of the cloud (e.g.Pineau des Forets et al. 1993;Bachiller & Pérez Gutiérrez 1997;Codella & Bachiller 1999;Wakelam et al. 2004Wakelam et al. , 2005;;Esplugues et al. 2013).However, the sulphur lines observed in B 335 present very narrow line widths (Figs.A.1-A.21), which is more consistent with the emission arising from the ambient quiescent cloud.Only for H 2 S, formed through the hydrogenation of sticked atomic sulphur on grains (Hatchell et al. 1998;Garrod et al. 2007;Esplugues et al. 2014;Oba et al. 2019) and thought to be the main sulphur reservoir in ices (Vidal et al. 2017;Navarro-Almaida et al. 2020), it is necessary to also consider a wider component (with ∆v=2.5 km s −1 ) to fit its only detected line, since the line profile shows the presence of faint wings.This could indicate that part of the H 2 S emission arises from the outflow, however, the line width of this component is much smaller (∼2.5 km s −1 ) than the ones observed for 12 CO (∼10-20 km s −1 ) in the outflow (Yen et al. 2010).Therefore, the wider H 2 S emission may be likely related to infall motions.
On the other hand, there are some species (CCS, CCCS, and H 2 CS) whose low transition lines present a double peak emission.The peak with the lowest intensity is located for all the cases at lower velocity (v LSR =7.39-8.1 km s −1 , blue peak) than for the more intense peak (v LSR =8.80-8.97km s −1 , red peak).Line widths are slightly larger (∆v 0.83-1.63km s −1 ) in the blue peaks than in the red ones (∆v 0.55-0.83km s −1 ).The comparison of these profiles with those from the less abundant isotopologues shows that these two-peak profiles are due to selfabsorption due to the large optical depths.All the detected lines present 0.005<T MB <1.5 K.

Rotational diagrams
For each molecule, we have computed a representative rotational temperature (T rot ) and a column density (N) by constructing a rotational diagram, assuming a single rotational temperature for all energy levels (Goldsmith & Langer 1999).Rotational diagrams represent a technique to analyse cloud properties from molecular line emission assuming local thermodynamic equilibrium (LTE).The standard relation for the rotational diagram analysis is ln( with N u /g u given by where N u is the column density of the upper level in the optically thin limit, N is the total column density, g u is the statistical weight of the upper state of each level, Z is the partition function evaluated at a rotational temperature T rot , E u /k is the energy of the upper level of the transition, ν ul is the frequency of the u→l transition, T MB dv is the velocity-integrated line intensity corrected from beam efficiency, and f is the beam filling factor. Assuming that the emission source has a 2D Gaussian shape, f is equal to f bf =θ 2 s /(θ 2 s + θ 2 b ) being θ b the HPBW of the telescope in arcsec and θ s the diameter of the Gaussian source in arcsec.In order to derive the values of f in each case, we have used the HPBW values listed in Table 1, and we have assumed a source size of 10 ′′ according to interferometric maps of C 17 O emission in B 335 (Cabedo et al. 2021).
The resulting rotational diagrams are shown in Fig. 1, while the column densities and rotational temperatures obtained are listed in Table 2.The uncertainties shown in Table 2 indicate the uncertainty obtained in the least squares fit of the rotational diagrams.These uncertainties include the 10% error in the observed line intensities due to calibration.The uncertainty obtained in the determination of the line parameters with the Gaussian fitting programme is included within the error bars at each point of the rotational diagram.We did not applied the rotational diagram method to NS since its four observed lines are close to the detection limit and/or blended with other species, which would introduce a high uncertainty in the results.The analysis of NS and others molecules with only one emission line detected (e.g.SO 2 ) has been done only with a LVG code (see next Section) instead.
The obtained rotational temperatures are 15 K.In the case of CCS, the molecule with the largest number of detected transitions, we observe a trend variation (Fig. 1) between the low and the high transition energies.In particular, the transitions with E u <40 K are well adjusted for a T rot ∼5 K, while transitions with E u >40 K are better adjusted for a T rot ∼15 K.This suggests that CCS arises from different regions characterised by distinct excitation temperatures.The rotational temperatures for the rest of molecules with several detected transitions (e.g.C 3 S, OCS, SO, H 2 CS) have 4.5 T rot 7.4 K.These low temperatures together  Numbers in parentheses are obtained from its corresponding 34 S isotopologue considering 32 S/ 34 S=24.4 (Mauersberger et al. 2004).
with their narrow line widths are consistent with emission arising from a cold and extended part of the dense core.Gas density at these scales is probably not too high as also indicated by the low T rot obtained for most of the S-molecules.All this suggests that the gas is sub-thermally excited with a moderate or low (<10 5 cm −3 ) density, in agreement with Frerking et al. (1987).In addition, the low number of transitions detected for several species (CS, 13 CS, C 34 S, C 33 S, HCS + , and HC 34 S + ) and the optically thick emission of some transitions of CCS, CCCS, and H 2 CS, also indicate that the rotational diagram method not suitable to derive column densities.

Column densities and abundances
We have therefore calculated molecular column densities for the detected S-species using a large velocity gradient (LVG) approximation.In particular, we have derived them by fitting the emission line profiles using the LVG code MADEX (Cernicharo 2012), which assumes that the radiative coupling between two relatively close points is negligible, and the excitation problem is local.The LVG models are based on the Goldreich & Kwan (1974) formalism.MADEX also includes corrections for the beam dilution of each line, depending on the different beam sizes at different frequencies.The final considered fit is the one that  reproduces more line profiles better from all the observed transitions.To carry out the fits, we assume uniform physical conditions (kinetic temperature, density, line width, radial velocity).
In particular, we have assumed a gas density of 5×10 4 cm −3 according to results from Sect.4.2 and to results from Cabedo et al. (2023), who deduced a density of ∼10 4 cm −3 for the outer re-gions of B 335 based on interferometric observations.We have also assumed a source size of 10 ′′ according to interferometric maps (Cabedo et al. 2021), and a kinetic temperature of the gas T K ∼15 K estimated by Shirley et al. (2011) from the formula for dust temperature considering only central heating by the B 335 protostar, and that dust and gas are coupled.
In the case of H 2 CS and H 2 S, we have used the code RADEX (van der Tak et al. 2007) since there are not available collisional rates in MADEX for these two molecules.The results are shown in Table 2.For those species for which one transition is observed ( 34 SO, SO 2 , H 2 S, H 2 34 S, HSCN, and HNCS), we have derived their column densities using only the LVG approximation and not the rotational diagram technique described in Section 4.2.For NS, since there are not available collisional rates, we have considered LTE approximation.
Column densities derived using the rotational diagram technique and LVG calculations agree within a factor 3. Sources of uncertainty are the low angular resolution of the used telescopes (with beam sizes between ∼36 ′′ and ∼54 ′′ for the 40m-Yebes telescope, and ∼14 ′′ and ∼29 ′′ for the 30m-IRAM telescope depending on frequency), which implies that the emission from the inner region of the B 335 is blended with the outer envelope, the possible overlap with the emission from other species (e.g. the case of H 2 34 S), the limited number of detected transitions (one transition in many cases, such as for H 2 S, H 2 34 S, SO 2 , HSCN, HNCS, and 34 SO), the lack of collisional rates for some species (e.g.NS), and the assumed source size derived from C 17 O interferometric observations.This is one of the most important uncertainties due to the lack of interferometric observations of S-species in B 335, since the emission from the different S-molecules may arise from distinct regions with slightly different sizes leading to a significant impact on the fits.In order to derive the uncertainty associated with the LVG results, we have run the LVG code varying the considered values for density, temperature, and source size.In particular, we have run models with density 5×10 4 cm −3 and 2.5×10 4 cm −3 (variation of a factor 2), temperatures of 15 K and 12 K, and source sizes of 10 ′′ and 8 ′′ (reduction of 20%).We found variations in the S-column densities 10% between these two T K which is comparable to the calibration errors, while the variations when changing the density and the source size were 20% and ∼30%, respectively.Taking all this into account, we estimate the uncertainties in the column densities derived with LVG approximation to be of the order of ∼50%.
The total column density of the sulphur-bearing species observed in B 335 is 6.2×10 14 cm −2 , with the highest contribution coming from CS (calculated from C 34 S considering 32 S/ 34 S=24.4 since CS is optically thick, Mauersberger et al. 2004), H 2 S, SO, and CCS.The lowest values are found for NS and HSCN.We have also observed several isotopologues relative to 13 C, 34 S, and 33 S. In particular, from CS, we obtain 34 S/ 33 S=10.7 (which is slightly higher than the solar abundance value 34 S/ 33 S=5.6-6.3,Anders & Grevesse 1989; Mauersberger et al. 2004), and 12 C/ 13 C=46.8,which is similar to that found in dark clouds by Cernicharo & Guelin (1987).
Table 2 also shows the fractional abundances relative to H 2 (with N H 2 obtained from Cabedo et al. 2021) of the sulphur molecules detected in the range λ = 2, 3, and 7 mm.Considering the detected S-molecules, we have obtained a total molecular sulphur abundance of 2.0×10 −8 , which is alike to the one found in similar sources, such as L 1544 (1.1×10 −8 , Vastel et al. 2018), L 483 (9.1×10 −9 , Agúndez et al. 2019), and B1-b (1.2×10 −8 , Fuente et al. 2016).The highest S-abundances are >10 −9 (CS, H 2 S, CCS, and SO), with also several relatively abundant (>10 −10 ) species, such as OCS and H 2 CS.In general, B 335 presents a great variety of S-bearing molecules, but less than other similar sources.In particular, although the total sulphur abundance in B 335 is similar to other sources, such as L 483, we detect a significant lower number of sulphur molecules in B 335 than in L 483.We are also missing some other S-  species already reported towards L 483, including SO + , HCS, and CH 3 SH, for which we have derived upper limits for their column densities of N(SO + )<1.5×10 11 cm −2 , N(HCS)<3×10 12 cm −2 , and N(CH 3 SH)<2×10 12 cm −2 .To derive these upper limits, once the physical parameters are fixed, we vary (for each species) the column density until the model fit reaches the observed intensity peak of any of the observed lines.We do not allow the model fit to be greater than any observed line.

Chemical comparison
Even though the total sulphur abundance in B 335 is similar to that from other Class 0 sources, the abundances of some specific molecules significantly vary among them.Figure 2 shows a comparison between abundances of different sulphur species in three Class 0 sources, B 335, L 483, and B1-b.The first thing that stands out is that sulphur carbon chains (CS, CCS, and CCCS) in B 335 are about one order of magnitude higher than the ones observed in L 483 and B1-b.Although high abundances of carbon-chain molecules have already been found in other cold clouds, such as TMC1 (Agúndez & Wakelam 2013), it is not a universal characteristic of this type of early regions, since many of them present different degrees of carbon chain richness (e.g.Suzuki et al. 1992;Hirota et al. 2009).This therefore highlights the nature of B 335 as a source especially rich in sulphur carbon chains.
Unlike carbon chains, sulphur molecules containing oxygen in B 335 have a similar (OCS) or even lower (SO and SO 2 ) abundance by up to one order of magnitude than those found in L 483 and B1-b (see also Fig. 3).These oxygen-sulphur molecules are usually detected in cold dense regions with abundances between 10 −10 -10 −8 , while in hot cores and hot corinos they are significantly enhanced (e.g.Tercero et al. 2010;Esplugues et al. 2013).In particular, SO 2 is a better tracer of warm gas than SO (Esplugues et al. 2013), while SO is a wellknown outflow tracer (e.g.Codella & Scappini 2003;Lee et al. 2010;Tafalla et al. 2010).SO seems to be more enhanced than SO 2 in shocks for timescales ∼10 4 yr (e.g.Codella & Bachiller 1999;Viti et al. 2004;Jiménez-Serra et al. 2005)  abundance found in B 335 with respect to the other two sources L 483 and B1-b suggests that its chemistry may not be significantly affected by shocks associated with the bipolar outflow (e.g.Hirano et al. 1988Hirano et al. , 1992;;Stutz et al. 2008;Bjerkeli et al. 2019).
Other sulphur molecules detected in B 335 with low abundances (Fig. 3) are those containing nitrogen, that is NS and HSCN (4×10 −11 and 3×10 −11 , respectively).HSCN, unlike HNCS which is the most stable isomer, is metastable and only detected on Earth under specific conditions (through UVphotolysis of HNCS (Wierzejewska & Mielke 2001) or through a low-pressure discharge (Brünken et al. 2009)).In contrast, HSCN is more easily observed in the interstellar medium.For instance, in Sgr B2 (where HSCN was first time detected) the ratio HNCS/HSCN∼3 (Halfen et al. 2009) and in TMC1 this ratio is ∼1 (Adande et al. 2010), while in L 1544 Vastel et al. (2018) only detect HSCN, but not HNCS.This shows the important role of non-equilibrium chemistry in dense clouds.In B 335, we only detect one transition of HNCS with a signal to noise ratio σ∼3.1, so we have considered its column density (N HNCS ≤4.7×10 12 cm −2 ) as an upper limit.This leads to a ratio HNCS/HSCN<4.7 consistent with the one found in Sgr B2.

Sulphur ions
Apart from neutral sulphur species, we have also observed the presence of some ions in B 335.Molecular ions are commonly observed in dense prestellar or protostellar regions (e.g.Caselli et al. 1998;Hogerheijde et al. 1998;Caselli et al. 2002;van der Tak et al. 2005).To a lesser extent, ions (such as HCO + and N 2 H + ) have also been observed in protostellar shocks (e.g.Bachiller & Pérez Gutiérrez 1997;Hogerheijde et al. 1998;Caselli et al. 2002;van der Tak et al. 2005).The presence of ions in shocked regions is a balance between their enhancement by the sputtering of dust grain mantles, which ejects molecules into the gas phase that undergo reactions forming ions, and their destruction by electronic recombination (e.g.Neufeld & Dalgarno 1989;Viti et al. 2002).For this reason, those ions that usually decrease in shocks (such as N 2 H + ) may indicate the presence of a pre-shock chemistry (Codella et al. 2013), while the ions that are enhanced by shock chemistry will be effective tracers of shocked gas.This is the case of HOCO + and SO + (Podio et al. 2014) since their abundances are predicted to be very low in the quiescent gas, but significantly enhanced in shocks following the release of CO 2 and S-species from dust grain mantles (Minh et al. 1991;Turner 1992Turner , 1994;;Deguchi et al. 2006).
Among the sulphur ions, in B 335 we do not detect SO + (strong tracer of shocked gas which forms from S + and OH Neufeld & Dalgarno 1989), but we detect HCS + and HC 34 S + with abundances of 4.7×10 −10 and 1.9×10 −11 , respectively.The value we have obtained for HCS + is similar to the one observed in the bow shock L 1157-B1 (Podio et al. 2014) and one order of magnitude higher than the ones obtained by Agúndez et al. (2019) and Fuente et al. (2016) in the Class 0 objects L 483 and B1-b, respectively, as also found for CS.The abundances of HCS + and CS are strongly related, since this ion forms through the reaction of CS with HCO + , H + 3 , H 3 O + (Millar et al. 1985).In addition, large abundances of CS and HCS + may also indicate that OCS is one of the main sulphur carriers on dust grains upon sputtering from the grain mantles according to previous theoretical (Podio et al. 2014) and observational (Wakelam et al. 2005;Codella et al. 2005) results.
Regarding the ratio R=HCS + /CS, we obtain from observations a value of 0.043.Thaddeus et al. (1981) determined an abundance ratio R∼0.01-0.03 in ten molecular sources located in Sgr B2 and Ori A, while Millar (1983) found that R could be as large as 0.1 in cold dark clouds.This result was also verified by Irvine et al. (1983).A similar R value to that from B 335 is found as well in the pre-stellar core L 1544 (0.03, Vastel et al. 2018), while in the protostellar sources L 483 and B1-b the HCS + /CS ratio is 0.017 and 0.008, respectively (Agúndez et al. 2019;Fuente et al. 2016).According to Clary et al. (1985), rate coefficients for ion-CS reactions rapidly increase at low temperatures suggesting that the sources B 335 and L 1544 are characterised by lower temperatures than L 483 and B1-b, which can be due to an earlier evolutionary stage of B 335 and L 1544 compared to L 483 and B1-b.

Time evolution of fractional abundances
Figures 4-7 show the evolution of several S-bearing species along 10 Myr for different physical conditions of the cosmicray (CR) ionisation rate (ζ H 2 ), gas temperature (T g ), hydrogen number density (n H ), and initial sulphur abundance (S + init ).Given that the CR ionisation rate is still uncertain and spreads over a range of values (van der Tak & van Dishoeck (2000) determined as ζ H 2 ∼3×10 −17 s −1 in dense clouds, Indriolo & McCall (2012) obtained a range (1.7-10.6)×−16 s −1 in a sample of diffuse clouds, and Neufeld & Wolfire (2017) derived a ζ H 2 of the order of a few 10 −16 s −1 in the Galactic disk), we have considered two values for the cosmic ionisation rate ζ H 2 =1.3×10 −16 s −1 and 1.3×10 −17 s −1 .For the gas temperature and density, given the results obtained in Sects.4.2 and 4.3, we have run models considering T g =7 K and 15 K, and n H =2×10 4 cm −3 and 2×10 5 cm −3 .For the case of the initial sulphur abundance, we have considered the solar elemental sulphur fractional abundance (S + init =1.5×10 −5 ) and also one factor of S depletion (1.5×10 −6 ) since to reproduce observations in hot corinos one needs to assume a significant sulphur depletion of at least one order of magnitude lower than the solar elemental sulphur abundance as previously stated in Sect. 1.
We have obtained these models using the Nautilus timedependent chemical code (Ruaud et al. 2016).Nautilus is a three-phase model in which gas, surface icy mantle, bulk icy Fig. 4. Evolution of fractional abundances of CS, CCS, C 3 S, OCS, SO, SO 2 , HCS + , H 2 CS, and NS as a function of time for an initial sulphur abundance S + init =1.5×10 −6 , a hydrogen number density n H =2×10 4 cm −3 , T gas =7 K, and two values of CR ionisation rate ζ=1.3×10 −17 s −1 (solid line) and ζ=1.3×10 −16 s −1 (dashed line).mantle, and their interactions are considered.Nautilus solves the kinetic equations for both the gas phase species and the surface species of interstellar dust grains and computes the evolution with time of chemical abundances for a given physical structure.The chemical network is based on the KInetic Database for Astrochemistry (KIDA 4 ).The used version here is detailed in Wakelam et al. (2021).In particular, this Nautilus version considers gas-phase processes including neutral-neutral and ionneutral reactions, direct cosmic ray ionisation or dissociation, ionisation or dissociation by UV photons, ionisation or dissociation produced by photons induced by cosmic-ray interactions with the medium (Prasad & Tarafdar 1983), and electronic recombinations.For species on the surfaces, there is a distinction between species in the most external layers (surface species) and species below these layers (mantle species).The species are adsorbed on the surface and become the mantle during the construction of the ices.Similarly, when the species desorb, only the species from the surface can desorb, but surface species are gradually replaced by the mantle species.Regarding surface reactions, the model includes thermal desorption, photodesorption, chemical desorption, and cosmic-ray heating.See 4 http://kida.obs.u-bordeaux1.fr/Wakelam et al. (2021) for extensive details about the surface parameters considered in each surface reaction.In all models, we have adopted the initial abundances shown in Table 3 and a visual extinction A V =15 mag 5 .
In general, we observe that the abundances of all the species considered in Figs.4-7 increase with time up to t∼5×10 5 yr, moment at which the abundances of various species decrease abruptly.This is the case of CS, SO, and SO 2 , whose abundances decrease by about 2-4 orders of magnitude between t∼5×10 5 yr and t∼5×10 6 yr since they react with other species to form more complex molecules.In particular, for t>10 5 yr, CS SO, and SO 2 are mainly destroyed by reacting with H + 3 and forming HCS + , HSO + , and HSO + 2 , respectively.For the abundances of the rest of sulphur species, this decrease after reaching the maximum value is lower and not greater than two orders of magnitude.From these figures, we also deduce the formation time scales 5 Derived from results in Table 2 and the expression where N i is the column density of the species i, X i its abundance, A V the visual extinction, and 1.6×10 21 the hydrogen column density at 1 mag of extinction (Bohlin et al. 1978). G.
of S-species.In particular, we have observed that CS is one of the molecules reaching first its maximum abundance value in all the models at t∼5×10 4 yr, followed by HCS + with its maximum abundance reached at t∼10 5 yr.Other molecules, such as H 2 CS, can be considered as late molecules, since their maximum abundance values are found when the chemical evolution exceeds one million years.
Regarding the influence of the physical conditions (cosmicray ionisation rate, temperature, density, and initial sulphur abundance) on the S-abundances, the parameter with the highest impact is the hydrogen nuclei number density, n H , since varying n H by one order of magnitude leads to abundance variations of up to three orders of magnitude for CS, CCS, OCS, HCS + , SO, and SO 2 given a specific evolutionary time (see Fig. 6).In general terms, we observe in Fig. 6 that the chemical evolution of each considered species is accelerated when the density increases since collisions are more frequent than in a low density regime.This results in the abundance peak of the species being reached at earlier times for larger densities, although the chemical behaviour along time is roughly maintained.In particular, for t∼10 4 yr, we obtain that SO 2 is mainly formed by the reaction O+SO for n H =2×10 4 and 2×10 5 cm −3 , while for t∼10 5 yr (when the SO 2 peak is reached) this molecule mainly forms by the reaction O+SO for n H =2×10 4 cm −3 and by OH+SO for n H =2×10 5 cm −3 .The most responsible species of destroying SO 2 for both densities is C for t∼10 4 yr, and H + 3 for t≥10 5 yr.Regarding CS, for t∼10 4 yr, it is mainly formed by the reaction HCS + +e − for the low density value, while CS is mainly formed through the reaction between atomic S and C 4 for n H =2×10 5 cm −3 .For t∼10 5 yr (when the CS abundance decreases for any model), CS is mainly destroyed by its reaction with H + 3 and with HCO + .Another parameter with a big impact on S-abundances is the cosmic-ray ionisation rate, ζ, since the variation from ζ=1.3×10 −17 to 1.3×10 −16 s −1 changes by about two orders of magnitude the abundances of all the sulphur species considered in the sample (Fig. 4).For the case of carbon-chains (CCS, C 3 S), we observe that the larger ζ, the larger the S-abundances for t 10 5 yr.Similar results are found for the species SO, SO 2 , NS, H 2 CS, and OCS.For t∼10 4 yr, we find that CCS is mainly formed by C+HCS for the low ζ model, while it is mainly formed through HC 3 S + +e − in a more efficient way for a model with larger ζ (1.3×10 −16 s −1 ).In the case of SO (one of the most affected molecules by ζ between 10 4 t 10 6 yr), it is mostly formed by the reaction O+HS in a model with low ζ, while the reaction between S+OH becomes more important to form SO 2 when increasing ζ by one order of magnitude.
On the other hand, changing the temperature (from T gas =7 K to 15 K) and especially the initial sulphur abundance (from S + init =1.5×10 −5 to 1.5×10 −6 ) in the models leads to abundance variations smaller than two orders of magnitude for t>10 4 yr.In particular, we have found that an increase of the S + init value by one order of magnitude also enhances the S-abundances by one order of magnitude during the entire duration of core evolution.
For the particular case of HCS + and CS, the main physical parameters significantly influencing their abundances are the density (the lower density, the higher the HCS + and CS abundances for t>10 4 yr) and the cosmic-ray ionisation rate (the higher the ζ, the higher the HCS + and CS abundances for t>5×10 5 yr).By contrast, other physical parameters, such as the gas temperature, barely affect the HCS + and CS abundances for t<10 6 yr.

Comparison with observations
We now compare the different models considered in the previous section with observational results to derive the current evolutionary stage of B335.show the ratio between the abundances obtained from the models (X model ) and from the observations (X obs ) for several sulphur-bearing species.In particular, we have considered those S-species for which we observe more than one emission line and those that are not close to the limit of detection.The results are shown for five specific evolution times (t=10 4 , 5×10 4 , 10 5 , 5×10 5 , and 10 6 yr).In order to find the cases with the closest values be-Fig.8. Evolution of fractional abundances of CS, CCS, C 3 S, OCS, SO, SO 2 , HCS + , H 2 CS, and NS as a function of time for an initial sulphur abundance S + init =1.5×10 −6 , a hydrogen number density n H =2×10 4 cm −3 , T gas =7 K, and two values of CR ionisation rate ζ=1.3×10 −17 s −1 (solid line) and ζ=1.3×10 −16 s −1 (dashed line).tween the theoretical and observational results, we have considered for each model the largest number of S-species for which -0.5≤log(X model /X obs )≤0.5 in order to only consider discrepancies between model and observational results smaller than a factor of ∼3.Results show that the evolutionary stages that better reproduce the observed abundances of the largest number of Sspecies in B 335 are those between t=10 4 -10 5 yr.This result obtained from a chemical approach agrees with the results obtained from Evans et al. (2015) who derived an age of ∼5×10 4 yr for B 335 through a dynamical perspective, by applying an inside-out collapse model that is compared with ALMA observations.Altogether this reveals, therefore, the particularly early evolutionary stage of B 335, which is comparable to the ages of pre-stellar condensations (e.g.Caselli & Ceccarelli 2012).As B 335 is an isolated globule, this may be due to the fact that the condensation has formed from very little dense gas (∼10 3 cm −3 ) without a dense pre-phase between the formation of the globule and the beginning of the collapse.This gives B 335 a different chemistry compared to other young protostars that have form in dense molecular clouds, such as L 483 and B1-b.Another plausible scenario is that material from the diffuse cloud surrounding B 335 may currently be accreting onto the protostellar envelope, thus rejuvinating its chemistry.This accretion has been found (e.g.Pineda et al. 2020) in another Class 0 source embedded in the Perseus Molecular Cloud Complex. Interferoetric observations are requested to distinguish between both scenarios.From Figs.A.22-A.25, we have also deduced that chemical models with values of n H =2×10 4 cm −3 , T g =15 K, S + init =1.5×10 −6 , and ζ=1.3×10 −17 s −1 best reproduce the observations of the largest number of S-species for t=10 4 -10 5 yr.This value of ζ agrees with results from Cabedo et al. (2023) who found a cosmic-ray ionisation rate 10 −16 s −1 in the outer regions of B 335, especially for an envelope radii>500 au.A sulphur depletion factor of 10 with respect to the sulphur cosmic elemental abundance is also consistent with results from Fuente et al. (2023) in low-mass star-forming regions, such as Taurus and Perseus.

Sulphur chemical ratios
We have considered several sulphur-bearing molecular ratios and have also analysed their evolution with time.Model results obtained with the Nautilus code are shown in Figures 8-11.We have also changed some of the physical conditions (cosmic-ray Fig. 11.Evolution of fractional abundances of CS, CCS, C 3 S, OCS, SO, SO 2 , HCS + , H 2 CS, and NS as a function of time for an initial sulphur abundance S + init =1.5×10 −6 , a hydrogen number density n H =2×10 4 cm −3 , T gas =7 K, and two values of CR ionisation rate ζ=1.3×10 −17 s −1 (solid line) and ζ=1.3×10 −16 s −1 (dashed line).ionisation rate, temperature, density, and initial sulphur abundance) in order to study their impact on the evolution of these ratios.We first notice that while there are some ratios, such as HCS + /CS, that present small variations of up 1 order of magnitude between t=10 4 -10 6 yr, there are other ratios that significantly change with time.In particular, the SO 2 /C 2 S, SO/CS, SO 2 /CS, NS/CS, H 2 S/CS, and H 2 CS/CS ratios show variations of up to ∼5 orders of magnitude in that time range (10 4 -10 6 yr).
The physical parameter affecting sulphur ratios the most is the density, leading to differences of up to ∼4-5 orders of magnitude especially between 5×10 4 and 10 6 yr.The cosmic-ray ionisation rate is the next physical parameter with the greatest effect on the sulphur ratios with differences of about one order of magnitude in their values.In particular, the SO 2 /CS, SO/CS, SO 2 /C 2 S, and C 2 S/CS ratios are among those mostly influenced by ζ, while the SO/SO 2 ratio is one of the least affected.The SO/SO 2 ratio is, however, more sensitive to the density variation.We also observe a large density impact on the NS/CS ratio, finding that the higher the density, the larger the NS/CS for t∼10 5 -10 6 yr, with NS being mostly formed during that time range through the neutral-neutral reaction N+HS and, to a lesser extent, through the dissociative recombination of HNS + .On the other hand, varying the gas temperature from 7 to 15 K (Fig. 9) mostly affects the H 2 S/SO ratio between 10 4 and 10 6 yr, and the SO/CS, SO 2 /CS, and H 2 S/CS ratios at late times (t 10 6 yr).Regarding the only ion detected in B 335, its HCS + /CS ratio slightly increases when the initial sulphur abundance is decreased, and also when the cosmic-ray ionisation rate or the density increase.
The use of molecular ratios as indicators of chemical evolution of young stellar sources has been previously approached from both observational and theoretical perspectives (e.g.Pineau des Forets et al. 1993;Charnley 1997;Hatchell et al. 1998;Codella & Bachiller 1999;Gibb et al. 2000;Boogert et al. 2000;Wakelam et al. 2004).Our models  show that SO and SO 2 constantly increase with time reaching both their peaks at t 5×10 5 yr, which prevents the use of their ratio as a chemical clock.Nevertheless, show that there are other sulphur ratios that could be used as evolutionary tracers since they maintain a constant trend over time.This is for instance the case of the SO 2 /C 2 S and SO/CS ratios, whose values constantly increase during the first ∼2×10 5 yr.In order to compare these theoretical results with observations, we have considered a sample of sources that includes starless cores (TMC1 and L 1544 considered as starless and prestellar cores, respectively, see for example Crapsi et al. 2005;Schnee et al. 2007;Lefloch et al. 2018;Agúndez et al. 2019;Cernicharo et al. 2021) and Class 0 objects (B 335, L 483, and B1-b).Observational results are shown in Fig. 12, where we have indeed found a clear trend where the SO 2 /C 2 S and SO/CS ratios increase with the age of the objects, suggesting that these sulphur ratios could be used as chemical clocks.We have also included in Fig. 12 observational results for the HCS + /CS ratio. I this case, we have obtained a negative trend in which the value of HCS + /CS decreases with the age of the objects.
Comparing the observational HCS + /CS, SO 2 /C 2 S, and CO/CS ratios for B 335 with the theoretical results (Figs. 8-11), we observe that, for the temperature and density values used in Sect.4.3 to derive observational results, the three ratios would be reproduced at 5×10 4 t 10 5 yr when considering a high cosmicray ionisation rate (ζ=1.3×10−16 s −1 ) and a S depletion of at least a factor 10 with respect to the solar elemental sulphur abundance.This depletion is also found by many authors (e.g.Wakelam et al. 2004;Crockett et al. 2014;Vastel et al. 2018;Bulut et al. 2021;Navarro-Almaida et al. 2021;Esplugues et al. 2022;Hily-Blant et al. 2022;Fuente et al. 2023) in other starforming regions.

Summary and conclusions
We have carried out a comprehensive observational and theoretical study of sulphur in the Class 0 object B 335.This source is an isolated dense Bok globule located at a distance of ∼100 pc.We have detected 20 different sulphur-bearing species (including isotopologues) using Yebes 40m and IRAM 30m telescope observations.In order to derive column densities, we have applied the method of rotational diagrams.For this purpose, we have also used non-LTE radiative transfer codes given the low number of transitions detected for some species and the optically thick emission of some lines.Our results show low rotational temperatures ( 15 K), suggesting that the gas is sub-thermally excited, and thus revealing a relatively low (<10 5 cm −3 ) gas density.These low temperatures, the narrow line widths, and the relatively low density are also consistent with emission arising from a cold and extended region, suggesting that the detected S-emission does not trace the hot corino, but its cold envelope.
From the column densities of all the observed sulphur species, we have derived a total molecular sulphur abundance of 2.0×10 −8 , which is similar to that found in the envelopes of other Class 0 objects.Nevertheless, we have also found that the abundances of some specific S-molecules vary significantly across the sample of Class 0 objects (B 335, L 483, and B1-b) considered in this paper.In particular, we have derived high abundances of sulphur carbon-chain molecules in B 335 compared to those found in other Class 0 objects, and also compared with those abundances of S-molecules containing oxygen or nitrogen.This points out the nature of B 335 as a source especially rich in sulphur carbon chains, which could be due to the fact that B 335 is an isolated source compared to other young protostars that have formed in dense molecular clouds, such as L 483 and B1-b, or to accretion of material from the diffuse cloud surrounding B 335 onto the protostellar envelope, thus rejuvenating its chemistry.The comparison between envelopes of different Class 0 objects has also allowed to deduce a low SO abundance in B 335 with respect to L 483 and B1-b.SO is a well-known outflow tracer, which is thought to be enhanced in shocks.Its low abundance in B 335 together with the non-detection of other shock tracers, such SO + , suggests that its chemistry in B 335 is not particularly influenced by shocks.
In this work, we have also used the time-dependent chemical code Nautilus to study the time evolution of fractional Sabundances, as well as the influence of several physical and chemical parameters (cosmic-ray ionisation rate, density, gas temperature, and initial sulphur abundance) on these abundances.The comparison between chemical model results and observations shows that the evolutionary stages that better reproduce the observed abundances of the largest number of S-species in B 335 correspond to t=10 4 -10 5 yr, in agreement with previous age estimations obtained by applying theoretical methods with the use of a dynamical model.The results from Nautilus have also revealed that the physical parameters with the highest influence on the S-abundances are the density and the cosmic-ray ionisation rate whose variations by one order of magnitude lead to abundance variations of several orders of magnitude for different sulphur-bearing species such as SO, SO 2 , CS, and CCCS.We have also used Nautilus to explore the behaviour of different S-molecular ratios with time, finding that there are some sulphur ratios that could be used as chemical evolution tracers since they maintain a constant trend over time.A comparison of the observational results of S-ratios in a sample of objects (including starless and prestellar cores, and Class 0 objects) characterised by different evolutionary stages confirms a clear trend where the SO 2 /C 2 S and SO/CS ratios increase with the age of the objects, while HCS + /CS decreases with the evolutionary stage.We, therefore, conclude that these S-ratios could be used as good chemical evolutionary indicators tracing the prestellar to protostellar transition.A.22.Ratios between the abundances obtained from the models, X model , and the observations, X obs , for different sulphur species at specific evolution times (t=10 4 , 5×10 4 , 10 5 , 5×10 5 , and 10 6 yr), when varying the cosmic-ray rate (ζ).A.23. Ratios between the abundances obtained from the models, X model , and the observations, X obs , for different sulphur species at specific evolution times (t=10 4 , 5×10 4 , 10 5 , 5×10 5 , and 10 6 yr), when varying the gas temperature.A.24. Ratios between the abundances obtained from the models, X model , and the observations, X obs , for different sulphur species at specific evolution times (t=10 4 , 5×10 4 , 10 5 , 5×10 5 , and 10 6 yr), when varying the hydrogen number density (n H ).  A.25. Ratios between the abundances obtained from the models, X model , and the observations, X obs , for different sulphur species at specific evolution times (t=10 4 , 5×10 4 , 10 5 , 5×10 5 , and 10 6 yr), when varying the initial sulphur abundance (S + init ).

Fig. 1 .
Fig. 1.Rotational diagrams of the detected sulphur-bearing molecules in B 335. Fitted values of the rotational temperature, T rot , column density, N, and their respective uncertainties are also indicated for each molecule.

Fig. 3 .
Fig. 3. Fractional abundance (X) ratios between B 335 and L 483, and B 335 and B1-b for all the detected sulphur-bearing species in B 335.

Fig. A. 1 .
Fig. A.1.Observed lines of CS in B 335 (black) with rms>3σ.The red dashed line indicates the systemic velocity of the source v LSR =8.3 km s −1 .

Fig. A. 6 .
Fig. A.6.Observed lines of CC 34 S in B 335 with rms>3σ.The red dashed line indicates the systemic velocity of the source v LSR =8.3 km s −1 .

Fig. A. 9 .
Fig. A.9. Observed lines of OCS in B 335 with rms>3σ.The red dashed line indicates the systemic velocity of the source v LSR =8.3 km s −1 .

Fig. A. 10 .
Fig. A.10. Observed lines of HCS + in B 335 with rms>3σ.The red dashed line indicates the systemic velocity of the source v LSR =8.3 km s −1 .

Fig. A. 11 .
Fig. A.11. Observed lines of HC 34 S + in B 335 with rms>3σ.The red dashed line indicates the systemic velocity of the source v LSR =8.3 km s −1 .

Fig. A. 14 .
Fig. A.14. Observed lines of o-H 2 CS in B 335 with rms>3σ.The red dashed line indicates the systemic velocity of the source v LSR =8.3 km s −1 .

Fig. A. 20 .
Fig. A.20. Observed line of HSCN in B 335 with rms>3σ.The red dashed line indicates the systemic velocity of the source v LSR =8.3 km s −1 .

Table 1 .
Yebes 40 m and IRAM 30 m telescope efficiency data along the covered frequency range.

Table 2 .
Column densities (N), rotational temperatures, and abundances (X) relative to H 2 of the sulphur-bearing molecules detected in B 335.
Esplugues et al.:Extremely young sulphur chemistry in the isolated Class 0 object B 335 G.
. The low SO G. Esplugues et al.: Extremely young sulphur chemistry in the isolated Class 0 object B 335

Table 3 .
Esplugues et al.:Extremely young sulphur chemistry in the isolated Class 0 object B 335 Abundances with respect to total hydrogen nuclei considered in the chemical code Nautilus.