© ESO 2018

1 Introduction

Filamentary structures play an important role in the process of star formation. Observations at different wavelengths based on various tracers have revealed that filaments are ubiquitous in the interstellar medium (e.g., Schneider & Elmegreen 1979; Molinari et al. 2010; André et al. 2010; Schisano et al. 2014; Ragan et al. 2014; Li et al. 2016). Filaments are seen in quiescent and star-forming clouds, in which a significant fraction of prestellar cores are located (André et al. 2010). Filamentary structures have wide ranges of masses (~1–10⁵ M_⊙) and lengths (~0.1–100 pc; e.g., Bally et al. 1987; Jackson et al. 2010; Arzoumanian et al. 2011; Hernandez et al. 2012; Hacar et al. 2013; Kirk et al. 2013; Palmeirim et al. 2013; Kainulainen et al. 2013, 2017; Beuther et al. 2015; Li et al. 2016; Abreu-Vicente et al. 2016; Zucker et al. 2018).

The processes of filament formation and filament fragmentation to star-forming cores are not well understood. Because of the wide range of filament size scales and masses these processes might also differ among filaments. High-resolution magnetohydrodynamical simulations of molecular cloud evolution and filament formation show subsonic motions in the inner dense regions of filaments, but the surrounding low density gas is supersonic (Padoan et al. 2001; Federrath 2016). Additionally, accretion flows along and radially onto the filament have been seen in observations and simulations (Schneider et al. 2010; Peretto et al. 2013, 2014; Henshaw et al. 2014; Smith et al. 2015). Therefore, the formation and evolution of filaments is a highly dynamical process and to constrain it is essential to study their kinematics.

Studies of filaments have targeted mainly sources in nearby star-forming regions, for example Orion, Musca and Taurus (Bally et al. 1987; Takahashi et al. 2013; Hacar et al. 2016; Kainulainen et al. 2015, 2017), where high resolution data (~ 0.01 pc, 0.1 km s⁻¹) reveals sub-structures like fibers (Hacar et al. 2013, 2018), or prominent mid-infrared extinction structures, like “Nessie” and infrared dark clouds like G11.11−0.12 (Johnstone et al. 2003; Pillai et al. 2006b; Schneider et al. 2010; Jackson et al. 2010; Kainulainen et al. 2013; Henshaw et al. 2014; Mattern et al. 2018). Detailed studies of these filaments led us to recognize their important role in star formation, and their internal structure, but studies of small samples do not allow to draw general conclusions. In particular the filaments towards the more distant, typically high-mass star forming regions have not yet been systematically studied. Therefore, it is necessary to study a large unbiased sample of filaments. Such studies have recently become feasible because of modern multiwavelength surveys, which cover the Galactic plane at high resolution and sensitivity.

Several catalogs of filamentary structures have been conducted in the last years, which can be divided in two groups. The filaments in the catalogs of Schisano et al. (2014), Koch & Rosolowsky (2015), and Li et al. (2016) were identified from continuum data and therefore, miss the kinematic information, and might be affected by line-of-sight projection effects. The catalogs of Ragan et al. (2014), Zucker et al. (2015), Abreu-Vicente et al. (2016), and Wang et al. (2015, 2016) concentrate on the longest filamentary structures in the Galaxy. While the identification methods and criteria vary in these studies, all filaments are tested for a velocity coherent behaviour.

In this study, we have targeted the largest catalog of filamentary structures published so far (Li et al. 2016), which is based on the ATLASGAL survey at 870 μm (Schuller et al. 2009). As these structures were identified in continuum dust emission data, the scope of this work is to use the SEDIGISM data (Schuller et al. 2017) to assess their velocity structure. Because of the large number of targets, it is necessary to perform the analysis in a fully-automated way, which are also presented in this work.

In this paper, we will refer to the structures identified by Li et al. (2016) as filament candidates. After the analysis of their velocity structure we will refer to the velocity coherent structures in the filament candidates as filaments, where one filament candidate can consist of multiple filaments. Some of these filaments may not meet the definitions of a filament, as they seem to be composed of a chain of dense clumps, or a dense clump with an elongated low column density environment. However, since filaments fragment, these structures could represent a late phase of evolution and should not be ignored.

The structure of the paper is as follows: Sect. 2 introduces the survey data used in this study and the targeted catalog of filament candidates. The methods used to separate the velocity components of a given filament candidate and to derive its filament parameters are described in Sect. 3. In Sect. 4 we present the resulting statistics of the velocity separation and the interpretation of the kinematics. We then discuss in Sect. 5 the dependency of the filament mass with increasing radius, and the origin of the correlation found between the line-mass (mass per unit length) and velocity dispersion of the filaments. Finally, we summarize our results in Sect. 6.

2 Data and filament sample

2.1 Survey data

Within this paper, we will make use of three surveys: ATLASGAL (APEX Telescope Large Area Survey of the Galaxy, Schuller et al. 2009), ATLASGAL+Planck (ATLASGAL combined with Planck, Csengeri et al. 2016) and SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic InterStellar Medium, Schuller et al. 2017).

The ATLASGAL survey was conducted with the Large APEX Bolometer Camera (LABOCA) at 870 μm between 2007 and 2010 at the Atacama Pathfinder Experiment (APEX) telescope (Güsten et al. 2006) located on the Chajnantor plateau in Chile. The resolution of the survey is 19.2′′ (6.0′′ per pixel) with a 1σ RMS noise in the range of 40–70 mJybeam⁻¹. It covers the inner Galactic plane between − 80° ≤ ℓ ≤ 60° and |b|≤ 1.5°. It is sensitive to the cold dust, and it traces mainly the high molecular hydrogen column density regions () of the ISM.

As the ATLASGAL data is missing the large scale low column density emission due to sky noise subtraction, Csengeri et al. (2016) combined the survey with the data observed by the HFI instrument at 353 GHz (850 μm) with a resolution of4.8′ on board the Planck satellite (Lamarre et al. 2010; Planck Collaboration I 2014). The combined ATLASGAL+Planck survey is sensitive to a wide range of spatial scales at a resolution of 21′′ covering thesame region as the original ATLASGAL data on the same pixel grid.

The SEDIGISM survey (Schuller et al. 2017) covers the inner Galactic plane between − 60°≤ ℓ ≤ 18° and |b | ≤ 0.5°, which was observed from 2013 to 2016 with the SHeFI heterodyne receiver (Vassilev et al. 2008) at the APEX telescope. The prime targets of the survey are the ¹³CO(2 – 1) and C¹⁸O(2 – 1) molecular lines. The average root-mean-square (RMS) noise of the survey is 0.9 K (T_MB) at a velocity resolution of0.25 km s⁻¹, an FWHM beam size of 30′′, and a pixel-size of 9.5′′. For this analysis we use the first data release (DR1, Schuller et al., in prep.).

2.2 The ATLASGAL sample of filaments

Based on ATLASGAL, Li et al. (2016) produced a catalog of filament candidates, which is the base for this study. The filaments were identified in the ATLASGAL only maps, after they were smoothed to a spatial resolution of 42′′. The source extraction was performed with the DisPerSE (Discrete Persistent Structures Extractor, Sousbie 2011) algorithm, which is optimized for the identification of large spatially coherent structures, and has been successfully used to trace filaments in previous studies (e.g., Hill et al. 2011; Arzoumanian et al. 2011). Because of the limited sensitivity and resolution (minimal mean column density ), the resulting catalog is unlikely to be complete, however, as it covers a large fraction of the Galactic plane it is likely to include the full range of sizes and masses of filamentary type structures.

Not all of the identified structures are filamentary, but they cover a range of morphologies and complexity from roundish clumps to large web-like structures. Therefore, the identified structures were categorized by Li et al. (2016) through visual inspection into six groups: unresolved clumps, marginally resolved elongated structures, filaments, networks of filaments, complexes, and unclassified structures. Here a filament was defined as single elongated linear structure with relatively few branches, an intensity clearly above the surrounding medium and an aspect ratio of at least 3, that is clearly resolved across its length and width. The high column densities found in the Galactic center region lead to a higher probability of identifying more complex structures. Therefore, the number of filament candidates in the catalog is lower towards the Galactic center. This classification resulted in a catalog of 517 filament candidates, providing the starting point of this study. For more details about the filament identification see Li et al. (2016). In the following, we will refer to the structures of the catalog as filament candidates, as they have been identified only in position–position space, which leaves the possibility of line-of-sight projection effects. One of the objectives of this study is to investigate their velocity coherence.

3 The automated filament analysis

The SEDIGISM survey covers the Galactic plane between − 60° ≤ ℓ ≤ 18° and |b | ≤ 0.5°, which is only a part of the ATLASGAL survey. Therefore, we analyze the 283 filament candidates in the area covered by all three surveys described in Sect. 2.1. This corresponds to ~ 55% of the total number of filaments, and can therefore be considered representative of such structures in the inner Galactic plane.

Because of the large number of filament candidates, it is necessary to use an automated approach to analyze them. However, as the sample is distributed over a large range of Galactic longitudes, it is unlikely to find homogeneous conditions in their surrounding material. Therefore, we choose a robust and efficient method to analyze the data in a systematic way, which leads to the following decisions for the analysis. We use the calculation of moments instead of multi-Gaussian fitting to identify the kinematics of the filaments. Also, we do not truncate or alter the skeletons of the filament candidates to fit the identified filaments more accurately, but rather neglect parts where we do not detect molecular gas. Therefore, there are two sets of pixels for a filament candidate used in the analysis: one that describes only the skeleton for the calculation of the kinematics and one that includes also the surrounding area within a dilation box with diameter of three beams used for the structure correlation. This approach results in larger uncertainties in the derived properties, but the homogeneous method enables the finding of correlations in the large scale properties of the filaments, which is the aim of this work.

From Li et al. (2016) we have a set of positions defining the skeleton of each filament, which trace the highest ATLASGAL intensities, that form the backbone of the structure. For each candidate we extract the data around the skeleton from the surveys using a rectangular box that is 5 ′ larger on each side than the extrema of the skeleton points. This showed to be sufficient for the analysis of the most nearby < 2 kpc filaments, where the angular extend of the structure is the largest. We now describe the analysis performed on every filament candidate.

3.1 The filament skeleton

For the analysis of each filament candidate we first have to transfer the skeleton coordinates (Fig. 1) onto the SEDIGISM grid. To do so, we check whether all positions of the skeleton are covered by the SEDIGISM observations, and remove the positions if they are not covered. This allows us to continue with structures that are partially truncated by the data limits. Then we overlay the skeleton coordinates on the pixel grid of the molecular line data. We mark the pixels within a dilation box around eachpixel that covers a position of the DisPerSE skeleton as part of the new pixel skeleton. The size of the dilation box is set to be larger than the maximum distance between two neighboring skeleton points. Here a width of 1 pixel (9.5′′) is sufficient. As the resulting skeleton mask might have a width larger than one pixel, we use the thinning algorithm of Gonzalez & Woods (1992) to truncate the pixel skeleton. The result is a “chain” of pixels which might have several branches (see Fig. 2).

Fig. 1 ATLASGAL contours (0.05, 0.1, 0.2, 0.5, 1.0, 2.0 Jy beam−1) and skeleton derived by DisPerSE for the filament candidate G333.297+00.073, overlaid on an infrared three color image of the field (red: MIPSGAL 24 μm; green: GLIMPSE 8.0 μm; blue GLIMPSE 3.6 μm.

Fig. 2 Skeleton of the filament candidate G333.297+00.073 derived by DisPerSE on top of the ATLASGAL grayscale contour map. The black contour indicates the dilation box used for the correlation in Sect. 3.3.

3.2 Identification of velocity components

In order to investigate whether the filament candidates form a single structure in velocity, spectroscopic observations are indispensable. From a study of filaments in the SEDIGISM first science field (seven filament candidates in 1.5 deg², Schuller et al. 2017), we know that one can observe line emission at very different velocities towards one continuum structure due to projection effects through the Galactic plane. Therefore, we average all spectra located on the skeletonand identify the velocity ranges that show emission peaks in this spectrum. To identify the velocities we smooth the average spectrum with a Gaussian kernel with a dispersion of four channels (= 1 km s⁻¹) to reduce the noise. In case the signal-to-noise ratio (S/N) is low, peak intensity ≤ 5σ, we double the kernel width. We then define the velocity range of each spectral component in the averaged spectrum as that between which the emission attains more than the 1σ noise level. This leaves us with a minimum separation limit of δv_min = 2.5 km s⁻¹, which is described in detail later on. Furthermore, we only consider components with an S/N ≥ 5 in their integrated intensity of the original data (e.g., see Fig. 3).

We then define the ends of the velocity components where the emission peaks of the average spectrum exceeds the 1σ noise level, so the emission peak is not likely to be truncated, and accept only components with an integrated intensity S∕N ≥ 5 (e.g., Fig. 3). The above procedure is only applied to the ¹³CO data because of their higher S/N. The C¹⁸O emission lines are narrowerthan the ¹³CO data lines and so we use the same velocity ranges to calculate the moments in the C¹⁸ O data. Thereafter, we calculate the zeroth, first, and second “order” moment of each velocity component, that indicate the integrated intensity, peak velocity, and velocity dispersion, respectively, for both molecules. This gives us a first impression of the kinematics of the filament.

Separating the velocity components of a filament along the line-of-sight is a crucial part of this work. Therefore, the technique of identifying the velocity range of emission needs to be tested in a systematic way. We created a simulated data cube with an RMS noise per position of 1 K, typical for SEDIGISM, and include two filaments at the same 2D location, with the same velocity dispersion and intensity, but with different peak velocities. For the emission of the filaments, we assume a Gaussian line profile. We then vary the peak-to-peak velocity, the linewidth, and the intensity and analyze these cubes in exactly the same way as the observed data. From this modeling, we find for filaments with an S∕N > 4, which is typical for our ¹³CO data (channel width 0.25 km s⁻¹), that the minimal separated peak-to-peak velocity depends linearly on the velocity dispersion σ_v like, (1)

shown in Fig. 4. Emission lines with a velocity dispersion of σ_v <0.75 km s⁻¹ (3 channels) are not identified as an emission line. For filaments with low intensities (S∕N ≤ 4) and for different intensities these limits must be degraded by Δσ_v = 0.25 km s⁻¹. We note, as the identification is done on the average spectrum over the filament skeleton, the velocity dispersion of an emission line can be larger than the intrinsic value because of velocity gradients along the skeleton. As a result, we are not resolving the kinematic substructure, like fibers (δv ≈ 1.0 km s⁻¹, Hacar et al. 2013), but can determine the large scale kinematics of the filament.

From previous studies (Schuller et al. 2017) we know that ¹³CO is likely to be optically thick towards the densest regions, which might lead to effects like self-absorption, and could affect the separation method. As the abundance of C¹⁸O is lower by ¹³CO ∕C ¹⁸ O = 8.3 (Miettinen 2012), it is likely to be optically thin over the whole filament. During visual comparison of the ¹³CO and C¹⁸ O spectra of all 604 ATLASGAL clumps within the filaments (position and velocity) optically thick ¹³CO emission is seen in 76 clumps, corresponding to 13%. However, these clumps do not show a significant effect on the average spectrum over the whole skeleton. Hence, our method is unlikely to separate velocity components because of self-absorption features. Additionally, we conclude that the contributionof the dense core emission is small when compared to the emission integrated up over the whole filament.

In the next step we use the same method as for the average spectrum on every pixel, hence spectrum, along the skeleton within the identified velocity ranges. In this way we identify which part of the skeleton is detected in the different velocity components. Velocity components in which less than ten positions of the skeleton are detected in ¹³CO are discarded, as these barely deviate from the noise, and we ensure a minimal elongated shape for all correlated structures (aspect ratio of 3 assuming a width of one beam). Additionally, we are able to detect multiple velocity components within the previously identified velocity range towards individual pixels. In the case where several velocity components are found, we only keep the calculated moments for the one with highest intensity. This is done for the ¹³CO and C¹⁸ O data. The separation of these subcomponents is limited by the smoothed velocity resolution. Also, we calculate the zeroth, first, and second “order” moments of skeleton pixels in each detected velocity component. With the first order measurements we potentially trace velocity gradients along the spine, but this is beyond the scope of this paper. However, the second order measurements, hence the velocity dispersion, of one pixel includes only the velocity gradient within one beam, which can be considered to be small. Therefore, the average over the skeleton pixels is a better estimate of the velocity dispersion of a filament than the value derived from the average spectrum and used in the further analysis. To check the results of these calculations we plot the derived moments on top of the position-velocity diagram of the filament candidate skeleton (Fig. 5). Additionally, we integrate over the velocity ranges of each velocity component; see Fig. 6.

Fig. 3 Average 13CO (top panel) and C18O (bottom panel) spectrum over the skeleton of filament candidate G333.297+00.073 (see Figs. 1 and 2). The red lines mark the identified emission intervals named by letters.

Fig. 4 Component separation limits dependent on the velocity dispersion of the two velocity components. The crosses indicate the modeled data points

Fig. 5 Position–velocity plot of the intensity along the skeleton of the filament candidate G333.297+00.073. The white stripes indicate the beginning/end of a skeleton branch. The first five branches show the longest connection through the skeleton from higher to lower galactic longitude and the last four are the branches to the side in the same direction. The horizontal red lines mark the identified emission intervals shown in Fig. 3 with intervals a, b, d, c from top to bottom. The jagged black and green lines mark the per pixel measured peak velocity and the 1σ interval of the detected emission peak.

3.3 Gas-dust correlation

We further estimate how much each velocity component contributes to the overall dust emission from a given filament candidate. We smooth and re-grid the ATLASGAL maps, which were used for the candidate identification, to the resolution (30.1′′) and pixel-grid (9.5′′) of the SEDIGISM cubes. For comparison with the ¹³CO integrated emission we restrict the maps (ATLASGAL and integrated ¹³CO intensities per velocity component) to an area within a dilation box around the skeleton with a width of 3 beams (9 pixels), which covers the emission seen in ATLASGAL. We scale the ATLASGAL intensities with the minimum and maximum value to an interval of [0–1], using (2)

where are the scaled ATLASGALintensities, and I^dust the original ones, and min() and max() describe the minimal and maximal value of the pixel intensities within the dilation box. As the ATLASGAL emission traces only the small scale high density gas, the minimun value is typically around 0 Jy beam⁻¹. For the ¹³CO data we define the intensity integrated over the velocity range of component i as , where T_mb (v) is the main beam temperature, and v_i is the velocity interval of one component as defined before. Also, for the correlation we do not apply any threshold. We use for the scaling the maximum and minimum value within the dilation box of the sum over the integrated intensity maps of all velocity components, (3)

where are the scaled ¹³CO intensities of the velocity component i, and I^gas(v_i) are the original ones. Assuming a constant gas-to-dust ratio, the ATLASGAL maps should correlate with the molecular line emission integrated over all velocity components. We perform pixel-to-pixel correlation between the scaled ATLASGAL maps and the scaled integrated ¹³CO maps of one velocity component using the same dilation box mask. Therefore, in cases of multiple components per candidate we will not finda one-to-one correlation, but we identify which velocity component shows the filamentary behaviour seen in dust emission. Additionally, noise in the observations and effects like CO depletion will affect the correlation plot. See Fig. 7 and Appendix A, for examples.

We analyze these correlation diagrams as follows. We calculate the standard deviation of the total (scaled) integrated intensities, , presenting an upper limit of the noise in the gas emission, and the standard deviation of the correlation, . To estimate the intensity contribution of one velocity component, we perform a linear fit with a slope of 1 on the data points with . We then derive the percentage of data points, which are above the σ_gas noise, p_gas, which are within ± σ_cor from the linear fit (red in Fig. 7), p_cor, and which meet both conditions, p_{cor, gas}. We then usethese three parameters to characterize the different velocity components with the limiting values shown in Table 1.

The limiting values for the characterization were obtained from the manual analysis of a test sample of filament candidates. This introduces a slight bias to the general analysis. The uncertainty of the characterization is difficult to determine, as the filament definition given by Li et al. (2016) can not be applied in a systematic way. Therefore, to get an objective, reproducible result, we decided to use this quantitative characterization over a visual approach as used by Li et al. (2016). A later visual inspection and different characterization will be possible, as all velocity components of all filament candidates, which are not uncorrelated, are handled in the same way in the subsequent analysis. Therefore, the characterization bias alters only the statistics of the characterization itself. However, the reliability of a filamentary shape decreases from fully correlated over partially correlated to diffuse structures, where partially correlated structures might not be continuous in position-position space, and diffuse structures might not show clearly enhanced emission from the surrounding. For simplicity we will still refer to all correlated structures as filaments. In total, we identified 422 filaments within the 283 filament candidates. More statistics of the characterization will be presented in Sects. 4.1 and 4.6.

Fig. 6 Integrated 13CO(2 – 1) emission of the four velocity components of the filament candidate G333.297+00.073. The intensity was integrated over the velocity ranges − 55.0 to − 40.5 km s−1 (top left panel), − 74.0 to − 66.0 km s−1 (top right panel), − 96.0 to − 86.0 km s−1 (bottom left panel) and − 82.5 to − 76 km s−1 (bottom right panel). The beamsize is indicated by the circle in the lower left and the contours show the ATLASGAL [0.1, 0.2, 0.5, 1.0, 2.0] Jy beam−1 levels. The letters in the upper right refer to the marked intervals in Fig. 3.

3.4 Thermal and non-thermal motions

The unresolved kinematic motions in a molecular cloud can be estimated by the observed linewidth. The total velocity dispersion can be separated into a thermal and a non-thermal component. The thermal motions depend on the observed molecule and gas temperature. For the molecular gas temperature we assume a typical value of T = 15 K (Pillai et al. 2006a; Wang et al. 2012, 2014), which is also in agreement with measured temperatures of ATLASGAL clumps (Urquhart et al. 2018). The non-thermal motions describe statistical motions of the gas, which are independent of the kinetic temperature of the the gas. The non-thermal velocity dispersion, σ_nt, is given by (4)

where σ_x is the measured second moment of ¹³CO or C¹⁸O, k_B the Boltzmann constant, m_p is the proton mass, and m_x is the molecular weight of the observed gas, here and . To derive the non-thermal velocity dispersion of the filament, σ_nt, we average the measurements of each pixel along the skeleton, where we choose C¹⁸ O if it is detected and ¹³CO otherwise. Therefore, this value is independent of velocity gradients along the skeleton, neglecting gradient effects within the beam, and it is less likely to be influenced by optical depth effects.

The thermal motion of the interstellar medium is given by the sound speed, (5)

where μ = 2.8 is the mean molecular weight of the mean free particle (Kauffmann et al. 2008), and other parameters as previously defined. The total velocity dispersion of a filament is given by .

3.5 Mass and length of filaments

To calculate physical parameters of the 422 filaments it is crucial to estimate the distance towards them. However, estimating distances towards extended and diffuse structures, like these filaments, is difficult. Especially, solving the ambiguity of kinematic distances. Therefore, we use a method similar to that discussed by Li et al. (2016), but including additional measurements.

As a first step we identify all ATLASGAL clumps (Contreras et al. 2013b; Urquhart et al. 2014) associated with the filaments. The distances towards these clumps have been estimated in Urquhart et al. (2018). As these estimates are based on kinematic distances, we must exclude the Galactic center region (|ℓ| < 5°), because of the large uncertainties. For filaments with an associated clump within the defined limits of the filament in position–position-velocity space) we simply assume the same distance. This provides distances for 222 filaments. In a second step we use friends-to-friends analysis to find adjacent clumps and adopt their distances for the filaments. This adds distances for another 114 filaments, but note the larger uncertainty for the distanceestimate. For the friends-of-friends analysis we allow sources with a spatial offset of at most 10′ (with 90% closer than 5′) and a kinematic offset of at most 10 km s⁻¹ (with 90% closer than 4 km s⁻¹). In total, we are able to assign a distance to 336 out of 422 filaments, including diffuse, partially correlated, and fully correlated ones. Additionally, we tested these estimates to be in agreement with (one of) the kinematic distances.

After we obtained the distances towards the filaments we can calculate their physical length. Here we take all pixels along the skeleton into account, towards which ¹³CO was detected for the single filament. This allows us to get accurate measurements of the angular length for complex structures or partially correlated structures, including the detected branches, by adding up the length over each relevant pixel. However, because of a possible inclination of the structures, the derived physical measurements represent lower limits to their true length.

For calculating the area and the mass of the filaments we again use a dilation box. With the distances in hand we are able to use a physical box-diameter. Here we take the box-diameter as a free parameter and in Sect. 5.1 we discuss the dependency between the filament mass and box-diameter. For calculating the filament mass we assume its diameter to be on the order of the star-forming size scale. For a first order approximation of the star-forming size scale we use the previously measured velocity dispersion within a filament, σ_v, and assume a star formation time of T_sf = 2 Myr (Evans et al. 2009). The size scale is then given by s_sf = σ_v ⋅ T_sf. As box-diameters are limited to discrete multiples of pixels, we interpolate linearly the values from measurements of the next bigger and smaller box size.

We estimate the area covered by the filament, by summing over all pixels within the dilation box of the integrated intensity maps, within which we detect ¹³CO emission (S∕N > 5). The same positions are used to estimate the filament mass. Here we follow two different approaches: first we use the integrated ¹³CO emission, , in combination with the ¹³CO X-factor, derived by Schuller et al. (2017). This has the advantage of tracing only the emission within the specific velocity component. The molecular hydrogen column density, N_i(H₂), in pixel i was then calculated by . We then computed the mass using the equation (6)

where N_i(H₂) is the H₂ column density computed for pixel i, A_i its area, μ =2.8 the mean molecular weight per H₂ molecule, and m_p the proton mass.

Second, we estimate the mass from dust emission maps of different surveys (ATLASGAL, ATLASGAL+Planck) using basic assumptions like a gas-dust ratio of R = 100, and a dust temperature of T_D = 15 K (Urquhart et al. 2018). The mass of the filament candidate is then computed through the equation (7)

where S_ν is the integrated flux density at the frequency ν of the used survey, d is the distance towards the structure, B_ν (T_D) is the Planck function at the given dust temperature, and κ_ν is the dust absorption coefficient, which is κ_870μm = 1.85 cm² g⁻¹ for the ATLASGAL emission (Ossenkopf & Henning 1994; Schuller et al. 2009).

However, because of the contribution of the Planck data, the ATLASGAL+Planck survey traces not only the filament and the low column density gas around the filament, but also the diffuse Galactic dust emission, which ideally should be removed. To do so, we exclude the filament area of the maps using the inverse of the filament dilation box used previously. The remaining pixels should be dominated by non-filament emission. However, as the box has a width of three beams, we find a few cases in the most nearby (< 2kpc) filaments where the emission extends clearly beyond the mask. Therefore, we use the 20th percentile value of the non filament pixels as estimate of the diffuse Galactic dust emission. We then correct the ATLASGAL+Planck maps for the diffuse emission and estimate the masses as shown previously. In Sect. 4.7 we discuss the differences of the these three mass estimates based on dust continuum emission and compare them to the ¹³CO emission estimate.

Fig. 7 Gas-dust correlation plot of the brightest velocity component (“a” in the spectrum shown in Fig. 3) of the filament candidate G333.297+00.073. The blue line gives the one-to-one correlation. The green area indicates values above the σgas limit (pgas = 0.51). The red line shows the fitting result, and the area within the dashed red lines marks the ± σcor surrounding (pcor = 0.80). pcor, gas = 0.70 is estimated from the overlap of these areas.

4 Results

4.1 Final catalog

Using the methods described in the previous section we derived a large set of filament parameters. With these parameters we created a catalog of velocity coherent structures. The derived parameters of the catalog are shown in Table 2, and the complete catalog is shown in Tables 5 and 6. However, as several parameters are distance dependent, they cannot be calculated for the whole filament catalog. The following description of the derived parameters includes only the structures with a distance estimate. The catalog contains all 422 filaments of the 283 observed filament candidates, which show some correlation with the dust emission. In Table 3 we show the statistics of the characterization and the subsample for which we have distance estimates. The names of the structures are based on the initial filament candidate name from Li et al. (2016) and are extended by an integer starting from 0, indicating the velocity component.

4.2 Detection of filaments in ¹³CO and C¹⁸O

Out of the 283 ATLASGAL filament candidates within the SEDIGISM survey we detect correlated ¹³CO emission for 260 filament candidates, which then show 422 velocity coherent (continuous kinematic structure, which cannot be resolved into separate components) filaments. We do not find a correlated ¹³CO velocity component for 23 filament candidates, which is partially because of the sensitivity of the SEDIGISM survey, and partially because the candidates result from line-of-sight alignments of diffuse gas clouds. About 20% of the detected filaments show ¹³CO emission at every position of the skeleton andfor about 60% we detected ¹³CO emission over at least half of the length ofthe skeleton. About 32% of the ¹³CO detected filaments show no detection of C¹⁸O on the skeleton, about 13% have C¹⁸O detected over at least half of the skeleton, and for no filament C¹⁸ O is detected over its entire length; see Fig. 8.

This difference in the detection rate is very likely due to the different abundances of the molecules (¹³CO∕C¹⁸O = 8.3; Miettinen 2012). The C¹⁸O line is expected to be weaker, resulting in a lower S/N and the observed lower detection rate.

Fig. 8 Cumulative histogram of the percentage of filament candidate skeletons detected in 13CO (blue) and C18O (green).

Fig. 9 Filaments with distance estimated plotted onto an artist’s impression of the Milky Way (Robert Hurt), the size is indicating the length of the filament, and the color indicates the 13CO velocity. Theyellow lines mark the range |ℓ|≤ 5°, where distances are uncertain and the blue lines mark the SEDIGISM survey limits.

4.3 Galactic distribution

Using the distance estimates we can derive the Galactocentric coordinates and plot the positions onto a face-on artist’s impression of the MilkyWay (Fig. 9). We find that a large fraction of the filaments are likely to be associated with the near Scutum-Centaurus arm. We also find some filaments located in the near Sagittarius arm, the near and far 3-kpc arm, and the near Norma arm, but also in some inter-arm regions. We note that we do not have distance estimates for the Galactic center region (|ℓ| ≤ 5°).

Histograms representing positions of the detected filaments with Galactic longitude and latitude are shown in Figs. 10 and 11, respectively. The distribution with Galactic longitude shows a peak around l = −21° with a strong decrease towards the outer Galaxy (only 14 structures for l < −45°), and a decrease towards the Galactic center. As the filament candidates were identified in the ATLASGAL maps that trace only high columndensity dust emission, it is more unlikely to find filaments towards the outer Galaxy, which contains fewer dense molecular cloud regions. However, the number of filaments is also suppressed in the direction of the Galactic center, where identification is difficult because many structures along the line-of-sight are confused, and were categorized as networks, complexes or unclassified, such as the Galactic center region (Li et al. 2016). Nevertheless, comparing the distribution of filaments to the distribution of ATLASGAL clumps presented by Beuther et al. (2012, Fig. 3) and Csengeri et al. (2014, Fig. 16) shows similarities for the location of peaks, indicating a possible correlation with active star-forming regions.

The distribution of the detected filaments with Galactic latitude (Fig. 11) shows a broad, almost flat behaviour similar to the findings of Li et al. (2016). However, we find the peak and mean (⟨b⟩ = +0.02°) of the distribution aligned with the Galactic mid-plane, which is in contrast to the general finding of more sources for b < 0.0° than for b > 0.0° (Beuther et al. 2012; Li et al. 2016). We note that our sample is not identical with that of Li et al. (2016), as we use only a sub-sample of their filament candidates and split some of these candidates in different velocity components, hence filaments.

Fig. 10 Distribution of the filament positions in Galactic longitude. The orange hatch marks the filaments with distance estimate.

4.4 Distributions of velocity dispersion, mass, length and distance

In the following, we give a short overview on the most interesting measured properties of the filaments, which are the non-thermal velocity dispersion along the skeleton, the mass derived from the ¹³CO emission, and the length of the filament, which we define as the sum over the detected parts of the skeleton.

For the calculation of the total velocity dispersion we assumed an isothermal medium of 15 K, see Sect. 3.4. The distribution of the resulting values is shown in Fig. 12. We find values reaching from about 0.5 to 2.5 km s⁻¹ with a relatively flat center between 0.8 and 1.4 km s⁻¹, and a mean of about 1.17 km s⁻¹. Concentrating on the 180 fully correlated, hence, the most reliable filaments, we find a similar distribution with the mean at about 1.20 km s⁻¹. In general, these values are higher than what Arzoumanian et al. (2013) find in nearby filaments (σ_v ≈ 0.3 km s⁻¹), but in agreement with studies of similar (more distant and more massive) objects like the DR21 filament (Schneider et al. 2010).

For the logarithmic distribution of the calculated masses, Fig. 13, we find a flat part between 1800 M_⊙ and 18 000 M_⊙ with a mean mass of 8600 M_⊙. Again the distribution of the fully correlated filaments (151 with distance estimate) is similar to the complete distribution with a mean of 11 000 M_⊙. For comparison we also show the mass ranges covered by other filament catalogs (Ragan et al. 2014; Wang et al. 2015; Zucker et al. 2015; Abreu-Vicente et al. 2016) and the study of Contreras et al. (2013a). These studies report filaments with similar or higher masses. However, several of these studies have tried to identify the largest structures in the Galaxy and therefore, are biased to larger structures. As a result, some filaments mentioned in this study are only parts of structures in the other catalogs. Also, we find filaments that are almost identical in several catalogs, like G11.046-00.069_2 (Snake) or G332.370-00.080_1.

We also find overlap between the catalogs for the lengths of the filaments (Fig. 14), where the shortest filaments of the other studies are as long as the mean of our sample (10.3 pc all, 11.1 pc fully correlated). In general, we cover the range from 2 to 100 pc with a peak around 8 pc.

Most filaments are found within 5 kpc from the Sun (Fig. 15), which is also the area where the other surveys found the long filaments. This area includes parts of the nearby Sagitarius and Scutum-Centaurus spiral arms. Another peak in the distance distribution is found at around 10 kpc, which is about the distance of the connection point of the Galactic bar with the Perseus spiral arm (see also Fig. 9).

Plottingthe filament lengths against the estimated distances separated by the categories (Fig. 16), we find that especially the fully correlated filaments follow the distance distribution, while the others are more equally distributed. Also, we find no correlation between the longest filaments and the distance. This results in a larger scatter of lengths for a given distance. The shortest filaments reproduce our minimal length criteria of at least 10 pixels.

Fig. 11 Distribution of the filament positions in Galactic latitude. The orange hatch marks the filaments with distance estimate.

Fig. 12 Distribution of the measured total velocity dispersion of all filaments (blue). The fully correlated filaments are marked by the orange-colored hatched area. The vertical lines indicate the mean value for the complete (black) and the sub-sample (orange).

Fig. 13 Distribution of the measured mass based on 13CO of all filaments (blue). The fully correlated filaments are marked by the orange hatch. The vertical lines indicate the mean value for the complete (black) and the sub-sample (orange). The horizontal lines mark the mass ranges measured by the studies mentioned above the lines.

Fig. 14 Distribution of the measured length over the detected skeleton of all filaments (blue). The fully correlated filaments are marked by the orange hatch. The vertical lines indicate the mean value for the complete (black) and the sub-sample (orange). The horizontal lines mark the length ranges measured by the studies mentioned above the lines.

4.5 ¹³CO – C¹⁸O velocity comparison

As mentioned before, C¹⁸O is less abundant than ¹³CO and traces mainly the bright, dense parts of the filaments, where ¹³CO is likely to be optically thick. However, to combine the kinematics of the two lines we need to be sure that both trace the same gas. Therefore, Fig. 17 shows the distribution of the absolute difference between the ¹³CO and C¹⁸ O peak velocities of each filament derived from the average spectrum along the full skeleton, which is supposed to be zero if both istotopologues trace the same gas. The logarithmic distribution shows a plateau between 0.17 and 1.0 km s⁻¹ and decreases steeply to both sides. Additionally, we find the largest difference in filaments with an S/N of the C¹⁸ O average spectra of S∕N < 5. In general, low S/N filaments peak at higher velocity differences (red).

We compare the observed distribution to a model distribution given by the mean velocity dispersion along the filament (see Sect. 4.4 and Fig. 12). Given the wide distribution of velocity dispersions this model gives only a first order impression of the expected distribution. We model the absolute difference between two velocities drawn from two Gaussian distributions. The dispersion of the differences is then given by . We generate an artificial difference distributions, using 10 000 draws to avoid statistical noise, bin the absolute values of the sample like the observed differences, and scale the height by 0.0373 to get a comparable total number of filaments as our sample. The resulting distribution (orange hatched in Fig. 17) does not agree with the observed one, as it is shifted to larger differences.

To further investigate this trend, we reduced the dispersion of the underlying velocity distribution until we found a distribution that matched the differences (black hatched) area. Its velocity dispersion is σ_v (model) = 0.35 km s⁻¹, which is about times the channel width (0.25 km s⁻¹). We speculate therefore, that this distribution is likely to arise from the sampling of the spectra.

However, we also see some filaments that show a larger difference between the ¹³CO and C¹⁸ O peak velocities than can be expected by the channel-width introduced sampling issues. For these filaments we speculate that they show a gradient along the skeleton and are only partially detected in C¹⁸ O. To investigate this, we plot the velocity difference against the ¹³CO velocity dispersion of the average spectrum (Fig. 18), as gradients along the skeleton result in a higher velocity dispersion. Additionally, filaments for which less than 10 % of the skeletons are detected in C¹⁸O are marked in red. We find that almost all filaments fall below the one-to-one line and that all velocity differences are smaller than . We also seethat 31 out of 43 filaments with a velocity difference larger than 1 kms⁻¹ show low C¹⁸ O detection rates.

In summary we rule out systematic differences between the kinematic of the isotopologues. The observed differences are based on observational limitations, like the velocity resolution and sensitivity.

Fig. 15 Distribution of the estimated distances of all filaments (blue). The fully correlated filaments are marked by the orange hatch. The vertical lines indicate the mean value for the complete (black) and the sub-sample (orange).

Fig. 16 Filament length plotted against the estimated distance. The three filament categories are indicated by blue, green, and red for diffuse component, partially correlated, and fully correlated, respectively. The black line shows the criteria of a minimum length of 10 pixels.

Fig. 17 Histogram of the absolute velocity difference. Sources with S∕N < 5 are shown in red. On top the model difference distributions given by a underlying Gaussian velocity distribution with a dispersion of1.17 km s−1 (orange) and0.35 km s−1 (black) are shown. The green dashed line indicates the velocity channel size of the data.

4.6 Multiplicity in velocity space

Filamentarystructures are often identified in continuum emission maps. But it is unknown whether these structures are actual continuous filaments or only an effect of line-of-sight projection of multiple velocity components. We address this question with our data.

The 260 detected filament candidates split up in 422 velocity coherent filaments in total. Kinematic subcomponents are identified in single spectra for 14 of the filaments, but will not be discussed any further as more detailed studies will be needed. Analysis of the velocity components shows that about 58% of the filament candidates exhibit one velocity component. Another significant fraction of the filaments, 27 and 12%, have 2 or 3 components, respectively. Six filaments have 4 or more velocity components with a maximum of 7 components, seen in only one filament (Fig. 19).

The categorization of the velocity components shows that a filament candidate can have several velocity components even in the case of one component being fully correlated. This is shown in Fig. 20. However, a filament candidate with a single component does not necessarily have a fully correlated structure. In general, we find that filament candidates with fewer velocity components are more likely to have a fully correlated component, and candidates with an increasing number of velocity components are more likely to have partially correlated and diffuse components.

Many structures are identified on continuum data that does not provide information about the velocity coherence. Therefore, we test whether there exists a correlation between the intensities of 2D-data and the number of velocity components. With the known multiplicity of the filament candidates, we are able to show that filament candidates with several velocity components tend to be brighter (see Fig. 21). We do so as follows: we derive the mean and maximum intensity of the ATLASGAL dust emission and the ¹³CO emission integrated over all velocity components. Because the statistical scatter of the intensity values shows a non-Gaussian distribution, we take the median and the 90th percentile value of the intensity distributions of each filament category (i.e. separated according to the number of components) as a qualitative measure. We also estimate theuncertainty of the median using a bootstrapping method. However, only the sample sizes of the filament categories with 1, 2, and 3 components are sufficient to use the bootstrapping method. In this method we draw new, random samples of intensities from among the observed values. We then calculate the median of these new, simulated samples. The resulting distribution of the median values then estimates the sampling function of the observed median and is used to estimate the uncertainty using the standard deviation. We find, that the medians of the ¹³CO peak intensities increase outside their uncertainties as the number of velocity components increases. The same increase is also seen for the 90th percentile values (Fig. 21). Our data suggests a similar increase for the ATLASGAL peak intensities, but a flat behaviour is also consistent with the data. We could not find such a behaviour for the mean intensities of the filaments (not shown in figures).

For filament candidates with multiple components we investigate whether two physically separate filaments can be located within the same spiral arm. Therefore, we show in Fig. 22 a histogram of the absolute peak-to-peak velocity difference (blue). The bins up to δv = 10 km s⁻¹ are likely tobe incomplete because of the component separation limit, as shown before (Fig. 4, Eq. (1)). We compare the distribution of the observed velocity differences with model distributions (hatched) of expected velocity difference from a spiral arm. We assume velocity dispersions of and following Reid et al. (2016) and Caldu-Primo et al. (2013). As we measure the absolute difference between two velocities drawn from a Gaussian distribution, the dispersion of the differences is given by . We sample the difference distributions with 10, 000 draws to avoid statistical noise, bin the absolute values of the sample like the observed differences, and scale by 0.016 to get a comparable total number of filaments. We find, that the observed and the model distribution of are similar for δv ≤ 30 km s⁻¹, but we see more observed filaments for larger velocity separations. The model distribution for does not describe the observedone. Therefore, we can conclude that a large fraction of separated filaments might be located in the same spiral arm with a velocity dispersion of , but we also see filaments from different Galactic structures along the line-of-sight. However, because of the kinematic distance ambiguity filaments located in different spiral arms can have similar line-of-sight velocities at specific Galactic Longitudes.

Fig. 18 Absolute velocity difference against the velocity dispersion of the filament, where filaments with a C18 O detection rate below 10% are indicated in red. The blue lines show the one-to-one (1σ, solid) and two-to-one (2σ, dashed) relations.

Fig. 19 Histogram of the number of velocity components per filament.

Fig. 20 Histogram of the number of velocity components for fully correlated (blue), partially correlated (green), and diffuse (red)structures.

4.7 Comparison of masses derived from gas and dust

Calculating the masses of the filaments is an important part of the analysis. However, it comes with some difficulties. Commonly, dust emission or dust extinction is used to calculate the mass of objects. However, because of the line-of-sight projection several filaments that may appear at the same position, and cannot be separated in the continuum data, this is not applicable here. Therefore, we need to use the CO emission to disentangle the projected emission from different structures in velocity space. For the mass estimate we then use the emission integrated over the velocity range of the filament. Specifically, we use the ¹³CO emission as it has a higher S/N and traces the lower column density gas around the skeleton, and calculate the mass like described in Sect. 3.5.

We first have to test whether this X-factor is a good approximation for the whole survey area. To do so we take a sample of filaments, which show only one velocity component and are fully correlated with the ATLASGAL dust emission. We calculate the masses for this sample using the integrated ¹³CO emission with Eq. (6) and using the ATLASGAL and ATLASGAL+Planck dust emission with Eq. (7). For all threedata-sets we use the same mask around the skeleton. The comparison of the resulting masses is shown in Fig. 23.

We find that masses derived from ¹³CO are systematically larger than the masses derived from ATLASGAL, but systematically smaller than the masses derived from ATLASGAL+Planck. This behaviour is expected as the ATLASGAL maps are sensitive to the small scale (2.5 ′), high columndensity dust structures and extended emission from the diffuse surrounding gas is filtered out because of the observing technique (see Schuller et al. 2009). Therefore, Csengeri et al. (2016) combined the ATLASGAL data with the Planck data, which traces the diffuse, large scale structures, but does not resolve the small scales because of the low resolution (4.8 ′). However, the combined data also traces the dust emission along the line-of-sight, i.e foreground and background. Thus, masses derived from ATLASGAL are likely to be underestimated while masses derived from ATLASGAL+Planck are likely to be overestimated.

As shown in Sect. 3.5, we corrected the ATLASGAL+Planck data for the line-of-sight emission towards every filament and used this data to derive another mass estimate. On average we find a mean Galactic emission of 0.52 Jy beam⁻¹ (beam size of 21′′). These corrected dust masses are in agreement with the ¹³CO derived masses within a factor of 2. Therefore, we conclude that the ¹³CO X-factor derived from the SEDIGISM science demonstration field (Schuller et al. 2017) is a good approximation for the whole survey area.

Fig. 21 Peak integrated 13CO intensity (top panel) and peak ATLASGAL dust intensity (bottom panel) of a filament against the number of velocity components. The blue circles mark the median, and the red crosses mark the 90th percentile for each number of components. The error-bars show the uncertainty of the median derived by bootstrapping.

Fig. 22 Histogram of the absolute difference in velocity between the neighboring velocity components of a filament. Upper right panel: complete distribution, main panel: only the lower velocity separations. The black and orange hatched distributions indicate the models for spiral arm velocity dispersions of 5 km s −1 and 10 km s−1. The dashed red line indicates the average velocity separation limit of 3.5 km s−1 (see Eq. (1) and mean velocity dispersion of 1.17 km s−1).

Fig. 23 Mass per fully correlated filament derived from dust versus the mass derived from integrated 13CO using an factor from Schuller et al. (2017). The black line indicates the one-to-one ratio.

5 Discussion

5.1 Radial filament profiles

Nearby (<500 pc), low line-mass (<100 M_⊙ pc⁻¹) filaments have been found to have an FWHM size on the order of 0.1 pc (Arzoumanian et al. 2011). The corrected ATLASGAL+Planck and ¹³CO data trace the wide range of column densities that is needed to study the filament profile. To ensure that we are looking only at true filaments we restrict our sample to the 151 fully correlated filaments with a distance estimate, but for completeness we show the results of the other filaments in Appendix B. However, measuring the filament profile is challenging as most filaments are not homogeneous, linear structures, but show branches and varying central densities.

Therefore, we do not extract the radial column density distribution directly from the data, but estimate the mass of the filaments within filament masks with increasing diameter, s_box, using the same equations and assumptions as in Sect. 3.5. The mass, M(R), is then given by (8)

where l is the length of the filament, Σ(r) is the column density of the gas at distance r from the skeleton (to both sides, with the skeleton at r = 0 pc), and R = 0.5 s_box the maximum radius. We normalize the values with the mass from a box-diameter of s_max = 4 pc, where the typically found radial profile is almost flat (Arzoumanian et al. 2011). The smallest box-diameter is given by the pixel size. The resulting mass curves of the ¹³CO emission areshown in Fig. 24, and of the continuum emission and less correlated filaments in Appendix B. As the physical resolution is changing with the distance towards the source, we group the filaments in four distance intervals d_i (d₁ < 2 kpc < d₂ < 4 kpc < d₃ < 8 kpc < d₄) and average the mass curves within these intervals (see Fig. 26).

The profiles of filamentary structures are found to be well described by a Plummer-like density distribution (Nutter et al. 2008; Arzoumanian et al. 2011; Contreras et al. 2013a), which is given by (9)

where ρ_c is the central density of the filament, and R_flat the characteristic radius of the flat inner part. The column density profile of the filament (Arzoumanian et al. 2011; Panopoulou et al. 2014)then is described by (10)

where Σ = N(H₂)μm_p, μ and m_p as previously defined, and , a finite constant for p > 1. Other studies (e.g., Arzoumanian et al. 2011) have shown that the inner part of the radial profile of a filament can also be well described by a Gaussian column density distribution. These two models are shown in Fig. 25, where they are normalized to an integrated intensity of 1.

The mass within a box around the theoretical filament is then given by Eq. (8). We plot the measured mass curves as well as the theoretical ones as M(R)∕M(2 pc) and test different values for the exponent, p, and the inner radius, R_flat, of the Plummer-like distribution and for the dispersion, w, of the Gaussian (see Fig. 26).

A detailed analysis of the density structure of filaments is beyond the scope of this paper. Still, we perform a rough visual comparison between the observed radial column density profiles and modeled ones. We find that the Plummer-like distribution is in agreement with the average profile of the ¹³CO and dust observationsfor p ≈ 1.5 ± 0.5 and R_flat ≈ 0.1 pc. Also, the Gaussian column density distribution with a dispersion of w = 1 pc describes the observation within the uncertainties. However, the two fitting models lead to different FWHMs for the filaments. While for a Plummer-like function the (FWHM_2.0 ≈ 0.17 pc, FWHM_1.5 ≈ 0.39 pc; Heitsch 2013b), for the Gaussian 2.36 pc.

One possible interpretation is that the Gaussian traces only the low column density surrounding of the filament, but not the dense inner part, hence the actual filament. From previous studies (see Arzoumanian et al. 2011; Panopoulou et al. 2014) and the Plummer-like function we see measurements of the FWHM between 0.1 and 0.6 pc. The physical beam size of the SEDIGISM data at a distance of 2 kpc is about 0.3 pc and therefore, we are at the resolution limit for the dense filament spine, but note that the mass curves (integral over the radial profile) have a dependency on R_flat in the Plummer-like case. However, we find that small changes of R_flat do not significantly change the agreement with the observation. Therefore, we conclude that the mass of the filament is dominated by the low column density gas surrounding the spine, and that the resolution of the SEDIGISM data is not sufficient to meaningfully fit the inner spine with a Gaussian radial profile.

The exponent of the average density profile is p ≈ 1.5, which is in agreement with the lower limit found by Arzoumanian et al. (2011). The single profiles scatter between p ≈ 1.0 and p ≈ 2.0, with the scatter decreasing with more distant filaments most likely because of the smaller sample. Also, we tested the effect of the beam size to the theoretical radial profiles by convolving the profile with a Gaussian beam. The resulting theoretical mass curves are shallower with increasing distance, but not significantly, given the scatter of the single observed mass curves. The study of Arzoumanian et al. (2011) analyzes prominent filaments in nearby molecular clouds using dust continuum emission. This selection of prominent filaments might give a bias towards higher exponents. Theoretically, an isolated, isothermal, cylindrical filament in hydrostatic equilibrium would be expected to have an exponent of p = 4 (Ostriker 1964). However, this exponent is typically not found in observations and models (Juvela et al. 2012; Kainulainen et al. 2015). Low resolution and S/N data explain only partially the observed exponents. Therefore, observations suggest that filaments are embedded in a surrounding molecular cloud (Fischera & Martin 2012), not isothermal (Recchi et al. 2013) and/or not in hydrostatic equilibrium (Heitsch 2013a,b).

It is important to mention that for an exponent of p < 2 mathematically the mass diverges with an increasing radius. This can be seen in Fig. 25. As a result the mass, M, and therefore, also the line-mass (mass per unit length), m = M∕l, are not well-defined. However, filaments are not isolated structures, but surrounded by low density gas, which sets boundary conditions that are not considered in this model. As mentioned before, we decided to use a radius dependent on the velocity dispersion of the filament to estimate the mass of the filaments in this study.

Fig. 24 Fraction of the filament mass derived from 13CO emission dependent on the box-diameter of the mask separated by distances. Top left panel: d < 2 kpc, top right panel: 2 kpc < d < 4 kpc, bottom left panel: 4 kpc < d < 8 kpc, bottom right panel: d > 8 kpc. One curve describes one fully correlated filament at its distance estimate. The gray lines indicate the physical beamsize at distances of 2, 3, 6, and 8 kpc. The black lines show the integrated theoretical radial profiles, which describe a Plummer-like distribution p = 1.5 (dashed) or p = 2.0 (dash-dotted), and a Gaussian distribution with a dispersion of w = 1.0 (dotted).

Fig. 25 Theoretical filament profiles normalized to an integrated intensity of 1, which describe a Plummer-like function with Rflat = 0.1 pc and p = 1.5 (dashed) or p = 2.0 (dash-dotted), and a Gaussian with a dispersion of w = 1.0 pc (dotted). The radial integration of these profiles is shown in Figs. 24 and 26.

Fig. 26 Average fraction of the filament mass derived from 13CO emission dependent on the box-diameter of the mask. The color indicates the distance of the filament with d1 orange, d2 red, d3 blue, and d4 green. The errorbars indicate the dispersion of the measured mass fraction and box-diameter. The black lines are same as Fig. 24.

5.2 Stability against collapse

Thermally supercritical filaments are commonly seen as star formation sites. Therefore, they need to build the connection between the diffuse gas of the molecular cloud and the dense gas in the star-forming cores. Inutsuka & Miyama (1992) showed that isothermal, infinitely long, self-gravitating cylinders will collapse radially if their line-mass (mass per unit length) exceed a critical value, and fragment along the axis in the sub-critical and equilibrium case. The critical line-mass is given by (11)

(Ostriker 1964), where G is the gravitational constant and c_s is the sound speed of the medium, which is dependent on the gas temperature T (Eq. (5)). Assuming a typical gas temperature of T = 15 K the critical line-mass is m_crit = 20 M_⊙ pc⁻¹.

Based on our observations and analysis we can estimate the line-mass for all the filaments with a distance estimate by m_obs = M∕l, where M is the mass estimated from the ¹³CO emission, and l is the length along the velocity coherent skeleton. Because of the separation of the velocity components this length does no longer securely describe the linear shape (with small branches) of the original filament candidate sample, especially for the not fully correlated filaments. Therefore, we concentrate this discussion on the fully correlated filaments, but also perform the calculations for the others.

The line-masses we observe with our resolution are significantly above the critical thermal value (see Fig. 27). This leaves us with two possible conclusions: either the filaments are collapsing radially or they have a supporting mechanism additional to the thermal pressure. Moreover, we find that the linewidth of the molecular gas is significantly larger than the sound speed, c_s = 0.21 km s⁻¹. This increased linewidth can support both theories, as it can be interpreted as either structured motions, like collapse, or turbulent motions within the gas.

Assuming that non-thermal motions contribute to the supporting mechanism, Eq. (11) can be modified to (12)

(Fiege & Pudritz 2000), where σ_nt is the non-thermal velocity dispersion of the filament, and m_crit,tot is the critical, total (thermal and non-thermal) line-mass. After determining the velocity dispersion for all filaments, we can calculate the critical non-thermal line-mass and compare it with the observed one (Fig. 27).

The uncertainty of the critical line-mass is given by the observed velocity dispersion and therefore depends on the velocity resolution and the quality of the signal. The main contributions for the uncertainties of the mass estimates are the X-factor (factor of 0.5–2), optically thick ¹³CO emission (factor of 1–2), and the distance. As the length is also dependent on the distance estimate, the line-mass is only linearly dependent on the distance, which adds another factor of 0.8–1.2 to the uncertainty. Additionally, the length is measured as projection on the sky and therefore, the observed line-mass is an upper limit considering possible inclinations. The typical uncertainty is given by the black cross in Fig. 27. Additionally, it needs to be mentioned that based on the resolution of our data we are only able to derive global parameters. Higher resolution data (spatial and kinematic) could reveal substructures, which might lead to different results.

We find that the critical, non-thermal line-mass is, within the uncertainties, in agreement with the observed line-mass. Therefore, Eq. (12) seems to describe a common relation between the observed linewidth and line-mass in the form . The sound speed, c_s, depends only on the temperature of the ISM, which can be assumed to be about constant. Hence, the line-mass is proportional to the non-thermal motion. We also find that partially correlated filaments and diffuse components follow the same relation as the fully correlated filaments, but with a slightly wider spread.

We now want to investigate where this relation comes from. As we discussed before, one explanation might be infall motion. Inutsuka & Miyama (1992) showed that infinitely long, isothermal filaments with a line-mass above the critical value collapse radially. Heitsch et al. (2009) and Heitsch (2013b) determined the accretion velocity profile, v(R), for gas in a steady-state free-fall onto the filament axis as (13)

where G is the gravitational constant, m is the line-mass of the filament, R_ref is the limiting, outer radius of the filament, and R is the radial position of the gas.

We use this radial velocity distribution, v(R), to estimate the signal which would be observed from a collapsing filament, similar to Heitsch (2013a). First, we derive the line-of-sight velocities, v_lsr, across the filament for an observer looking edge on (see Fig. 28), (14)

where R is the radial distance to the center, and x is the position in the x-axis direction of the Cartesian coordinate system. Second, we draw for each position in the filament 50 values from a Gaussian distribution centered on the derived velocities with a thermal velocity dispersion of c_s = 0.21 km s⁻¹. Third, we plot a weighted histogram of the velocities with bins identical to the SEDIGISM channel width of 0.25 km s⁻¹, where we use the density at the position of the filament as weight. The density is given by a Plummer-like distribution (see Sect. 5.1). From the weighted histogram we calculate the standard deviation, hence the theoretically observed velocity dispersion.

Within this template we vary the line-mass, m, the exponent of the density profile, p, and the area of the filament we observe (complete, middle, edge) to study their effect on theoretical signal. We choose the outer radius R_ref = 2.0 pc and the characteristic radius of the density distribution R₀ = 0.1 pc, see Sect. 5.1. The results are shown in Table 4 and Fig. 29.

We find that the infall motions indeed show the observed relation between line-mass and velocity dispersion, (Table 4, models a–c). However, the theoretical velocity dispersion is higher by a factor of 2 than the observed values. Additionally, we did not take turbulent motions into account, which would be caused by the infall (Heitsch 2013b), and give a widertheoretical signal. To fit the observations, the collapse needs to be slowed down, which can be caused by the turbulence created by the collapse. Finding a combination of collapse and turbulent motions that would reproduce the observed kinematics is beyond the scope of this study.

Anotherway to identify ongoing collapse could be an analysis of the shape of the emission lines (Schneider et al. 2010). The theoretical line profiles show double-peaked shapes towards the center of the filament and for filaments with a steep density profile. However, a comparison of the theoretical line shapes with the observed ones is difficult as several other effects can alter the shape, like spatially unresolved motions of the filament within the beam, or self-absorption, also see Heitsch (2013a) for different inclination angles. In a few filaments we find some evidence for a double-peaked velocity profile (see Sect. 3.2). But a more detailed analysis of, preferably, geometrically simple filaments with higher spatial and kinematic resolution would be necessary to address the effects of complex kinematics.

Rapid radial collapse would lead to extremely narrow filaments, which rarely have been observed up to now (e.g., Stutz & Gould 2016). However, Heitsch (2013b) shows that the fragmentation timescales are shorter than the collapse timescales, which is supported by the fact that fragmentation is seen in almost all filaments (Jackson et al. 2010; Kainulainen et al. 2013, 2017; Takahashi et al. 2013; Wang et al. 2014; Beuther et al. 2015; Teixeira et al. 2016). Our finding of slowed collapse would increase the difference between the timescales, which is still in agreement with the observations. Also, simulations of filament evolution by Chira et al. (2018) show a start of fragmentation before the filaments reach the critical mass for gravitational collapse.

However, we also cannot rule out the possibility that the observed velocity dispersion is mainly created by turbulent motions. These turbulent motions are discussed to increase the internal pressure and support the filament against gravitational collapse. Therefore, this theory is also in agreement with the observations.

Fig. 27 Critical, non-thermal line-mass derived from the velocity dispersion versus observed line-mass using integrated 13CO. The fully correlated filaments are indicated in blue, and the other filaments of this study in a shaded green and red. The black solid line shows the one-to-one correlation, the gray lines indicate a factor 2 uncertainty, and the dashed line shows the critical thermal line-mass.

Fig. 28 Line-of-sight velocity distribution across a slice of a filament base on Eq. (13) using a line-mass of m = 500 M⊙ pc−1 and Rref = 2.0 pc. The observer is located on the right side, observing the whole slice, the inner part (dashed lines), and outer part (dotted lines), which where analyzed separately.

6 Conclusions

In this study we studied spectral line emission from 283 filament candidates detected with ATLASGAL continuum dust emission from the catalog of Li et al. (2016) in the SEDIGISM ¹³CO and C¹⁸ O survey. As these candidates can be the result of line-of-sight projection of multiple structures, we tested the candidates for coherence in velocity space and derived the mass, size, and a collection of kinematic properties. To do so we developed an automated analysis tool that finds the different velocity components of a candidate, if existing, separates them and checks for correlation with the original ATLASGAL emission. We found 422 velocity-coherent filaments that correlate completely or partially with the original candidate. For these filaments we find the following:

• Two-thirds of the filament candidates are single velocity-coherent structures. The other candidates are line-of-sight projections of mainly two and three velocity components, and up to one candidate with seven velocity components. Also, we found a possible indication for a correlation between the maximum intensity within the filament candidate and the number of velocity components for the integrated ¹³CO and ATLASGAL dust emission, but a flat behaviour is within the uncertainties.

• Comparing the kinematics of the filaments seen in ¹³CO and C¹⁸O, we could show that both isotopologues trace the same gas. Differences found in the comparison could be identified as biases arising from low S/N C¹⁸O data.

• The filament profiles are on average in agreement with a Plummer-like density distribution with an exponent of p ≈ 1.5 ± 0.5. The inner radius cannot be constrained exactly because of the limited resolution of the data. This low exponent indicates that filaments are typically located within larger molecular clouds, and therefore, the outer radius of a filament cannot be well defined. For the mass estimates we chose a radius which includes the gas that can take part in star formation within the next 2 Myr.

• The observed line-mass of the filaments is in agreement with the critical non-thermal line-mass and significantly higher than the critical thermal line-mass. However, we do not know the source of the observed velocity dispersion. Comparing the relation we find between velocity dispersion and line-mass with the theoretical infall velocity profile based on Heitsch (2013b) generally does not reveal evidence for free-fall collapse. However, radial infall of the gas onto the skeleton can possibly explain the relation.

In this study we analyzed the kinematics of 283 filament candidates, finding 180 reliable velocity coherent filaments, 151 with distance estimates between 1 kpc and 13 kpc, and 242 other velocity coherent filamentary structures in the line-of-sight of the candidates, leading to the largest statistics of filament parameters so far. However, due to the spatial resolution of 30′′ and velocity resolution of0.25 km s⁻¹, the derived parameters generally only describe global behaviour of the filaments. As the evolution and fragmentation of filaments is a hierarchical process it will be necessary to also study the smaller scales. High resolution observations to recover the small scales (< 0.1 pc) are essential, and this study can be the starting point for the selection of a representative sample for such higher resolution follow-ups.

Fig. 29 Theoretical molecular line profiles with varying parameters, see Table 4.

Acknowledgments

We want to thank the anonymous referee for constructive and very detailed comments. Also, we thank Guang-XingLi for supporting this study by offering additional data and information on the filament candidates. Also, we thank Prof. KostasTassis, Dr. Dario Colombo, Carsten König, and the rest of the SEDIGISM consortium for support and productive discussions. M.M. is supported for this research through a stipend from the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities of Bonn and Cologne. T.Cs. acknowledges support from the Deutsche Forschungsgemeinschaft, DFG via the SPP (priority program) 1573 “Physics of the ISM”. H.B. acknowledges support from the European Research Council under the Horizon 2020 Framework Program via the ERC Consolidator Grant CSF-648505. L.B. acknowledges support by CONICYT project Basal AFB-170002. This work of J.K. has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 639459 (PROMISE).

Appendix A Gas dust correlation examples

In Sect. 3.3 we are describing the method of comparing the ATLASGAL dust intensities with the integrated ¹³CO intensity maps to identify the spatially correlated velocity components. Based on this correlation, we assign one of four categories to the structure, which are uncorrelated, completely correlated, partially correlated, and diffuse component. To visualize this categorization, we show one example in Sect. 3.3 (Figs. 6, and 7). To give a complete picture, we present here the correlation plots of the other components of the filament candidate “G333.297+00.073” (Fig. A.2, for intensity maps see Fig. 6), and ¹³CO spectra (Fig. A.1), integrated intensity maps and the corresponding correlation plots for the filament candidates “G339.116-00.405” (Fig. A.3) and “G346.293+00.109” (Figs. A.4 and A.5).

Fig. A.1 Average 13CO spectrum over the skeleton of filament candidate “G339.116-00.405” (top panel) and “G346.293+00.109” (bottom panel). The red lines mark the identified emission intervals named by letters.

Fig. A.2 Gas-dust correlation plots of the velocity components “b”,“c”, and “d” of the filament candidate G333.297+00.073 (see Figs. 3, 5, and 6). The blue line gives the one-to-one correlation. The green area indicates values above the σgas limit. The red line shows the fitting result, and the area within the dashed red lines marks the ± σcor surrounding. pcor,gas is estimated from the overlap of these areas.

Fig. A.3 Integrated 13CO intensity map (left panel) and gas-dust correlation plot (right panel, description as in Fig. A.2) of the single velocity component filament candidate “G339.116-00.405”. Because of the good correlation between the only identified velocity component with the ATLASGAL emission, it was categorized as fully correlated filament and can be found as “G339.116-00.405_0” in the final catalog.

Fig. A.4 Integrated 13CO intensity map (left panel) and gas-dust correlation plot (right panel, description as in Fig. A.2) of the first velocity component of the filament candidate “G346.293+00.109”. The partially correlated filament is listed as “G346.293+00.109_0”. The other three velocity components of the candidate are shown in Fig. A.5.

Fig. A.5 Integrated 13CO intensity maps (left panels) and gas-dust correlation plots (right panels, description as in Fig. A.2) of the second (top panels), third (middle panels), and fourth (bottom panels) velocity component of the filament candidate “G346.293+00.109”. The second diffuse component is listed as “G346.293+00.109_1” and the third partially correlated component is listed as “G346.293+00.109_2”. The uncorrelated component is not listed in the catalog. The first velocity components of the candidate is shown in Fig. A.4.

Appendix B Filament profiles

In Sect. 5.1 we showed how the gas mass of correlated filaments is increasing with increasing radius. To be complete, we also show the mass curves of the correlated filaments based on the corrected ATLASGAL+Planck data (Fig. B.1), and the mass curves of the partially correlated (blue, Fig. B.2) and diffuse component filaments (red, Fig. B.2) based on the integrated ¹³CO observations.

Fig. B.1 Fraction of the filament mass derived from corrected ATLASGAL+Planck dust continuum emission dependent on the box-diameter of the mask separated with distances. Top left panel: d < 2 kpc, top right panel: 2 kpc < d < 4 kpc, bottom left: 4 kpc < d < 8 kpc, bottom right panel: d > 8 kpc. The gray lines indicate the physical beam size at these distances. The black lines show the theoretical profiles, which describe a Plummer-like p = 1.5 (dashed) or p = 2.0 (dash-dotted), and a Gaussian with w = 1.0 (dotted).

The mass curves of of the correlated filaments are in mostly in agreement with ¹³CO observations. However, for the most nearby filaments (<2 kpc), and especially in the continuum data, we find profiles indicating a slope of p < 0. This can be explained by line-of-sight confusion within the large boxes around the filament skeleton. As the dust continuum data traces all emission along the line-of-site, it is possible that strong emission, which is not related to the filament but located nearby, is taken into account for the mass estimate for larger radii. Therefore, the masses will be overestimated. This effect is more likely for nearby filaments, as larger angular sizes are taken into account for the same physical size. The gas mass curves of the partially correlated and diffuse component filaments show on average similar results as the correlated filaments, but with a larger scatter, as the skeleton no longer necessarily represents the shape of the structure.

Fig. B.2 Fraction of the filament mass for partially correlated (blue) and diffuse component (red) filaments derived from 13CO emission dependent on the box-diameter of the mask separated with distances. Top left panel: d < 2 kpc, top right panel: 2 kpc < d < 4 kpc, bottom left panel: 4 kpc < d < 8 kpc, bottom right panel: d > 8 kpc. The vertical gray lines indicate the physical beam size at these distances. The black lines show the theoretical profiles, which describe a Plummer-like p = 1.5 (dashed) or p = 2.0 (dash-dotted), and a Gaussian with w = 1.0 (dotted).

References

All Tables

All Figures

	Fig. 1 ATLASGAL contours (0.05, 0.1, 0.2, 0.5, 1.0, 2.0 Jy beam−1) and skeleton derived by DisPerSE for the filament candidate G333.297+00.073, overlaid on an infrared three color image of the field (red: MIPSGAL 24 μm; green: GLIMPSE 8.0 μm; blue GLIMPSE 3.6 μm.
In the text

	Fig. 2 Skeleton of the filament candidate G333.297+00.073 derived by DisPerSE on top of the ATLASGAL grayscale contour map. The black contour indicates the dilation box used for the correlation in Sect. 3.3.
In the text

	Fig. 3 Average 13CO (top panel) and C18O (bottom panel) spectrum over the skeleton of filament candidate G333.297+00.073 (see Figs. 1 and 2). The red lines mark the identified emission intervals named by letters.
In the text

	Fig. 4 Component separation limits dependent on the velocity dispersion of the two velocity components. The crosses indicate the modeled data points
In the text

Fig. 5 Position–velocity plot of the intensity along the skeleton of the filament candidate G333.297+00.073. The white stripes indicate the beginning/end of a skeleton branch. The first five branches show the longest connection through the skeleton from higher to lower galactic longitude and the last four are the branches to the side in the same direction. The horizontal red lines mark the identified emission intervals shown in Fig. 3 with intervals a, b, d, c from top to bottom. The jagged black and green lines mark the per pixel measured peak velocity and the 1σ interval of the detected emission peak.

In the text

Fig. 6 Integrated 13CO(2 – 1) emission of the four velocity components of the filament candidate G333.297+00.073. The intensity was integrated over the velocity ranges − 55.0 to − 40.5 km s−1 (top left panel), − 74.0 to − 66.0 km s−1 (top right panel), − 96.0 to − 86.0 km s−1 (bottom left panel) and − 82.5 to − 76 km s−1 (bottom right panel). The beamsize is indicated by the circle in the lower left and the contours show the ATLASGAL [0.1, 0.2, 0.5, 1.0, 2.0] Jy beam−1 levels. The letters in the upper right refer to the marked intervals in Fig. 3.

In the text

Fig. 7 Gas-dust correlation plot of the brightest velocity component (“a” in the spectrum shown in Fig. 3) of the filament candidate G333.297+00.073. The blue line gives the one-to-one correlation. The green area indicates values above the σgas limit (pgas = 0.51). The red line shows the fitting result, and the area within the dashed red lines marks the ± σcor surrounding (pcor = 0.80). pcor, gas = 0.70 is estimated from the overlap of these areas.

In the text

	Fig. 8 Cumulative histogram of the percentage of filament candidate skeletons detected in 13CO (blue) and C18O (green).
In the text

	Fig. 9 Filaments with distance estimated plotted onto an artist’s impression of the Milky Way (Robert Hurt), the size is indicating the length of the filament, and the color indicates the 13CO velocity. Theyellow lines mark the range |ℓ|≤ 5°, where distances are uncertain and the blue lines mark the SEDIGISM survey limits.
In the text

	Fig. 10 Distribution of the filament positions in Galactic longitude. The orange hatch marks the filaments with distance estimate.
In the text

	Fig. 11 Distribution of the filament positions in Galactic latitude. The orange hatch marks the filaments with distance estimate.
In the text

	Fig. 12 Distribution of the measured total velocity dispersion of all filaments (blue). The fully correlated filaments are marked by the orange-colored hatched area. The vertical lines indicate the mean value for the complete (black) and the sub-sample (orange).
In the text

	Fig. 13 Distribution of the measured mass based on 13CO of all filaments (blue). The fully correlated filaments are marked by the orange hatch. The vertical lines indicate the mean value for the complete (black) and the sub-sample (orange). The horizontal lines mark the mass ranges measured by the studies mentioned above the lines.
In the text

	Fig. 14 Distribution of the measured length over the detected skeleton of all filaments (blue). The fully correlated filaments are marked by the orange hatch. The vertical lines indicate the mean value for the complete (black) and the sub-sample (orange). The horizontal lines mark the length ranges measured by the studies mentioned above the lines.
In the text

	Fig. 15 Distribution of the estimated distances of all filaments (blue). The fully correlated filaments are marked by the orange hatch. The vertical lines indicate the mean value for the complete (black) and the sub-sample (orange).
In the text

	Fig. 16 Filament length plotted against the estimated distance. The three filament categories are indicated by blue, green, and red for diffuse component, partially correlated, and fully correlated, respectively. The black line shows the criteria of a minimum length of 10 pixels.
In the text

	Fig. 17 Histogram of the absolute velocity difference. Sources with S∕N < 5 are shown in red. On top the model difference distributions given by a underlying Gaussian velocity distribution with a dispersion of1.17 km s−1 (orange) and0.35 km s−1 (black) are shown. The green dashed line indicates the velocity channel size of the data.
In the text

	Fig. 18 Absolute velocity difference against the velocity dispersion of the filament, where filaments with a C18 O detection rate below 10% are indicated in red. The blue lines show the one-to-one (1σ, solid) and two-to-one (2σ, dashed) relations.
In the text

	Fig. 19 Histogram of the number of velocity components per filament.
In the text

	Fig. 20 Histogram of the number of velocity components for fully correlated (blue), partially correlated (green), and diffuse (red)structures.
In the text

	Fig. 21 Peak integrated 13CO intensity (top panel) and peak ATLASGAL dust intensity (bottom panel) of a filament against the number of velocity components. The blue circles mark the median, and the red crosses mark the 90th percentile for each number of components. The error-bars show the uncertainty of the median derived by bootstrapping.
In the text

Fig. 22 Histogram of the absolute difference in velocity between the neighboring velocity components of a filament. Upper right panel: complete distribution, main panel: only the lower velocity separations. The black and orange hatched distributions indicate the models for spiral arm velocity dispersions of 5 km s −1 and 10 km s−1. The dashed red line indicates the average velocity separation limit of 3.5 km s−1 (see Eq. (1) and mean velocity dispersion of 1.17 km s−1).

In the text

	Fig. 23 Mass per fully correlated filament derived from dust versus the mass derived from integrated 13CO using an factor from Schuller et al. (2017). The black line indicates the one-to-one ratio.
In the text

Fig. 24 Fraction of the filament mass derived from 13CO emission dependent on the box-diameter of the mask separated by distances. Top left panel: d < 2 kpc, top right panel: 2 kpc < d < 4 kpc, bottom left panel: 4 kpc < d < 8 kpc, bottom right panel: d > 8 kpc. One curve describes one fully correlated filament at its distance estimate. The gray lines indicate the physical beamsize at distances of 2, 3, 6, and 8 kpc. The black lines show the integrated theoretical radial profiles, which describe a Plummer-like distribution p = 1.5 (dashed) or p = 2.0 (dash-dotted), and a Gaussian distribution with a dispersion of w = 1.0 (dotted).

In the text

	Fig. 25 Theoretical filament profiles normalized to an integrated intensity of 1, which describe a Plummer-like function with Rflat = 0.1 pc and p = 1.5 (dashed) or p = 2.0 (dash-dotted), and a Gaussian with a dispersion of w = 1.0 pc (dotted). The radial integration of these profiles is shown in Figs. 24 and 26.
In the text

	Fig. 26 Average fraction of the filament mass derived from 13CO emission dependent on the box-diameter of the mask. The color indicates the distance of the filament with d1 orange, d2 red, d3 blue, and d4 green. The errorbars indicate the dispersion of the measured mass fraction and box-diameter. The black lines are same as Fig. 24.
In the text

Fig. 27 Critical, non-thermal line-mass derived from the velocity dispersion versus observed line-mass using integrated 13CO. The fully correlated filaments are indicated in blue, and the other filaments of this study in a shaded green and red. The black solid line shows the one-to-one correlation, the gray lines indicate a factor 2 uncertainty, and the dashed line shows the critical thermal line-mass.

In the text

	Fig. 28 Line-of-sight velocity distribution across a slice of a filament base on Eq. (13) using a line-mass of m = 500 M⊙ pc−1 and Rref = 2.0 pc. The observer is located on the right side, observing the whole slice, the inner part (dashed lines), and outer part (dotted lines), which where analyzed separately.
In the text

	Fig. 29 Theoretical molecular line profiles with varying parameters, see Table 4.
In the text

	Fig. A.1 Average 13CO spectrum over the skeleton of filament candidate “G339.116-00.405” (top panel) and “G346.293+00.109” (bottom panel). The red lines mark the identified emission intervals named by letters.
In the text

Fig. A.2 Gas-dust correlation plots of the velocity components “b”,“c”, and “d” of the filament candidate G333.297+00.073 (see Figs. 3, 5, and 6). The blue line gives the one-to-one correlation. The green area indicates values above the σgas limit. The red line shows the fitting result, and the area within the dashed red lines marks the ± σcor surrounding. pcor,gas is estimated from the overlap of these areas.

In the text

Fig. A.3 Integrated 13CO intensity map (left panel) and gas-dust correlation plot (right panel, description as in Fig. A.2) of the single velocity component filament candidate “G339.116-00.405”. Because of the good correlation between the only identified velocity component with the ATLASGAL emission, it was categorized as fully correlated filament and can be found as “G339.116-00.405_0” in the final catalog.

In the text

	Fig. A.4 Integrated 13CO intensity map (left panel) and gas-dust correlation plot (right panel, description as in Fig. A.2) of the first velocity component of the filament candidate “G346.293+00.109”. The partially correlated filament is listed as “G346.293+00.109_0”. The other three velocity components of the candidate are shown in Fig. A.5.
In the text

Fig. A.5 Integrated 13CO intensity maps (left panels) and gas-dust correlation plots (right panels, description as in Fig. A.2) of the second (top panels), third (middle panels), and fourth (bottom panels) velocity component of the filament candidate “G346.293+00.109”. The second diffuse component is listed as “G346.293+00.109_1” and the third partially correlated component is listed as “G346.293+00.109_2”. The uncorrelated component is not listed in the catalog. The first velocity components of the candidate is shown in Fig. A.4.

In the text

Fig. B.1 Fraction of the filament mass derived from corrected ATLASGAL+Planck dust continuum emission dependent on the box-diameter of the mask separated with distances. Top left panel: d < 2 kpc, top right panel: 2 kpc < d < 4 kpc, bottom left: 4 kpc < d < 8 kpc, bottom right panel: d > 8 kpc. The gray lines indicate the physical beam size at these distances. The black lines show the theoretical profiles, which describe a Plummer-like p = 1.5 (dashed) or p = 2.0 (dash-dotted), and a Gaussian with w = 1.0 (dotted).

In the text

Fig. B.2 Fraction of the filament mass for partially correlated (blue) and diffuse component (red) filaments derived from 13CO emission dependent on the box-diameter of the mask separated with distances. Top left panel: d < 2 kpc, top right panel: 2 kpc < d < 4 kpc, bottom left panel: 4 kpc < d < 8 kpc, bottom right panel: d > 8 kpc. The vertical gray lines indicate the physical beam size at these distances. The black lines show the theoretical profiles, which describe a Plummer-like p = 1.5 (dashed) or p = 2.0 (dash-dotted), and a Gaussian with w = 1.0 (dotted).

In the text