Clumpiness in Star Forming Regions

In this article, some aspects of the clumpy nature of molecular clouds are reviewed. In particular the observational evidence for small-scale structures both in low and high mass star forming regions will be discussed. I will review some examples of `clumpiness' such as: i) the molecular clumps ahead of HH objects and how the study of the physical and chemical nature of these clumps is important for the understanding of the clumpiness of the Interstellar Medium; and ii)hot cores and their use as a tool to study the early phases of massive star formation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Assessing the performance of molecular gas clump identification algorithms

The detection of clumps(cores) in molecular clouds is an important issue in sub-millimetre astronomy. However, the completeness of the identification and the accuracy of the returned parameters of the automated clump identification algorithms are still not clear. In this work, we test the performance and bias of the GaussClumps, ClumpFind, FellWalker, Reinhold, and Dendrograms algorithms in identifying simulated clumps. By designing the simulated clumps with various sizes, peak brightness, and crowdedness, we investigate the characteristics of the algorithms and their performance. In the aspect of detection completeness, FellWalker, Dendrograms, and GaussClumps are the first, second, and third best algorithms, respectively. The numbers of correct identifications of the six algorithms gradually increase as the size and signal-to-noise ratio(SNRs) of the simulated clumps increase and they decrease as the crowdedness increases. In the aspect of the accuracy of retrieved parameters, FellWalker and Dendrograms exhibit better performance than the other algorithms. The average deviations in clump parameters for all algorithms gradually increase as the size and SNR of clumps increase. Most of the algorithms except FellWalker exhibit significant deviation in extracting the total flux of clumps. Taken together, FellWalker, GaussClumps,and Dendrograms exhibit the best performance in detection completeness and extracting parameters. The deviation in virial parameter for the six algorithms is relatively low. When applying the six algorithms to the clump identification for the Rosette molecular cloud, ClumpFind1994, ClumpFind2006, GaussClumps,FellWalker, and Reinhold exhibit performance that is consistent with the results from the simulated test.     

Effects of clumping on temperature – I. Externally heated clouds

We present a study of radiative transfer in dusty, clumpy star-forming regions. A series of self-consistent, 3D, continuum radiative transfer models are constructed for a grid of models parametrized by central luminosity, filling factor, clump radius and face-averaged optical depth. The temperature distribution within the clouds is studied as a function of this parametrization. Among our results, we find that: (i) the effective optical depth in clumpy regions is less than in equivalent homogeneous regions of the same average optical depth, leading to a deeper penetration of heating radiation in clumpy clouds, and temperatures higher by over 60 per cent; (ii) penetration of radiation is driven by the fraction of open sky (FOS) – which is a measure of the fraction of solid angle along which no clumps exist; (iii) FOS increases as clump radius increases and as filling factor decreases; (iv) for values of   FOS >0.6–0.8  the sky is sufficiently 'open' that the temperature distribution is relatively insensitive to FOS; (v) the physical process by which radiation penetrates is preferentially through streaming of radiation between clumps as opposed to diffusion through clumps; (vi) filling factor always dominates the determination of the temperature distribution for large optical depths, and for small clump radii at smaller optical depths; (vii) at lower face-averaged optical depths, the temperature distribution is most sensitive to filling factors of 1–10 per cent, in accordance with many observations; (viii) direct shadowing by clumps can be important for distances approximately one clump radius behind a clump.     

Effects of clumping on the observed properties of dusty galaxies

We present Monte Carlo radiative-transfer simulations for spiral galaxies modelled as a stellar disc and a two-phase clumpy dust distribution. We divide the volume occupied by the dust into a three-dimensional grid and assign each cell a clump or smooth medium status. Cell dimension, clump dust mass and spatial distribution are derived from the observed properties of giant molecular clouds and molecular gas in the Galaxy. We produce models for several values of the optical depth and fraction of the interstellar medium residing in clumps. As a general result, clumpy models are less opaque than the corresponding homogeneous models. For the adopted parameters, the increase in the fraction of energy that escapes the disc is moderate, resulting in surface-brightness profiles that are less than one magnitude brighter than those of the homogeneous models. The effects of clumping are larger for edge-on views of the disc. This is in contrast with previous preliminary results for clumping in the literature. We show how differences arise from the different parametrization and clump distribution adopted. We also consider models in which a fraction of the stellar radiation is emitted within the clumps. In this case, galaxies are less transparent than in the case when only dust is clumped. The opacity can be even higher than in the homogeneous case, depending on the fraction of embedded stellar emission. We point out the implications of the results for the determination of the opacity and dust mass of spiral galaxies.     

Disentangling Saturn's F Ring. I. Clump orbits and lifetimes

A comprehensive analysis of the Voyager images reveals the kinematics and lifetimes of clumps in the F Ring. At any given time, the ring has 2-3 major clumps, each several times brighter than the typical ring, plus numerous smaller features. A total of 34 individual clumps have been tracked over periods of 1-7 weeks. The clumps orbited at measurably different rates, implying a 100-km range of semimajor axes centered on 140,220 km. Most are centered around the nominal mean motion of the ring's core, but a few outliers may be associated with a different strand, or with no strand at all. Most clumps change very little over the ∼30 days that they can be detected; however, no clump persisted for the nine-month interval between the two Voyager encounters. The brightest Voyager 2 clump is unusual in that it travels at a rapid mean motion and seems to be associated with the formation of several other clumps.     

The migrating embryo model for disk evolution

Abstract– A new view of disk evolution is emerging from self‐consistent numerical simulation modeling of the formation of circumstellar disks from the direct collapse of prestellar cloud cores. This has implications for many aspects of star and planet formation, including the growth of dust and high‐temperature processing of materials. A defining result is that the early evolution of a disk is crucially affected by the continuing mass loading from the core envelope, and is driven into recurrent phases of gravitational instability. Nonlinear spiral arms formed during these episodes fragment to form gaseous clumps in the disk. These clumps generally migrate inward due to gravitational torques arising from their interaction with a trailing spiral arm. Occasionally, a clump can open up a gap in the disk and settle into a stable orbit, revealing a direct pathway to the formation of companion stars, brown dwarfs, or giant planets. At other times, when multiple clumps are present, a low mass clump may even be ejected from the system, providing a pathway to the formation of free‐floating brown dwarfs and giant planets in addition to low mass stars. Finally, it has been suggested that the inward migration of gaseous clumps can provide the proper conditions for the transport of high‐temperature processed solids from the outer disk to the inner disk, and even possibly accelerate the formation of terrestrial planets in the inner disk. All of these features arising from clump formation and migration can be tied together conceptually in a migrating embryo model for disk evolution that can complement the well‐known core accretion model for planet formation.     

Application of a clumpy core model to NGC 2071

The multi-transitional observations of CS molecules towards the NGC 2071 core have been re-analysed by using a tri-dimensional Monte Carlo radiative transfer code. Better agreement with the observations is made by an introduction of clumpiness to this model than by smoothly varying density to the 1D microturbulent one. The best-fitting model shows that, when a unique density is assumed for clumps, the volume filling factor of the clumps varies as r ⁻² with an average of ∼5 per cent over the entire core, and that the H₂ number density and the CS abundance of the clump relative to H₂ are ∼ 2 × 10⁶ cm⁻³ and ∼ 6 × 10⁻¹⁰, respectively. The radial density gradient ∝ r ⁻² obtained from our clumpy core model is steeper than that (∝ r ^−1.3) obtained from the microturbulent model. Since all clumps are subject to random bulk motions in this 3D clumpy macroturbulent model, synthesized line profiles do not show self-absorption dips even for opaque transitions and the resulting linewidth is in good accordance with the observations.     

Physical processes in the interstellar matter determining parameters of radio recombination lines in meter and decameter ranges

We give a review of problems connected with the interpretation of meter and decameter carbon radiolines. The lines are formed inside clumps of molecular clouds in layers with a column density N ≈ 6 · 10²¹ cm⁻². These clumps are very typical structures. The distribution of physical parameters (number density, temperature, etc.) inside the clumps is poorly known. The most difficult and important question is the penetration of subcosmic rays into the clumps. Observations show that the ionization rate is ζ = (1–7) · 10⁻¹⁷ s⁻¹ inside molecular clouds and significantly greater in the diffuse gas. Long-wave radio recombination lines can probably be used for the analysis of the distribution of subcosmic rays inside molecular clouds. The interpretation is complicated by the influence of low-temperature dielectron recombination and poorl known variations of carbon depletion in the clumps.     

Molecular line mapping of the giant molecular cloud associated with RCW 106 – II. Column density and dynamical state of the clumps

We present a fully sampled C¹⁸O (1–0) map towards the southern giant molecular cloud (GMC) associated with the H  ii region RCW 106, and use it in combination with previous ¹³CO (1–0) mapping to estimate the gas column density as a function of position and velocity. We find localized regions of significant ¹³CO optical depth in the northern part of the cloud, with several of the high-opacity clouds in this region likely associated with a limb-brightened shell around the H  ii region G333.6−0.2. Optical depth corrections broaden the distribution of column densities in the cloud, yielding a lognormal distribution as predicted by simulations of turbulence. Decomposing the ¹³CO and C¹⁸O data cubes into clumps, we find relatively weak correlations between size and linewidth, and a more sensitive dependence of luminosity on size than would be predicted by a constant average column density. The clump mass spectrum has a slope near −1.7, consistent with previous studies. The most massive clumps appear to have gravitational binding energies well in excess of virial equilibrium; we discuss possible explanations, which include magnetic support and neglect of time-varying surface terms in the virial theorem. Unlike molecular clouds as a whole, the clumps within the RCW 106 GMC, while elongated, appear to show random orientations with respect to the Galactic plane.     

The effect of magnetic field morphology on the structure of massive IRDC clumps

Infrared dark clouds (IRDCs) have dense elongated clumps and filaments with the favorable viewing condition of being on the near-side of a bright mid-infrared background. The clumps usually have multiple cores around the center. In this work, we study the effect of magnetic field morphology on the structure of massive IRDC clumps. To achieve this goal, we consider an axisymmetric isothermal oblate IRDC clump, embedded into a constant external magnetic field. We assume a polynomial function for the magnetic field morphology inside the clump. We use the numerical iterative methods to solve the equations: the successive over-relaxation method to find the magnetic and gravitational fluxes, and then the biconjugate gradient method to find the optimized values of mass and current densities. The results show that the IRDC clump will be very elongated along the perpendicular direction of the external magnetic field lines. Also, the assumption of choosing of a polynomial function for the magnetic field morphology leads to the formation of dense regions around the center. The greater the density of the central region, the larger the density of these dense regions and the closer to the center. The presence of these dense regions can lead to the formation of cores at these points.     

The dynamical evolution of substructure

The evolution of substructure embedded in non-dissipative dark haloes is studied through N -body simulations of isolated systems, both in and out of initial equilibrium, complementing cosmological simulations of the growth of structure. We determine by both analytic calculations and direct analysis of the N -body simulations the relative importance of various dynamical processes acting on the clumps, such as the removal of material by global tides, clump–clump heating, clump–clump merging and dynamical friction. The ratio of the internal clump velocity dispersion to that of the dark halo is an important parameter; as this ratio approaches a value of unity, heating by close encounters between clumps becomes less important, while the other dynamical processes continue to increase in importance. Our comparison between merging and disruption processes implies that spiral galaxies cannot be formed in a protosystem that contains a few large clumps, but can be formed through the accretion of many small clumps; elliptical galaxies form in a more clumpy environment than do spiral galaxies. Our results support the idea that the central cusp in the density profiles of dark haloes is the consequence of self-limiting merging of small, dense haloes. This implies that the collapse of a system of clumps/substructure is not sufficient to form a cD galaxy, with an extended envelope; plausibly, subsequent accretion of large galaxies is required. The post-collapse system is in general triaxial, with rounder systems resulting from fewer, but more massive, clumps. Persistent streams of material from disrupted clumps can be found in the outer regions of the final system, and at an overdensity of around 0.75, can cover 10 to 30 per cent of the sky.     

Fragmentation in molecular clouds and its connection to the IMF

We present an analysis of star-forming gas cores in a smooth particle hydrodynamics simulation of a giant molecular cloud. We identify cores using their deep potential wells. This yields a smoother distribution with clearer boundaries than density. Additionally, this gives an indication of future collapse, as bound potential cores (p-cores) represent the earliest stages of fragmentation in molecular clouds. We find that the mass function of the p-cores resembles the stellar initial mass function and the observed clump mass function, although p-core masses  (∼0.7 M_⊙)  are smaller than typical density clumps. The bound p-cores are generally subsonic, have internal substructure and are only quasi-spherical. We see no evidence of massive bound cores supported by turbulence. We trace the evolution of the p-cores forward in time, and investigate the connection between the original p-core mass and the stellar mass that formed from it. We find that there is a poor correlation, with considerable scatter suggesting accretion on to the core is dependent on more factors than just the initial core mass. During the accretion process the p-cores accrete from beyond the region first bound, highlighting the importance of the core environment to its subsequent evolution.     

A possible EGRET source with a two-sided radio jet structure on a parsec scale – 3EG J0808+5114

Recent observations have revealed that damped Lyα clouds (DLAs) host star formation activity. In order to examine if such star formation activity can be triggered by ionization fronts, we perform high-resolution hydrodynamics and radiative transfer simulations of the effect of radiative feedback from propagating ionization fronts on high-density clumps. We examine two sources of ultraviolet (UV) radiation field to which high-redshift ( z ∼ 3) galaxies could be exposed: one corresponding to the UV radiation originating from stars within the DLA, itself, and the other corresponding to the UV background radiation. We find that, for larger clouds, the propagating I-fronts created by local stellar sources can trigger cooling instability and collapse of significant part, up to 85 per cent, of the cloud, creating conditions for star formation in a time-scale of a few Myr. The passage of the I-front also triggers collapse of smaller clumps (with radii below ∼4 pc), but in these cases the resulting cold and dense gas does not reach conditions conducive to star formation. Assuming that 85 per cent of the gas initially in the clump is converted into stars, we obtain a star formation rate of  ∼0.25 M_⊙ yr⁻¹ kpc⁻²  . This is somewhat higher than the value derived from recent observations. On the other hand, the background UV radiation which has harder spectrum fails to trigger cooling and collapse. Instead, the hard photons which have long mean free-path heat the dense clumps, which as a result expand and essentially dissolve in the ambient medium. Therefore, the star formation activity in DLAs is strongly regulated by the radiative feedback, both from the external UV background and internal stellar sources and we predict quiescent evolution of DLAs (not starburst-like evolution).     

Formation of Expanding Envelopes in Wolf-Rayet Stars

Some features of the envelopes of WR stars have not yet been explained in terms of the standard stellar wind. One of these features is the presence of gas condensations (clumps) in the envelopes. In this paper the possible reasons for clump formation are examined along with their role in the structure of envelopes. Clumps can be formed in the transitional zone between a star and its envelope because of instabilities in the gas. This zone lies much deeper than the "optical radus" of the star, so it cannot be observed. A clump expands as it moves under the action of radiation pressure from the star and its density decreases at the same time. The clump mixes with the surrounding gas if its mass is low. Large clumps can reach the visible region of the envelope.     

Molecular Clouds: Observation to Experiment

Our ongoing investigation of how ‘Pillars’ and other structure form in molecular clouds irradiated by ultraviolet (UV) stars has revealed that the Rayleigh–Taylor instability is strongly suppressed by recombination in the photoevaporated outflow, that clumps and filaments may be key, that the evolution of structure is well-modeled by compressible hydrodynamics, and that directionality of the UV radiation may have significant effects. We discuss a generic, flexible set of laboratory experiments that can test these results.     

Experiments to Study Radiation Transport in Clumpy Media

Clumpiness of the interstellar medium may play an important role in the transfer of infrared continuum radiation in star forming regions (Boisse, 1990). For example, in homogeneous models, C II emission should be confined to the cloud edge (Viala, 1986). However, in star formation regions (such as M17SW, M17 and W51), it is observed to extend deep into the molecular cloud (Stutzki et al., 1988; Keene et al., 1985). One plausible interpretation of these observations is that, due to their clumpiness, the clouds are penetrated by UV radiation far deeper than expected from simple homogeneous models. The interaction of H II regions around young massive stars with a clumpy medium is another area of interest. Molecular clouds are well established to be clumpy on length scales down to the limits of observational resolution. Clumps can act as localized reservoirs of gas which can be injected into the surroundings by photoionization and/or hydrodynamic ablation (Dyson et al., 1995; Mathis et al., 1998). The calculation of radiation transport in hot, clumpy materials is a challenging problem. Approximate, statistical treatments of this problem have been developed by several workers, but their application has not been tested in detail. We describe laboratory experiments, using the Omega laser to test modelling of radiation transport through clumpy media in the form of inhomogeneous plasmas.     

一种针对MWISP项目分子云团块的3D CNN证认方法

分子云团块是恒星的诞生地. 分子团块的普查和其性质的全面研究将有助于了解恒星的形成乃至星系和宇宙的演化过程. 随着银河画卷计划(MWISP)项目的深入进行, 这类研究方案变得切实可行. 但是项目产生的分子云观测数据是海量的, 因此迫切需要一种能够自动识别和证认分子团块的方法. 目前应用广泛的3维分子云数据处理方法有很多, 典型的包括GaussClumps、ClumpFind、FellWalker、Reinhold等, 但都需要输入多个参数来控制它们的性能, 并且进行反复的参数优化和目测才能得到比较满意的结果. 对于大规模的观测数据, 利用现有方法进行分子团块的证认将是一项耗时耗力的任务. 为了克服传统分子云团块检测算法的局限性, 人工智能(AI)的方法将提供一个很好的解决方案. 提出了一种3D CNN (Convolutional Neural Network)方法, 它可以自动处理3D分子谱线数据, 整个过程分为检出和验证两个步骤. 首先, 通过设置较低阈值使用ClumpFind以检出候选对象, 然后通过训练好的3D CNN模型进行验证. 利用仿真数据所做的一系列的实验结果表明, 该方法的综合表现优于4种传统方法. 将该方法应用于实际的MWISP数据表明, 3D CNN方法的性能也令人满意.     

Very high and ultra-high-energy gamma-rays from giant molecular clouds

Very high and ultra-high-energy gamma-ray fluxes expected from giant molecular clouds, and the prospects of detection of these fluxes by means of air Cherenkov telescopes and EAS arrays of new generation are discussed.     

Investigating thermal evolution of the self-gravitating one dimensional molecular cloud by smoothed particle hydrodynamics

The heating of the ion-neutral (or ambipolar) diffusion may affect the thermal phases of the molecular clouds. We present an investigation on the effect of this heating mechanism in the thermal instability of the molecular clouds. A weakly ionized one-dimensional slab geometry, which is allowed for self-gravity and ambipolar diffusion, is chosen to study its thermal phases. We use the thermodynamic evolution of the slab to obtain the regions where slab cloud becomes thermally unstable. We investigate this evolution using the model of ambipolar diffusion with two-fluid smoothed particle hydrodynamics, as outlined by Hosking and Whitworth. Firstly, some parts of the technique are improved to test the pioneer works on behavior of the ambipolar diffusion in an isothermal self-gravitating slab. Afterwards, the improved two-fluid technique is used for thermal evolution of the slab. The results show that the thermal instability may persist inhomogeneities with a large density contrast at the intermediate parts of the cloud. We suggest that this feature may be responsible for the planet formation in the intermediate regions of a collapsing molecular cloud and/or may also be relevant to the formation of star forming dense cores in the clumps.     

Equivalent width, shape and proper motion of the iron fluorescent line emission from molecular clouds as an indicator of the illuminating source X-ray flux history

Observations of the diffuse emission in the 8–22 keV energy range, elongated parallel to the Galactic plane, and detection of the strong 6.4-keV fluorescent line with ∼ 1 keV equivalent width from some giant molecular clouds (e.g. Sgr B2) in the Galactic Centre region suggest that the neutral matter of these clouds is (or was) illuminated by powerful X-ray radiation, which gave rise to the reprocessed radiation. The source of this radiation remains unknown. A transient source close to the Sgr B2 cloud, or a short outburst of the X-ray emission from a supermassive black hole at the Galactic Centre are the two prime candidates under consideration. We argue that a new generation of X-ray telescopes combining very high sensitivity and excellent energy and angular resolutions would be able to discriminate between these two possibilities when studying time-dependent changes of the morphology of the surface brightness distribution, the equivalent width and the shape of the fluorescent line in Sgr B2 and other molecular clouds in the region. We note also that detection of broad and complex structures near the 6.4-keV line in the spectra of distant AGNs, which are X-ray weak now, may prove the presence of violent activity in the central engines of these objects in the past. Accurate measurements of the line shape may provide information on the time elapsed since the outburst. Proper motion (super- or subluminal) of the fluorescent radiation wave front can give additional information on the location of the source. Observations of the described effects can provide unique information on the matter distribution inside Sgr B2 and other giant molecular clouds.