Free Access
Issue
A&A
Volume 597, January 2017
Article Number A74
Number of page(s) 16
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201628608
Published online 03 January 2017

© ESO, 2017

1. Introduction

Pre-stellar cores within molecular clouds constitute the stage previous to the star formation process. They provide the very initial conditions in which stars are formed. They are typically cold (10 K) with temperature gradients decreasing towards the center (e.g., Ruoskanen et al. 2011; Wilcock et al. 2012; Launhardt et al. 2013), and densities >104 cm-3, with a chemistry affected by molecular freeze-out onto dust grains (e.g., Tafalla et al. 2002; Marsh et al. 2014). Their physical structure can be well explained through a Bonnor-Ebert profile (Bonnor 1956; Ebert 1955), that is, the profile corresponding to an isothermal gas sphere in hydrostatic equilibrium (e.g., Evans et al. 2001; Kandori et al. 2005; Roy et al. 2014). However, their density and temperature gradient structures are still a matter of debate (e.g., Sipilä et al. 2011; Hardegree-Ullman et al. 2013) and the implications of heavy molecular depletion on other physical properties, such as dust grain properties and polarization, are not clear.

FeSt 1-457 (Feitzinger & Stuewe 1984), also known as Core 109 in the Lombardi et al. (2006) catalog, is a prototypical starless core, as previous studies have shown that this core is quiescent with no IRAS nor Spitzer Space Telescope point sources associated (Forbrich et al. 2009, 2015; Ascenso et al. 2013). Indeed, the internal luminosity, estimated using the non-detection from the 70 μm Herschel1 image (3 × rms = 0.02 Jy pixel-1) and using the relation between the 70 μm flux and the internal luminosity (Dunham et al. 2008), is Lint< 0.004L. It is also the most chemically evolved starless core among those isolated, embedded in the most pristine part of the Pipe nebula, the bowl, at a distance of 145 pc (Alves & Franco 2007). Frau et al. (2015) suggest the presence of two filaments with north-south and east-west directions colliding along the northwest-southeast direction at the bowl region, just where FeSt 1-457 is located.

Frau et al. (2012a) carried out a 15 GHz bandwidth molecular survey at 3 mm towards FeSt 1-457 with the Institute of Millimeter Radioastronomy (IRAM) 30 m telescope. They find a rich chemistry (specially compared with most of the other Pipe starless cores): many lines are detected from early- (C3H2, HCN, CS, C34S and CN) and late-time (N2H+, N2D+ and DCO+) molecules, suggesting that FeSt 1-457 is an evolved starless core. Aguti et al. (2007) observed the 3 mm lines of N2H+, HCO+, CS and C18O with the IRAM 30 m. They find evidence of CO and HCO+ depletion in the center of the core. However, the N2H+ (10) may be undepleted to depths of 40 mag of visual extinction.

FeSt 1-457 has a mass of ~4 M (Frau et al. 2010; Román-Zúñiga et al. 2010), quasi spherical and compact structure and shows signs of gravitational instability (Kandori et al. 2005; Frau et al. 2010). Aguti et al. (2007) find that FeSt 1-457 might be pulsating, based on expansions of the outer layers. However, their Jeans mass measurement is compatible with the mass of the core and they propose a quasi-stable state near hydrodynamic equilibrium.

The core is embedded in a magnetised medium (Alves et al. 2008, 2014; Franco et al. 2010), which suggests that the magnetic field could be a source of external support. Alves et al. (2014), in a multi scale polarization study towards FeSt 1-457, find a polarization “hole” for radii 55 arcsec which they propose to be the result of loss of grain alignment with the magnetic field due to the lack of an internal source of radiation.

In this work we study the starless core FeSt 1-457 through N2H+ (10) observations obtained with the IRAM 30 m and Plateau de Bure Interferometer (PdBI) telescopes2. N2H+ is known to be a very good tracer of the densest inner parts of dense cores as it is not expected to be depleted until densities >105 cm-3 (e.g., Bergin & Tafalla 2007). In addition, we present mapping information of 15 molecular lines detected at 3 mm within a frequency range spanning 30 GHz obtained with the IRAM 30 m which allowed us also to perform a detailed study of the abundance radial profiles and to relate them to the known grain growth and depolarization in the center.

2. Observations

2.1. IRAM 30 m

2.1.1. Line observations

Table 1

FTS molecular transitions.

The observations were carried out on July 2012 at the IRAM 30 m telescope in Granada, Spain. We used the Eight MIxer Receiver (EMIR) in the E090 configuration. The backend used was the FTS spectrometer with a channel resolution of 195 kHz (0.6 km s-1 at 3 mm), which provided a total bandwidth of 30 GHz covering the frequency ranges from 78.6 to 86.4 GHz, from 87.6 to 102.0 GHz and from 103.2 to 111.0 GHz (see Table 1). In addition, the VESPA (Versatile Spectrometer Array) correlator was used simultaneously with the FTS in order to observe the N2H+ (10) line with a better channel resolution, 20 kHz (corresponding to 0.063 km s-1), with a total bandwidth of 40 MHz. The observations were done using the frequency-switching mode with a frequency throw of ~7 MHz. The maps were done in on-the-fly (OTF) mode covering a region of 4.5 arcmin × 4.5 arcmin centered at RA(J2000)=17h35m47.s70\hbox{$)=17^{\rm h} 35^{\rm m} 47\fs70$}, Dec(J2000\hbox{$)=-25^{\circ}32'52\farcs9$}. In the OTF mode, the maps were done with scans either along the right ascension or either along the declination (with Nyquist separation between scans). This allows us better sampling on the observed region. The maps were combined and regrided to have a final pixel size of 6 arcsec. Pointing and focusing were done following the recommended procedure at the telescope, which should give a pointing accuracy of ~2 arcsec. Absolute calibration accuracy with EMIR is better than 10%. In this paper we present the spectra in main beam brightness temperature scale, which was obtained by using the efficiencies reported on the IRAM website. The weather conditions were good with system temperatures of 125–192 K and zenith opacities around 0.07. The reduction of the final maps and most of the figures have been done with the gildas3 software. Given the highly distorted baselines generated in the frequency switching mode, the spectral baselines were subtracted only over a small frequency range around a given line. The typical rms noise per channel for the VESPA spectral map is ΔTmb ≃ 0.18 K. The half power beam width (HPBW) for each frequency is given in Table 1.

2.1.2. Continuum observations

FeSt 1-457 was observed on April and May 2009 and on January 2010 with the 117-receiver Max-Planck Millimetre Bolometer (MAMBO-II) of the IRAM 30 m telescope (previously described in Frau et al. 2010).

2.2. PdBI

The PdBI observations of N2H+ (10) were performed on several runs in October and November 2011, and March 2012, in the C array configuration, which provides a baseline range from 24 m to 176 m (7.5–55 kλ at 3.2 mm). The observations were done in mosaic mode with an hexagonal pattern of 12 pointings and with a separation of ~18 arcsec. This allows to uniformly cover an area of ~102 × 102 arcsec2 centered on RA(J2000)=17h35m47.s804\hbox{$)=17^{\rm h}35^{\rm m}47\fs804$}, Dec(J2000\hbox{$)=-25^{\circ}33'01\farcs10$}. The primary beam at 3.2 mm is ~54 arcsec. The quasars J1733130 and J1743169 were used as gain calibrators and 3C 273, J2015+371, 3C 454.3 and J1751+096 as bandpass calibrators. Absolute flux calibration was performed on MWC349 and on the suitable solar system object available at the time of the observations. The absolute flux accuracy is ~10%. The N2H+ (10) line was observed in a spectral window of 20 MHz with 512 channels, which provides a spectral resolution of 0.039 MHz (0.125 km s-1). The system temperatures and the precipitable water vapor values were between ~100–300 K and ~2–6 mm respectively. The data were reduced using clic and mapping packages of the gildas software using the standard procedures. The rms noise per channel achieved was 25 mJy beam-1. The synthesized beam was 10.6 × 5.7 arcsec2, with a PA of 11.94°.

2.3. PdBI+IRAM 30 m

We combined the N2H+ single dish data taken with the IRAM 30 m telescope (with an rms noise per channel of ~180 mK) with the PdBI data (rms noise per channel of ~53 mK) to add the short-spacing information not detected by the millimetre interferometer. We used the short-spacing processing widget of the mapping package of the gildas software which produces visibilities from the single dish observations with the PdBI spectral resolution and merges them with the PdBI data. After trying different single dish weight factors, we found the best value to be 0.1. This value was the best compromise between adding excessive noise from the PdBI data due to the weak signal and giving to much weight to the single dish data not being able to see the contribution of the interferometer data.

When creating the dirty map we used a uv-taper of 65 m (20.3 kλ) to avoid the noise contribution of the longer uv-distances, recovering the highest signal-to-noise ratio. For the cleaning process we used the Steer, Dewdney and Ito (SDI) algorithm (Steer et al. 1984) to improve some strips that appeared using the Högbom algorithm due to the extended emission. The velocity resolution of the channel maps was 0.12 km s-1 as in the PdBI observations and the cutoff value to stop the cleaning was 1.5 times the rms, 0.06 Jy beam-1. The synthesized beam was 10.78 arcsec × 6.11 arcsec, with a PA of 11°.

2.4. Extinction data

We have made use of the high angular resolution extinction maps by Román-Zúñiga et al. (2009, 2010). These maps were constructed from a concerted deep near-infrared imaging survey using several telescopes (ESO-VLT, ESO-NTT and CAHA 3.5 m) as well as the 2MASS data. These maps have a resolution three times higher than the previous extinction map of this cloud by Lombardi et al. (2006), allowing to resolve structures down to ~2000 au.

3. Results

3.1. N2H+ (10)

3.1.1. VESPA IRAM 30 m

thumbnail Fig. 1

N2H+ (10) IRAM 30 m VESPA data. a) Moment-zero map (using the three central hyperfine lines) of the IRAM 30 m data. Contours are 35, 45, ..., 95 percent of the value of the peak, 1.9 K km s-1. b) Color scale: column density map. Black contours: IRAM 30 m MAMBO-II 1.2 mm dust continuum emission at an angular resolution of 21 arcsec. Contour levels are 30, 40, 50, 60, 70 and 80 percent of the peak value, 122 mJy beam-1. c) Velocity map and d) line width map resulting from the hyperfine fit. Numbers in panel c) indicate the position of the fit spectra shown in Fig. 2 and Table 2. The HPBW of 27.8 arcsec is shown at bottom left corner of panel c).

In Fig. 1a we present the moment-zero map (i.e., velocity integrated intensity map) of the N2H+ line emission observed with the IRAM 30 m telescope. The moment-zero emission was computed in the 4–8 km s-1 velocity range, which covers the three central hyperfine lines. The structure is compact with a size of ~2 arcmin (~19 000 au) and the intensity increases towards the center of the map up to 2 K km s-1.

In order to infer the physical and chemical properties in FeSt 1-457 from the N2H+ (10) IRAM 30 m data, we used the HfS tool4 (Estalella 2016) to fit the hyperfine structure of our spectra. With the HfS_fit_cube procedure we could fit all the spectra of our data cube pixel-by-pixel. The fit was done only for the spectra with a peak intensity larger than five times the rms noise. We obtained the physical parameters vLSR, Δv and τ and derived the excitation temperature Tex and the column density N following Caselli et al. (2002).

thumbnail Fig. 2

Upper pannels: HfS fit of the hyperfine structure of N2H+ (10) IRAM 30 m spectral data at different positions of the core (labeled in Fig. 1c). Black and red lines are the observational data and the fit, respectively. Lower panels: ARTIST simulated N2H+ (10) spectral data (see Sect. 4.2). Black and blue lines are the observational and simulated data, respectively. We note that the simulated core from ARTIST is spherically symmetric. The discrepancies between the observational and the simulated spectra in panel (4) may be due to the slightly elongated shape of the core towards the west.

Table 2

HfS fit parameters at different positions of FeSt 1-457.

The procedure used for fitting the hyperfine structure of N2H+ (1–0) assumes local thermodynamic equilibrium (LTE) for all the seven hyperfine lines of the transition. However, it has been found that the hyperfine lines F1, F = 1,0 → 1,1 and 1,2 → 1,2 can be out of LTE. This is the case of the quiescent core L1512, where these excitation anomalies are attributed to collisional pumping, if the collisional coefficients of the hyperfines were different, or to radiative trapping effects depending on the different optical depth of the hyperfines (Caselli et al. 1995; Daniel et al. 2006). Thus, in order to avoid this problem, we performed the fit of the hyperfine structure of N2H+ (1–0) using only the five hyperfine components without excitation anomalies. In Fig. 2 and Table 2 we show the results from the fit for six selected positions (shown in Fig. 1c) of FeSt 1-457. As can be seen in Fig. 2, although there are some slight discrepancies in the amplitudes of the two hyperfine lines out of LTE (the two with the highest velocity), the synthetic spectra (assuming LTE) fit very well the observed spectra.

In Fig. 1b we show the resulting N2H+ column density map overlapped with the 1.2 mm dust continuum emission. The values range from 3 × 1012 cm-2 at the edges to 9.0 × 1012 cm-2 towards the center, with hints of an arc-like structure near the continuum peak marginally resolved (see Sect. 3.1.3 below).

The velocity map shown in Fig. 1c presents a very smooth gradient in the northeast-southwest direction. The velocity ranges from 5.72 ± 0.01 km s-1 at the northeast to 5.90 ± 0.01 km s-1 at the southwest. Excluding the southwesternmost region, the velocity gradient along the northeast-southwest direction is 1.8 km s-1 pc-1. Aguti et al. (2007) find a similar velocity pattern in the C18O (10) and (21) channel maps.

The typical line widths (see Fig. 1d) are 0.17–0.23 (±0.03 km s-1), except in the southwesternmost region, where it raises up to 0.30 ± 0.03 km s-1 (also in agreement with the results from Aguti et al. 2007). We note that the eastern side of the core shows the lowest line widths. The thermal line width for N2H+ has a value of 0.12 km s-1, for a kinetic temperature of 10 K. Thus, the three-dimensional non thermal dispersion (σ3D,nth) ranges between 0.09 and 0.21 km s-1. The isothermal sound speed value cs is 0.19 km s-1, which gives a Mach number (defined as σ3D,nth/cs) between 0.46 to 1.1 (following Palau et al. 2015). These values indicate that the non thermal motions seem to be subsonic (as reported by Frau et al. 2010, 2012b; Alves et al. 2014) and the thermal pressure seems to be dominant over the turbulent pressure in the core (in agreement with Aguti et al. 2007) except at the southwestern corner where the ratio seems to be more balanced.

3.1.2. PdBI data

thumbnail Fig. 3

PdBI channel map of N2H+ (10). Contours are 6, 3, 3, 6, 9, 12 times the rms level of the map, 25 mJy beam-1 (~53 mK). The synthesized beam located at the bottom right corner of the third panel is 10.6 × 5.74 arcsec2, with a PA of 11.94°. We note that the rms noise increases at the edge of the images. Last panel shows the PdBI spectrum (red line) at the position of the cross. The IRAM 30 m N2H+ spectrum (black) has been divided by a factor of five to compare the hyperfine structure with the PdBI spectrum.

In Fig. 3 we show a section of the channel map showing the N2H+ (10) emission detected with the PdB interferometer. The weak emission (6σ) is better detected in the second panel corresponding to the strongest line of the three central hyperfine lines, at a velocity of 5.80 km s-1. The last panel of Fig. 3 shows the spectra (in red) at the position of the 6σ detection. The IRAM 30 m spectra (in black) is also shown as comparison, and we can see that the hyperfine structure is also seen in the PdBI spectrum indicating that indeed there is a real N2H+ detection from the PdBI. The emission filtered out by the PdBI with respect to the IRAM 30 m is estimated to be ~98%, from a comparison of the PdBI spectra (convolved with a Gaussian to obtain the IRAM 30 m beam size) and the IRAM 30 m spectra at the peak position. Given the smallest antenna baseline of PdBI observations, the largest angular scale detectable by the PdBI observations is ~12 arcsec or ~1700 au (following appendix of Palau et al. 2010). Thus, about 98% of the emission detected with the IRAM 30 m comes from structures at scales 1700 au.

3.1.3. IRAM 30 m+PdBI

thumbnail Fig. 4

Color scale: moment-zero map of the IRAM 30 m and PdBI combined data. Contours are 57, 67, 77, 88, 97, 99.5, 100 percent of the peak, 1 Jy beam-1 km s-1. The synthesized beam located at bottom right is 10.8 × 6.11 arcsec2, with a PA of 10.79°. Black contours: MAMBO-II 1.2 mm dust continuum emission at an angular resolution of 15 arcsec. Contours are 59, 79, 89, 99 percent of the peak, 57 mJy beam-1.

In Fig. 4 we show the moment-zero map of the merged IRAM 30 m and PdBI data (see Sect. 2.3), overlapped with IRAM 30 m MAMBO-II map of the dust continuum emission at 1.2 mm. The peak is 1 Jy beam-1 km s-1 and the arc structure hinted at the column-density single-dish map (panel b of Fig. 1) is now clearly seen.

3.2. Chemical survey (FTS)

thumbnail Fig. 5

Red contours: averaged intensity maps of the observed molecular transition indicated in each panel. Contours in a), j) and l) are 10, 20, 30, ... times the rms of the maps, 8, 6 and 9 mK, respectively. In panels m) and n) contours are 20, 23, 26, ... and −3, 3, 6, 9, ... times the rms, 53 mK, respectively. In panels b) and h), k) and p) countours are −3, 3, 6, 9, 12, 15, 20, ... times the rms, 21 mK and 16 mK, respectively. In the rest of the maps contours are −3, 3, 6, 9, 12, 15, 20, 25, 30, 35, ... times the rms, 10 mK. Panel n) shows a fainter component of the 13CO molecular line. The beam size of each molecular line is shown at the bottom right corner of each panel. Gray scale: IRAM 30 m MAMBO-II map of the dust continuum emission at 1.2 mm at an angular resolution of 21 arcsec.

Table 3

Chemical differentiation.

In this section, we describe the maps of 15 molecular lines detected with the IRAM 30 m FTS spectrometer towards FeSt 1-457 (see Fig. 5 and Table 1). To generate the maps we performed a Gaussian smoothing of three pixels and we used the map_aver task of greg software to average the intensity of only the channels with detected emission to obtain a better signal-to-noise ratio and to better characterize the structure of the faintest lines. As presented in previous works (Frau et al. 2010, 2012b) we observed a rich chemistry with emission from early- (e.g., C3H2, HCN, CS) and late-time molecules (e.g., N2H+). For the nitrogen-bearing molecules (e.g., HNC and N2H+; panels a to f) the gas emission is very compact and the peak coincides well with the peak of the 1.2 mm dust continuum emission. A similar behavior is seen in C3H2 (panels g and h). However, for C4H and molecules containing oxygen in our sample (e.g., HCO+, CH3OH; panels i to k) the emission is shifted ~44′′ towards the southwest with respect to the dust. An extreme case of this would be the C18O and 13CO molecules, which present strong extended emission peaking to the southwest, at ~120′′ and ~60′′ of the dust continuum, respectively (panels l and m). For the molecules containing sulfur (CS and SO, panels o and p) the emission is located to the east with respect to the dust peak, ~35′′ and ~25′′, respectively. The peak offsets of the molecular lines with respect to the 1.2 mm dust continuum peak are comparable to the beam sizes of the Gaussian smoothed integrated intensity maps (~45′′; see Table 3) and therefore significant. Panel n shows the emission of a second (and fainter) velocity component of the 13CO emission at 3 km s-1.

Regarding the intensity of the lines, apart from the 13CO and C18O which do not trace the core well, the brightest lines are HNC, HCO+ and N2H+. We found the weaker lines to be C4H, NH2D and H13CN. This is in agreement with the previous observations by Aguti et al. (2007). However, they find the morphology of N2H+ more compact compared to our data. This difference could be due to the better S/N of our observations which allowed us to detect weaker extended emission. As the FTS backend has poor spectral resolution, the specific values of the line intensities cannot be compared with other works due to the spectral dilution. This does not affect the relative positions as the lines are narrow throughout the core. We note that the VESPA and FTS N2H+ maps present very similar structure (see Figs. 1a and 5e).

4. Analysis

4.1. The Bonnor-Ebert sphere model fit

A suitable model to fit the Pipe nebula starless cores is the Bonnor-Ebert sphere (Bonnor 1956; Ebert 1955). It describes a self-graviting, pressure-confined, isothermal gas sphere in hydrostatic equilibrium. The density profile can be derived by solving the Lane-Emden equation 1ξ2ddξ(ξ2dφdξ)=eφ,\begin{equation} \frac{1}{\xi^2}\frac{{\rm d}}{{\rm d}\xi}\left(\xi^2\frac{{\rm d}\phi}{{\rm d}\xi}\right)={\rm e}^{-\phi}, \label{eq_BE} \end{equation}(1)where ξ is the non-dimensional radius, ξ=r4πGρcs,\begin{equation} \xi=r\frac{\sqrt{4 \pi G \rho}}{c_{\rm s}}, \label{eq_xi} \end{equation}(2)and φ the logarithm of the density contrast, φ=ln(ρcρ)·\begin{equation} \phi=\ln{\left(\frac{\rho_{\rm c}}{\rho}\right)}\cdot \label{eq_denCont} \end{equation}(3)r, ρ, and ρc are the radius, volume density and central volume density, respectively. The sound speed, cs, is defined as cs=kT/μ\hbox{$c_{\rm s}=\sqrt{kT/\mu}$}. G, k, T, and μ are the gravitational constant, the Boltzmann constant, the temperature, and the mean molecular mass, assumed to be 2.33. Imposing boundary conditions at the core center, forcing the density to be ρc and its derivative to be 0 (i.e., ρ = ρc and dρ/ dr = 0 at r = 0), Eq. (1)can be solved numerically. The last ingredient left is the outer confining pressure, Pout, exerted at the outer radius, Rout, from which one can derive the ξmax parameter (ξmax=(Rout/cs)4πGρ\hbox{$\xi_{\rm max}=(R_{\rm out}/c_{\rm s})\sqrt{4 \pi G \rho}$}) that uniquely characterizes the Bonnor-Ebert solution. The critical value is ξmax = 6.5, which corresponds to (ρc/ρ)max = 14.1. For values larger than ξmax the equilibrium is unstable to gravitational collapse.

thumbnail Fig. 6

Upper panel: intensity profile from the 1.2 mm dust continuum emission map (Frau et al. 2010, 2012b). Black line and dots: observed values with vertical bars depicting the ± 1-σ range. Green (long-dashed) line: isothermal Bonnor-Ebert fit. Red (solid) line: fit Bonnor-Ebert profile adding the temperature profile shown in the upper panel of Fig. 7. The best fitting central density is ρc = 2 × 105 cm-3. The fit core boundary radius, Rout = 20 000 au is marked by a vertical dashed line. Lower panel: AV profile from the dust extinction map (Román-Zúñiga et al. 2009, 2010). Black, green (long-dashed) and red (solid) lines as in upper panel. The fit core boundary radius, Rout = 20 000 au, and AV value of the surrounding medium, AVbg=9.5\hbox{$A_\mathrm{V}^{\rm bg}=9.5$} mag, are labeled and marked by a vertical dashed line and horizontal dotted line, respectively.

In order to fit the observational data of FeSt 1-457, synthetic 3D cores were generated with density profiles obeying the numerical integration of a Bonnor-Ebert sphere (Eqs. (1)–(3)) with a given set of parameters. Then, the emission was integrated along the line-of-sight to generate 2D synthetic maps following a similar recipe as in Frau et al. (2011). We took into account the presence of the molecular cloud surrounding the core by adding a term in the extinction data, Avbg\hbox{$A_{\rm v}^{\rm bg}$}, which is constant along the extinction map. Since the IRAM 30 m bolometer filters out the extended emission, we assumed that the molecular cloud does not contribute to the 1.2 continuum emission. Two maps were generated, a visual extinction map and a 1.2 mm dust continuum emission map, convolved to the respective observational beams (20′′ and 21.5′′ from Román-Zúñiga et al. 2010; Frau et al. 2010, 2012b, respectively). A dust opacity of κ250 GHz = 0.007 cm2 g-1 was adopted for dust grains with thin ice mantles and for gas density <106 cm-3 (Ossenkopf & Henning 1994). Finally, radial profiles of the observed and synthetic maps were obtained using the ellint task from the Miriad software (Sault et al. 1995). This task averages the intensity in concentric annuli. The quality of the fit was assessed by computing the χ2 residual fit of the two radial profiles. Figure 6 shows the comparison between the observational (black lines and dots) and synthetic Bonnor-Ebert emission radial profiles (green (long-dashed) lines).

Table 4

Bonnor-Ebert fit parameters of FeSt 1-457.

The best fit was obtained for the values shown in Table 4. The first row lists the fit parameters of the Bonnor-Ebert profile (ρc, Rout, and T), the physical parameters derived (ξmax, ρc/ρR, Pout, and mass), the fit background visual extinction arising from the surrounding ambient medium (AVbg\hbox{$A_\mathrm{V}^{\rm bg}$}), the gravitational stability status, and an estimate of the age (Aikawa et al. 2005). There is no degeneracy between the fit parameters ρc, Rout, and T as the visual extinction map does not depend on temperature and ρc and Rout modify the height and width of the radial profile, respectively. The temperature only modifies the 1.2 mm dust continuum emission map. We found that the temperature derived from the Bonnor-Ebert fit (10.5 ± 0.1 K) is in fair agreement with the temperature derived from ammonia observations (9.5 ± 0.1 K; Rathborne et al. 2008). The value obtained for the derived parameter ξmax was 10.5 ± 0.7. The error was estimated through error propagation analysis, being ρc and Rout the main contributors to its uncertainty (see Table 4). This value is significantly higher than the critical value (6.5), thus the core is gravitationally unstable, in accordance with Kandori et al. (2005) results.

FeSt 1-457 is observed with an excellent signal-to-noise to test the Bonnor-Ebert scenario. By using the isothermal Bonnor-Ebert sphere we found a good fit for the AV map but non-negligible differences in the dust emission profile (see green (long-dashed) line in Fig. 6). This can be interpreted as a departure from the isothermal approximation because the extinction map (AV), contrary to the dust continuum emission, is only sensitive to the column density and does not depend on temperature. Previous studies in starless cores (e.g., Crapsi et al. 2007) show that those with nH2> 105 cm-3 may have an inner temperature gradient. These cores are colder at the center, reaching temperatures as low as 6 K, than at the outer shells, where they may be heated by the interstellar radiation field (e.g., Marsh et al. 2014). We tested different temperature profiles to check whether the fits improved. We found that a temperature profile following: T(r)=Tout(ToutTin)(1+(r/rin)2),\begin{equation} T(r)= T_{\mathrm{out}}-\frac{(T_{\mathrm{out}}-T_{\mathrm{in}})}{(1+(r/r_{\mathrm{in}})^2)}, \end{equation}(4)(rin = 2500 au) with Tin = 6 K at the center of the core and increasing until Tout = 12 K at a radius of 9000 au fit much better the data. The second row of Table 4 shows the best fit parameters of the Bonnor-Ebert+temperature gradient model. Figure 7 shows the Bonnor-Ebert density and temperature profiles that appear to fit best the data. In addition, the red line in Fig. 6 shows the predicted 1.2 mm dust intensity and extinction radial profile from this Bonnor-Ebert+temperature gradient model. Forbrich et al. (2015) also find a clear decrease of the effective dust temperature towards the center of FeSt 1-457 from Herschel observations. However, they find the lowest value at the center to be 13.5 K, higher than kinetic gas temperature derived from ammonia (9.5 K; Rathborne et al. 2008) and higher than the temperature at the center of the core adopted by us. They argue this difference could be related to the effect of averaging the emission from Herschel along the line of sight which would take into account dust from the cloud at higher temperatures. We note that the temperature profile we found fits well the FeSt 1-457 core is similar to the one found for the L1544 core (Crapsi et al. 2007), a core with significantly higher density, but also to core 16 from Marsh et al. (2014), which has a density similar to that of FeSt1-457.

thumbnail Fig. 7

Temperature profile (upper pannel) and Bonnor-Ebert density profile (lower panel) used in this work (see Sect. 4.1).

Even though we did not modify the Bonnor-Ebert density profile (lower panel of Fig. 7), which assumes a constant temperature, we finally considered the temperature profile shown in the upper panel of Fig. 7 to generate the synthetic maps. This apparent inconsistency is not important in our case. Sipilä et al. (2011) show that the isothermal and non-isothermal Bonnor-Ebert density profiles are fairly similar. In addition, as the gas pressure is proportional to the temperature, reducing the temperature by half, reduces the pressure by half as well, which makes the core even more unstable against gravity. Therefore, this model is good enough to say that FeSt 1-457 has a density profile of a very young core, it is gravitationally unstable and seems to be colder at the center than in the outer parts.

4.2. ARTIST

Once the physical model of temperature and density of the core was determined (see previous section), we estimated the abundances for the molecular lines detected with the IRAM 30 m observations.

We used the Adaptable Radiative Transfer Innovations for Submillimeter Telescopes (ARTIST; Brinch & Hogerheijde 2010; Padovani et al. 2011, 2012; Jørgensen et al. 2014) software to generate a synthetic spectral cube of the molecular lines detected. We assume that the source is spherically symmetric with the density and temperature radial profiles derived from the Bonnor-Ebert model derived in the previous section.

For the simulations we assumed the dust temperature equal to the gas temperature. We assumed no velocity field in the core. The velocity gradient in the core is very small (0.05 km s-1 beam-1 at the center of the core), so this does not affect the results. Initially, the only free parameter used in ARTIST was the abundance. First, we assumed a constant value for the abundance. If the model did not fit well the observed intensity profile, we used a step function with two or three different abundances to obtain more accurate fits. We systematically tried different values for the radius at which there was a change in abundance, until obtaining the smallest χ2 parameter. One or two step functions (at a radius of 55 arcsec (~8000 au) and 100 arcsec (14 500 au)) for the abundance was necessary to obtain good fits for ten out of the 12 molecules, in which a constant abundance value did not fit well the observed data. We made use of the Leiden Atomic and Molecular Database (LAMDA)5 to load the collisional coefficients of the observed molecules (Green & Thaddeus 1974; Green & Chapman 1978; Flower 1999; Chandra & Kegel 2000; Daniel et al. 2005; Lique et al. 2006a,b; Rabli & Flower 2010; Yang et al. 2010). The spectral cube generated with ARTIST was convolved with a Gaussian to mimic the angular resolution of each molecular line observed. The output from ARTIST is given in units of brightness temperature.

In order to compare both the observed and the synthetic molecular line maps, we generated line intensity profiles for both data sets in the same way using the Miriad software. We first computed integrated intensity maps (moment-zero) using the moment task over a velocity range that included all the emission. Then, we used the ellint task for the integrated intensity maps to compute the line radial profiles with annuli of 10 arcsec width, which is a value between one third and one fifth of the angular resolution of the observed molecular lines.

Figures 8 and 9 show the annuli averaged integrated intensity (moment-zero) vs. the annuli radius for the observed (open circles) and ARTIST produced data.

thumbnail Fig. 8

N2H+ ARTIST abundance estimate. Annuli averaged integrated intensity vs. annuli radius. Black line: IRAM 30 m VESPA data. Blue dashed line: ARTIST radial integrated intensity profile for a constant abundance of 2.5 × 10-10. Red line: ARTIST radial integrated intensity profile with a step function abundance of 2.5 × 10-10 and 1.5 × 10-10 for smaller and larger radius than ~8000 au (55 arcsec; indicated by the brown vertical dotted line), respectively.

thumbnail Fig. 9

FTS molecular lines abundance fits. Annuli averaged integrated intensity vs. annuli radius. Black and red lines: FTS data profiles. Orange and violet lines: ARTIST data profiles with estimated abundances. Brown vertical dashed lines indicate the radius where there is a change in abundance (see Table 5).

4.2.1. N2H+ abundance fit

To obtain the moment-zero map of N2H+ (see Fig. 1a), we only included the three central hyperfine lines to avoid the anomalous hyperfine lines (see Sect. 3.1.1) and we obtained a good fit of the data. The high spectral resolution ARTIST fit is shown in the lower panels of Fig. 2 (we note that ARTIST handles line overlapping: Brinch & Hogerheijde 2010). We first used a constant abundance (with respect to H2) to fit the radial profile. A value of 2.5 × 10-10 fit well the data except for the outer radius where the model gave slightly higher values (blue dashed line of Fig. 8). Thus, we used a step function varying the radius at which there was a change in abundance systematically in steps of 10 arcsec, and a better fit was obtained with values of 2.5 × 10-10 and 1.5 × 10-10 for smaller and larger radius than 55 arcsec (~8000 au), respectively (red line of Fig. 8). The error of the radius where the abundance changes is ±15 arcsec, estimated from a χ2 analysis (e.g., Gibb et al. 2000). We also estimated the abundance uncertainty for the N2H+ line to be less than ~15%. This uncertainty was estimated from a χ2 analysis using different values of the abundance at the same range of radii of the core.

4.2.2. FTS molecular lines abundance fits

We followed the same procedure used for N2H+ to fit the other molecular lines detected with the FTS spectrometer, except for C4H for which the collisional coefficients file was not available in the LAMDA catalog. H13CN and C3H2 2(0, 2)1(1, 1) were not fit as they present weak emission and we had better signal to noise in the similar molecule and transition HCN and C3H2 2(1, 2)1(0, 1), respectively. The observed intensity profiles are presented in Fig. 9 (open circles). Our ARTIST modeling indicated that some molecular lines are well fit with just one abundance value. This is the case of HC3N and C3H2 2(1, 2) (see panels c and e of Fig. 9). HNC and NH2D are well fit with a one step function, one value at the inner part of the core until 55 arcsec (~8000 au) radius and a smaller value at larger radii (see panels a and d of Fig. 9). However, our ARTIST model revealed that most of the lines needed 3 different abundance values changing at 55 and 100 arcsec (~8000 and 14 500 au) radius to obtain a good fit (radii where we assumed a change in abundance are marked in Fig. 9 by vertical dashed lines).

We estimated the abundance uncertainty for the CH3OH (one of the molecules with the lowest signal-to-noise ratio) and HCN (one of the molecules with the highest signal-to-noise ratio) following the same procedure as with the N2H+ giving similar results.

Table 5

Estimated molecular abundances using ARTIST.

In all cases where we introduced a two step function profile in the abundance, values are higher in the outer part of the core and decrease towards the center. We consider a slightly different analysis for the abundances of C18O and 13CO lines. As their emission is very asymmetric presenting their peak at the southwest part with respect to the dust continuum (see panels l and m in Fig. 5), we divided the image in two sections, one including the strongest emission (southwest) and the other including the faintest emission to the northeastern part, which resembles more the emission pattern of the other species. Thus this part should trace more reliably the behavior of this molecule in the core. We performed the same procedure for the two sections separately. Both molecular lines show very flat profiles fit with a two step function for the abundance as the emission is very extended (see panels j and k of Fig. 9). Table 5 summarizes the abundances for each molecule. We note that by using the main beam brightness temperature scale the intensity of the lines with extended emission, namely C18O and specially 13CO, may be overestimated. This may result in some overestimation of their abundances but not of the observed relative drop within the core.

4.3. Estimation of the magnetic field strength

We made use of the polarization data presented in Alves et al. (2014) to estimate the magnetic field strength on the plane of the sky. The observations were performed at the APEX 12 m telescope at 345 GHz obtaining a FWHM of ~20′′. The magnetic field lines are relatively uniform with a weighted mean position angle of 130° (from north to east direction).

Using the Chandrasekhar-Fermi technique we are able to estimate the magnetic field strength. We used the Ostriker et al. (2001) approximation for estimates of the plane-of-sky field strength under strong-field cases (δφ ≤ 25°). The projected field strength (Bpos) is given by: Bpos=8.5nH2/(106cm-3)Δv/(kms-1)-1δφdeg-1mG,\begin{equation} B_\mathrm{pos}=8.5\frac{\sqrt{n_{\mathrm{H_2}}/(10^6\mathrm{~cm}^{-3})}\Delta v/(\mathrm{km~s}^{-1})^{-1}}{\delta\phi \mathrm{~deg}^{-1}} \mathrm{~mG}, \end{equation}(5)where the number density of molecular hydrogen nH2 is taken as 2.33mHnH, Δv is the observed line width and δφ is the dispersion of the polarization angles.

The observed dispersion δφobs, is produced by contributions from the measurement uncertainty of the polarization angle σPA and the intrinsic dispersion δφint. Taking into account every pair of polarization angles (PA) and their uncertainties σPA according to the equation δφint2=i,j(PAiPAj)2(σPAi2+σPAj2)N,\begin{equation} \delta\phi_{\mathrm{int}}^2=\frac{\sum\limits_{i,j}\left(\mathrm{PA}_i-\mathrm{PA}_j\right)^2-\left(\sigma_{\mathrm{PA}_i}^2+\sigma_{\mathrm{PA}_j}^2\right)}{N}, \end{equation}(6)where N is the number of data pairs, we obtained \hbox{$\delta\phi_{\mathrm{int}}=9\fdg22$}.

Frau et al. (2010) obtain nH2 = 4 × 105 cm-3, however, we adopted nH2 = 6 × 104 cm-3, as Alves et al. (2014) find depolarization at higher densities. From our N2H+ (1–0) observations, Δv ~ 0.2 km s-1. Thus, we obtained Bpos ~ 45.2 μG and assumed to be uncertain by a factor of two.

5. Discussion

5.1. Molecular depletion

After our analysis with ARTIST (see Sect. 4.2), we found that molecular abundances are well fit with three values, presenting a decrease of abundance towards the inner region of the core in most of the molecules detected (seven out of 12). The simplest explanation for the decrease of molecular abundance is that the molecules are depleted onto the dust grains once the density of the core increases and the temperature decreases. Some of these molecules are strongly depleted in the core interior with a decrease of abundances by at least one or two orders of magnitude. The strongest depletion was found in 13CO and C18O with a difference in abundance of two orders of magnitude, from ~10-6 to ~10-8 and ~10-7 to ~10-9, respectively. Aguti et al. (2007) find CO is depleted only by a factor of five but by averaging the abundance along the line of sight. CS and SO also present strong depletion with a decrease of abundance of a factor of ~30 (it is known that sulfurated molecules tend to deplete quite fast in chemical models of dense cores, e.g., Ruffle et al. 1999). Similar results are found in Tafalla et al. (2002, 2004) in several starless cores. We also found a decrease in abundance of one order of magnitude in HCN and HCO+. The HCN shows emission peaks at both sides of the continuum peak in an east-west strip. The decrease of emission around the continuum peak is probably an indication of heavily depletion at that position. Curiously, the shape of the HCN and HNC emissions are quite different close to the continuum peak. Loison et al. (2014) find that the HCN/HNC ratio tends to decrease when carbon atoms or CO are depleted, due to isomerisation of HCN into HNC. In fact, we see an increase of abundance of a factor of two of HNC at the inner region of the core. Thus, the combination of depletion of HCN plus the conversion of HCN into HNC could be enough to explain the observed difference in emission.

Our best studied molecule, N2H+, is well fit with two abundance values with a small increase towards the interior of the core (see Sect. 4.2.1). The abundance value ~10-10 is in agreement with Aguti et al. (2007) results. In Sect. 3.1.1, we showed that the N2H+ column density follows an arc-structure around the 1.2 mm continuum peak (see Fig. 1b) which is clearly seen using the combined emission of IRAM 30 m and PdBI (see Fig. 4). This might not be seen in the ARTIST analysis because we averaged the emission in concentric annuli and the feature is only present at a very small scale of 5 arcsec or 725 au, suggesting that N2H+ might be frozen-out onto the dust grains at this scale, which corresponds to a density of ~2 × 105 cm-3. Aguti et al. (2007) also suggest posible depletion of N2H+ in FeSt 1-457 at depths >40 mag of visual extinction, which corresponds to radii 3350 au and densities 1.3 × 105 cm-3. Our results are consistent with previous observations towards starless cores, where N2H+ is depleted at densities >5 × 105 cm-3 (e.g., Belloche & André 2004; Pagani et al. 2007). The arc-like structure is not likely due to opacity effects as in the calculation of the column density we already took opacity into account. Higher angular resolution and with even better sensitivity observations are needed to discard other explanations for the arc structure seen in N2H+ (such as a real physical structure).

NH2D appears to have a significant abundance only in the inner region of the core, at radius 55 arcsec (8000 au). Previous observational studies have shown evidence of NH2D abundance increase towards starless cores (e.g., Busquet et al. 2010; Masqué et al. 2013). Deuterated ammonia is expected to be abundant in the coldest and densest parts of starless cores, where the depletion of CO enhances the deuterium fractionation through a combination of gas phase (e.g., Roberts et al. 2004; Hatchell et al. 2013) and surface reactions (e.g., Fedoseev et al. 2015; Fontani et al. 2015).

Finally, HC3N and C3H2 are well fit with a constant abundance value throughout the core. Tafalla et al. (2006) study the molecular abundances of two starless cores (L1498 and L1517B) and find depletion of one order of magnitude of HC3N towards the center of L1498 and of only a factor of two in L1517B. However, they find that HC3N seems to survive in the gas phase to higher densities than other species as CO and SO. They also find a similar behavior in L1498 for C3H2 and suggest that the chemistry of C3H2 and HC3N could be related. In our case, we can only see some depletion effect in the molecular profiles of HC3N and C3H2 in the innermost position closer to the center of the core (20 arcsec, see panels c and e of Fig. 9).

In summary, for most of the molecules detected in our sample (HCN, HCO+, CH3OH, CS, SO, 13CO and C18O), we found a clear decrease in their abundance towards the inner region of the core.

5.2. Chemical differentiation in the core’s outskirts

As we have seen in Sect. 3.2, FeSt 1-457 presents a clear chemical differentiation for nitrogen, oxygen and sulfurated molecules. Their emission appears with different offsets with respect to the peak of the 1.2 mm dust continuum emission (see Sect. 3.2 and Table 3).

We computed the CS/C34S, SO/34SO and 13CO/C18O line intensity ratios from the integrated spectra of an area of 30 arcsec radius around the position where the main isotopologue line is brighter. Because these ratios are done for the same rotational transition, they depend only on the optical depth and abundance ratio of the isotopologues (we adopted an abundance ratio of 5.6 for 13CO/C18O and 24 for 32S/34S: Wilson & Matteucci 1992; Chin et al. 1995). The optical depth for CS (21), SO 3(2)2(1) and 13CO (10) transitions is 1.54, 2.11 and 2.14, respectively. We expect this value to decrease as we go away from the peak. Thus, the observed different morphologies are not likely due to opacity. In addition, the critical densities and energies of the upper level of N2H+ (10), CS (21) and HCO+ (10) are similar, ~105 cm-3 and 4–7 K, respectively (from LAMDA database; see Table 1). Therefore, we can exclude different excitation conditions as well. This suggests that different types of molecules are tracing different chemical conditions in the core.

The characteristic asymmetries in the emission of C4H, HCO+ and CH3OH, peaking (south)west of the dust peak, could be evidence of the penetration of UV radiation into FeSt 1-457.

Moreover, H13CN and HC3N also show a slight displacement of their emission peak towards the southwest. This may be a similar situation to the emission enhancement of these molecules in PDRs (photo-dominated regions) and in molecular condensations ahead of Herbig-Haro objects (e.g., Hollenbach & Tielens 1997; Girart et al. 1998, 2002; Viti et al. 2006). In the second case, chemical models show that this enhancement is due to the UV radiation arising from the Herbig-Haro objects (Viti & Williams 1999; Viti et al. 2003). Viti et al. (2003) find that the UV radiation affects the chemistry enhancing HCO+ and CH3OH for AV< 4 mag, meaning that the enhanced emission should trace the most outer layer of the core. Additionally, PDR models show that C4H reach higher abundances at somewhat less deep layers (lower AV values) than HCN. Indeed, C4H has been found to peak close to the edge of PDRs as its production is enhanced by the UV radiation field in a relatively narrow range of visual extinction values (e.g., Pety et al. 2005; Morata & Herbst 2008; Rimmer et al. 2012). Moreover, in some PDR models (e.g., Morata & Herbst 2008; Rimmer et al. 2012), HC3N also has an abundance peak at similar layers to HCO+.

FeSt 1-457 is, in projection, embedded in a diffuse molecular cloud with average AVbg~9.5\hbox{$A_{\rm V}^{\rm bg} \sim9.5$} mag. Thus, the enhanced C4H, HCO+ and CH3OH emission could only happen if the core is near the border of the molecular cloud, such that one side of the core is exposed to the interstellar UV radiation without much absorption. We do not expect the penetration of the UV photons to be too large, but it is probably strong enough to produce some differences in the distribution of the molecular emission. Infrared extinction maps show that FeSt 1-457 is surrounded by gas with relatively low extinction, AV = 5–10 mag (Román-Zúñiga et al. 2010). If the medium around FeSt 1-457 is not homogeneous, some UV radiation from a nearby OB star or the interstellar radiation field could be penetrating quite deep into the (western side of the) core. Indeed, Gritschneder & Lin (2012) proposed that the Pipe nebula might be an H II region shell swept by the nearby θ Ophiuchi (a B2 IV β Cephei star). θ Ophiuchi lies in a direction approximately west-northwest of FeSt 1-457 and might be the source of the UV radiation required to explain the morphology of the emission of C4H, HCO+ and CH3OH. Future works must be done to confirm this possible scenario.

thumbnail Fig. 10

Polarization maps. Red contours: averaged intensity maps of a) C4H, b) HCO+, c) CS and d) SO. Contours in panels a) and c) are −3, 3, 6, 9, 12, 15, 20, 25, 30, ... times the rms of the maps, 10 mK. Contours in panel b) and d) are 10, 20, 30, ... and −3, 3, 6, 9, 12, 15 times the rms of the maps, 6, and 16 mK. The beam size of each molecular line is shown at the bottom right corner of each panel. Gray scale: IRAM 30 m MAMBO-II map of the dust continuum emission at 1.2 mm. Blue segments: position angle of the linearly polarized data from Alves et al. (2014). The length of the segments is proportional to the polarization fraction. Black arrows in panels a) and c) indicate the 1-D slice shown in Fig. 11.

Regarding the asymmetries of CS and SO emission (see panels o and p in Fig. 5, respectively), it has been studied that in PDRs (see Morata & Herbst 2008) CS has a relatively shallow abundance closer to the edge of the PDR (around or inwards the C4H peak) and a higher abundance in inner layers, while SO has high abundances in even deeper layers into the PDR. The combination of depletion and the influence of the UV field could explain the shape of the emission of these maps. If the eastern side of the core is more shielded from the radiation, these molecules will reach higher abundances and have their emission peaks there, while they will have lower abundances on the western side of the core (more markedly for SO), and they will be depleted towards the center of the core. Finally, N2H+ and deuterated molecules, from both the observational and chemical point of view, are supposed to be abundant in the inner and most shielded parts of the cloud (e.g., Roberts et al. 2003; Crapsi et al. 2005), which is what our emission maps show (see Figs. 5 and 9).

5.3. A correlation between depletion, grain growth and polarization

As we have seen in Sect. 4.2, and discussed in Sect. 5.1, molecular abundances are well fit with two or three values, presenting a decrease of abundance towards the inner region of the core in most of the molecules detected. This change in the chemical properties suggests FeSt 1-457 experiences a variation in the physical properties of its inner region.

It is interesting to note that Forbrich et al. (2015) find evidence of grain growth at the inner region of the core from Herschel-derived dust opacity results and near-infrared extinction measurements. Depletion of the molecular gas onto the dust grains due to a decrease of temperature and increase of density in the center of the core could make easier for the grains to stick together and contribute to the grain growth seen by Forbrich et al. (2015). This is supported by the fact that we found a depletion radius very similar to the radius where Forbrich et al. find grain growth. It has been studied that water ice on the surface of dust particles enhance grain-grain adhesion favoring grain growth (e.g., Gundlach & Blum 2015; Wang et al. 2005). Furthermore, Liu et al. (2013) find that CO depletion positively correlates with a decrease of the dust emissivity index, which could be related to the size of the dust grains (e.g., Lommen et al. 2007; Ricci et al. 2012; Piétu et al. 2014). This also agrees with our results of low abundance of CO at the center of the core (see Fig. 5) and the grain growth seen by Forbrich et al. (2015).

In a different work, Alves et al. (2014) present a study of the polarization properties in FeSt 1-457 at different wavelengths and find that the emission is more polarized to the western side of the core (see Alves et al. 2014, Fig. 1). It is globally accepted that non-spherical dust grains are aligned perpendicularly with respect to magnetic field lines producing linearly polarized thermal continuum emission (Davis & Greenstein 1951; Hildebrand 1988). One of the most important theories for the dust grain alignment is via radiative torques (Dolginov & Mitrofanov 1976; Draine & Weingartner 1996, 1997; Lazarian & Hoang 2007). Thus, the stronger polarized emission to the west also supports the idea of a UV radiation field penetrating into the core from the west which would be affecting not only the polarization properties but also the chemistry in the core (see previous section).

thumbnail Fig. 11

Intensity from CS (red), C4H (black) and polarization intensity (blue) along a 1-D slice from the southwest part of the core towards the east passing through the C4H and CS peaks (see panels a) and c) of Fig. 10). Both C4H and the polarization intensity increase at the same positions in the core along the slice. CS shows an anti-correlation with the polarization intensity.

In Fig. 10 we show the emission of four molecules peaking to the southwest (C4H and HCO+) and to the east (CS and SO) of the dust peak, overlaid with the position angles of the polarization data from Alves et al. (2014). The C4H and HCO+ emission seem to correlate with the stronger polarized emission towards the west. On the other hand, the emission peaks of CS and SO seem to present an anti-correlation with the polarized emission. To study further these possible correlations, we performed a 1D slice across the core passing through the C4H and CS peaks (black arrows in Fig. 10). In Fig. 11 we show the polarization intensity (blue) and the emission of C4H (black) and CS (red) along the slice. The polarization intensity presents a similar distribution and seem to correlate well with the C4H emission and anti-correlate with the CS emission. To quantify these correlations we performed the statistical tests Kendall’s τ and Spearman’s ρ6 and found a good correlation (τ = 0.647, ρ = 0.804) between C4H emission and polarization intensity and anti-correlation values (τ = −0.062, ρ = −0.058) between CS and polarization intensity. More observations are necessary to further study these results.

On the other hand, Alves et al. (2014) also find a steep decrease in the polarization efficiency at visual extinctions (AV) larger than 30 mag. Alves et al. (2014) propose the lack of an internal source of radiation results in a loss of grain alignment with the magnetic field, which is consistent with the theory of dust grain alignment via radiative torques. Our work suggests that depolarization could be produced as a combination of lack of internal source and dust grain growth possibly affected by molecular depletion. Due to the “grain clustering” yielding to larger grains, the elongated shape of the grains could be distorted. The paramagnetic nuclei of the grains would be surrounded by a thick layer of ice (unpolarized) and the alignment of the grains could be harder to achieve resulting in a depolarized inner region as seen by Alves et al. (2014).

Therefore, we conclude that depletion, grain growth and depolarization seem to be correlated in FeSt 1-457.

5.4. Magnetic properties of FeSt 1-457

Previous studies in star forming cores have shown that magnetic fields can play an important role in their evolution (e.g., Girart et al. 2006).

Magnetic forces could prevent a cloud from collapsing if its mass-to-flux ratio is below the critical value, M/ΦB<(M/ΦB)cr,\begin{equation} M/\Phi_B<(M/\Phi_B)_{\mathrm{cr}}, \end{equation}(7)where M is its mass and ΦB = πBR2 is the magnetic flux through it, derived from the mean magnetic field and the radius of the cloud cross-section perpendicular to the magnetic field. The critical mass-to-flux ratio is (M/ ΦB)cr = 1/(2πG1/2) (Nakano & Nakamura 1978), where G is the gravitational constant. We used the ratio in terms of the column density following Eq. (7) of Pillai et al. (2015)(M/ΦB)(M/ΦB)cr=0.76(NH21023cm-2)(Btot1000μG)-1·\begin{equation} \frac{(M/\Phi_B)}{(M/\Phi_B)_{\mathrm{cr}}}=0.76 \left(\frac{N_{\rm H_2}}{10^{23}~\mathrm{cm}^{-2}}\right)\left(\frac{B_{\mathrm{tot}}}{1000~\mu G}\right)^{-1}\cdot \end{equation}(8)Frau et al. (2010) obtain a core+cloud column density NH2 = 47.6 × 1021 cm-2. From our Bonnor-Ebert model fit (see Sect. 4.1) we found the background visual extinction arising from the surrounding medium (AVbg\hbox{$A_\mathrm{V}^{\rm bg}$}) to be 9.5 mag which corresponds to a column density of 8.9 × 1021 cm-2. Thus, we adopted a column density for the core of 38.7 × 1021 cm-2. On the other hand, we assumed Bpos ~ Btot ~ 45.2 μG as Alves et al. (2014) find the magnetic field direction mainly onto the plane-of-sky. We obtained a ratio of 6.5 (magnetically supercritical) which implies that gravity dominates over the magnetic support. Our Bonnor-Ebert model fit showed that the core is gravitationally unstable (see Sect. 4.1). Finally, we can see that the magnetic field is not enough to avoid the collapse.

6. Conclusions

We have used the IRAM 30 m telescope and the PdBI to study the chemical and physical properties of the starless core FeSt 1-457 (Core 109) in the Pipe nebula. Our main conclusions are as follows:

  • 1.

    We fit the hyperfine structure of the N2H+ (10) IRAM 30 m data. This allows us to measure with high precision the velocity field and line widths in the core. The N2H+ emission shows a clear southwest-northeast velocity gradient of 1.78 km s-1 pc-1, and larger line widths (0.3 km s-1) at the southwest part of the map. The typical line widths are 0.170.23 km s-1. The column density map presents an arc-structure around the 1.2 mm dust continuum peak.

  • 2.

    Combining both IRAM 30 m and PdBI N2H+ (10) data we can resolve the arc-like structure hinted at the single dish map.

  • 3.

    The core presents a rich chemistry with emission from early (C3H2, HCN, CS) and late-time molecules (e.g., N2H+), with a clear chemical spatial differentiation for nitrogen (centrally peaked), oxygen (peaking to the southwest) and sulfurated molecules (peaking to the east).

  • 4.

    The chemical difference in the core could be due to an external UV radiation field penetrating into the core from the (south)west which could be also affecting the polarization properties. This implies that the core should be close to the edge of the molecular cloud.

  • 5.

    FeSt 1-457 is well fit with a Bonnor-Ebert sphere model and introducing a temperature gradient decreasing towards the center from 12 to 6 K. We found the core is gravitationally unstable and that the magnetic field is not enough to stop the collapse.

  • 6.

    We have analyzed the abundances of the molecular lines using the ARTIST software. For most of the molecules detected (HCN, HCO+, CH3OH, CS, SO, 13CO and C18O), abundances are best fit with three values, presenting a clear decrease of abundance of at least one or two orders of magnitude towards the inner region of the core. This is consistent with chemical models of starless cores, that show that there is a significant depletion of molecules onto the dust grains. On the other hand, N2H+, HC3N and C3H2 are well fit with a constant abundance throughout the core and the abundance of HNC and NH2D increases in the inner region of the core (at radius 55 arcsec).

Finally, we have seen that depletion of molecules onto the dust grains, grain growth and depolarization take place at the inner region of the core. This strongly suggests that these properties could be correlated in FeSt 1-457.


1

Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.

2

IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain).

Acknowledgments

We would like to thank the referee for her/his useful comments. We also thank Phil Myers and Jan Forbrich for thoughtful discussions. C.J., J.M.G. and R.E. acknowledge support from MICINN (Spain) AYA2014-57369-C3 grant. C.J. and R.E. also acknowledge MDM-2014-0369 of ICCUB (Unidad de Excelencia “María de Maeztu”). J.M.G. also acknowledges the support from the MECD (Spain) PRX15/00435 travel grant and from the SI CGPS award “Magnetic Fields and Massive Star Formation”. A.P. acknowledges financial support from UNAM-DGAPA-PAPIIT IA102815 grant, Mexico. O.M. is supported by the MOST (Ministry of Science and Technology, Taiwan) ALMA-T grant MOST 103-2119-M-001-010-MY to the Institute of Astronomy & Astrophysics, Academia Sinica.

References

  1. Aguti, E. D., Lada, C. J., Bergin, E. A., Alves, J. F., & Birkinshaw, M. 2007, ApJ, 665, 457 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aikawa, Y., Herbst, E., Roberts, H., & Caselli, P. 2005, ApJ, 620, 330 [NASA ADS] [CrossRef] [Google Scholar]
  3. Alves, F. O., & Franco, G. A. P. 2007, A&A, 470, 597 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Alves, F. O., Franco, G. A. P., & Girart, J. M. 2008, A&A, 486, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Alves, F. O., Frau, P., Girart, J. M., et al. 2014, A&A, 569, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Ascenso, J., Lada, C. J., Alves, J., Román-Zúñiga, C. G., & Lombardi, M. 2013, A&A, 549, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Belloche, A., & André, P. 2004, A&A, 419, L35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bergin, E. A., & Tafalla, M. 2007, ARA&A, 45, 339 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bonnor, W. B. 1956, MNRAS, 116, 351 [CrossRef] [Google Scholar]
  10. Brinch, C., & Hogerheijde, M. R. 2010, A&A, 523, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Busquet, G., Palau, A., Estalella, R., et al. 2010, A&A, 517, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Caselli, P., Myers, P. C., & Thaddeus, P. 1995, ApJ, 455, L77 [NASA ADS] [CrossRef] [Google Scholar]
  13. Caselli, P., Walmsley, C. M., Zucconi, A., et al. 2002, ApJ, 565, 344 [NASA ADS] [CrossRef] [Google Scholar]
  14. Chandra, S., & Kegel, W. H. 2000, A&AS, 142, 113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Chin, Y.-N., Henkel, C., Whiteoak, J. B., Langer, N., & Churchwell, E. B. 1995, VizieR Online Data Catalog: J/A+A/305/960 [Google Scholar]
  16. Crapsi, A., Caselli, P., Walmsley, C. M., et al. 2005, ApJ, 619, 379 [Google Scholar]
  17. Crapsi, A., Caselli, P., Walmsley, M. C., & Tafalla, M. 2007, A&A, 470, 221 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Daniel, F., Dubernet, M.-L., Meuwly, M., Cernicharo, J., & Pagani, L. 2005, MNRAS, 363, 1083 [NASA ADS] [CrossRef] [Google Scholar]
  19. Daniel, F., Cernicharo, J., & Dubernet, M.-L. 2006, ApJ, 648, 461 [NASA ADS] [CrossRef] [Google Scholar]
  20. Davis, Jr., L., & Greenstein, J. L. 1951, ApJ, 114, 206 [Google Scholar]
  21. Dolginov, A. Z., & Mitrofanov, I. G. 1976, Ap&SS, 43, 291 [NASA ADS] [CrossRef] [Google Scholar]
  22. Draine, B. T., & Weingartner, J. C. 1996, ApJ, 470, 551 [NASA ADS] [CrossRef] [Google Scholar]
  23. Draine, B. T., & Weingartner, J. C. 1997, ApJ, 480, 633 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dunham, M. M., Crapsi, A., Evans, II, N. J., et al. 2008, ApJS, 179, 249 [NASA ADS] [CrossRef] [Google Scholar]
  25. Ebert, R. 1955, ZAp, 37, 217 [Google Scholar]
  26. Estalella, R. 2016, PASP, in press [arXiv:1608.04088] [Google Scholar]
  27. Evans, II, N. J., Rawlings, J. M. C., Shirley, Y. L., & Mundy, L. G. 2001, ApJ, 557, 193 [NASA ADS] [CrossRef] [Google Scholar]
  28. Fedoseev, G., Ioppolo, S., & Linnartz, H. 2015, MNRAS, 446, 449 [NASA ADS] [CrossRef] [Google Scholar]
  29. Feitzinger, J. V., & Stuewe, J. A. 1984, A&AS, 58, 365 [NASA ADS] [Google Scholar]
  30. Flower, D. R. 1999, MNRAS, 305, 651 [NASA ADS] [CrossRef] [Google Scholar]
  31. Fontani, F., Busquet, G., Palau, A., et al. 2015, A&A, 575, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Forbrich, J., Lada, C. J., Muench, A. A., Alves, J., & Lombardi, M. 2009, ApJ, 704, 292 [NASA ADS] [CrossRef] [Google Scholar]
  33. Forbrich, J., Lada, C. J., Lombardi, M., Román-Zúñiga, C., & Alves, J. 2015, A&A, 580, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Franco, G. A. P., Alves, F. O., & Girart, J. M. 2010, ApJ, 723, 146 [Google Scholar]
  35. Frau, P., Girart, J. M., Beltrán, M. T., et al. 2010, ApJ, 723, 1665 [NASA ADS] [CrossRef] [Google Scholar]
  36. Frau, P., Galli, D., & Girart, J. M. 2011, A&A, 535, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Frau, P., Girart, J. M., & Beltrán, M. T. 2012a, A&A, 537, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Frau, P., Girart, J. M., Beltrán, M. T., et al. 2012b, ApJ, 759, 3 [NASA ADS] [CrossRef] [Google Scholar]
  39. Frau, P., Girart, J. M., Alves, F. O., et al. 2015, A&A, 574, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Gibb, E., Nummelin, A., Irvine, W. M., Whittet, D. C. B., & Bergman, P. 2000, ApJ, 545, 309 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  41. Girart, J. M., Estalella, R., & Ho, P. T. P. 1998, ApJ, 495, L59 [NASA ADS] [CrossRef] [Google Scholar]
  42. Girart, J. M., Viti, S., Williams, D. A., Estalella, R., & Ho, P. T. P. 2002, A&A, 388, 1004 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Girart, J. M., Rao, R., & Marrone, D. P. 2006, Science, 313, 812 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  44. Green, S., & Chapman, S. 1978, ApJS, 37, 169 [NASA ADS] [CrossRef] [Google Scholar]
  45. Green, S., & Thaddeus, P. 1974, ApJ, 191, 653 [NASA ADS] [CrossRef] [Google Scholar]
  46. Gritschneder, M., & Lin, D. N. C. 2012, ApJ, 754, L13 [NASA ADS] [CrossRef] [Google Scholar]
  47. Gundlach, B., & Blum, J. 2015, ApJ, 798, 34 [NASA ADS] [CrossRef] [Google Scholar]
  48. Hardegree-Ullman, E., Harju, J., Juvela, M., et al. 2013, ApJ, 763, 45 [NASA ADS] [CrossRef] [Google Scholar]
  49. Hatchell, J., Wilson, T., Drabek, E., et al. 2013, MNRAS, 429, L10 [NASA ADS] [CrossRef] [Google Scholar]
  50. Hildebrand, R. H. 1988, QJRAS, 29, 327 [NASA ADS] [Google Scholar]
  51. Hollenbach, D. J., & Tielens, A. G. G. M. 1997, ARA&A, 35, 179 [NASA ADS] [CrossRef] [Google Scholar]
  52. Jørgensen, J., Brinch, C., Girart, J. M., et al. 2014, ARTIST: Adaptable Radiative Transfer Innovations for Submillimeter Telescopes, Astrophysics Source Code Library [record ascl:1402.014] [Google Scholar]
  53. Kandori, R., Nakajima, Y., Tamura, M., et al. 2005, AJ, 130, 2166 [NASA ADS] [CrossRef] [Google Scholar]
  54. Launhardt, R., Stutz, A. M., Schmiedeke, A., et al. 2013, A&A, 551, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Lazarian, A., & Hoang, T. 2007, MNRAS, 378, 910 [Google Scholar]
  56. Lique, F., Dubernet, M.-L., Spielfiedel, A., & Feautrier, N. 2006a, A&A, 450, 399 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Lique, F., Spielfiedel, A., & Cernicharo, J. 2006b, A&A, 451, 1125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Liu, T., Wu, Y., & Zhang, H. 2013, ApJ, 775, L2 [NASA ADS] [CrossRef] [Google Scholar]
  59. Loison, J.-C., Wakelam, V., & Hickson, K. M. 2014, MNRAS, 443, 398 [Google Scholar]
  60. Lombardi, M., Alves, J., & Lada, C. J. 2006, A&A, 454, 781 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Lommen, D., Wright, C. M., Maddison, S. T., et al. 2007, A&A, 462, 211 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Marsh, K. A., Griffin, M. J., Palmeirim, P., et al. 2014, MNRAS, 439, 3683 [NASA ADS] [CrossRef] [Google Scholar]
  63. Masqué, J. M., Girart, J. M., Anglada, G., et al. 2013, ApJ, 776, 28 [NASA ADS] [CrossRef] [Google Scholar]
  64. Morata, O., & Herbst, E. 2008, MNRAS, 390, 1549 [NASA ADS] [Google Scholar]
  65. Nakano, T., & Nakamura, T. 1978, PASJ, 30, 671 [NASA ADS] [Google Scholar]
  66. Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943 [NASA ADS] [Google Scholar]
  67. Ostriker, E. C., Stone, J. M., & Gammie, C. F. 2001, ApJ, 546, 980 [NASA ADS] [CrossRef] [Google Scholar]
  68. Padovani, M., Jørgensen, J. K., Bertoldi, F., et al. 2011, in Computational Star Formation, eds. J. Alves, B. G. Elmegreen, J. M. Girart, & V. Trimble, IAU Symp., 270, 451 [Google Scholar]
  69. Padovani, M., Brinch, C., Girart, J. M., et al. 2012, A&A, 543, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Pagani, L., Bacmann, A., Cabrit, S., & Vastel, C. 2007, A&A, 467, 179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Palau, A., Sánchez-Monge, Á., Busquet, G., et al. 2010, A&A, 510, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Palau, A., Ballesteros-Paredes, J., Vázquez-Semadeni, E., et al. 2015, MNRAS, 453, 3785 [NASA ADS] [CrossRef] [Google Scholar]
  73. Pety, J., Teyssier, D., Fossé, D., et al. 2005, A&A, 435, 885 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Piétu, V., Guilloteau, S., Di Folco, E., Dutrey, A., & Boehler, Y. 2014, A&A, 564, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Pillai, T., Kauffmann, J., Tan, J. C., et al. 2015, ApJ, 799, 74 [NASA ADS] [CrossRef] [Google Scholar]
  76. Rabli, D., & Flower, D. R. 2010, MNRAS, 406, 95 [NASA ADS] [CrossRef] [Google Scholar]
  77. Rathborne, J. M., Lada, C. J., Muench, A. A., Alves, J. F., & Lombardi, M. 2008, ApJS, 174, 396 [NASA ADS] [CrossRef] [Google Scholar]
  78. Reid, M. J., Schneps, M. H., Moran, J. M., et al. 1988, ApJ, 330, 809 [NASA ADS] [CrossRef] [Google Scholar]
  79. Ricci, L., Trotta, F., Testi, L., et al. 2012, A&A, 540, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Rimmer, P. B., Herbst, E., Morata, O., & Roueff, E. 2012, A&A, 537, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Roberts, H., Herbst, E., & Millar, T. J. 2003, ApJ, 591, L41 [NASA ADS] [CrossRef] [Google Scholar]
  82. Roberts, H., Herbst, E., & Millar, T. J. 2004, A&A, 424, 905 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Román-Zúñiga, C. G., Lada, C. J., & Alves, J. F. 2009, ApJ, 704, 183 [NASA ADS] [CrossRef] [Google Scholar]
  84. Román-Zúñiga, C. G., Alves, J. F., Lada, C. J., & Lombardi, M. 2010, ApJ, 725, 2232 [NASA ADS] [CrossRef] [Google Scholar]
  85. Roy, A., André, P., Palmeirim, P., et al. 2014, A&A, 562, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Ruffle, D. P., Hartquist, T. W., Caselli, P., & Williams, D. A. 1999, MNRAS, 306, 691 [NASA ADS] [CrossRef] [Google Scholar]
  87. Ruoskanen, J., Harju, J., Juvela, M., et al. 2011, A&A, 534, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Data Analysis Software and Systems IV, eds. R. A. Shaw, H. E. Payne, & J. J. E. Hayes, ASP Conf. Ser., 77, 433 [Google Scholar]
  89. Sipilä, O., Harju, J., & Juvela, M. 2011, A&A, 535, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Steer, D. G., Dewdney, P. E., & Ito, M. R. 1984, A&A, 137, 159 [NASA ADS] [Google Scholar]
  91. Tafalla, M., Myers, P. C., Caselli, P., Walmsley, C. M., & Comito, C. 2002, ApJ, 569, 815 [NASA ADS] [CrossRef] [Google Scholar]
  92. Tafalla, M., Myers, P. C., Caselli, P., & Walmsley, C. M. 2004, A&A, 416, 191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Tafalla, M., Santiago-García, J., Myers, P. C., et al. 2006, A&A, 455, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Viti, S., & Williams, D. A. 1999, MNRAS, 310, 517 [NASA ADS] [CrossRef] [Google Scholar]
  95. Viti, S., Girart, J. M., Garrod, R., Williams, D. A., & Estalella, R. 2003, A&A, 399, 187 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Viti, S., Girart, J. M., & Hatchell, J. 2006, A&A, 449, 1089 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  97. Wang, H., Bell, R. C., Iedema, M. J., Tsekouras, A. A., & Cowin, J. P. 2005, ApJ, 620, 1027 [NASA ADS] [CrossRef] [Google Scholar]
  98. Wilcock, L. A., Ward-Thompson, D., Kirk, J. M., et al. 2012, MNRAS, 424, 716 [NASA ADS] [CrossRef] [Google Scholar]
  99. Wilson, T. L., & Matteucci, F. 1992, A&ARv, 4, 1 [NASA ADS] [CrossRef] [Google Scholar]
  100. Yang, B., Stancil, P. C., Balakrishnan, N., & Forrey, R. C. 2010, ApJ, 718, 1062 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

FTS molecular transitions.

Table 2

HfS fit parameters at different positions of FeSt 1-457.

Table 3

Chemical differentiation.

Table 4

Bonnor-Ebert fit parameters of FeSt 1-457.

Table 5

Estimated molecular abundances using ARTIST.

All Figures

thumbnail Fig. 1

N2H+ (10) IRAM 30 m VESPA data. a) Moment-zero map (using the three central hyperfine lines) of the IRAM 30 m data. Contours are 35, 45, ..., 95 percent of the value of the peak, 1.9 K km s-1. b) Color scale: column density map. Black contours: IRAM 30 m MAMBO-II 1.2 mm dust continuum emission at an angular resolution of 21 arcsec. Contour levels are 30, 40, 50, 60, 70 and 80 percent of the peak value, 122 mJy beam-1. c) Velocity map and d) line width map resulting from the hyperfine fit. Numbers in panel c) indicate the position of the fit spectra shown in Fig. 2 and Table 2. The HPBW of 27.8 arcsec is shown at bottom left corner of panel c).

In the text
thumbnail Fig. 2

Upper pannels: HfS fit of the hyperfine structure of N2H+ (10) IRAM 30 m spectral data at different positions of the core (labeled in Fig. 1c). Black and red lines are the observational data and the fit, respectively. Lower panels: ARTIST simulated N2H+ (10) spectral data (see Sect. 4.2). Black and blue lines are the observational and simulated data, respectively. We note that the simulated core from ARTIST is spherically symmetric. The discrepancies between the observational and the simulated spectra in panel (4) may be due to the slightly elongated shape of the core towards the west.

In the text
thumbnail Fig. 3

PdBI channel map of N2H+ (10). Contours are 6, 3, 3, 6, 9, 12 times the rms level of the map, 25 mJy beam-1 (~53 mK). The synthesized beam located at the bottom right corner of the third panel is 10.6 × 5.74 arcsec2, with a PA of 11.94°. We note that the rms noise increases at the edge of the images. Last panel shows the PdBI spectrum (red line) at the position of the cross. The IRAM 30 m N2H+ spectrum (black) has been divided by a factor of five to compare the hyperfine structure with the PdBI spectrum.

In the text
thumbnail Fig. 4

Color scale: moment-zero map of the IRAM 30 m and PdBI combined data. Contours are 57, 67, 77, 88, 97, 99.5, 100 percent of the peak, 1 Jy beam-1 km s-1. The synthesized beam located at bottom right is 10.8 × 6.11 arcsec2, with a PA of 10.79°. Black contours: MAMBO-II 1.2 mm dust continuum emission at an angular resolution of 15 arcsec. Contours are 59, 79, 89, 99 percent of the peak, 57 mJy beam-1.

In the text
thumbnail Fig. 5

Red contours: averaged intensity maps of the observed molecular transition indicated in each panel. Contours in a), j) and l) are 10, 20, 30, ... times the rms of the maps, 8, 6 and 9 mK, respectively. In panels m) and n) contours are 20, 23, 26, ... and −3, 3, 6, 9, ... times the rms, 53 mK, respectively. In panels b) and h), k) and p) countours are −3, 3, 6, 9, 12, 15, 20, ... times the rms, 21 mK and 16 mK, respectively. In the rest of the maps contours are −3, 3, 6, 9, 12, 15, 20, 25, 30, 35, ... times the rms, 10 mK. Panel n) shows a fainter component of the 13CO molecular line. The beam size of each molecular line is shown at the bottom right corner of each panel. Gray scale: IRAM 30 m MAMBO-II map of the dust continuum emission at 1.2 mm at an angular resolution of 21 arcsec.

In the text
thumbnail Fig. 6

Upper panel: intensity profile from the 1.2 mm dust continuum emission map (Frau et al. 2010, 2012b). Black line and dots: observed values with vertical bars depicting the ± 1-σ range. Green (long-dashed) line: isothermal Bonnor-Ebert fit. Red (solid) line: fit Bonnor-Ebert profile adding the temperature profile shown in the upper panel of Fig. 7. The best fitting central density is ρc = 2 × 105 cm-3. The fit core boundary radius, Rout = 20 000 au is marked by a vertical dashed line. Lower panel: AV profile from the dust extinction map (Román-Zúñiga et al. 2009, 2010). Black, green (long-dashed) and red (solid) lines as in upper panel. The fit core boundary radius, Rout = 20 000 au, and AV value of the surrounding medium, AVbg=9.5\hbox{$A_\mathrm{V}^{\rm bg}=9.5$} mag, are labeled and marked by a vertical dashed line and horizontal dotted line, respectively.

In the text
thumbnail Fig. 7

Temperature profile (upper pannel) and Bonnor-Ebert density profile (lower panel) used in this work (see Sect. 4.1).

In the text
thumbnail Fig. 8

N2H+ ARTIST abundance estimate. Annuli averaged integrated intensity vs. annuli radius. Black line: IRAM 30 m VESPA data. Blue dashed line: ARTIST radial integrated intensity profile for a constant abundance of 2.5 × 10-10. Red line: ARTIST radial integrated intensity profile with a step function abundance of 2.5 × 10-10 and 1.5 × 10-10 for smaller and larger radius than ~8000 au (55 arcsec; indicated by the brown vertical dotted line), respectively.

In the text
thumbnail Fig. 9

FTS molecular lines abundance fits. Annuli averaged integrated intensity vs. annuli radius. Black and red lines: FTS data profiles. Orange and violet lines: ARTIST data profiles with estimated abundances. Brown vertical dashed lines indicate the radius where there is a change in abundance (see Table 5).

In the text
thumbnail Fig. 10

Polarization maps. Red contours: averaged intensity maps of a) C4H, b) HCO+, c) CS and d) SO. Contours in panels a) and c) are −3, 3, 6, 9, 12, 15, 20, 25, 30, ... times the rms of the maps, 10 mK. Contours in panel b) and d) are 10, 20, 30, ... and −3, 3, 6, 9, 12, 15 times the rms of the maps, 6, and 16 mK. The beam size of each molecular line is shown at the bottom right corner of each panel. Gray scale: IRAM 30 m MAMBO-II map of the dust continuum emission at 1.2 mm. Blue segments: position angle of the linearly polarized data from Alves et al. (2014). The length of the segments is proportional to the polarization fraction. Black arrows in panels a) and c) indicate the 1-D slice shown in Fig. 11.

In the text
thumbnail Fig. 11

Intensity from CS (red), C4H (black) and polarization intensity (blue) along a 1-D slice from the southwest part of the core towards the east passing through the C4H and CS peaks (see panels a) and c) of Fig. 10). Both C4H and the polarization intensity increase at the same positions in the core along the slice. CS shows an anti-correlation with the polarization intensity.

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.