Surveys of clumps, cores, and condensations in Cygnus X: Temperature and nonthermal velocity dispersion revealed by VLA NH3 observations

The physical properties, evolution, and fragmentation of massive dense cores (MDCs, $\sim$ 0.1 pc) are fundamental pieces in our understanding of high-mass star formation. We aim to characterize the temperature, velocity dispersion, and fragmentation of the MDCs in the Cygnus X giant molecular cloud and to investigate the stability and dynamics of these cores. We present the Karl G. Jansky Very Large Array (VLA) observations of the NH$_3$ (J,K) = (1,1) and (2,2) inversion lines towards 35 MDCs in Cygnus X, from which we calculated the temperature and velocity dispersion. We extracted 202 fragments ($\sim$ 0.02 pc) from the NH$_3$ (1,1) moment-0 maps with the GAUSSCLUMPS algorithm. We analyzed the stability of the MDCs and their NH$_3$ fragments through evaluating the corresponding kinetic, gravitational potential, and magnetic energies and the virial parameters. The MDCs in Cygnus X have a typical mean kinetic temperature T$_K$ of $\sim$ 20 K. Our virial analysis shows that many MDCs are in subvirialized states, indicating that the kinetic energy is insufficient to support these MDCs against their gravity. The calculated nonthermal velocity dispersions of most MDCs are at transonic to mildly supersonic levels, and the bulk motions make only a minor contribution to the velocity dispersion. Regarding the NH$_3$ fragments, with T$_K$ $\sim$ 19 K, their nonthermal velocity dispersions are mostly trans-sonic to subsonic. Unless there is a strong magnetic field, most NH$_3$ fragments are probably not in virialized states. We also find that most of the NH$_3$ fragments are dynamically quiescent, while only a few are active due to star formation activity.


Introduction
High-mass stars (> 8 M ⊙ ) play a vital role in shaping their host galaxies via their intense feedback effects in the form of powerful outflows, strong stellar winds, and violent UV radiation (Motte et al. 2018).However, the formation of high-mass stars is still an open question.Unlike the star formation model built in low-mass stars (Shu et al. 1993), accretion onto high-mass protostars is difficult to model because of the radiation pressure problem (Wolfire & Cassinelli 1987) and the violent wind (Krumholz 2015).In observations, it is difficult to catch the high-mass stellar embryos as their birthplaces are usually far away and deeply embedded (Zinnecker & Yorke 2007).The small populations and short timescales of critical evolutionary phases mean that the sample of high-mass stellar embryos is even more insufficient (Zinnecker & Yorke 2007).
Several candidate theories of high-mass star formation have been proposed, such as competitive accretion in protostellar ⋆ Table of the physical parameters of NH 3 fragments is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/* after publishing.clusters (Bonnell et al. 1997(Bonnell et al. , 2001)), mergers of protostars (Bonnell et al. 1998), turbulent core accretion (McKee & Tan 2003), and inertial inflow (Padoan et al. 2020).To verify these theories, it is necessary to observationally investigate the initial conditions of high-mass star formation.Observations with a range of tracers in a variety of wavelengths and resolutions have been carried out towards massive dense cores (MDCs), the possible birth places of high-mass stars (e.g., Beuther et al. 2002;Motte et al. 2007;Lu et al. 2014;Rathborne et al. 2016;Cao et al. 2019;Zhang et al. 2019).By characterizing the physical conditions of highmass star formation with observations, we can test and provide practical constraints to the corresponding theories.
The CENSUS project (Surveys of Clumps, CorEs, and CoN-denSations in CygnUS X, PI: Keping Qiu) is dedicated to a systematic study of the 0.01-10 pc hierarchical molecular structures and the high-mass star formation process in the molecular complex Cygnus X.In order to investigate the temperature structures and dynamic conditions of MDCs in our project, we make use of the NH 3 (1,1) and (2,2) inversion emission lines, which trace gas with critical density of ∼ 10 4 cm −3 and excitation temperature of ∼ 10-30 K (Ho & Townes 1983;Mangum & Shirley 2015).These lines are widely used for studying the initial conditions of Article number, page 1 of 37 arXiv:2403.03845v1[astro-ph.GA] 6 Mar 2024 A&A proofs: manuscript no.VLA_NH3 star formation.Previous NH 3 surveys find a mean temperature of ≲ 15 K in less-evolved structures and ∼ 20 K in protostellar cores with significant star-forming activity (e.g., Pillai et al. 2006;Wu et al. 2006;Chira et al. 2013;Kerr et al. 2019).The velocity dispersion of these lines is found to be dominated by nonthermal components in most cases (e.g., Lu et al. 2014;Billington et al. 2019).Furthermore, analyses of the dynamical stability of starforming structures at various scales and evolutionary stages can be done with these lines via the virial analysis (e.g., Jijina et al. 1999;Zhang et al. 2020).Therefore, a high-resolution unbiased NH 3 survey towards the MDCs in the Cygnus X cloud complex offers a remarkable observational sample that can be used to unravel the initial conditions of high-mass star formation.
The Cygnus X cloud, located at 1.4 kpc from the Sun (Rygl et al. 2012), is one of the richest and most active molecular and HII complexes in the Galaxy.It is massive (∼ 4 × 10 6 M ⊙ ) and large (extends over ∼ 100 pc in diameter) (Leung & Thaddeus 1992), and consists of substantial HII regions (Wendker et al. 1991) and OB associations (Uyanıker et al. 2001).All these features make Cygnus X an ideal object for studying high-mass star formation.Several surveys and numerous case studies towards Cygnus X have been carried out over recent decades (e.g., Motte et al. 2007;Gottschalk et al. 2012;Takekoshi et al. 2019).
The CENSUS project has published several works on Cygnus X: Cao et al. (2019) generated a complete sample of MDCs in Cygnus X and performed an analysis of their far-infrared emission; Wang et al. (2022) provide radio continuum views of the MDCs; Pan et al. (2023) provide a survey of the possible circumstellar disks in Cygnux X; Yang et al. (2023) analyzed the highvelocity outflow using Submillimeter Array (SMA) observations of .In the CENSUS project, we adopt the terminology that clumps, cores, and condensations correspond to cloud structures of a physical full width at half maximum (FWHM) size of ∼ 1, ∼ 0.1, and ∼ 0.01 pc, respectively, which is also widely used in other studies (e.g., Pokhrel et al. 2018;Zhang et al. 2019).
As a part of CENSUS, we report an unbiased NH 3 survey towards 37 fields in Cygnus X observed with the Karl G. Jansky Very Large Array (VLA) with a spatial resolution of ∼ 0.02 pc.This paper is organized as follows.Section 2 outlines details of the sample and observational setups.In Sect.3, we describe the methods we used to analyze the NH 3 data.The derived physical parameters and their distributions and correlations are presented in Sect. 4. The NH 3 morphology, a virial analysis, the dynamical states of MDCs, and a comparison between dust and gas temperatures are discussed in Sect. 5. Finally, the study is summarized in Sect.6, where we outline our conclusions.

Sample selection
The interferometer observations of the CENSUS project (PI: Keping Qiu, also see Cao et al. 2019) were designed to target most MDCs in Cygnus X (Wang et al. 2022).We first selected 42 sources with estimated masses of greater than 30 M ⊙ from Motte et al. (2007), though one source (S26) was later found to be outside of Cygnus X (Rygl et al. 2012) and another (N58) should have a much lower mass (see Table 1 for more details).In this work, we exclude 10 sources in the DR 21 region, which were observed by the VLA in a mosaic setup (Project ID: 14A-241; PI: Qizhou Zhang) and will be presented in a separate work.The other 32 sources, as well as an additional source (NW12, see Table 1 for more details), were observed with the VLA in 29 fields.In the present work, we also include VLA observations of another eight fields covering bright JCMT SCUBA-2 sources (Cao et al. 2019) located in the OB2 association region1 .Here, VLA observations of a total of 37 fields are therefore presented (Table 1).As part of the CENSUS project, Cao et al. (2021) constructed a column density map of the entire Cygnus X complex based on SED fitting of multiband Herschel observations and extracted more than 8000 cores with the getsources algorithm.In the following analysis, we adopt the parameters of dense cores from Cao et al. (2021).
The VLA data were calibrated using the Common Astronomy Software Applications (CASA) package (McMullin et al. 2007).The flux calibrator and bandpass calibrator were both set as the quasar 1331+305 (3C286).Gain calibration was performed with observations of J2007+4029.The calibrated data were imaged using the CASA task clean, gridded with a pixel size of 1 ′′ .The rms of these fields is ≲ 0.013 Jy beam −1 in a 0.1 km s −1 channel width.The coordinates of the field centers and noise levels are also listed in Table 1.The NH 3 (1,1) and (2,2) emission lines are the main interests of this paper.The rest frequencies of the (1,1) and (2,2) emission lines are 23.6945 and 23.72263 GHz, while upper energy levels are 23.4 and 64.9 K, respectively.Observations of the other spectral lines in 17A-107 are beyond the scope of this paper, and are listed in Appendix A.

MIR archive data from the Spitzer Space Telescope
We use the Spitzer IRAC 3.6, 4.5, and 8.0 µm data from the Spitzer Legacy Survey of the Cygnus X Complex3 (Kraemer et al. 2010) in this paper, which are downloaded from the NASA/IPAC Infrared Science Archive4 .All these data cover the entire Cygnus X region, with a pixel size of 2 ′′ and 1σ rms noise levels of 25.1, 27.1, and 47.9 MJy•sr −1 , respectively.

Methods
We checked detections in all data cubes before fitting the NH 3 lines.In order to determine whether there is detection in a spec-trum, the NH 3 (1,1) data cubes are smoothed to a velocity resolution of 0.7 km s −1 , or seven spectral channels in our observations, to improve signal-to-noise ratios (S/Ns).The criterion of significant detection within a beam is set as S/N ≥ 5 in the smoothed NH 3 (1,1) data.The corresponding original spectra were fitted in the following processes if their smoothed spectra exceeded the threshold.The detections and noise levels in all the fields are reported in Table 2.The detections in fields are marked as marginal if the emissions are only detected in one or two beams.We adopt the distance of 1.4 kpc (Rygl et al. 2012) for all the fields except for Field 21 (also known as AFGL2591), of which the distance is 3.3 kpc (Rygl et al. 2012).

Ammonia line fitting
Following the theoretical framework in Mangum & Shirley (2015), we built spectral models of the NH 3 (1,1) and (2,2) emission with the centroid velocity, line width, excitation temperature, NH 3 column density in the upper transition state, and beam-filling factor as free parameters and with the following assumptions: (1) the same beam-filling factor and excitation temperature for all hyperfine transitions; (2) local thermodynamic equilibrium (LTE).The NH 3 (1,1) and (2,2) data cubes were fitted simultaneously using the model to derive the aforementioned physical parameters.
We can get the relationship between the brightness temperature and excitation temperature as where f is the beam-filling factor and T bg is the background temperature.We take the temperature of cosmic microwave background (CMB) as T bg , which is 2.73K.The optical depth τ ν can be written as where N u is the column density of molecules in the upper transition state.Here, we assume a Gaussian profile, where ν 0 is the central frequency of the emission.The Einstein coefficient A ul is where µ lu and R s are the permanent dipole moment of NH 3 and the relative strengths of the hyperfine structures, respectively.In our model, R s is set as 1/2, 5/36, and 1/9 for the main, inner satellite, and outer satellite components of NH 3 (1,1) emission, respectively.As for the main and inner satellite components of NH 3 (2,2) emission, R s is set as 0.796 and 0.052, respectively.We can derive the total NH 3 column density N tot from the relationship between N u and N tot (Boltzmann distribution), where the rotation partition function Q rot ≈ 168.7 T 3 ex /B 2 C (B = 298.117GHz, C = 186.726GHz), and g u is the total degeneracy for the upper energy level, which can be found in Mangum & Shirley (2015).Using these formulas, we built a blended spectral model rather than the all-hyperfine model after comparing their fitted results.This model consists of five Gaussian components for NH 3 (1,1) emission and three Gaussian components for NH 3 (2,2) emission.We do not include the two outer satellite components in our model as most NH 3 (2,2) lines do not show remarkable satellite components in our data.The excitation temperature and the column density are derived only when the S/Ns of the (2,2) moment-0 maps are greater than 2 in order to avoid large uncertainties.

Temperature threshold and kinetic temperature
Regarding the NH 3 inversion emission lines, the excitation temperature refers to the rotation temperature of NH 3 , T rot .We set the upper limit and lower limit of T rot in the fitting processes to 50 K and 8 K, respectively.The lower limit refers to the typical temperature of interstellar molecular clouds without starforming activity, that is, ∼ 10K, and the upper bound corresponds to the fact that the NH 3 (1,1) and (2,2) lines cannot be a good temperature probe when the T rot of the gas exceeds ∼ 30 K.
To obtain the kinetic temperature, T K , from T rot , we employ the empirical approximation provided by Ott et al. (2011), T K = 6.06e 0.061(T rot /K) . (6) The temperature mentioned hereafter refers to the kinetic temperature.Tafalla et al. (2004) and Swift et al. (2005) also derived T rot -T K relations.Comparing these results, we find that the empirical approximation from Ott et al. (2011) applies for a wider range of T rot (T rot ≲ 500 K), while the ranges are 5 K < T rot < 20 K and Article number, page 3 of 37 A&A proofs: manuscript no.VLA_NH3   20:38:36.4 42:37:34.8 3.3×2.6; -29.0 Notes. (a) MDCs from Cao et al. (2021) covered by the field. (b) MDCs in Motte et al. (2007) covered by the field, extracted from 1.2 mm continuum data from IRAM 30-m telescope (with effective angular resolution of ∼ 11 ′′ and spatial resolution of ∼ 0.09 pc). (c) The beams are shown as "major axes × minor axes; position angles". (d) Dense core 774 and 742 have mass of 25.0 and 19.4 M ⊙ , slightly lower than 30 M ⊙ .However, both of them have remarkable NH 3 detections in our observations.We therefore took them into the following analysis and do not distinguish them from other MDCs. (e) The mass of N58 could be highly overestimated in Motte et al. (2007) due to the significant contribution of free-free emission in this core.N58 is therefore excluded from the sample in Wang et al. (2022).However, MDC 532 is included in this work as the SED-fitted mass of MDC 532 is 68.6 M ⊙ , which is still larger than 30 M ⊙ . (f) NW12 is an interesting isolated core in its parental clump.It is covered by VLA 14A-107 though its mass is lower than 30 M ⊙ .The mass of MDC 839 extracted by Cao et al. (2021) is also lower than 30 M ⊙ . (g) S26, also known as AFGL 2591, is excluded from our sample due to its different distance from Cygnux X region, of namely 3.3 kpc (Rygl et al. 2012).
T K ≪ 41 K for the relations in Tafalla et al. (2004) and Swift et al. (2005), respectively.In our results, a remarkable number of detections are hotter than 20 K, which is why we choose the relation from Ott et al. (2011).

Intrinsic line width and velocity dispersion
The line width derived from our fitting is blended by several hyperfine lines, because we do not conduct the full-hyperfine fitting.To derive the intrinsic line width that truly reflects the gas motion, we firstly subtract the channel width (0.1 km s −1 ) from the FWHM line widths derived from the line fitting quadrati-cally.We then apply the following empirical relation of Barranco & Goodman (1998) to eliminate the effect of hyperfineline blending: where C n (τ 0 ) is given as Here, τ 0 is the total optical depth of NH 3 (1,1) emission.We take the optical depth of the NH 3 (1,1) main line, τ 11,m , as the approximation of τ 0 .The conversion coefficient, K nm , is listed in  (a) √ for remarkable detection, * for marginal detection, and -for nondetection.× shows the lack of related cubes as we only obtain NH 3 (1,1) and (2,2) cubes for Fields 1, 2, and 37. (b) Nondetection cases are shown as - (c) The number of H 2 O maser spots.
3. As a polyfit result, this relation has its limitations for large FWHM line widths, which will make ∆v int larger than ∆v blended .In this case, we use ∆v blended as an approximation of ∆v int ; that is, ∆v int ≈ ∆v blended .Furthermore, we also tested the full-hyperfine fitting in Fields 3 and 5 and compared the fitted line width with that derived from our fitting results.In Field 3, the mean and median differences between the line widths from these two fitting methods are 0.034 and 0.027 km s −1 (5.4% and 6.7% for relative deviations, respectively).Regarding Field 5, these values are 0.062 (2.4%) and 0.008 (1.5%) km s −1 , respectively.The empirical relation is therefore reliable and the differences are totally negligible.
The intrinsic line width can be further decomposed into , where σ v t,NH 3 is the thermal velocity dispersion of NH 3 and σ v,NT is the nonthermal component.We can calculate the thermal velocity dispersion via σ v t,m = 0.288 km s −1 •(m/m H ) −1/2 •(T K /10K) 1/2 .Thus, we can derive the nonthermal velocity dispersion σ v,NT .Taking the molec-  (Kauffmann et al. 2008), we can calculate the thermal velocity dispersion of the whole molecular gas, σ v t,ISM .
We then derive the total velocity dispersion for the whole molecular gas from Article number, page 5 of 37 A&A proofs: manuscript no.VLA_NH3 We refer to σ v when we discuss velocity dispersion in the following context.

Uncertainties of the results
The uncertainties of the derived physical parameters could originate from both the observational process and the data reduction, including the missing flux, the calibration errors, the systematic errors of the calculating algorithms we use, and so on.The uncertainty of kinetic temperatures comes from two processes: the calculation of T rot and the transformation from T rot to T K .Here, T rot is mainly determined by the ratio of the main components of the NH 3 (1,1) and (2,2) inversion lines.In this process, the uncertainty of T rot is ∼ 10.3/S/N K, where S/N is the signal-tonoise ratio of the NH 3 (2,2) inversion line (Li et al. 2003).The S/Ns of our results range from ∼ 2 to ≳ 15, corresponding to the rms errors from ∼ 5.2 to ≲ 0.69 K.For the transformation process, we expect an uncertainty of ∼ 4 K according to Ott et al. (2011).
Regarding the velocity dispersion, there are three physical parameters involved in its calculation: the blended line width, the optical depth, and the kinetic temperature.The optical depth has little effect on the calculation of intrinsic line width (Barranco & Goodman 1998).According to Lu et al. (2014), the uncertainty of the blended line width can be written as where F is the flux intensity.δF/F is estimated to be ∼ 30% in our data, indicating an uncertainty on the blended line width of ∼ 15%.The empirical relation between ∆v int and ∆v blended would increase this uncertainty to about 30%.Coupled with the uncertainty of kinetic temperatures, the uncertainty of the velocity dispersion could be ≳ 30%.

Moment-0 paps and fragment extraction
In order to show the distributions of the spectral lines, we produce the moment-0 maps of all fields by integrating their intensities in the velocity ranges of the NH 3 (1,1) main components.
The NH 3 (1,1) and (2,2) moment-0 maps of Field 3 are shown as an example in Fig. 2. The NH 3 (1,1) and (2,2) moment-0 maps are plotted from 3σ and 2σ levels, respectively.Moment-0 maps of other fields and the other emission lines are presented in Appendix A.
As mentioned in Sect. 1, the NH 3 (1,1) moment-0 maps could roughly illustrate the distribution of cold dense gas in our MDCs.We extracted 202 fragments from the NH 3 (1,1) moment-0 maps using the Starlink GAUSSCLUMPS algorithm (Berry et al. 2007;Currie et al. 2014).The field 21 was excluded from the extraction as it is not at the distance of Cygnus X.The results of the identified NH 3 fragments in Field 3 are also plotted in Fig. 2 with blue ellipses.In addition, the sizes of NH 3 fragments are defined as the deconvolved radius, following where a and b are the semi-major and semi-minor axes of NH 3 fragments.
We also show the positions of the condensations identified in the SMA 1.3 mm continuum data from Cao et al. (2021).There are 148 condensations overlapping with the fields in this paper.The majority of them (103 out of 148) are distributed in the identified NH 3 fragments, indicating that the NH 3 fragments can probe the density structures well in most cases.As shown in Fig. 3, the (1,1) main components of most NH 3 fragments have moderate optical depths, and the mean and median values are 2.2 and 1.9.Under such a moderate optical depth, the NH 3 (1,1) lines might show different morphologies from the density profiles of the dust continuum emission, which could account for the condensations uncorrelated to the NH 3 fragments.In conclusion, the majority of the NH 3 fragments can be considered as density structures.

Maps of the derived physical parameters
The line-fitting process of the NH 3 (1,1) and (2,2) emission lines provides the spatial distributions of the kinetic temperature, centroid velocity, velocity dispersion, optical depth of the NH 3 (1,1) main component, NH 3 column density, and the beam-filling fac-Article number, page 6 of 37 tor in all the fields, which allow us to investigate the temperature and dynamic conditions of the MDCs in Cygnus X and their fragments.We present maps of these derived parameters in Appendix B, and the maps of Field 3 in Fig. 4 as an example.Ranges of the derived physical parameters are listed by each field in Table 4. Using these maps, we obtain the average kinetic temperature, NH 3 column density, velocity dispersion, and its nonthermal component for all the NH 3 fragments.The contribution from the velocity gradients within the NH 3 fragments, which is usually considered as the bulk motion, is taken into account when calculating the velocity dispersion, σ v , and its nonthermal component, σ v,NT .We also fit the velocity gradients (following the method in e.g., Goodman et al. 1993;Wu et al. 2018) and deduct them from the nonthermal velocity dispersions, which are provided as σ v,ng .The detailed calculation can be found in Appendix E. A brief statistic analysis of these physical parameters is reported in Table 5 and the histograms of their distributions are shown in Fig. 6.The detailed physical parameters of the NH 3 fragments are available in our online materials at the CDS.
Statistically, the NH 3 fragments have a mean deconvolved radius of 0.02 pc thanks to the high spatial resolution of our observations.Recently, Zhang et al. (2019) obtained the same value for the typical radius of condensations from their PdBI 1.3 mm continuum observations.In the other two papers of CENSUS (Pan et al. 2023;Yang et al. 2023), the properties of condensations in the Cygnus X regions are explored in more detail using submillimeter continuum and spectral observations from SMA.
The temperatures of NH 3 fragments range from ∼ 11.0 K to 34.1 K, with mean and median values of 18.6 K and 17.5 K, respectively.The mean temperature of the condensations in Zhang et al. (2019), which are selected from several massive protocluster clumps, is 28.3 K, which is higher than the value found in the present study.Lu et al. (2014) derived a mean temperature of ∼ 18.3 K for NH 3 dense cores (with a typical size of ∼ 0.08 pc) in 62 high-mass star-forming regions, which is consistent with our results.Temperature is one of the commonly used indicators of the evolutionary stage of molecular structures.For example, Sánchez-Monge et al. ( 2013) calculated the mean tem-peratures of the quiescent starless cores and the protostellar cores (with typical sizes of ∼ 0.05 pc) as 18.8 K and 28.8 K, respectively.To verify this relation in our data, we matched the NH 3 fragments with the condensations from Cao et al. (2021) and the IR sources from the Cygnus X Archive catalog (the Spitzer Cygnus X Legacy Survey, Kraemer et al. 2010).The temperature distribution of the NH 3 fragments coincident with the continuum condensations (70 out of 202) is shown in Fig. 5, indicating that these fragments are more likely to be density structures.These condensation-related fragments are divided into two groups: those associated and not associated with the IR sources (44 and 26 in number).The mean temperature of the former group is 20.4±4.5 K, slightly higher than that of the latter group of 19.1±4.4K. Also, their median temperature values are 19.4K and 17.3 K.In general, those with IR sources are considered to be more evolved than those without IR sources.We performed a two-sample Kolmogorov-Smirnov test on the temperature distribution of these two groups and the p-value is calculated as 0.18.The p-value is higher than the threshold 0.05, indicating the two groups are not distinguishable.As shown in Fig. 5, there are also considerable overlaps in temperatures between the two groups, suggesting that there is no significant difference in internal heating for the sources in the two groups.Instead, we find that the morphology of the NH 3 fragments is more correlated with the evolutionary stage; we discuss this topic in Sect.5.1.
The mean values of σ v , σ v,NT and σ v,ng are 0.46±0.20 km s −1 , 0.37±0.23 km s −1 , and 0.34±0.21km s −1 , respectively.The typical velocity dispersion of the condensations in Zhang et al. ( 2019) is ∼ 0.9 km s −1 .Considering that their temperature is also higher than that of our NH3 fragments, the larger velocity dispersion could also result from them being a more evolved sample.In contrast, our results are comparable to those of quiescent starless cores (0.

Correlation analysis
In this section, we present analyses of the correlations between physical parameters of the NH 3 fragments.We find one possible correlation between the nonthermal velocity dispersion and the kinetic temperature in our results, which are presented in Fig. 7.
Figures 7 (a) and (b) show σ v,NT and σ v,ng versus kinetic temperature.Though the distributions show large scatters, both of them show weak positive correlations: σ v,NT = 0.012T 1.13±0.17K with a correlation coefficient of r = 0.43, and σ v,ng = 0.007T 1.28±0.17K with r = 0.47.Similar correlations have been reported in previous NH 3 and H 2 CO observations, such as T K ∝ ∆v 0.53±0.14 in Wu et al. (2006), T rot ∝ ∆v 0.65 in Sánchez-Monge et al. (2013), ∆v ∝ T rot 1.16±0.29 in Lu et al. (2014), and in Tang et al. (2018).This relation is interpreted by Wu et al. (2006) and Tang et al. (2018) as the turbulent heating effect, which converts the turbulent energy into the internal energy of gas (Guesten et al. 1985).In particular, Tang et al. (2018) obtain a positive correlation between temperature and nonthermal velocity dispersion with the APEX H 2 CO observations.Using the equations of turbulent temperature from Ao et al. (2013), Tang et al. (2018) suggest that the turbulent heating could play Article number, page 7 of 37   Notes. (a) The nonthermal velocity dispersion (b) Also known as AFGL 2591 (c) Derived from the upper limit of T rot , 50 K. (d) Field 37 is fitted with two velocity components an important role in widely existing massive star formation regions on scales of ∼ 0.06-2 pc.However, such an effect usually dominates in the central molecular zone (Ginsburg et al. 2016;Immer et al. 2016) or other regions with strong nonthermal motions, which is not the case in our sample (e.g., ∼ 0.6-4 km s −1 in Tang et al. 2018, much larger than our result of ∼ 0.4 km s −1 ).Using the same method as Ao et al. (2013), we find a turbulent temperature of only ∼ 8 K, which is comparable to the background temperature (∼ 10 K) of cold molecular gas.Therefore, the turbulent heating does not seem to be remarkable in our results.Lu et al. (2014) linked temperature with velocity dispersion using the black-body radiation law (L ∼ T 4 ), the massluminosity relation (L ∼ M α ), and the virial equilibrium approximation (M ∼ σ 2 ), which gives an index varying from 0.67 to 2. Nevertheless, the assumptions of black-body radiation and virial equilibrium are questionable in molecular cores, and the massluminosity relation is very uncertain.In our results, the relation is considered as a natural consequence of the star-forming activity.An NH 3 fragment with stronger star-forming activity tends to be hotter and to have stronger infall and shearing motions that contribute to its nonthermal velocity dispersion.Conversely, an NH 3 fragment with weaker star-forming activity would be cooler and would have a smaller nonthermal velocity dispersion.This scenario is also consistent with the fact that our derived temperatures and velocity dispersions are higher at the positions of water masers and radio sources from Wang et al. (2022), where strong star-forming activities reside.
We also plot σ v,NT and σ v,ng versus the radius of fragments in panels (c) and (d).Such a relation is known as one of Larson's laws (Larson 1981), and is found in a number of works (e.g., Jijina et al. 1999;Wu et al. 2006;Lu et al. 2014).However, this relation is not found in our results.The main reason for the absence of these correlations could be the narrow range of the radius in our sample, which is within an order of magnitude.On the other hand, our mean value of σ v,ng 0.34 km s −1 coincides with the prediction of Larson's laws of 0.32 km s −1 at the scale of 0.02 pc.

Comparison between ammonia observations and H 2 column density
The NH 3 (1,1) emissions in our data show diverse morphologies in the moment-0 maps, such as filaments (Field 28), core-like structures (Field 6), hub-like structures (Field 13), and dispersed features (Field 17).The morphology of the NH 3 (1,1) emission can be affected by protostellar feedback during core evolution.
During the star formation process, increasing feedback from protostars can distort their envelopes (Motte et al. 2018) and heat the surrounding gas.The NH 3 (1,1) line gets much weaker and becomes undetectable when the gas temperature rises beyond 100 K (Ho & Townes 1983).Therefore, the morphology of the NH 3 in MDCs could evolve from a concentrated distribution at the peak of the H 2 density profile (hereafter, the M1 stage) to a dispersed shape around the protostar (hereafter, M2), and finally become undetectable (hereafter, M3).In order to construct an evolutionary sequence of the MDCs in Cygnus X, we compared the morphology of the NH 3 (1,1) emission to that of the H 2 column density (from the column density map, hereafter N-map, derived from spectral energy distribution (SED) fitting in Cao et al. 2019) of the MDCs, and classified the MDCs with infrared emissions.We present our findings in the following subsection.
As shown in Fig. 8 and Appendix C, the NH 3 (1,1) emission is strongly correlated with dense regions of the N-map in most fields, indicating that the NH 3 (1,1) line emission generally traces the dense structures in molecular gas.Nevertheless, the difference in resolution (∼ 3 ′′ .1 for NH 3 data and ∼ 20 ′′ for the N-map) leads to a more complicated morphology of NH 3 (1,1) emission than that of the N-map for our MDCs.We also find that many of the peaks of the NH 3 (1,1) emission show slight offsets from the peaks in the N-map, which could result from the depletion of NH 3 around protostars (Tobin et al. 2011;Sewiło et al. 2017) and from the difference in resolution between the N-map and the NH 3 cubes.
In light of their morphologies, we find that most (27 out of 35) MDCs are at the M1 stage.This result is consistent with our expectations, because the majority of our MDCs are not very evolved.As for the MDCs at the M2 stage (8 out of 35), we find evident star-forming activity (radio continuum emission, water masers, or bright sources on Spitzer images) in all of them, which corroborates the mentioned mechanism for distorted NH 3 (1,1) morphology.The MDCs at the M3 stage are not included in Table 6 because of their nondetection with the NH 3 (1,1) line.Figures 8 and 9 present some example MDCs in each stage.

A possible evolutionary sequence in our sample
The evolution of MDCs is a fundamental piece in our understanding of high-mass star formation.There have been many studies focusing on this topic, both observational and theoretical (e.g., Chambers et al. 2009;Battersby et al. 2014b;Molinari et al. 2019).Based on theoretical frameworks and observations, Article number, page 9 of 37 a simple four-stage evolutionary scenario of MDCs is proposed as CDMC→HDMC→DAMS→FIMS by Zinnecker & Yorke (2007), where CDMC, HDMC, DAMS and FIMS stand for cold dense massive core, hot dense massive core, disk-accreting main sequence star, and final main sequence star, respectively.Recently, Motte et al. (2018) updated the star-formation pictures in different scales: (1) at the scale of 1-10 pc, the gas is fed into clumps by the inflow on ridges and hubs; (2) at the scale of ∼ 0.1 pc, the beginning of star-forming activity will make starless MDCs grow into protostellar MDCs; (3) the stellar embryos in high-mass protostars grow from low-mass ones into high-mass ones, which transform the high-mass protostars from IR-quiet into IR-bright at ∼ 0.02 pc scales; and (4) the HII region phase at 0.01-10 pc scales is the end of this scenario.
Although the evolution of high-mass protostars cannot be properly determined based on their infrared SEDs -as is possible for the low-mass ones-due to their complicated environments (Andre et al. 1993;Marseille et al. 2010), the infrared properties can still serve as an indicator of the evolution of embedded high-mass protostars.In order to obtain a better estimation of evolutionary stage, we adopt the infrared classification of Wang et al. (2022) 5 .The MDCs are divided into three groups: 8 IR-bright MDCs, 25 IR-quiet MDCs, and 2 starless MDCs.The IR-bright MDCs correspond to MDCs with possible embedded high-mass protostar(s), while the IR-quiet ones only harbor lowmass protostars.Regarding the MDCs without any star-forming activity according to the criterion of Wang et al. (2022), we take them as starless MDC candidates.To illustrate the evolutionary stages of MDCs, we combine the data at hand, which include NH 3 emission data, the N-map, water masers, radio continuum sources, and mid-infrared pseudo three-color maps (8, 4.5, and 3.6 µm emission colored in red, green, and blue, respectively), into a composite image for each field (see Figures 8 and 9).
Figure 8 shows three MDCs at the M1 stage.The MDC 774 and 1179 are the only two starless MDCs in our sample.In the IR-quiet MDC 1225, we find two radio sources and two water masers, indicating star-forming activities and the existence of protostar(s).As the protostar or protostars grow, stronger feedback distorts the morphology of NH 3 , which could be the situation shown in Fig. 9.The MDC 725 and 1454 are considered to be at stage M2.The MDC 725 is an IR-bright MDC with a radio source; that is, the high-mass stellar embryo or embryos have formed in it, while the MDC 1454 is an IR-quiet MDC with NH 3 (1,1) emission distorted near the density peak.The MDC 1201 and 675 illustrate the latest stages in our sample, where most of the NH 3 is destroyed and NH 3 (1,1) emission is almost undetected within them.The MDC 1201 and 675 are IR-bright and IR-quiet, respectively.In conclusion, the morphology of the NH 3 (1,1) emission could be a possible indicator for the evolutionary stage of MDCs.

Virial analysis
The virial parameter, which is the ratio between the kinetic energy and the gravitational energy, is an important indicator of the dynamical conditions of molecular structures, especially for those in early evolutionary stages.Typically, the turbulent core accretion model requires a virialized MDC as the birth place of the high-mass star (McKee & Tan 2002).Previous studies in which a virial analysis was carried out using different gas tracers at different size scales found that the distribution of virial parameters varies greatly with the sizes of molecular structures (e.g., Larson 1981;Heyer et al. 2001;Chen et al. 2019).The virial parameters of molecular clouds (∼ 10 pc) are mostly ≫ 2, which is the critical value of an isothermal sphere in hydrostatic equilibrium embedded in a pressurized medium (Bonnor 1956;Ebert 1955;Larson 1981), while virial (1 ≤ α vir ≤ 2) and subvirial structures (α vir < 1) are found in the fragments (∼ 0.01-1 pc) within molecular clouds (e.g., Pillai et al. 2011;Li et al. 2013).
In this section, we perform a virial analysis of MDCs and try to understand their stability.

Calculation of virial parameters
As defined in Bertoldi & McKee (1992), the virial parameter, α vir , is in the form of where σ v is the velocity dispersion and R is the radius of the structure.In this work, the total velocity dispersion of the whole molecular gas, σ v , is obtained from the NH 3 results, and the mass of MDCs is provided by Cao et al. (2021).We therefore calculate the virial parameters for MDCs that have both NH 3 (1,1) and Article number, page 11 of 37 (2,2) detections.It is worth noting that we include the contributions from both turbulence and kinematics in the calculation of σ v .Detailed calculation processes are listed in Appendix E.
In many MDCs, the detected NH3 emission only covers part of the spatial extent of MDCs (such as the MDC 725 in Fig. 9).For these cases, we take the sizes of their detected parts, R sel , and calculate the corresponding mass, M sel , by multiplying their areas by the column densities taken from the N-map.Using R sel and M sel , the virial parameters are calculated for the detected parts of the MDCs, which are called "partial" virial parameters in our analysis.The partial virial parameters are considered to be an upper limit value for their host MDCs.Such overestimation comes from two aspects, R/M and σ v .First, we compare R sel /M sel of the selected parts and R total /M total of their host MDCs. Figure 10 shows that all the ratios between R sel /M sel and R total /M total are larger than 1, suggesting R/M is overestimated in the partial virial parameters.Second, as shown in Appendices A and B, the velocity dispersion is usually enhanced where the integral NH 3 (1,1) emission is strong in most sources.The undetected parts should therefore have lower velocity dispersion than the detected parts, indicating the velocity dispersion is also overestimated in the partial virial parameters.

Distribution of virial parameters
Through the above calculations, we determined the virial parameters for 35 MDCs, including 14 full ones and 21 partial ones, as listed in Table 6.The mean and median values of the virial parameter are 1.0 and 0.9, respectively, with a standard deviation of 1.0.The majority of MDCs (23 out of 35) have a virial parameter of smaller than 1.Moreover, most full virial parameters (11 out of 14) are smaller than 1, with only one exceeding 2. Furthermore, the partial virial parameters, mostly distributed around 1, are considered as upper limits in our analysis, indicating that the majority of MDCs could be subvirialized.This result is consistent with previous findings that many MDCs are moderately to strongly subvirial (e.g., Pillai et al. 2011;Wienen et al. 2012;Li et al. 2013;Lu et al. 2014;Keown et al. 2017;Kerr et al. 2019;Traficante et al. 2020;Zhang et al. 2015;Ohashi et al. 2016;Li et al. 2020).Recently, Singh et al. (2021) proposed three possible sources of bias that could result in a false subvirial result when calculating virial parameters: neglected bulk motions, an overestimated density profile correction factor, and mass contribution from the fore-and/or background.We tried to exclude these sources of bias in our calculation: (1) The bulk motions are considered in the total kinetic energy of MDCs; (2) the density profile correction factor is taken as 1 as defined in Singh et al. (2021), which could result in an overestimated virial parameter rather than an underestimated one; and (3) the mass of our MDCs is derived from the SED fitting results of Cao et al. (2019), where the fore-and background are removed in the fitting processes.Our subvirial results should therefore not be affected by these sources of bias and we consider them to be trustworthy.
As a result, our MDCs appear to be mostly subvirialized, which is not expected based on the turbulence core accretion model (McKee & Tan 2002).This suggests that MDCs could undergo dynamical collapse in nonequilibrium in the early stages of star formation unless there is additional support.Regarding high-mass star formation, the subvirial prestellar cores, according to the simulation work by Rosen et al. (2019), are more likely to form massive stars than virialized ones due to their higher accretion rates and lower levels of fragmentation.Additionally, an anti-correlation between the virial parameter and mass has been found in many works using only a single gas tracer (e.g., Lada Article number, page 12 of 37 As mentioned above, the subvirial MDCs would collapse if there were no additional support.Previous studies suggested that the magnetic field could provide significant support in MDCs (e.g., Frau et al. 2014;Hull & Zhang 2019;Liu et al. 2020Liu et al. , 2022)).Following the methods in Liu et al. (2020), the total virial mass -which includes the contribution of magnetic field-can be calculated as where M k = α vir M is the kinetic virial mass, and M B here is We then calculated the total virial parameter, α k+B = M k+B /M, by assuming a typical magnetic field strength of ∼ 0.5 mG (Hezareh et al. 2014).The results are shown in Fig. 12 for comparison with the original results in Fig. 11.
Article number, page 13 of 37  11.Virial parameters vs. mass of the 35 MDCs in Cygnus X.In the legends, "B", "Q", and "S" refer to different infrared classes, IR-bright, IR-quiet, and starless, respectively."F" and "N" are the same as those in Table 6.The red line shows the least-linear-square fitting result.The gray dashed lines show where virial parameters are equal to 1 and 2. The mean and median values of the total virial parameter are 2.2 and 2.1, with a standard deviation of 1.2.The majority of total virial parameters are approximately 2 (as shown in Kauffmann et al. 2013, the critical virial parameter is ∼ 2 for Bonnor-Ebert spheres), indicating that most MDCs are expected to be in stability with a typical magnetic field strength of 0.5 mG.In addition, if there is a strong magnetic field strength up to 1 mG, the total virial parameter would increase by ∼ 1-2 in value.In addition, Liu et al. (2022) find an anti-correlation between the magnetic virial parameter (α B = M B /M) and the column density of H 2 .Considering the typical column density of our MDCs, of namely ∼ 10 23 cm −3 , the magnetic virial parameters are ∼ 0.5 in Liu et al. (2022), indicating that the magnetic field could provide remarkable support against the gravity.In our results, the nonthermal velocity dispersion of MDCs is lower than that predicted by McKee & Tan (2002), ∼ 1.65(m/30M ⊙ ) 1/4 Σ 1/4 km s −1 .Subsequent studies using this model (Krumholz & McKee 2005, e.g.,) also proposed significant support from the magnetic field against the gravity at the scale of MDCs.The Cygnus X North and South regions are also covered by the lower resolution NH 3 survey (KEYSTONE) presented in Keown et al. (2019).The virial parameter of the NH 3 cores presented by these authors (≳ 0.1 pc, referring to the NH 3 "leaves" in that paper) show higher values than our results.A significant part of the NH 3 cores are at supervirial states in Keown et al. (2019).Therefore, the parent structures of our NH 3 fragments have larger virial parameters, indicating that the accretion of "cores" may not be dominated by the gravity.The gas might be fed to the "core" via the filament or the initial inflow suggested by Motte et al. (2018) and Padoan et al. (2020).

Dynamical states in star-forming regions
We produced nonthermal velocity dispersion maps of each field and mark the regions where the nonthermal velocity dispersion is higher than the sound speed with green contours (Fig. 13).We find the dynamical states vary among different MDCs.As shown in Fig. 13, the MDC 248 is dominated by the trans-sonic to supersonic nonthermal velocity dispersion (M NT ≥ 1), while the nonthermal velocity dispersion in MDC 892 is almost entirely subsonic.The majority of MDCs consist of both supersonic and subsonic regions, such as MDCs 274, 584, and 774.Maps of other fields are reported in Appendix D.
In Sect.4.2, we provide the σ v,NT and σ v,ng of NH 3 fragments, where the resolved bulk motions are included and excluded, respectively.We also calculate σ v,NT and σ v,ng for MDCs.Then, using σ v,NT and σ v,ng , we obtain the Mach numbers, M NT and M ng , for both NH 3 fragments and MDCs (Fig. 14).For MDCs, M NT has mean and median values of 2.2 and 2.3, with a standard deviation of 0.8, while the mean and median values of M ng are 1.9 and 1.9, with a standard deviation of 0.7.As for NH 3 fragments, the majority of both M NT and M ng are smaller than 2 (trans-sonic to subsonic).The mean and median values of M NT are 1.4 and 1.2, with a standard deviation of 0.8, and the mean and median values of M ng are 1.3 and 1.1, with a standard deviation of 0.7.These distributions of NH 3 fragments are similar to previous results (e.g., Sokolov et al. 2018;Yue et al. 2021;Li et al. 2020) The differences between M NT and M ng of MDCs are more remarkable than those of NH 3 fragments, suggesting that the bulk motions have a larger contribution to the nonthermal velocity dispersions of MDCs than those of NH 3 fragments.In contrast, the M ng and M NT of NH 3 fragments are close.Moreover, both panels show a decreasing trend in Mach number from MDCs to NH 3 fragments.
Recently, high-spatial-resolution observations revealed similar transonic and subsonic nonthermal velocity dispersion in both low-mass and high-mass star-forming regions, such as L1517 (Hacar & Tafalla 2011), NGC 1333 (Hacar et al. 2017), the Orion integral filament (Hacar et al. 2018), OMC 2/3 (Yue et al. 2021), and NGC 6334S (Li et al. 2020).Comparing N 2 H + J=1-0 observations from the Nobeyama 45m telescope and ALMA, Yue et al. (2021) propose the nonthermal motions would transit from supersonic motions to subsonic motions at a scale of ∼ 0.05 pc.In our results, the transition occurs at a scale of NH 3 fragments, ∼ 0.04 pc, which is consistent with the result of Yue et al. (2021).These latter authors interpret this transition as a resolution-dependent bias, where more bulk motions are involved in larger beams and give a greater nonthermal velocity dispersion.The large-scale dispersion measured from Nobeyama observations, σ ls , is decomposed in Yue et al. (2021) into smallscale dispersion measured from ALMA observations, σ ss , bulk motion probed by dispersion between the peak velocity of each ALMA beam, σ bm , and residual dispersion, σ rd , with the relation σ 2 ls =σ 2 ss +σ 2 bm +σ 2 rd .Comparing M ng of MDCs with that of NH 3 fragments as shown in the panel (b) of Fig. 14, where the bulk motions are deducted, the decreasing trend of M ng from MDCs to NH 3 fragments still remains.The differences between Article number, page 15 of 37 The ellipses are the same as those in Fig. 2 the M ng of MDCs and that of NH 3 fragments could refer to σ rd .A possible explanation for σ rd is turbulent vortex, as put forward by Yue et al. (2021).In most cases, the decrease in nonthermal velocity dispersion from large scale to small scale is also thought to be the dissipation of turbulence, as suggested by the relation between velocity dispersion and size at lower spatial resolutions (∼ 0.1 -100 pc) in Larson (1981).Considering that the typical size of NH 3 fragments (∼ 0.04 pc) is larger than the typical Kolmogorov scale in molecular clouds (∼ 10 −5 -10 −4 pc) (McKee & Ostriker 2007;Kritsuk et al. 2011;Qian et al. 2018), the turbulent dissipation could be a possible cause for the decreasing trend of M ng .Furthermore, the lower Mach numbers imply that magnetic fields can play a more important role if the fragments are in a virialized state.However, as shown in Liu et al. (2022), α B would decrease with increasing column density.Considering that α B is about 0.5 in our MDCs, the gravity could be more dominant than the magnetic field in these fragments (usually the peaks on the column density map).As estimated by Sokolov et al. (2018), magnetic fields with strengths over ∼ 1 mG are needed to achieve virialization in these fragments.However, such a strong magnetic fields have only been detected in a small number of starforming structures (e.g., Curran & Chrysostomou 2007;Fish & Sjouwerman 2007;Palau et al. 2021).Therefore, the fragments of MDCs are probably not in virialized states.

Gas and dust temperature
In general, there are two ways to trace the temperature of molecular clouds: from the dust continuum emission or from spectral line emission of a certain molecule.As the temperatures of MDCs are also provided in Cao et al. (2021) (T dust column in Table 6), which were calculated from the SED fitting of the dust continuum data, we compared these temperatures with the temperatures derived from the NH 3 data; we present our findings below.
In order to compare the NH 3 temperature and the dust temperature in MDCs, the NH 3 data cubes are convolved to the resolution of temperature map in Cao et al. (2019) and the NH 3 temperatures of MDCs are rederived (T gas column in Table 6).The correlation between them is plotted in Fig. 15.The gray line illustrates where T gas = T dust and the red line presents the linear regression fitting result, with the correlation coefficient r = 0.66.According to this relation, the gas temperature is mostly higher than the dust temperature when T gas ≳ 20 K.In the scatter plot, the turning point is also around 20 K. We find that T dust is almost equal to T gas when T gas ≲ 20 K and deviates from T gas when T gas ≳ 20 K.
The gas and dust are expected to be well coupled -that is, they can be characterized by the same temperature-within a molecular core with a density of ≳ 10 4.5 cm −3 (Goldsmith 2001;Young et al. 2004).This coupling has been confirmed by many works (e.g., Dunham et al. 2010;Merello et al. 2019).However, although all the MDCs in our sample have densities higher than 10 5 cm −3 , T gas and T dust are not thermally well coupled.Similarly, in high-mass star-forming regions, T gas is usually measured to be higher than T dust , which is the case for example for the quiescent clumps in Battersby et al. (2014a), the S 140 in Koumpia et al. (2015), and the W40 in Tursun et al. (2020).
The origin of the excess gas temperature is a complex open question.In the photon-dominated region S 140, Koumpia et al. (2015) found an excess of ∼ 5-15 K in their deeply embedded regions, where T dust and T gas are expected to be well coupled in the KOSMA-τ PDR model (Röllig et al. 2013).Koumpia et al. (2015) interpret these differences as the result of inhomogeneous clumpy structures, which could lead to a deeper penetration of the UV radiation and heating of the gas in the embedded cores.In addition, differences between T gas and T dust can also arise from different layers traced by the gas and the dust (Battersby et al. 2014a).Such differences arise from both different excitation conditions and different angular resolutions.The NH 3 (1,1) and (2,2) lines are usually excited in the dense gas, while the dust continuum emission arises from a range of densities.In the present work, the resolutions of the NH 3 data and the temperature map from Cao et al. (2019) are quite different.Only a few MDCs have sufficient temperature pixels as shown in Fig. 4 and  Appendix B. The spatially extended NH 3 emission, which usually has lower temperatures, is missed due to the filtering effect of VLA.In our results, T gas deviates from T dust at T gas ≳ 20 K (Fig. 15).The MDCs with T gas ≳ 20 K should have internal heating from the star-forming activity, but this heating is diluted by the large beam in the dust temperature data.Therefore, on average, we obtain a higher T gas than T dust in our results.

Summary
The initial conditions of the birth places of high-mass stars are of utmost importance in the study of high-mass star formation.We present VLA NH 3 (1,1) and ( 2 Cygnus X at a distance of 1.4 kpc.As a part of the CENSUS project, we provide the temperature and dynamical information derived from the simultaneous line fitting of the NH 3 (1,1) and (2,2) inversion emission lines.We used the MDCs catalog from Cao et al. (2021) in our analysis.We compared the morphology between the H 2 column density map from Cao et al. (2019) and NH 3 (1,1) emission within each MDC and provide a possible indicator for estimating their evolutionary stages.We also performed a virial analysis for each MDC to evaluate the stability of the MDCs.The dynamical states of the MDCs are discussed based on their nonthermal velocity dispersions.We also compare the gas temperature with the dust temperature of MDCs.In addition, we identified 202 fragments from the NH 3 (1,1) moment-0 maps using the GAUSSCLUMPS algorithm.We derived the physical parameters of these NH 3 fragments and conducted a correlation analysis among them.We also analyzed the dynamical states of the NH 3 fragments and compared them with those of the MDCs.The main results of this paper can be summarized as follows: 1.The NH 3 fragments have a mean deconvolved radius of 0.02 pc, a mean velocity dispersion of 0.46 km s −1 , a mean nonthermal velocity dispersion of 0.37 km s −1 , a mean kinetic temperature of 18.6 K, and a mean NH 3 column density of 3.9 × 10 15 cm −2 .Comparison with the temperature of structures in previous works indicates that the NH 3 fragments are mostly early evolved structures.2. We find a correlation between nonthermal velocity dispersion and kinetic temperature.This correlation is likely the result of different levels of star-forming activity.A molecular structure with more active star-forming activity is usually hotter and with larger velocity dispersion.3. We perform a morphological analysis of the NH 3 (1,1) emission and column density map of MDCs.We propose a possible evolutionary sequence of MDCs based on the NH 3 (1,1) emission: in the early, relatively little-evolved MDCs, the NH 3 emission is concentrated around density peaks; the NH 3 emission is subsequently dispersed by the star-forming activity, after which, in the evolved MDC, the NH 3 emission becomes undetected.4. We perform a virial analysis of all MDCs.The majority of the MDCs of our sample are subvirialized, indicating that the kinetic energy cannot support the MDCs alone in our sample.Considering the support from magnetic fields, MDCs could attain their stability states if there were a typical magnetic field strength of ∼ 0.5 mG. 5. We investigate the dynamic states of MDCs and the NH 3 fragments with their nonthermal Mach numbers.The Mach numbers show a decreasing trend from MDCs to NH 3 fragments, which could be an indication of the dissipation of turbulence.When the resolved bulk motions are excluded from the nonthermal velocity dispersion, we find the NH 3 fragments are mostly subsonic to trans-sonic, while the majority of MDCs are transonic.6.The gas temperature derived from NH 3 (1,1) and (2,2) lines is compared with the dust temperature calculated by SED fitting from Cao et al. (2019).The two are well coupled when the gas temperature is ≲ 20 K, while the gas temperature is Article number, page 17 of 37 mostly larger than the dust temperature when the gas temperature is ≳ 20 K.

Fig. 2 .
Fig. 2. (a) Moment-0 map of the NH 3 (1,1) line for Field 3. The contour is plotted from 3σ level and in steps of 2σ.The red ellipses show the MDCs identified from the H 2 column density map (resolution ∼ 20 ′′ ) by Cao et al. (2021), and the blue ellipses illustrate the fragments extracted using the GAUSSCLUMPS algorithm.The magenta boxes show the condensations identified from the SMA 1.3 mm continuum.(b) Moment-0 map of the NH 3 (2,2) line for Field 3. The contour is plotted from 2σ level and in steps of 2σ

Fig. 3 .
Fig. 3. Histogram of the average optical depth of the main component of the NH 3 (1,1) line for identified NH 3 fragments.

Fig. 4 .
Fig. 4. Physical parameter maps of Field 3. (a) Centroid velocity map.(b) Velocity dispersion map.(c) Kinetic temperature map.The ellipses are the same as those in Fig. 2Table 4. Ranges of derived physical properties in each field

Fig. 6 .
Fig. 6.Distributions of the physical parameters of NH 3 fragments.(a) Histogram of kinetic temperatures T K .(b) Histogram of total velocity dispersions σ v (including bulk motions and using the thermal component of the entirety of the molecular gas).(c) Histogram of deconvolved radii of NH 3 fragments R dec .(d) Histogram of NH 3 column density N(NH 3 ).(e) Histogram of nonthermal velocity dispersion with bulk motions σ v,NT .(f) Histogram of nonthermal velocity dispersion without bulk motions σ v,ng .

Fig. 7 .
Fig. 7. Correlations between physical parameters of NH 3 fragments.Panels (a) and (b) plot nonthermal velocity dispersions including or excluding the bulk motions versus kinetic temperature.The red and gray lines show the fitted result of the linear regression and the thermal velocity dispersion.There are positive correlations in both of them.Panels (c) and (d) plot nonthermal velocity dispersions including or excluding the bulk motions versus deconvolved radius.The gray line shows Larson's law.

Fig. 8 .
Fig. 8. MDCs at the M1 stage (i.e., the NH 3 (1,1) emission is roughly concentrated around the column density peaks).The background is pseudo three-color images of each field, with Spitzer 8, 4.5, and 3.6µm emission colored in red, green, and blue, respectively.The contours of H 2 column density and NH 3 (1,1) integrated emission are over-plotted with green and orange, respectively.Cyan and magenta crosses mark the positions of any H 2 O masers or radio continuum sources, respectively.The starless, IR-quiet, and IR-bright MDCs are shown as pink dashed, red dashed, and red solid ellipses, respectively.More figures are showed in Appendix C.

XFig. 9 .Fig. 10 .
Fig. 9. MDCs at the M2 stage (MDC 725 and 1454) and M3 stage (MDC 1201 and 675).Marks are the same as those in Fig. 8.For MDC 1201, there is nearly no NH 3 (1,1) emission in it.As for MDC 675, NH 3 (1,1) emission is only detected in a few isolated beams around the density peak.We therefore consider both of them to be MDCs at the M3 stage.

Fig. 12 .
Fig. 12.Total virial parameter v.s.mass.The legends are the same as those in Fig. 11.

Fig. 13 .
Fig. 13.Nonthermal velocity dispersion maps of Fields 3, 5, and 26.The green contours show the regions with Mach numbers of greater than 1.The ellipses are the same as those in Fig.2

Fig. 14 .Fig. 15 .
Fig. 14.Normalized density distribution of M NT and M ng in NH 3 fragments and MDCs.The contribution of velocity gradient is included in (a) and excluded in (b) . A.1.Moment-0 maps with remarkable detection of NH 3 (J,K) = (1,1), (2,2), (3,3), (4,4), and (5,5) inversion lines among our observations.The (1,1) maps are plotted from a 3 rms level and in steps of 2σ, while others are plotted from the 2 rms level and in steps of 2σ.The IRbright, IR-quiet, and starless MDCs are shown as red solid, red dashed, and pink dashed ellipses, respectively.The smaller ellipses present the fragments extracted using the GAUSSCLUMPS algorithm.The blue ones are included in the parameter analysis while the green ones are not.Appendix B: Distributions of kinetic temperature, centroid velocity and velocity dispersionFig.B.1.Distributions of kinetic temperature (left panel), centroid velocity (middle panel), and velocity dispersion (right panel) derived from the simultaneous fitting of NH 3 (1,1) and (2,2) inversion lines.The ellipses are the same as those in Fig. A.

Table 1 .
Observational fields and related MDCs

Table 2 .
Detections in each field Notes.

Table 3 .
Values of K nm

Table 5 .
Statistics of physical parameters of NH 3 fragments Article number, page 8 of 37 X.Zhang et al.: NH 3 observations of Cygnus X Fig. 5. Histogram of average kinetic temperatures of the NH 3 fragments coinciding with the condensations of 1.3mm continuum emission.Those matched with IR sources are colored in red, while the rest are colored in cyan.

Table 6 .
Properties of MDCs with NH 3 Detection MDC IR-class Morph a NH 3