Source author record

Guang-Xing Li

Guang-Xing Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2026arXiv

Scale-Aware Adversarial Analysis: A Diagnostic for Generative AI in Multiscale Complex Systems

Complex physical systems, from supersonic turbulence to the macroscopic structure of the universe, are governed by continuous multiscale dynamics. While modern machine learning architectures excel at mapping the high-dimensional observables of these systems, it remains unclear whether they internalize the governing physical laws or merely interpolate discrete statistical correlations. Standard Explainable AI (XAI) architectures, particularly perturbation-based and gradient-saliency methods, rely on pixel-wise perturbations, which generate unphysical artifacts and push inputs off the valid empirical distribution. To resolve this, we introduce a diagnostic framework driven by Constrained Diffusion Decomposition (CDD), a diffusion-based multiscale data decomposition algorithm that enables physically constrained data generation and model evaluation via scale-aware modifications. Applying this framework to a Denoising Diffusion Probabilistic Model (DDPM), we execute deterministic interventions directly within the continuous, CDD-based scale space. We demonstrate that under moderate physical perturbations, the unconstrained generative model exhibits localized structural freezing and non-linear instability rather than continuous PDE-like responses. The network fails to maintain cross-scale continuity, causing the generative trajectory to diverge when pushed into unseen physical states. By synthesizing a continuum of physically coherent states, this scale-informed methodology establishes a controlled test ground to evaluate algorithmic vulnerabilities, providing the rigorous physical constraints necessary for future architectures to respect the multiscale causality of the natural universe.

preprint2023arXiv

Effective Shielding of $\lesssim$ 10 GeV Cosmic Rays from Dense Molecular Clumps

The density of cosmic rays inside molecular clouds determines the ionization rate in the dense cores where stars form. It is also one of the drivers of astrochemistry leading to the creation of complex molecules. Through Fermi Large Area Telescope observations of nearby giant molecular clouds, we observed deficits (holes) in the gamma-ray residual map when modelling with the expected gamma-ray diffuse emission from uniform cosmic rays interacting with the molecular content. We propose that the deficit is due to the lack of penetration of the low-energy (sub-GeV to GeV) cosmic rays into denser regions or clumps. This differs from the prevailing view of fast cosmic ray transport in giant molecular clouds where the magnetic turbulence is suppressed by neutral-ion damping, as our results require a slow diffusion inside dense molecular clumps. Through modelling we find that while the shielding is negligible on the cloud scale, it becomes important in the denser, parsec-sized regions where the gravitational collapse is already at play, changing the initial condition of star formation and astrochemistry.

preprint2022arXiv

A massive Keplerian protostellar disk with flyby-induced spirals in the Central Molecular Zone

Accretion disks are an essential component in the paradigm of the formation of low-mass stars. Recent observations further identify disks surrounding low-mass pre-main-sequence stars perturbed by flybys. Whether disks around more massive stars evolve in a similar manner becomes an urgent question. We report the discovery of a Keplerian disk of a few solar masses surrounding a 32 solar-mass protostar in the Sagittarius C cloud around the Galactic Center. The disk is gravitationally stable with two embedded spirals. A combined analysis of analytical solutions and numerical simulations demonstrates that the most likely scenario to form the spirals is through external perturbations induced by a close flyby, and one such perturber with the expected parameters is identified. The massive, early O-type star embedded in this disk forms in a similar manner with respect to low-mass stars, in the sense of not only disk-mediated accretion, but also flyby-impacted disk evolution.

preprint2022arXiv

Density Exponent Analysis: Gravity-driven steepening of the density profiles of star-forming regions

The evolution of molecular interstellar clouds is a complex, multi-scale process. The power-law density exponent describes the steepness of density profiles, and it has been used to characterize the density structures of the clouds yet its usage is usually limited to spherically symmetric systems. Importing the Level-Set Method, we develop a new formalism that generates robust maps of a generalized density exponent $k_ρ$ at every location for complex density distributions. By applying it to high fidelity, high dynamical range map of the Perseus molecular cloud constructed using data from the Herschel and Planck satellites, we find that the density exponent exhibits a surprisingly wide range of variation ($-3.5 \lesssim k_ρ \lesssim -0.5$). Regions at later stages of gravitational collapse are associated with steeper density profiles. Inside a region, gas located in the vicinities of dense structures has very steep density profiles with $k_ρ \approx -3$, which forms because of depletion. This density exponent analysis reveals diverse density structures, forming a coherent picture that gravitational collapse leads to a continued steepening of the density profile. We expect our method to be effective in studying other power-law-like density structures, including granular materials and the Large-Scale Structure of the Universe.

preprint2022arXiv

Discovery of a coherent, wave-like velocity pattern for the Radcliffe Wave

Recently studies discovered that part of the Gould Belt belongs to a 2.7 kpc-long {coherent, thin} wave consisting of a chain of clouds, where a damped undulation pattern has been identified from the spatial arrangement of the clouds. We use the proper motions of Young Stellar Objects (YSOs) anchored inside the clouds to study the kinematic structure of the Radcliffe Wave in terms of $v_z$, and identify a damped, wave-like pattern from the $v_z$ space, which we call "velocity undulation". We propose a new formalism based on the Ensemble Empirical Mode Decomposition (EEMD) to determine the amplitude, period, and phase of the undulation pattern, and find that the spatial and the velocity undulation share an almost identical spatial frequency of about 1.5 kpc, and both are damped when measured from one side to the other. Measured for the first cycle, they exhibit a phase difference of around $2π/3$. The structure is oscillating around the midplane of the Milky Way disk with an amplitude of $\sim\,130\,\pm\,20\,\rm pc$. The vertical extent of the Radcliffe Wave exceeds the thickness of the molecular disk, suggesting that the undulation of the undulation signature might originate from a perturbation, e.g. the passage of a dwarf galaxy.

preprint2022arXiv

Kinematics of the molecular interstellar medium probed by Gaia: steep velocity dispersion-size relation, isotropic turbulence, and location-dependent energy dissipation

The evolution of the molecular interstellar medium is controlled by processes such as turbulence, gravity, stellar feedback, and Galactic shear. AL a part of the ISM-6D https://gxli.github.io/ISM-6D/ project, using Gaia astrometric measurements towards a sample of young stellar objects (YSOs), we study morphology and kinematic structure of the associated molecular gas. We identify 150 YSO associations with distance $d \lesssim 3 \;\rm kpc$. The YSO associations are elongated, with a median aspect ratio of 1.97, and are oriented parallel to the disk midplane, with a median angle of 30$^{\circ}$. The turbulence in the molecular clouds as probed by the YSOs is isotropic, and the velocity dispersions are related to the sizes by $σ_{v,{\rm 2D}} = 0.74\;(r/{\rm pc})^{0.67} \;({\rm km/s})\;$. The slope is on the steeper side, yet consistent with previous measurements. The energy dissipation rate of turbulence $\dotε = σ_{v,{\rm 3D}}^3 /L$ decreases with the Galactocentric distance, with a gradient of 0.2 $\rm dex \; kpc^{-1}$, which can be explained if turbulence is driven by cloud collisions. In this scenario, the clouds located in the inner Galaxy have higher chances to accrete smaller clouds and are more turbulent. Although the density structures of the complexes are anisotropic, the turbulence is consistent with being isotropic. If the alignment between density structures and the Galactic-disk mid-plane is due to shear, we expect $t_{\rm cloud} \gtrsim t_{\rm shear}\approx 30\; \rm Myr$. This cloud lifetime is longer than the turbulence crossing time, and a continuous energy injection is required to maintain the turbulence.

preprint2022arXiv

Multi-scale decomposition of astronomical maps -- a constrained diffusion method

We propose a new, efficient multi-scale method to decompose a map (or signal in general) into components maps that contain structures of different sizes. In the widely-used wave transform, artifacts containing negative values arise around regions with sharp transitions due to the application of band-limited filters. In our approach, the decomposition is achieved by solving a modified, non-linear version of the diffusion equation. This is inspired by the anisotropic diffusion methods, which establish the link between image filtering and partial differential equations. In our case, the artifact issue is addressed where the positivity of the decomposed images is guaranteed. Our new method is particularly suitable for signals which contain localized, non-linear features, as typical of astronomical observations. It can be used to study the multi-scale structures of astronomical maps quantitatively and should be useful in observation-related tasks such as background removal. We thus propose a new measure called the "scale spectrum", which describes how the image values distribute among different components in the scale space, to describe maps. The method allows for input arrays of an arbitrary number of dimensions, and a python3 implementation of the algorithms is included in the Appendix and available at https://gxli.github.io/Constrained-Diffusion-Decomposition/.

preprint2022arXiv

Network of Star Formation: Fragmentation controlled by scale-dependent turbulent pressure and accretion onto the massive cores revealed in the Cygnus-X GMC complex

Molecular clouds have complex density structures produced by processes including turbulence and gravity. We propose a triangulation-based method to dissect the density structure of a molecular cloud and study the interactions between dense cores and their environments. In our {approach}, a Delaunay triangulation is constructed, which consists of edges connecting these cores. Starting from this construction, we study the physical connections between neighboring dense cores and the ambient environment in a systematic fashion. We apply our method to the Cygnus-X massive GMC complex and find that the core separation is related to the mean surface density by $Σ_{\rm edge} \propto l_{\rm core }^{-0.28 }$, which can be explained by {fragmentation controlled by a scale-dependent turbulent pressure (where the pressure is a function of scale, e.g. $p\sim l^{2/3}$)}. We also find that the masses of low-mass cores ($M_{\rm core} < 10\, M_{\odot}$) are determined by fragmentation, whereas massive cores ($M_{\rm core} > 10\, M_{\odot}$) grow mostly through accretion. The transition from fragmentation to accretion coincides with the transition from a log-normal core mass function (CMF) to a power-law CMF. By constructing surface density profiles measured along edges that connect neighboring cores, we find evidence that the massive cores have accreted a significant fraction of gas from their surroundings and thus depleted the gas reservoir. Our analysis reveals a picture where cores form through fragmentation controlled by scale-dependent turbulent pressure support, followed by accretion onto the massive cores, {and the method can be applied to different regions to achieve deeper understandings in the future.

preprint2022arXiv

Weather Forecast of the Milky Way: Shear and Stellar feedback determine the lives of Galactic-scale filaments

The interstellar medium (ISM) is an inseparable part of the Milky Way ecosystem whose evolutionary history remains a challenging question. We trace the evolution of the molecular ISM using a sample of Young Stellar Objects (YSO) association --molecular cloud complex (YSO-MC complex). We derive their three-dimensional (3D) velocities by combining the Gaia astrometric measurements of the YSO associations and the CO observations of the associated molecular clouds. Based on the 3D velocities, we simulate the motions of the YSO-MC complexes in the Galactic potential and forecast the ISM evolution by tracing the motions of the individual complexes, and reveal the roles of shear and stellar feedback in determining ISM evolution: Galactic shear stretches Galactic-scale molecular cloud complexes, such as the G120 Complex, into Galactic-scale filaments, and it also contributes to the destruction of the filaments; while stellar feedback creates interconnected superbubbles whose expansion injects peculiar velocities into the ISM. The Galactic-scale molecular gas clumps are often precursors of the filaments and the Galactic-scale filaments are transient structures under a constant stretch by shear. This evolutionary sequence sets a foundation to interpret other gas structures. Animations are available at https://gxli.github.io/ISM-6D/movie.html.

preprint2022arXiv

Widespread subsonic turbulence in Ophiuchus North 1

Supersonic motions are common in molecular clouds. (Sub)sonic turbulence is usually detected toward dense cores and filaments. However, it remains unknown whether (sub)sonic motions at larger scales ($\gtrsim$1~pc) can be present in different environments or not. Located at a distance of about 110 pc, Ophiuchus North 1 (Oph N1) is one of the nearest molecular clouds that allows in-depth investigation of its turbulence properties by large-scale mapping observations of single-dish telescopes. We carried out the $^{12}$CO ($J=1-0$) and C$^{18}$O ($J=1-0$) imaging observations toward Oph N1 with the Purple Mountain Observatory 13.7 m telescope. The observations have an angular resolution of $\sim$55\arcsec (i.e., 0.03~pc). Most of the whole C$^{18}$O emitting regions have Mach numbers of $\lesssim$1, demonstrating the large-scale (sub)sonic turbulence across Oph N1. Based on the polarization measurements, we estimate the magnetic field strength of the plane-of-sky component to be $\gtrsim$9~$μ$G. We infer that Oph N1 is globally sub-Alfv{é}nic, and is supported against gravity mainly by the magnetic field. The steep velocity structure function can be caused by the expansion of the Sh~2-27 H{\scriptsize II} region or the dissipative range of incompressible turbulence. Our observations reveal a surprising case of clouds characterised by widespread subsonic turbulence and steep size-linewidth relationship. This cloud is magnetized where ion-neutral friction should play an important role.

preprint2021arXiv

The DR21(OH) Trident -- Resolving the Massive Ridge into Three Entangled Fibers As the Initial Condition of Cluster Formation

DR21(OH) ridge, the central part of a high-mass star and cluster forming hub-filament system, is resolved spatially and kinematically into three nearly parallel fibers (f1, f2, and f3) with a roughly north-south orientation, using the observations of molecular transitions of H$^{13}$CO$^+$ (1-0), N$_2$H$^+$ (1-0), and NH$_2$D (1$_{1,1}$-1$_{0,1}$) with the Combined Array for Research in Millimeter Astronomy. These fibers are all mildly supersonic ($σ_{\rm V}$ about 2 times the sound speed), having lengths around 2 pc and widths about 0.1 pc, and they entangle and conjoin in the south where the most active high-mass star formation takes place. They all have line masses 1 - 2 orders of magnitude higher than their low-mass counterparts and are gravitationally unstable both radially and axially. However, only f1 exhibits high-mass star formation all the way along the fiber, yet f2 and f3 show no signs of significant star formation in their northern parts. A large velocity gradient increasing from north to south is seen in f3, and can be well reproduced with a model of free-fall motion toward the most massive and active dense core in the region, which corroborates the global collapse of the ridge and suggests that the disruptive effects of the tidal forces may explain the inefficiency of star formation in f2 and f3. On larger scales, some of the lower-density, peripheral filaments are likely to be the outer extensions of the fibers, and provide hints on the origin of the ridge.

preprint2020arXiv

A mean density of $112\, M_{\odot}\,\rm pc^{-3}$ for Central Molecular Zone clumps -- Evidences for shear-enabled pressure equilibrium in the Galactic Center

We carry out a systematic study of the density structure of gas in the Central Molecular Zone (CMZ) in the Galactic center by extracting clumps from the APEX Telescope Large Area Survey of the Galaxy survey at 870 $μ$m. We find that the clumps follow a scaling of $m = ρ_0 r^3$ which corresponds to a characteristic density of $n_{\rm H_2} = 1.6 \times 10^3\,\rm cm^{-3}$ ($ρ_0 =112\;M_{\odot}\;\rm pc^{-3}$) with a variation of $\approx 0.5\,\rm dex$, where we assumed a gas-to-dust mass ratio of 100. This characteristic density can be interpreted as the result of thermal pressure equilibrium between the molecular gas and the warm ambient interstellar medium. Such an equilibrium can plausibly be established since shear has approximately the same strength as self-gravity. Our findings may explain the fact that star formation in the CMZ is highly inefficient compared to the rest of the Milky Way disk. We also identify a population of clumps whose densities are two orders of magnitudes higher in the vicinity of the Sgr B2 region, which we propose are produced by collisions between the clumps of lower densities. For these collisions to occur, processes such as compressive tides probably have created the appropriate condition by assembling the clumps together.

preprint2020arXiv

Edge collapse and subsequent longitudinal accretion in Filament S242

Filament S242 is 25 pc long with massive clumps and YSO clusters concentrated in its end regions; it is considered a good example of edge collapse. We mapped this filament in the $^{12}$CO(1-0) and $^{13}$CO(1-0) lines. A large-scale velocity gradient along filament S242 has been detected; the relative velocity between the two end-clumps is $\sim$ 3 km s$^{-1}$, indicating an approaching motion between them. These signatures are consistent with the filament S242 being formed through the collapse of a single elongated entity, where an effect known as "gravitational focusing" drives the ends of the filament to collapse (edge collapse). Based on this picture, we estimate a collapse timescale of $\sim$ 4.2 Myr, which is the time needed for a finite and elongated entity evolving to the observed filament S242. For the whole filament, we find that increases in surface densities lead to increases in velocity dispersion, which can be consistently explained as the result of self-gravity. We also calculated the contribution of longitudinal collapse to the observed velocity dispersion and found it to be the dominant effect in driving the gas motion near the end-clumps. We propose that our filament S242 is formed through a two-stage collapse model, where the edge collapse of a truncated filament is followed by a stage of longitudinal accretion toward the dense end-clumps.

preprint2020arXiv

Probing the initial conditions of high-mass star formation -- IV. Gas dynamics and NH$_2$D chemistry in high-mass precluster and protocluster clumps

The initial stage of star formation is a complex area study because of its high density and low temperature. Under such conditions, many molecules become depleted from the gas phase by freezing out onto dust grains. However, the deuterated species could remain gaseous and are thus ideal tracers. We investigate the gas dynamics and NH$_2$D chemistry in eight massive pre/protocluster clumps. We present NH$_2$D 1$_{11}$-1$_{01}$ (at 85.926 GHz), NH$_3$ (1, 1) and (2, 2) observations in the eight clumps using the PdBI and the VLA, respectively. We find that the distribution between deuterium fractionation and kinetic temperature shows a number density peak at around $T_{\rm kin}=16.1$ K, and the NH$_2$D cores are mainly located at a temperature range of 13.0 to 22.0 K. We detect seven instances of extremely high deuterium fractionation of $1.0 \leqslant D_{\rm frac} \leqslant 1.41$. We find that the NH$_2$D emission does not appear to coincide exactly with either dust continuum or NH$_3$ peak positions, but often surrounds the star-formation active regions. This suggests that the NH$_{2}$D has been destroyed by the central young stellar object (YSO) due to its heating. The detected NH$_2$D lines are very narrow with a median width of $\rm 0.98\pm0.02 km/s$. The extracted NH$_2$D cores are gravitationally bound ($α_{\rm vir} < 1$), are likely prestellar or starless, and can potentially form intermediate-mass or high-mass stars. Using NH$_3$ (1, 1) as a dynamical tracer, we find very complicated dynamical movement, which can be explained by a combined process with outflow, rotation, convergent flow, collision, large velocity gradient, and rotating toroids. High deuterium fractionation strongly depends on the temperature condition. NH$_2$D is a poor evolutionary indicator of high-mass star formation in evolved stages, but a useful tracer in the starless and prestellar cores.

preprint2016arXiv

ATLASGAL -- A Galaxy-wide sample of dense filamentary structures

[Abridged] Aims. We study the properties of filamentary structures from the ATLASGAL survey. Methods. We use the DisPerSE algorithm to identify spatially coherent structures located across the inner-Galaxy (300 < l < 60 and |b| < 1.5). Results. We have determined distances, masses and physical sizes for 241 of the filamentary structures. We find a median distance of 3.8 kpc, a mean mass of a few 10^3 m_sun, a mean length of ~6pc and a mass-to-length ratio of (M/L) ~200-2000M_sun/ pc. We also find that these filamentary structures are tightly correlated with the spiral arms in longitude and velocity, and that their semi-major axis is preferentially aligned parallel to the Galactic mid-plane and therefore with the direction of large-scale Galactic magnetic field. We find many examples where the dense filaments identified in ATLASGAL are associated with larger scale filamentary structures (~100 pc), and argue that this is likely to be common, and as such these may indicate a connection between large-scale Galactic dynamics and star formation. Conclusions. We have produced a large and Galaxy-wide catalogue of dense filamentary structures that are representative of a particular size and mass range not previously well studied in the literature. Analyses of the properties and distribution of these filaments reveals that they are correlated with the spiral arms and make a significant contribution to star formation in the Galaxy. Massive star formation is ongoing within ~20% of the filaments and is strongly correlated with the filaments with the largest mass-to- length ratios. The luminosity of the embedded sources has a similar distribution to the Galactic-wide samples of young massive stars and can therefore be considered to be representative.

preprint2016arXiv

Constructing multi-scale gravitational energy spectra from molecular cloud surface density PDF -- Interplay between turbulence and gravity

(Abridged) We derive an analytical formula which provides estimates on multiscale gravitational energy distribution using the observed surface density PDF. Our analytical formalism also enables one to convert the observed column density PDF into an estimated volume density PDF, and to obtain average radial density profile $ρ(r)$. For a region with $N_{\rm col} \sim N^{-γ_{\rm N}}$, the gravitational energy spectra is $E_{\rm p}(k)\sim k^{-4(1 - 1/γ_{\rm N})}$. We apply the formula to observations of molecular clouds, and find that a scaling index of $-2$ of the surface density PDF implies that $ρ\sim r^{-2}$ and $E_{\rm p}(k) \sim k^{-2}$. The results are valid from the cloud scale (a few parsec) to around $\sim 0.1 \;\rm pc$. Because of the resemblance the scaling index of the gravitational energy spectrum and the that of the kinetic energy power spectrum of the Burgers turbulence (where $E\sim k^{-2}$), our result indicates that gravity can act effectively against turbulence over a multitude of physical scales. This is the critical scaling index which divides molecular clouds into two categories: clouds like Orion and Ophiuchus have shallower power laws, and the amount of gravitational energy is too large for turbulence to be effective inside the cloud. Because gravity dominates, we call this type of cloud g-type clouds. On the other hand, clouds like the California molecular cloud and the Pipe nebula have steeper power laws, and turbulence can overcome gravity if it can cascade effectively from the large scale. We call this type of cloud t-type clouds. The analytical formula can be used to determine if gravity is dominating cloud evolution when the column density probability distribution function (PDF) can be reliably determined.

preprint2016arXiv

Criteria for gravitational instability and quasi-isolated gravitational collapse in turbulent medium

We study the evolution of structures in turbulent, self-gravitating media, and present an analytical criterion $M_{\rm crit} \approx ε_{\rm cascade}^{2/3} η^{-2/3} G^{-1} l^{5/3}$ (where $M_{\rm crit}$ is the critical mass, $l$ is the scale, $ε_{\rm cascade}\approx ησ_{\rm v}^3 /l $ is the turbulence energy dissipation rate of the ambient medium, $G$ is the gravitational constant, $σ_{\rm v}$ is the velocity dispersion, $l$ is the scale and $η\approx 0.2$ is an efficiency parameter) for an object to undergo quasi-isolated gravitational collapse. The criterion also defines the critical scale ($l_{\rm crit} \approx ε_{\rm cascade}^{1/2} η^{-1/2} G^{-3/4} ρ^{-3/4}$) for turbulent gravitational instability to develop. The analytical formalism explains the size dependence of the masses of the progenitors of star clusters ($M_{\rm cluster} \sim R_{\rm cluster}^{1.67}$) in our Galaxy.

preprint2016arXiv

Gravitational acceleration and edge effects in molecular clouds

Gravity plays important roles in the evolution of molecular clouds. We present an acceleration mapping method to estimate the acceleration induced by gravitational interactions in molecular clouds based on observational data. We find that the geometry of a region has a significant impact on the behavior of gravity. In the Pipe nebula which can be approximated as a gas filament, we find that gravitational acceleration can effectively compress the end of this filament, which may have triggered star formation. We identify this as the "gravitational focusing" effect proposed by Burkert & Hartman (2004). In the sheet-like IC348-B3 region, gravity can lead to collapse at its edge, while in the centrally condensed NGC1333 cluster-forming region gravity can drive accretion towards the center. In general, gravitational acceleration tends to be enhanced in the localized regions around the ends of the filaments and the edges of sheet-like structures. Neglecting magnetic fields, these "gravitational focusing" and "edge collapse" effects can promote the formation of dense gas in a timescale that is much shorter than the global dynamical time. Since the interstellar medium is in general structured, these edge effects should be prevalent.

preprint2016arXiv

N131: A dust bubble born from the disruption of a gas filament

OB type stars have strong ionizing radiation, and drive energetic winds. The ultraviolet (UV) radiation from ionizing stars may heat dust and ionize gas to sweep up an expanding bubble shell. This shell may be the result of feedback leading to a new generation of stars. N131 is an infrared dust bubble residing in a molecular filament. We study the formation and fragmentation of this bubble with multi-wavelength dust and gas observations. Towards the bubble N131, we analyzed archival multi-wavelength observations including 3.6, 4.5, 5.8, 8.0, 24, 70, 160, 250, 350, 500 $μ$m, 1.1 mm, and 21 cm. In addition, we performed new observations of CO (2-1), CO (1-0), and $^{13}$CO (1-0) with the IRAM 30-m telescope. Multi-wavelength dust and gas observations reveal a ringlike shell with compact fragments, two filamentary structures, and a secondary bubble N131-A. The bubble N131 is a rare object with a large hole at 24 $μ$m and 21 cm in the direction of its center. The dust and gas clumps are compact and might have been compressed at the inner edge of the ringlike shell, while they are extended and might be pre-existing at the outer edge. The column density, excitation temperature, and velocity show a potentially hierarchical distribution from the inner to outer edge of the ringlike shell. We also detected the front and back sides of the secondary bubble N131-A in the direction of its center. The derived Lyman-continuum ionizing photon flux within N131-A is equivalent to an O9.5 star. Based on the above, we suggest that the bubble N131 might be triggered by the strong stellar winds from a group of massive stars inside the bubble. We propose a scenario in which the bubble N131 forms from the disruption of a gas filament by expansion of HII region, strong stellar winds, and fragments under self-gravity.

preprint2016arXiv

Probing the multi-scale interplay between gravity and turbulence - Power-law like gravitational energy spectra of the Orion Complex

Gravity plays a determining role in the evolution of the molecular ISM. In \citet{2016arXiv160304342L}, we proposed a measure called gravitational energy spectrum to quantify the importance of gravity on multiple physical scales. In this work, using a wavelet-based decomposition technique, we derive the gravitational energy spectra of the Orion A and the Orion B molecular cloud from observational data. The gravitational energy spectra exhibit power-law-like behaviours. From a few pc down to $\sim 0.1 $ pc scale, the Orion A and Orion B molecular cloud have $E_{\rm p}(k)\sim k^{-1.88}$ and $E_{\rm p}(k)\sim k^{-2.09}$, respectively. These scaling exponents are close to the scaling exponents of the kinetic energy power spectrum of compressible turbulence (where $E\sim k^{-2}$), with a near-equipartition of turbulent versus gravitational energy on multiple scales. This provides a clear evidence that gravity is able to counteract effectively against turbulent motion for these length scales. The results confirm our earlier analytical estimates. For the Orion A molecular cloud, gravity inevitably dominates over turbulence inside the cloud. Our results provide a clear observational proof that gravity is playing a determining role in the evolution these molecular clouds from the cloud scale down to $\sim 0.1\;\rm pc$. { However, turbulence is likely to dominate in clouds like California. } The method is general and should be applicable to all the astrophysical problems where gravity plays a role.

preprint2016arXiv

Revealing a spiral-shaped molecular cloud in our galaxy - Cloud fragmentation under rotation and gravity

The dynamical processes that control star formation in molecular clouds are not well understood, and in particular, it is unclear if rotation plays a major role in cloud evolution. We investigate the importance of rotation in cloud evolution by studying the kinematic structure of a spiral-shaped Galactic molecular cloud G052.24+00.74. The cloud belongs to a large filament, and is stretching over ~ 100 pc above the Galactic disk midplane. The spiral-shaped morphology of the cloud suggests that the cloud is rotating. We have analysed the kinematic structure of the cloud, and study the fragmentation and star formation. We find that the cloud exhibits a regular velocity pattern along west-east direction - a velocity shift of ~ 10 km/s at a scale of ~ 30 pc. The kinematic structure of the cloud can be reasonably explained by a model that assumes rotational support. Similarly to our Galaxy, the cloud rotates with a prograde motion. We use the formalism of Toomre (1964) to study the cloud's stability, and find that it is unstable and should fragment. The separation of clumps can be consistently reproduced assuming gravitational instability, suggesting that fragmentation is determined by the interplay between rotation and gravity. Star formation occurs in massive, gravitational bound clumps.

preprint2015arXiv

G-virial: Gravity-based structure analysis of molecular clouds

We present the G-virial method (available at http://gxli.github.io/G-virial/) which aims to quantify (1) the importance of gravity in molecular clouds in the position-position-velocity (PPV) space, and (2) properties of the gas condensations in molecular clouds. Different from previous approaches that calculate the virial parameter for different regions, our new method takes gravitational interactions between all the voxels in 3D PPV data cubes into account, and generates maps of the importance of gravity. This map can be combined with the original data cube to derive relations such as the mass-radius relation. Our method is important for several reasons. First, it offers the the ability to quantify the centrally condensed structures in the 3D PPV data cubes, and enables us to compare them in an uniform framework. Second, it allows us to understand the importance of gravity at different locations in the data cube, and provides a global picture of gravity in clouds. Third, it offers a robust approach to decomposing the data into different regions which are gravitationally coherent. To demonstrate the application of our method we identified regions from the Perseus and Ophiuchus molecular clouds, and analyzed their properties. We found an increase in the importance of gravity towards the centers of the individual molecular condensations. We also quantified the properties of the regions in terms of mass-radius and mass-velocity relations. Through evaluating the virial parameters based on the G-virial, we found that all our regions are almost gravitationally bound. Cluster-forming regions appear are more centrally condensed.

preprint2013arXiv

A 500 pc filamentary gas wisp in the disk of the Milky Way

Star formation occurs in molecular gas. In previous studies, the structure of the molecular gas has been studied in terms of molecular clouds, but has been overlooked beyond the cloud scale. We present an observational study of the molecular gas at 49.5 degree <l<52.5 degree and -5.0 km/s <v_lsr <17.4 km/s. The molecular gas is found in the form of a huge (>= 500 pc) filamentary gas wisp. This has a large physical extent and a velocity dispersion of ~5 km/s. The eastern part of the filamentary gas wisp is located ~130 pc above the Galactic disk (which corresponds to 1.5-4 e-folding scale-heights), and the total mass of the gas wisp is >= 1 X 10^5 M_sun. It is composed of two molecular clouds and an expanding bubble. The velocity structure of the gas wisp can be explained as a smooth quiescent component disturbed by the expansion of a bubble. That the length of the gas wisp exceeds by much the thickness of the molecular disk of the Milky Way is consistent with the cloud-formation scenario in which the gas is cold prior to the formation of molecular clouds. Star formation in the filamentary gas wisp occurs at the edge of a bubble (G52L nebula), which is consistent with some models of triggered star formation.

preprint2013arXiv

Turbulent entrainment origin of protostellar outflows

Protostellar outflow is a prominent process that accompanies the formation of stars. It is generally agreed that wide-angled protostellar outflows come from the interaction between the wind from a forming star and the ambient gas. However, it is still unclear how the interaction takes place. In this work, we theoretically investigate the possibility that the outflow results from interaction between the wind and the ambient gas in the form of turbulent entrainment. In contrast to the previous models, turbulent motion of the ambient gas around the protostar is taken into account. In our model, the ram-pressure of the wind balances the turbulent ram-pressure of the ambient gas, and the outflow consists of the ambient gas entrained by the wind. The calculated outflow from our modelling exhibits a conical shape. The total mass of the outflow is determined by the turbulent velocity of the envelope as well as the outflow age, and the velocity of the outflow is several times higher than the velocity dispersion of the ambient gas. The outflow opening angle increases with the strength of the wind and decreases with the increasing ambient gas turbulence. The outflow exhibits a broad line width at every position. We propose that the turbulent entrainment process, which happens ubiquitously in nature, plays a universal role in shaping protostellar outflows.

preprint2010arXiv

Emergent Spectra From Disks Surrounding Kerr Black Holes: Effect of Photon Trapping and Disk Self-Shadowing

Based on a new estimation of their thickness, the global properties of relativistic slim accretion disks are investigated in this work. The resulting emergent spectra are calculated using the relativistic ray-tracing method, in which we neglect the self-irradiation of the accretion disk. The angular dependence of the disk luminosity, the effects of the heat advection and the disk thickness on the estimation of the black hole spin are discussed. Compare to the previous works, our improvements are that we use the self-consistent disk equations and we consider the disk self-shadowing effect. We find that at the moderate accretion rate, the radiation trapped in the outer region of the accretion disks will escape in the inner region of the accretion disk and contribute to the emergent spectra. At the high accretion rate, for the large inclination and large black hole spin, both the disk thickness and the heat advection have significant influence on the emergent spectra. Consequently, these effects will influence the measurement of the black hole spin based on the spectra fitting and influence the angular dependence of the luminosity. For the disks around Kerr black holes with $a=0.98$, if the disk inclination is greater than $60^\circ$, and their luminosity is beyond 0.2 Eddington luminosity, the spectral model which is based on the relativistic standard accretion disk is no longer applicable for the spectra fitting. We also confirm that the effect of the self-shadowing is significantly enhanced by the light-bending, which implies that the non-relativistic treatment of the self-shadowing is inaccurate. According to our results, the observed luminosity dependence of the measured spin suggests that the disk self-shadowing significantly shapes the spectra of GRS 1915+105, which might lead to the underestimation of the black hole spin for the high luminosity states.