Researcher profile

Shaobo Zhang

Shaobo Zhang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2025arXiv

Symmetry-Checking in Band Structure Calculations on a Noisy Quantum Computer

Band crossings in electronic band structures play an important role in determining the electronic, topological, and transport properties in solid-state systems, making them central to both condensed matter physics and materials science. The emergence of noisy intermediate-scale quantum (NISQ) processors has sparked great interest in developing quantum algorithms to compute band structure properties of materials. While significant research has been reported on computing ground state and excited state energy bands in the presence of noise that breaks the degeneracy, identifying the symmetry at crossing points using quantum computers is still an open question. In this work, we propose a method for identifying the symmetry of bands around crossings and anti-crossings in the band structure of bilayer graphene with two distinct configurations on a NISQ device. The method utilizes eigenstates at neighbouring $\mathbf{k}$ points on either side of the touching point to recover the local symmetry by implementing a character-checking quantum circuit that uses ancilla qubit measurements for a probabilistic test. We then evaluate the performance of our method under a depolarizing noise model, using four distinct matrix representations of symmetry operations to assess its robustness. Finally, we demonstrate the reliability of our method by correctly identifying the correct band crossings of AA-stacked bilayer graphene around $K$ point, using the character-checking circuit implemented on a noisy IBM quantum processor $ibm\_marrakesh$.

preprint2023arXiv

Magneto hydrodynamic simulations of the supernova remnant G1.9+0.3

The youngest Galactic supernova remnant G1.9+0.3 shows a discrete feature between its radio and X-ray morphologies. The observed radio morphology features a single maximum in the north, while the X-ray observation shows two opposite 'ears' on the east and west sides. Using 3D magneto hydrodynamical simulations, we investigate the formation of the discrete feature of the remnant. We have tested different parameters for better simulation and reproduced similar discrete features under an environment with density gradient and an environment with clump, which provides a possible explanation of the observation.

preprint2022arXiv

CO Emission Delineating the Interface between the Milky Way Nuclear Wind Cavity and the Gaseous Disk

Based on the MWISP survey, we study high-z CO emission toward the tangent points, in which the distances of the molecular clouds (MCs) are well determined. In the region of l=12-26 deg and |b|<5.1 deg, a total of 321 MCs with |z|> 110 pc are identified, of which nearly 30 extreme high-z MCs (EHMCs at |z|> 260 pc) are concentrated in a narrow region of R_GC=2.6-3.1 kpc. The EHMC concentrations, together with other high-z MCs at R_GC=2.3-2.6 kpc, constitute molecular crater-wall structures surrounding the edges of the HI voids that are physically associated with the Fermi bubbles. Intriguingly, some large high-z MCs, which lie in the crater walls above and below the Galactic plane, show cometary structures with the head toward the plane, favouring the scenario that the entrained molecular gas moves with the multi-phase flows from the plane to the high-z regions. We suggest that the Milky Way nuclear wind has a significant impact on the Galactic gaseous disk. The powerful nuclear wind at ~3-6 Myr ago is likely responsible for the observational features, (1) the enhanced CO gas lying in the edges of the HI voids, (2) the deficiency of atomic and molecular gas within R_GC<3 kpc, (3) the possible connection between the EHMC concentrations and the 3-kpc arm, and (4) the elongated high-z MCs with the tail pointing away from the Galactic plane.

preprint2022arXiv

CO(J = 1-0) Observations toward the Filamentary Cloud in the Galactic Region of $153.60^{\circ} \leqslant l \leqslant 156.50^{\circ}$ and $1.85^{\circ} \leqslant b \leqslant 3.50^{\circ}$

We present observations of $J$=1-0 transition lines of ${ }^{12} \mathrm{CO}$, ${ }^{13} \mathrm{CO}$, and $\mathrm{C}^{18} \mathrm{O}$ towards the Galactic region of $153.60^{\circ} \leqslant l \leqslant 156.50^{\circ}$ and $1.85^{\circ} \leqslant b \leqslant 3.50^{\circ}$, using the Purple Mountain Observatory (PMO) 13.7 m millimeter telescope. Based on the \tht data, one main filament and five sub-filaments are found together as a network structure in the velocity interval of $[-42.5, -30.0] \,\mathrm{km} \mathrm{\,s}^{-1}$. The kinematic distance of this molecular cloud (MC) is estimated to be $\sim4.5 \mathrm{\,kpc}$. The median length, width, excitation temperature, line mass of these filaments are $\sim49 \mathrm{\,pc}$, $\sim2.9 \mathrm{\,pc}$, $\sim8.9 \mathrm{\,K}$, and $\sim39 \,M_{\odot} \mathrm{pc}^{-1}$, respectively. The velocity structures along these filaments exhibit oscillatory patterns, which are likely caused by the fragmentation or accretion process along these filaments. The maximum accretion rate is estimated to be as high as $\sim700 \,M_{\odot} \mathrm{pc}^{-1}$. A total of $\sim162$ \tht clumps and $\sim 103$ young stellar objects (YSOs) are identified in this region. Most of the clumps are in gravitationally bound states. Three \hii regions (G154.359+2.606, SH2-211, SH2-212) are found to be located in the apexes of the filaments. Intense star forming activities are found along the entire filamentary cloud. The observed results may help us to better understand the link between filaments and massive star formation.

preprint2022arXiv

Comparative analysis of error mitigation techniques for variational quantum eigensolver implementations on IBM quantum system

Quantum computers are anticipated to transcend classical supercomputers for computationally intensive tasks by exploiting the principles of quantum mechanics. However, the capabilities of the current generation of quantum devices are limited due to noise or errors, and therefore implementation of error mitigation and/or correction techniques is pivotal to reliably process quantum algorithms. In this work, we have performed a comparative analysis of the error mitigation capability of the [[4,2,2]] quantum error-detecting code (QEC method), duplicate circuit technique, and the Bayesian read-out error mitigation (BREM) approach in the context of proof-of-concept implementations of variational quantum eigensolver (VQE) algorithm for determining the ground state energy of H$_2$ molecule. Based on experiments on IBM quantum device, our results show that the duplicate circuit approach performs superior to the QEC method in the presence of the hardware noise. A significant impact of cross-talk noise was observed when multiple mappings of the duplicate circuit and the QEC method were implemented simultaneously $-$ again the duplicate circuit approach overall performed better than the QEC method. To gain further insights into the performance of the studied error mitigation techniques, we also performed quantum simulations on IBM system with varying strengths of depolarising circuit noise and read-out errors which further supported the main finding of our work that the duplicate circuit offer superior performance towards mitigating of errors when compared to the QEC method. Our work reports a first assessment of the duplicate circuit approach for a quantum algorithm implementation and the documented evidence will pave the way for future scalable implementations of the duplicated circuit techniques for the error-mitigated practical applications of near-term quantum computers.

preprint2022arXiv

Dependence of Molecular Cloud Samples on Angular Resolution, Sensitivity, and Algorithms

In this work, we investigate the observational and algorithmic effects on molecular cloud samples identified from position-position-velocity (PPV) space. By smoothing and cutting off the high quality data of the Milky Way Imaging Scroll Painting (MWISP) survey, we extract various molecular cloud samples from those altered data with the DBSCAN (density-based spatial clustering of applications with noise) algorithm. Those molecular cloud samples are subsequently used to gauge the significance of sensitivity, angular/velocity resolution, and DBSCAN parameters. Two additional surveys, the FCRAO Outer Galaxy Survey (OGS) and the CfA-Chile 1.2 m complete CO (CfA-Chile) survey, are used to verify the MWISP results. We found that molecular cloud catalogs are not unique and the boundary and therefore the number shows strong variation with angular resolution and sensitivity. At low angular resolution (large beam sizes), molecular clouds merge together in PPV space, while low sensitivity (high cutoffs) misses small faint molecular clouds and takes bright parts of large molecular clouds as single ones. At high angular resolution and sensitivity, giant molecular clouds (GMCs) are resolved into individual clouds, and their diffuse components are also revealed. Consequently, GMCs are more appropriately interpreted as clusters or aggregates of molecular clouds, i.e., GMCs represent molecular cloud samples themselves.

preprint2022arXiv

Gas Column Density Distribution of Molecular Clouds in the Third Quadrant of the Milky Way

We have obtained column density maps for an unbiased sample of 120 molecular clouds in the third quadrant of the Milky Way mid-plane (b$\le |5|^{\circ}$) within the galactic longitude range from 195$^{\circ}$ to 225$^{\circ}$, using the high sensitivity $^{12}$CO and $^{13}$CO ($J=1-0$) data from the Milky Way Imaging Scroll Painting (MWISP) project. The probability density functions of the molecular hydrogen column density of the clouds, N-PDFs, are fitted with both log-normal (LN) function and log-normal plus power-law (LN+PL) function. The molecular clouds are classified into three categories according to their shapes of N-PDFs, i.e., LN, LN+PL, and UN (unclear), respectively. About 72\% of the molecular clouds fall into the LN category, while 18\% and 10\% into the LN+PL and UN categories, respectively. A power-law scaling relation, $σ_s\propto N_{H_2}^{0.44}$, exists between the width of the N-PDF, $σ_s$, and the average column density, $N_{H_2}$, of the molecular clouds. However, $σ_s$ shows no correlation with the mass of the clouds. A correlation is found between the dispersion of normalized column density, $σ_{N/\rm <N>}$, and the sonic Mach number, $\mathcal{M}$, of molecular clouds. Overall, as predicted by numerical simulations, the N-PDFs of the molecular clouds with active star formation activity tend to have N-PDFs with power-law high-density tails.

preprint2022arXiv

Molecular Gas Structures traced by $^{13}$CO Emission in the 18,190 $^{12}$CO Molecular Clouds from the MWISP Survey

After the morphological classification of the 18,190 $^{12}$CO molecular clouds, we further investigate the properties of their internal molecular gas structures traced by the $^{13}$CO($J=$ 1$-$0) line emissions. Using three different methods to extract the $^{13}$CO gas structures within each $^{12}$CO cloud, we find that $\sim$ 15$\%$ of $^{12}$CO clouds (2851) have $^{13}$CO gas structures and these $^{12}$CO clouds contribute about 93$\%$ of the total integrated flux of $^{12}$CO emission. In each of 2851 $^{12}$CO clouds with $^{13}$CO gas structures, the $^{13}$CO emission area generally does not exceed 70$\%$ of the $^{12}$CO emission area, and the $^{13}$CO integrated flux does not exceed 20$\%$ of the $^{12}$CO integrated flux. We reveal a strong correlation between the velocity-integrated intensities of $^{12}$CO lines and those of $^{13}$CO lines in both $^{12}$CO and $^{13}$CO emission regions. This indicates the H$_{2}$ column densities of molecular clouds are crucial for the $^{13}$CO lines emission. After linking the $^{13}$CO structure detection rates of the 18,190 $^{12}$CO molecular clouds to their morphologies, i.e. nonfilaments and filaments, we find that the $^{13}$CO gas structures are primarily detected in the $^{12}$CO clouds with filamentary morphologies. Moreover, these filaments tend to harbor more than one $^{13}$CO structure. That demonstrates filaments not only have larger spatial scales, but also have more molecular gas structures traced by $^{13}$CO lines, i.e. the local gas density enhancements. Our results favor the turbulent compression scenario for filament formation, in which dynamical compression of turbulent flows induces the local density enhancements. The nonfilaments tend to be in the low-pressure and quiescent turbulent environments of the diffuse interstellar medium.

preprint2022arXiv

Search for monopole-dipole interactions at the sub-millimeter range with a $^{129}$Xe-$^{131}$Xe-Rb comagnetometer

Monopole-dipole interactions involving scalar couplings between a spin and a massive particle violate both P- and T-symmetry, and can be mediated by axions. We use a $^{129}$Xe-$^{131}$Xe-Rb atomic cell comagnetometer to measure the ratio of precession frequencies between the two xenon isotopes, and search for changes of the ratio correlated with the distance between the atomic cell and a non-magnetic bismuth germanate (BGO) crystal. A modulated Rb polarization scheme is used to suppress systematic effects by two orders of magnitude. The null results of this search improve the upper limit on the coupling strength $g_{s}^Ng_{p}^{n}$ over the interaction range 0.11 - 0.55 mm, and by a maximum improvement factor of 30 at 0.24 mm. The corresponding propagator mass range of this new excluded region covers 0.36 - 1.80 meV.

preprint2021arXiv

Preliminary analysis on the noise characteristics of MWISP data

Noise is a significant part within a millimeter-wave molecular line datacube. Analyzing the noise improves our understanding of noise characteristics, and further contributes to scientific discoveries. We measure the noise level of a single datacube from MWISP and perform statistical analyses. We identified major factors which increase the noise level of a single datacube, including bad channels, edge effects, baseline distortion and line contamination. Cleaning algorithms are applied to remove or reduce these noise components. As a result, we obtained the cleaned datacube in which noise follows a positively skewed normal distribution. We further analyzed the noise structure distribution of a 3D mosaicked datacube in the range l = 40°.7 to 43°.3 and b = -2°.3 to 0°.3 and found that noise in the final mosaicked datacube is mainly characterized by noise fluctuation among the cells.

preprint2021arXiv

Testing the star formation scaling relations in the clumps of the North American and Pelican cloud complexes

The processes which regulate the star-formation within molecular clouds are still not well understood. Various star-formation scaling relations have been proposed to explain this issue by formulating a relation between star-formation rate surface density ($\rm Σ_{SFR}$) and the underlying gas surface density ($\rm Σ_{gas}$). In this work, we test various star formation scaling relations, such as Kennicutt-Schmidt relation, volumetric star-formation relation, orbital time model, crossing time model, and multi free-fall time scale model towards the North American and Pelican Nebulae complexes and in cold clumps associated with them. Measuring stellar-mass from young stellar objects and gaseous mass from CO measurements, we estimated mean $\rm Σ_{SFR}$, star formation rate per free-fall time, and star formation efficiency (SFE) for clumps to be 1.5 $\rm M{_\odot}~yr^{-1}~kpc^{-2}$, 0.009, 2.0$\%$, respectively, while for the entire NAN complex the values are 0.6 $\rm M{_\odot}~yr^{-1}~kpc^{-2}$, 0.0003, and 1.6$\%$, respectively. For clumps, we notice that the observed properties are in line with the correlation obtained between $\rm Σ_{SFR}$ and $\rm Σ_{gas}$, and between $\rm Σ_{SFR}$ and $\rm Σ_{gas}$ per free-fall time and orbital time for Galactic clouds. At the same time, we do not observe any correlation with $\rm Σ_{gas}$ per crossing time and multi free-fall time. Even though we see correlations in former cases, however, all models agree with each other within a factor of 0.5 dex, and discriminating between these models is not possible due to the current uncertainties in the input observables. We also test the variation of $\rm Σ_{SFR}$ versus the dense gas, but due to low statistics, a weak correlation is seen in our analysis.

preprint2020arXiv

In Search for Infall Motion in molecular clumps II: HCO+ (1-0) and HCN (1-0) Observations toward a Sub-sample of Infall Candidates

Gravitational accretion accumulates the original mass, and this process is crucial for us to understand the initial phases of star formation. Using the specific infall profiles in optically thick and thin lines, we searched the clumps with infall motion from the Milky Way Imaging Scroll Painting (MWISP) CO data in previous work. In this study, we selected 133 sources of them as a sub-sample for further research and identification. The excitation temperatures of these sources are between 7.0 and 38.5 K, while the H_2 column densities are between 10^21 and 10^23 cm^-2. We have observed optically thick lines HCO+ (1-0) and HCN (1-0) using the DLH 13.7-m telescope, and found 56 sources of them with blue profile and no red profile in these two lines, which are likely to have infall motions, with the detection rate of 42\%. It suggests that using CO data to restrict sample can effectively improve the infall detection rate. Among these confirmed infall sources, there are 43 associated with Class 0/I young stellar objects (YSOs), and 13 are not. These 13 sources are probably associated with the sources in earlier evolutionary stage. By comparison, the confirmed sources which are associated with Class 0/I YSOs have higher excitation temperatures and column densities, while the other sources are colder and have lower column densities. Most infall velocities of the sources we confirmed are between 10^-1 to 10^0 km s^-1, which is consistent with previous studies.

preprint2020arXiv

Local Molecular Gas toward the Aquila Rift Region

We present the results of a ~250 square degrees CO mapping (+26d<l<+50d and -5d<b<+5d) toward the Aquila Rift region at a spatial resolution of ~50&#34; and a grid spacing of 30&#34;. The high dynamic range CO maps with a spectral resolution of ~0.2km/s display highly structured molecular cloud (MC) morphologies with valuable velocity information, revealing complex spatial and dynamical features of the local molecular gas. In combination with the MWISP CO data and the Gaia DR2, distances of the main MC structures in the local ISM are well determined toward the Aquila Rift. We find that the total MC mass within 1 kpc is about >4.1x10^5 Msun in the whole region. In fact, the mass of the molecular gas is dominated by the W40 giant molecular cloud (GMC) at ~474 pc (~1.4x10^5 Msun) and the GMC complex G036.0+01.0 at ~560-670 pc (~2.0x10^5 Msun), while the MCs at ~220-260pc have gas masses of 10^2-10^3 Msun. Interestingly, an ~80pc long filamentary MC G044.0-02.5 at a distance of ~404 pc shows a systematic velocity gradient along and perpendicular to the major axis of the filament. The HI gas with the enhanced emission has the similar spatial morphologies and velocity features compared to the corresponding CO structure, indicating that the large-scale converging HI flows are probably responsible for the formation of the MC. Meanwhile, the long filamentary MC consists of many sub-filaments with the lengths ranging from ~0.5 pc to several pc, as well as prevalent networks of filaments in other large-scale local MCs.

preprint2020arXiv

Molecular Clouds Surrounding Supernova Remnant G43.9+1.6: Associated and Non-associated

Many supernova remnants (SNRs) are considered to evolve in molecular environments, but the associations between SNRs and molecular clouds (MCs) are often unclear. Being aware of such ambiguous case, we report our study on the molecular environment towards the SNR G43.9+1.6 by CO line observations. We investigated the correlations between the SNR and MCs at different velocities, and found two velocity components, i.e. $\sim$5 km s$^{-1}$ and $\sim$50 km s$^{-1}$ velocity components, showing spatial correlations with the remnant. However, no dynamical evidence of disturbance was found for the $\sim$5 km s$^{-1}$ velocity component. At the distance of the $\sim$5 km s$^{-1}$ velocity component, either near or far distance, the derived physical parameters are unreasonable too. We conclude that the SNR is not associated with the $\sim$5 km s$^{-1}$ velocity component, and their spatial correlation is just a chance correlation. For the $\sim$50 km s$^{-1}$ velocity component, dynamical evidence of disturbances, as well as the spatial correlation, indicate that it is associated with the SNR. We found that all the CO spectra extracted from the molecular clumps distributed along the border of the remnant are with broadened components presented, which can be fitted by Gaussian functions. By further analysis, we suggest that the SNR is at a near kinematic distance of about 3.1 kpc.

preprint2020arXiv

Searching for Molecular Outflows with Support Vector Machines: Dark Cloud Complex in Cygnus

We present a survey of molecular outflows across the dark cloud complex in the Cygnus region, based on 46.75 deg^2 field of CO isotopologues data from Milky Way Imaging Scroll Painting (MWISP) survey. A supervised machine learning algorithm, Support Vector Machine (SVM), is introduced to accelerate our visual assessment of outflow features in the data cube of 12CO and 13CO J = 1-0 emission. A total of 130 outflow candidates are identified, of which 77 show bipolar structures and 118 are new detections. Spatially, these outflows are located inside dense molecular clouds and some of them are found in clusters or in elongated linear structures tracing the underlying gas filament morphology. Along the line of sight, 97, 31, and 2 candidates reside in the Local, Perseus, and Outer arm, respectively. Young stellar objects as outflow drivers are found near most outflows, while 36 candidates show no associated source. The clusters of outflows that we detect are inhomogeneous in their properties; nevertheless, we show that the outflows cannot inject turbulent energy on cloud scales. Instead, at best, they are restricted to affecting the so called &#34;clump&#34; and &#34;core&#34; scales, and this only on short (~0.3 Myr) estimated timescales. Combined with outflow samples in the literature, our work shows a tight outflow mass-size correlation.