Catalog footprint

What is connected

113works
20topics
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

113 published item(s)

preprint2026arXiv

Evidence for a Damped Millisecond Quasi-Periodic Structure in a Fast Radio Burst

Fast radio bursts (FRBs) are millisecond-duration transients of unknown origin, likely associated with compact astrophysical objects. We report evidence for a damped millisecond quasi-periodic structure in a non-repeat FRB~20190122C. The burst consists of eight closely spaced radio pulses separated by $\sim$1 ms, with pulse amplitudes exhibiting an exponential decay starting from the brightest component. Combined Gaussian fitting and time-series analysis reveal a quasi-periodic oscillation (QPO) at $\sim$1 kHz. The observed QPO is consistent with damped magnetospheric oscillations. Assuming an Alfvén wave origin, we estimate a surface magnetic field of $\sim 10^{12}$ G and a characteristic spin period of $\sim$1 s, favoring a low-field magnetar or young neutron star scenario. The absence of frequency drift and the presence of exponential damping disfavor a merger-driven origin. These results suggest the first detection of an exponentially decaying QPO in any FRB, marking a rare detection of coherent oscillatory behavior in FRBs.

preprint2026arXiv

Large Language Models for Limited Noisy Data: A Gravitational Wave Identification Study

This work investigates whether large language models (LLMs) offer advantages over traditional neural networks for astronomical data processing, in regimes with non-Gaussian, non-stationary noise and limited labeled samples. Gravitational wave observations provide an suitable test case, using only 90 LIGO events, finetuned LLMs achieve 97.4\% accuracy for identifying signals. Further experiments show that, in contrast to traditional networks that rely on large simulated datasets, additional simulated samples do not improve LLM performance, while scaling studies reveal predictable gains with increasing model size and dataset size. These results indicate that LLMs can extract discriminative structure directly from observational data and provide an efficient assessment for gravitational wave identification. The same strategy may extend to other astronomical domains with similar noise properties, such as radio or pulsar observations.

preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

preprint2023arXiv

LGN-Net: Local-Global Normality Network for Video Anomaly Detection

Video anomaly detection (VAD) has been intensively studied for years because of its potential applications in intelligent video systems. Existing unsupervised VAD methods tend to learn normality from training sets consisting of only normal videos and regard instances deviating from such normality as anomalies. However, they often consider only local or global normality in the temporal dimension. Some of them focus on learning local spatiotemporal representations from consecutive frames to enhance the representation for normal events. But powerful representation allows these methods to represent some anomalies and causes miss detection. In contrast, the other methods are devoted to memorizing prototypical normal patterns of whole training videos to weaken the generalization for anomalies, which also restricts them from representing diverse normal patterns and causes false alarm. To this end, we propose a two-branch model, Local-Global Normality Network (LGN-Net), to simultaneously learn local and global normality. Specifically, one branch learns the evolution regularities of appearance and motion from consecutive frames as local normality utilizing a spatiotemporal prediction network, while the other branch memorizes prototype features of the whole videos as global normality by a memory module. LGN-Net achieves a balance of representing normal and abnormal instances by fusing local and global normality. In addition, the fused normality enables LGN-Net to generalize to various scenes more than exploiting single normality. Experiments demonstrate the effectiveness and superior performance of our method. The code is available online: https://github.com/Myzhao1999/LGN-Net.

preprint2022arXiv

Active Coding Piezoelectric Metasurfaces

The manipulation of acoustic waves plays an important role in a wide range of applications. Currently, acoustic wave manipulation typically relies on either acoustic metasurfaces or phased array transducers. The elements of metasurfaces are designed and optimized for a target frequency, which thus limits their bandwidth. Phased array transducers, suffering from high-cost and complex control circuits, are usually limited by the array size and the filling ratio of the control units. In this work, we introduce active coding piezoelectric metasurfaces; demonstrate commonly implemented acoustic wave manipulation functionalities such as beam steering, beam focusing and vortex beam focusing, acoustic tweezers; and eventually realize ultrasound imaging. The information coded on the piezoelectric metasurfaces herein is frequency independent and originates from the polarization directions, pointing either up or down, of the piezoelectric materials. Such a piezoelectric metasurface is driven by a single electrode and acts as a controllable active sound source, which combines the advantages of acoustic metasurfaces and phased array transducers while keeping the devices structurally simple and compact. Our coding piezoelectric metasurfaces can lead to potential technological innovations in underwater acoustic wave modulation, acoustic tweezers, biomedical imaging, industrial non-destructive testing and neural regulation.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): How do dense core properties affect the multiplicity of protostars?

During the transition phase from a prestellar to a protostellar cloud core, one or several protostars can form within a single gas core. The detailed physical processes of this transition, however, still remain unclear. We present 1.3 mm dust continuum and molecular line observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward 43 protostellar cores in the Orion Molecular Cloud Complex ($λ$ Orionis, Orion B, and Orion A) with an angular resolution of $\sim$ 0.35" ($\sim$ 140 au). In total, we detect 13 binary/multiple systems. We derive an overall multiplicity frequency (MF) of 28$\%$ $\pm$ 4$\%$ and a companion star fraction (CSF) of 51$\%$ $\pm$ 6$\%$, over a separation range of 300-8900 au. The median separation of companions is about 2100 au. The occurrence of stellar multiplicity may depend on the physical characteristics of the dense cores. Notably, those containing binary/multiple systems tend to show higher gas density and Mach number than cores forming single stars. The integral-shaped filament (ISF) of Orion A giant molecular cloud (GMC), which has the highest gas density and hosts high-mass star formation in its central region (the Orion Nebula cluster), shows the highest MF and CSF among the Orion GMCs. In contrast, the $λ$ Orionis Giant Molecular Cloud (GMC) has a lower MF and CSF than the Orion B and Orion A GMCs, indicating that feedback from HII regions may suppress the formation of multiple systems. We also find that the protostars comprising a binary/multiple system are usually at different evolutionary stages.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VIII. A search for hot cores by using C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH lines

Hot cores characterized by rich lines of complex organic molecules are considered as ideal sites for investigating the physical and chemical environments of massive star formation. We present a search for hot cores by using typical nitrogen- and oxygen-bearing complex organic molecules (C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH), based on ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS). The angular resolutions and line sensitivities of the ALMA observations are better than 2 arcsec and 10 mJy/beam, respectively. A total of 60 hot cores are identified with 45 being newly detected, in which the complex organic molecules have high gas temperatures ($>$ 100 K) and small source sizes ($<$ 0.1 pc). So far this is the largest sample of hot cores observed with similar angular resolution and spectral coverage. The observations have also shown nitrogen and oxygen differentiation in both line emission and gas distribution in 29 hot cores. Column densities of CH$_3$OH and CH$_3$OCHO increase as rotation temperatures rise. The column density of CH$_3$OCHO correlates tightly with that of CH$_3$OH. The pathways for production of different species are discussed. Based on the spatial position difference between hot cores and UC~H{\sc ii} regions, we conclude that 24 hot cores are externally heated while the other hot cores are internally heated. The observations presented here will potentially help establish a hot core template for studying massive star formation and astrochemistry.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation.

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Chemical variations across the TMC-1 boundary: molecular tracers from translucent phase to dense phase

We investigated the chemical evolutions of gas phase and grain surface species across the Taurus molecular cloud-1 (TMC-1) filament from translucent phase to dense phase. By comparing observations with modeling results from an up-to-date chemical network, we examined the conversion processes for the carbon-, oxygen-, nitrogen- and sulfur-bearing species, i.e.from their initial atomic form to their main molecular reservoir form both in the gas phase and on the grain surface. The conversion processes were found to depend on the species and A$_V$. The effect of initial carbon to oxygen elemental abundances ratio (C/O) by varying O on the chemistry was explored, and an initial carbon elemental abundance of 2.5 $\times$ 10$^{-4}$ and a C/O ratio of 0.5 could best reproduce the abundances of most observed molecules at TMC-1 CP, where more than 90 molecules have been identified. Based on the TMC-1 condition, we predicted a varied grain ice composition during the evolutions of molecular clouds, with H$_2$O ice as the dominant ice composition at A$_V$ $>$ 4 mag, CO$_2$ ice as the dominant ice composition at A$_V$ $<$ 4 mag, while CO ice severely decreased at A$_V$ around 4--5 mag.

preprint2022arXiv

Cold Gas in Massive Galaxies as A Critical Test of Black Hole Feedback Models

Black hole feedback has been widely implemented as the key recipe to quench star formation in massive galaxies in modern semi-analytic models and hydrodynamical simulations. As the theoretical details surrounding the accretion and feedback of black holes continue to be refined, various feedback models have been implemented across simulations, with notable differences in their outcomes. Yet, most of these simulations have successfully reproduced some observations, such as stellar mass function and star formation rate density in the local Universe. We use the recent observation on the change of neutral hydrogen gas mass (including both ${\rm H_2}$ and ${\rm HI}$) with star formation rate of massive central disc galaxies as a critical constraint of black hole feedback models across several simulations. We find that the predictions of IllustrisTNG agree with the observations much better than the other models tested in this work. This favors IllustrisTNG's treatment of active galactic nuclei - where kinetic winds are driven by black holes at low accretion rates - as more plausible amongst those we test. In turn, this also indirectly supports the idea that the massive central disc galaxy population in the local Universe was likely quenched by AGN feedback.

preprint2022arXiv

Compact and variable radio emission from an active galaxy with supersoft X-ray emission

RX J1301.9+2747 is a unique active galaxy with supersoft X-ray spectrum that lacks significant emission at energies above 2 keV. In addition, it is one of few galaxies displaying quasi-periodic X-ray eruptions that recur on a timescale of 13-20 ks. We present multi-epoch radio observations of RX J1301.9+2747 using GMRT, VLA and VLBA. The VLBA imaging at 1.6 GHz reveals a compact radio emission unresolved at a scale of <0.7 pc, with a brightness temperature of T_b>5x10^7 K. The radio emission is variable by more than a factor of 2.5 over a few days, based on the data taken from VLA monitoring campaigns. The short-term radio variability suggests that the radio emitting region has a size as small as 8x10^{-4} pc, resulting in an even higher brightness temperature of T_b ~10^{12} K. A similar limit on the source size can be obtained if the observed flux variability is not intrinsic and caused by the interstellar scintillation effect. The overall radio spectrum is steep with a time-averaged spectral index alpha=-0.78+/-0.03 between 0.89 GHz and 14 GHz. These observational properties rule out a thermal or star-formation origin of the radio emission, and appear to be consistent with the scenario of episodic jet ejections driven by magnetohydrodynamic process. Simultaneous radio and X-ray monitoring observations down to a cadence of hours are required to test whether the compact and variable radio emission is correlated with the quasi-periodic X-ray eruptions.

preprint2022arXiv

Deep simultaneous limits on optical emission from FRB 20190520B by 24.4 fps observations with Tomo-e Gozen

We conduct 24.4~fps optical observations of repeating Fast Radio Burst (FRB) 20190520B using Tomo-e Gozen, a high-speed CMOS camera mounted on the Kiso 105-cm Schmidt telescope, simultaneously with radio observations carried out using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We succeeded in the simultaneous optical observations of 11 radio bursts that FAST detected. However, no corresponding optical emission was found. The optical fluence limits as deep as 0.068 Jy ms are obtained for the individual bursts (0.029 Jy ms on the stacked data) corrected for the dust extinction in the Milky Way. The fluence limit is deeper than those obtained in the previous simultaneous observations for an optical emission with a duration $\gtrsim 0.1$ ms. Although the current limits on radio--optical spectral energy distribution (SED) of FRBs are not constraining, we show that SED models based on observed SEDs of radio variable objects such as optically detected pulsars, and a part of parameter spaces of theoretical models in which FRB optical emission is produced by inverse-Compton scattering in a pulsar magnetosphere or a strike of a magnetar blastwave into a hot wind bubble, can be ruled out once a similar fluence limit as in our observation is obtained for a bright FRB with a radio fluence $\gtrsim 5$ Jy ms.

preprint2022arXiv

Dependence of Chemical Abundance on the Cosmic Ray Ionization Rate in IC 348

Ions (e.g., H$_3^+$, H$_2$O$^+$) have been used extensively to quantify the cosmic-ray ionization rate (CRIR) in diffuse sightlines. However, measurements of CRIR in low-to-intermediate density gas environments are rare, especially when background stars are absent. In this work, we combine molecular line observations of CO, OH, CH, and HCO$^+$ in the star-forming cloud IC~348, and chemical models to constrain the value of CRIR and study the response of the chemical abundances distribution. The cloud boundary is found to have an $A_{\rm V}$ of approximately 4 mag. From the interior to the exterior of the cloud, the observed $^{13}$CO line intensities drop by an order of magnitude. The calculated average abundance of $^{12}$CO (assuming $^{12}$C/$^{13}$C = 65) is (1.2$\pm$0.9) $\times$10$^{-4}$, which increases by a factor of 6 from the interior to the outside regions. The average abundance of CH (3.3$\pm$0.7 $\times$ 10$^{-8}$) is in good agreement with previous findings in diffuse and translucent clouds ($A_{\rm V}$ $<$ 5 mag). However, we did not find a decline in CH abundance in regions of high extinction ($A_{\rm V}\simeq$8 mag) as previously reported in Taurus. By comparing the observed molecular abundances and chemical models, we find a decreasing trend of CRIR as $A_{\rm V}$ increases. The inferred CRIR of $ζ_{cr}$ = (4.7$\pm$1.5) $\times$ 10$^{-16}$ s$^{-1}$ at low $A_{\rm V}$ is consistent with H$^+_3$ measurements toward two nearby massive stars.

preprint2022arXiv

FAST observations of an extremely active episode of FRB 20201124A: II. Energy Distribution

We report the properties of more than 800 bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during an extremely active episode on UTC September 25-28, 2021 in a series of four papers. In this second paper of the series, we mainly focus on the energy distribution of the detected bursts. The event rate initially increased exponentially but the source activity stopped within 24 hours after the 4th day. The detection of 542 bursts in one hour during the fourth day marked the highest event rate detected from one single FRB source so far. The bursts have complex structures in the time-frequency space. We find a double-peak distribution of the waiting time, which can be modeled with two log-normal functions peaking at 51.22 ms and 10.05 s, respectively. Compared with the emission from a previous active episode of the source detected with FAST, the second distribution peak time is smaller, suggesting that this peak is defined by the activity level of the source. We calculate the isotropic energy of the bursts using both a partial bandwidth and a full bandwidth and find that the energy distribution is not significantly changed. We find that an exponentially connected broken-power-law function can fit the cumulative burst energy distribution well, with the lower and higher-energy indices being $-1.22\pm0.01$ and $-4.27\pm0.23$, respectively. Assuming a radio radiative efficiency of $η_r = 10^{-4}$, the total isotropic energy of the bursts released during the four days when the source was active is already $3.9\times10^{46}$ erg, exceeding $\sim 23\%$ of the available magnetar dipolar magnetic energy. This challenges the magnetar models invoking an inefficient radio emission (e.g. synchrotron maser models).

preprint2022arXiv

FAST observations of an extremely active episode of FRB 20201124A: III. Polarimetry

As the third paper in the multiple-part series, we report the statistical properties of radio bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio telescope (FAST) during an extremely active episode between the 25th and the 28th of September 2021 (UT). We focus on the polarisation properties of 536 bright bursts with $\mathrm{S/N}>50$. We found that the Faraday rotation measures (RMs) monotonically dropped from $-579 \ {\rm rad \ m^{-2}}$ to $-605 \ {\rm rad \ m^{-2}}$ in the 4-day window. The RM values were compatible with the values ($-300$ to $-900\ {\rm rad \ m^{-2}}$ ) reported 4 month ago (Xu et al. 2022). However, the RM evolution rate in the current observation window was at least an order of magnitude smaller than the one ($\sim 500\ {\rm rad \ m^{-2}\, day^{-1}}$) previously reported during the rapid RM-variation phase, but is still higher than the one ($\le 1\ {\rm rad \ m^{-2} day^{-1}}$ ) during the later RM no-evolution phase. The bursts of FRB 20201124A were highly polarised with the total degree of polarisation (circular plus linear) greater than 90% for more than 90\% of all bursts. The distribution of linear polarisation position angles (PAs), degree of linear polarisation ($L/I$), and degree of circular polarisation ($V/I$) can be characterised with unimodal distribution functions. During the observation window, the distributions became wider with time, i.e. with larger scatter, but the centroids of the distribution functions remained nearly constant. For individual bursts, significant PA variations (confidence level 5-$σ$) were observed in 33% of all bursts. The polarisation of single pulses seems to follow certain complex trajectories on the Poincaré sphere, which may shed light on the radiation mechanism at the source or the plasma properties along the path of FRB propagation.

preprint2022arXiv

FAST observations of an extremely active episode of FRB 20201124A: IV. Spin Period Search

We report the properties of more than 800 bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio telescope (FAST) during an extremely active episode on UTC September 25th-28th, 2021 in a series of four papers. In this fourth paper of the series, we present a systematic search of the spin period and linear acceleration of the source object from both 996 individual pulse peaks and the dedispersed time series. No credible spin period was found from this data set. We rule out the presence of significant periodicity in the range between 1 ms to 100 s with a pulse duty cycle $< 0.49\pm0.08$ (when the profile is defined by a von-Mises function, not a boxcar function) and linear acceleration up to $300$ m s$^{-2}$ in each of the four one-hour observing sessions, and up to $0.6$ m s$^{-2}$ in all 4 days. These searches contest theoretical scenarios involving a 1 ms to 100 s isolated magnetar/pulsar with surface magnetic field $<10^{15}$ G and a small duty cycle (such as in a polar-cap emission mode) or a pulsar with a companion star or black hole up to 100 M$_{\rm \odot}$ and $P_b>10$ hours. We also perform a periodicity search of the fine structures and identify 53 unrelated millisecond-timescale "periods" in multi-components with the highest significance of 3.9 $σ$. The "periods" recovered from the fine structures are neither consistent nor harmonically related. Thus they are not likely to come from a spin period. We caution against claiming spin periodicity with significance below $\sim$ 4 $σ$ with multi-components from one-off FRBs. We discuss the implications of our results and the possible connections between FRB multi-components and pulsar micro-structures.

preprint2022arXiv

Frequency-dependent polarization of repeating fast radio bursts-implications for their origin

The polarization of fast radio bursts (FRBs), bright astronomical transients, contains crucial information about their environments. We report polarization measurements of five repeating FRBs, the abundant signals of which enable wide-band observations with two telescopes. A clear trend of lower polarization at lower frequencies was found, which can be well characterized by a single parameter rotation-measure-scatter (σRM) and modeled by multi-path scatter. Sources with higher σRM have higher RM magnitude and scattering timescales. The two sources with the most substantial σRM, FRB 20121102A and FRB 20190520B, are associated with a compact persistent radio source. These properties indicate a complex environment near the repeating FRBs, such as a supernova remnant or a pulsar wind nebula, consistent with their arising from young populations.

preprint2022arXiv

J-comb: An image fusion algorithm to combine observations covering different spatial frequency ranges

Ground-based, high-resolution bolometric (sub)millimeter continuum mapping observations on spatially extended target sources are often subject to significant missing fluxes. This hampers accurate quantitative analyses. Missing flux can be recovered by fusing high-resolution images with observations that preserve extended structures. However, the commonly adopted image fusion approaches do not maintain the simplicity of the beam response function and do not try to elaborate the details of the yielded beam response functions. These make the comparison of the observations at multiple wavelengths not straightforward. We present a new algorithm, J-comb, which combines the high and low-resolution images linearly. By applying a taper function to the low-pass filtered image and combining it with a high-pass filtered image using proper weights, the beam response functions of our combined images are guaranteed to have near-Gaussian shapes. This makes it easy to convolve the observations at multiple wavelengths to share the same beam response functions. Moreover, we introduce a strategy to tackle the specific problem that the imaging at 850 um from the present-date ground-based bolometric instrument and that taken with the Planck satellite do not overlap in the Fourier domain. We benchmarked our method against two other widely-used image combination algorithms, CASA-feather and MIRIAD-immerge, with mock observations of star-forming molecular clouds. We demonstrate that the performance of the J-comb algorithm is superior to those of the other two algorithms. We applied the J-comb algorithm to real observational data of the Orion A star-forming region. We successfully produced dust temperature and column density maps with ~10" angular resolution, unveiling much greater details than the previous results.

preprint2022arXiv

Observation of Ultrafast Interfacial Exciton Formation and Recombination in Graphene/MoS2 Heterostructure

In this study,we combined time-resolved terahertz spectroscopy along with transient absorption spectroscopy to revisit the interlayer non-equilibrium carrier dynamics in largely lateral size Gr/MoS2 heterostructure fabricated with chemical vapor deposition method. Our experimental results reveal that, with photon-energy below the A-exciton of MoS2 monolayer, hot electrons transfer from graphene to MoS2 takes place in time scale of less than 0.5 ps, resulting in ultrafast formation of interfacial exciton in the heterostructure, subsequently, recombination relaxation of the interfacial exciton occurs in time scale of ~18 ps. A new model considering carrier heating and photogating effect in graphene is proposed to estimate the amount of carrier transfer in the heterostructure, which shows a good agreement with experimental result. Moreover, when the photon-energy is on-resonance with the A-exciton of MoS2, photogenerated holes in MoS2 are transferred to graphene layer within 0.5 ps, leading to the formation of interfacial exciton, the subsequent photoconductivity (PC) relaxation of graphene and bleaching recovery of A-exciton in MoS2 take place around ~10 ps time scale, ascribing to the interfacial exciton recombination. The faster recombination time of interfacial exciton with on-resonance excitation could come from the reduced interface barrier caused by bandgap renormalization effect. Our study provides deep insight into the understanding of interfacial charge transfer as well as the relaxation dynamics in graphene-based heterostructures, which are promising for the applications of graphene-based optoelectronic devices.

preprint2022arXiv

Physical Publicly Verifiable Randomness from Pulsars

We demonstrate how radio pulsars can be used as random number generators. Specifically, we focus on publicly verifiable randomness (PVR), in which the same sequence of trusted and verifiable random numbers is obtained by multiple parties. PVR is a critical building block for many processes and algorithms (including cryptography, scientific trials, electoral audits and international treaties). However, current approaches (based on number theory) may soon become vulnerable to quantum computers, motivating a growing demand for PVR based on natural physical phenomena. In this context, we explore pulsars as a potential physical PVR source. We first show that bit sequences extracted from the measured flux densities of a bright millisecond pulsar can pass standardised tests for randomness. We then quantify three illustrative methods of bit-extraction from pulsar flux density sequences, using simultaneous observations of a second pulsar carried out with the Parkes telescope in Australia and the Five-hundred-metre Aperture Spherical radio Telescope (FAST) in China, supported by numerical simulations. We demonstrate that the same bit sequence can indeed be obtained at both observatories, but the ubiquitous presence of radiometer noise needs to be accounted for when determining the expected bit error rate between two independent sequences. We discuss our results in the context of an imaginary use-case in which two mutually distrusting parties wish to obtain the same random bit sequence, exploring potential methods to mitigate against a malicious participant.

preprint2022arXiv

Radio detection of an elusive millisecond pulsar in the Globular Cluster NGC 6397

We report the discovery of a new 5.78 ms-period millisecond pulsar (MSP), PSR J1740-5340B (NGC 6397B), in an eclipsing binary system discovered with the Parkes radio telescope (now also known as Murriyang), Australia, and confirmed with the MeerKAT radio telescope in South Africa. The measured orbital period, 1.97 days, is the longest among all eclipsing binaries in globular clusters (GCs) and consistent with that of the coincident X-ray source U18, previously suggested to be a 'hidden MSP'. Our XMM-Newton observations during NGC 6397B's radio quiescent epochs detected no X-ray flares. NGC 6397B is either a transitional MSP or an eclipsing binary in its initial stage of mass transfer after the companion star left the main sequence. The discovery of NGC 6397B potentially reveals a subgroup of extremely faint and heavily obscured binary pulsars, thus providing a plausible explanation to the apparent dearth of binary neutron stars in core-collapsed GCs as well as a critical constraint on the evolution of GCs.

preprint2022arXiv

Revealing the dust grain polarization properties as a function of extinction and distance towards NGC 1893

Dust polarization observations at optical wavelengths help to understand the dust grain properties and trace the plane-of-the-sky component of the magnetic field. In this study, we make use of the $I$-band polarization data acquired from AIMPOL along with the distances ($d$) and extinction ($A_{V}$) data to study the variation of polarization fraction ($P$) as a function of $A_{V}$ and $d$ towards the star-forming region, NGC 1893. We employ a broken power-law fit and Bayesian analysis on extinction ($A_{V}$) versus polarization efficiency ($P$/$A_{V}$) and distance ($d$) versus rate of polarization ($P$/$d$). We find that $P$/$A_{V}$ shows a break at an extinction of $\sim$0.9 mag, whereas $P/d$ exhibits a break at a distance of $\sim$1.5 kpc. Based on these, we categorize the dust towards NGC 1893 into two populations: (i) foreground dust confined to $A_{V}$ $<$ $\sim$1 mag and distance up to $\sim$2 kpc and (ii) Perseus spiral arm dust towards NGC 1893 characterized with $A_{V}$ $>$ $\sim$1 mag and distance beyond $\sim$2 kpc. Foreground dust exhibits higher polarization efficiency but a lower polarization rate, whereas Perseus dust shows a lower polarization efficiency but a slightly higher polarization rate. Hence, we suggest that while polarization efficiency reveals the dust grain alignment, the rate of polarization infers about the distribution of dust grains towards NGC 1893. Further, we also shed a light on the spatial variation of intrinsic polarization and magnetic field orientation, and other parameters within the intra-cluster medium of NGC 1893.

preprint2022arXiv

Simulating high-time resolution radio-telescope observations

We describe a new software package for simulating channelised, high-time resolution data streams from radio telescopes. The software simulates data from the telescope and observing system taking into account the observation strategy, receiver system and digitisation. The signatures of pulsars, fast radio bursts and flare stars are modelled, including frequency-dependent effects such as scattering and scintillation. We also simulate more generic signals using spline curves and images. Models of radio frequency interference include signals from satellites, terrestrial transmitters and impulsive, broadband signals. The simulated signals can also be injected into real data sets. Uses of this software include the production of machine learning training data sets, development and testing of new algorithms to search for anomalous patterns and to characterise processing pipelines.

preprint2022arXiv

Temperature-driven Emergence of Negative Photoconductivity in Semimetal MoTe2 Film Probed with Terahertz Spectroscopy

Layered two-dimensional (2D) materials MoTe2 have been paid special attention due to the rich optoelectronic properties with various phases. The nonequilibrium carrier dynamics as well as its temperature dependence in MoTe2 are of prime importance, as it can shed light on understanding the anomalous optical response and potential applications in far infrared (IR) photodetection. Hereby, we employ time-resolved terahertz (THz) spectroscopy to study the temperature dependent nonequilibrium carrier dynamics in MoTe2 films. After photoexcitation of 1.59 eV, the 1T'-phase MoTe2 at high temperature behaves only THz positive photoconductivity (PPC) with relaxation time of less than 1 ps. In contrast, the Td-phase MoTe2 at low temperature shows ultrafast THz PPC initially followed by emerging THz negative photoconductivity (NPC), and the THz NPC signal relaxes to the equilibrium state in hundreds of ps time scale. Small polaron formation induced by hot carrier has been proposed to be ascribed to the THz NPC in the polar semimetal MoTe2 at low temperature. The polaron formation time after photoexcitation increases slightly with temperature, which is determined to be ~0.4 ps at 5 K and 0.5 ps at 100 K. Our experimental result demonstrates for the first time the dynamical formation of small poalron in MoTe2 Weyl semimetal, this is fundamental importance on the understanding the temperature dependent electron-phonon coupling and quantum phase transition, as well as the designing the MoTe2-based far IR photodetector.

preprint2022arXiv

Temporal Scattering, Depolarization, and Persistent Radio Emission from Magnetized Inhomogeneous Environments Near Repeating Fast Radio Burst Sources

Some repeating fast radio burst (FRB) sources exhibit complex polarization behaviors, including frequency-dependent depolarization, variation of rotation measure (RM), and oscillating spectral structures of polarized components. Very recently, Feng et al. (2022) reported that active repeaters exhibit conspicuous frequency-dependent depolarization and a strong correlation between RM scatter ($σ_{\rm RM}$) and the temporal scattering time ($τ_{\rm s}$), $σ_{\rm RM}\proptoτ_{\rm s}^{1.0\pm0.2}$, both of which can be well described by multi-path propagation through a magnetized inhomogeneous plasma screen. This observation strongly suggests that the temporal scattering and RM scatter originate from the same region. Besides, a particular finding of note in Feng et al. (2022) is that the FRBs with compact persistent radio sources (PRS) tend to have extreme $σ_{\rm RM}$. In this work, we focus on some theoretical predictions on the relations among temporal scattering, depolarization by RM scatter, and PRS contributed by the magnetized plasma environment close to a repeating FRB source. The behaviors of the RM scatter imply that the magnetized plasma environment is consistent with a supernova remnant or a pulsar wind nebula, and the predicted $σ_{\rm RM}$-$τ_{\rm s}$ relation is $σ_{\rm RM}\proptoτ_{\rm s}^{(0.54-0.83)}$ for different astrophysical scenarios. We further make a general discussion on PRS that does not depend on specific astrophysical scenarios. We show that the specific luminosity of a PRS should have a positive correlation with the RM contributed by the plasma screen. This is consistent with the observations of FRB 121102 and FRB 190520B.

preprint2022arXiv

The Distribution of UV Radiation Field in the Molecular Clouds of Gould Belt

The distribution of ultraviolet (UV) radiation field provides critical constraints on the physical environments of molecular clouds. Within 1 kpc of our solar system and fostering protostars of different masses, the giant molecular clouds in the Gould Belt present an excellent opportunity to resolve the UV field structure in star forming regions. We performed spectral energy distribution (SED) fitting of the archival data from the Herschel Gould Belt Survey (HGBS). Dust radiative transfer analysis with the DUSTY code were applied to 23 regions in 14 molecular complexes of the Gould Belt, resulting in the spatial distribution of radiation field in these regions. For 10 of 15 regions with independent measurements of star formation rate, their star formation rate and UV radiation intensity largely conform to a linear correlation found in previous studies.

preprint2022arXiv

The Potential of Detecting Radio-flaring Ultracool Dwarfs at L band in the FAST Drift-scan Survey

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) completed its commissioning and began the Commensal Radio Astronomy FasT Survey (CRAFTS), a multi-year survey to cover 60% of the sky, in 2020. We present predictions for the number of radio-flaring ultracool dwarfs (UCDs) that are likely to be detected by CRAFTS. Based on the observed flaring UCDs from a number of unbiased, targeted radio surveys in the literature, we derive a detection rate of $\ge$3%. Assuming a flat radio spectrum $νL _ν\propto ν^{β+1} $ with $β$ = -1.0 for UCD flares, we construct a flare luminosity function $d N/d L \propto L^{-1.96 \pm 0.45}$ (here $L=νL_ν$). CRAFTS is found to be sensitive enough for flares from UCDs up to $\sim$180 pc. Considering the Galactic thin disk, we carry out a 3D Monte Carlo simulation of the UCD population, which is then fed to mock CRAFTS observations. We estimate that $\sim$170 flaring UCDs would be detected through transient searches in circular polarization. Though only marginally sensitive to the scale height of UCDs, the results are very sensitive to the assumed spectral index $β$. For $β$ from 0 to -2.5, the number of expected detections increases dramatically from $\sim$20 to $\sim$3460. We also contemplate the strategies for following up candidates of flaring UCDs, and discuss the implications of survey results for improving our knowledge of UCD behavior at $L$ band and dynamos.

preprint2022arXiv

Topological superconductivity in a topological insulator

Topological superconductivity is an exotic quantum phenomenon with coupled nontrivial topological order and superconductivity together. A direct idea for producing topological superconductors is to create superconductivity based on the well recognized topological insulators. The topological insulating states in highly efficient thermoelectric materials Bi$_2$Te$_3$ and Bi$_2$Se$_3$ and their alloy Bi$_{2}$Te$_{3-x}$Se$_{x}$ have been established from angle-resolved photoemission and transport experiments. Superconductivity was also observed based on these popular topological insulators by the application of pressure chemical dopant and heterostructures. However, the experiments mainly focusing on Bi$_{2}$Se$_3$ doped by metals have not provided the consistent evidence to support the topological superconductivity. Here we carry out a systematic high-pressure study on a topological insulator Bi$_{2}$Te$_{2.7}$Se$_{0.3}$ to provide the convincing evidence for the expected topological superconductivity. Four phases with different structures are found upon compression. The topological surface state is identified in the entire initial phase, while superconductivity is found to coexist with such a state of the compressed material after its passing the electronic topological transition, followed by three other superconducting phases without topological character. For these superconducting phases, we observe that the upper critical field follows with the temperature in the critical exponent ${2/3}$ for the first one with the topological surface state and $1$ for the left. These observations support the realization of the topological superconductivity in the initial phase according to the theoretically proposed critical field measure. This work also points out a big pool and new direction for finding topological superconductors from topological thermoelectric materials.

preprint2022arXiv

Ultrafast Dynamics of Defect-Assisted Auger process in PdSe2 films: Synergistic Interaction Between Defect Trapping and Auger Effect

Strong Coulomb interactions in two-dimensional systems, together with quantum confinement, make many-body processes particularly effective for carrier dynamics, which plays a crucial role in determining carrier lifetime, photoconductivity, and emission yield of the materials. Hereby, by using optical pump and terahertz probe spectroscopy, we have investigated the photocarrier dynamics in the PdSe2 films with different thickness. The experimental results reveal that the photocarrier relaxation consists of two components: a fast component of 2.5 ps that shows the layer-thickness independence, and a slow component has typical lifetime of 7.3 ps decreasing with the layer thickness. Surprisingly, the relaxation times for both fast and slow components are exhibited both pump fluence and temperature independence, which suggests that synergistic interactions between defect trapping and Auger effect dominate the photocarrier dynamics in PdSe2 films. A model involving defect-assisted Auger process is proposed, which can reproduce the experimental results well. The fitting results reveal that the layer dependent lifetime is determined by the defect density rather than carrier occupancy rate after photoexcitation. Our results underscore the interplay between Auger process and defects in two-dimensional semiconductors.

preprint2021arXiv

A FAST Survey of HINSA in PGCCs Guided by HC3N

Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we search for HI narrow-line self-absorption (HINSA) features in twelve Planck Galactic cold clumps (PGCCs), one starless core L1521B and four star forming sources. Eight of the 12 PGCCs have emission of J=2-1 of cyanoacetylene (HC3N). With an improved HINSA extraction method more robust for weaker and blended features with high velocity resolution, the detection rates of HINSA in PGCCCs are high, at 92% overall (11/12) and 87% (7/8) among sources with HC3N J=2-1 emissions. Combining the data of molecular spectra and Planck continuum maps, we studied the morphologies, abundances and excitations of HI, CO and HC3N in PGCCs. The distribution of HINSA is similar to that of CO emission. HINSA tends to be not detected in regions associated with warm dust and background ionizing radiation, as well as regions associated with stellar objects. The abundances of HI in PGCCs are approximately 3E-4, and vary within a factor of ~3. The non-thermal velocity dispersions traced by C18O J=1-0 and HINSA are consistent with each other (0.1-0.4 km/s), larger than those of HC3N (~0.1 km/s). Carbon chain molecule abundant PGCCs provide a good sample to study HINSA.

preprint2021arXiv

A GPU based single-pulse search pipeline (GSP) with database and its application to the commensal radio astronomy FAST survey (CRAFTS)

We developed a GPU based single-pulse search pipeline (GSP) with candidate-archiving database. Largely based upon the infrastructure of Open source pulsar search and analysis toolkit (PRESTO), GSP implements GPU acceleration of the de-dispersion and integrates a candidate-archiving database. We applied GSP to the data streams from the commensal radio astronomy FAST survey (CRAFTS), which resulted in a quasi-real-time processing. The integrated candidate database facilitates synergistic usage of multiple machine-learning tools and thus improves efficient identification of radio pulsars such as rotating radio transients (RRATs) and Fast Radio Bursts (FRBs). We first tested GSP on pilot CRAFTS observations with the FAST Ultra-Wide Band (UWB) receiver. GSP detected all pulsars known from the the Parkes multibeam pulsar survey in the respective sky area covered by the FAST-UWB. GSP also discovered 13 new pulsars. We measured the computational efficiency of GSP to be ~120 times faster than the original PRESTO and ~60 times faster than a MPI-parallelized version of PRESTO.

preprint2021arXiv

ALMA observations of NGC 6334S. II. Subsonic and Transonic Narrow Filaments in a High-mass Star Formation Cloud

We present a study of narrow filaments toward a massive infrared dark cloud, NGC 6334S, using the Atacama Large Millimeter/submillimeter Array (ALMA). Thirteen gas filaments are identified using the H$^{13}$CO$^{+}$ line, while a single continuum filament is revealed by the continuum emission. The filaments present a compact radial distribution with a median filament width of $\sim$0.04 pc narrower than the previously proposed `quasi-universal' 0.1~pc filament width. The higher spatial resolution observations and higher-density gas tracer tend to identify even narrower and lower mass filaments. The filament widths are roughly twice the size of embedded cores. The gas filaments are largely supported by thermal motions. The nonthermal motions are predominantly subsonic and transonic in both identified gas filaments and embedded cores, which may imply that stars are likely born in environments of low turbulence. A fraction of embedded objects show a narrower velocity dispersion compared with their corresponding natal filaments, which may indicate that the turbulent dissipation is taking place in these embedded cores. The physical properties (mass, mass per unit length, gas kinematics, and width) of gas filaments are analogous to those of narrow filaments found in low- to high-mass star-forming regions. The more evolved sources are found to be farther away from the filaments, a situation that may have resulted from the relative motions between the YSOs and their natal filaments.

preprint2021arXiv

Molecular Cloud Cores with High Deuterium Fractions: Nobeyama Mapping Survey

We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$μ$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$μ$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion.

preprint2021arXiv

Studying Infall in Infrared Dark Clouds with Multiple HCO+ Transitions

We investigate the infall properties in a sample of 11 infrared dark clouds (IRDCs) showing blue-asymmetry signatures in HCO$^{+}$ J=1--0 line profiles. We used JCMT to conduct mapping observations in HCO$^{+}$ J=4--3 as well as single-pointing observations in HCO$^{+}$ J =3--2, towards 23 clumps in these IRDCs. We applied the HILL model to fit these observations and derived infall velocities in the range of 0.5-2.7 km s$^{-1}$, with a median value of 1.0 km s$^{-1}$, and obtained mass accretion rates of 0.5-14$\times$10$^{-3}$ Msun yr$^{-1}$. These values are comparable to those found in massive star forming clumps in later evolutionary stages. These IRDC clumps are more likely to form star clusters. HCO$^{+}$ J =3--2 and HCO$^{+}$ J =1--0 were shown to trace infall signatures well in these IRDCs with comparable inferred properties. HCO$^{+}$ J=4--3, on the other hand, exhibits infall signatures only in a few very massive clumps, due to smaller opacties. No obvious correlation for these clumps was found between infall velocity and the NH3/CCS ratio.

preprint2021arXiv

The first evidence for three-dimensional spin-velocity alignment in pulsars

More than 50 years after the discovery of pulsars and confirmation of their association with supernova explosions, the origin of the initial spin and velocity of pulsars remains largely a mystery. The typical space velocities of several hundred km/s have been attributed to "kicks" resulting from asymmetries either in the supernova ejecta or in the neutrino emission. Observations have shown a strong tendency for alignment of the pulsar space velocity and spin axis in young pulsars but, up to now, these comparisons have been restricted to two dimensions. We report here the first evidence for three-dimensional alignment between the spin and velocity vectors, largely based on observations made with the Five-hundred-meter Aperture Spherical radio Telescope of the pulsar PSR~J0538+2817 and its associated supernova remnant S147. Analysis of these and related observations has enabled us to determine the location of the pulsar within the supernova remnant and hence its radial velocity. Current simulations of supernova explosions have difficulty producing such three-dimensional alignment. Our results, which depend on the unprecedented sensitivity of the new observations, add another dimension to the intriguing correlation between pulsar spin-axis and birth-kick directions, thus deepening the mysteries surrounding the birth of neutron stars.

preprint2021arXiv

The Stellar "Snake" I: Whole Structure and Properties

To complement our previous discovery of the young snake-like structure in the solar neighborhood and reveal the structure's full extent, we build two samples of stars within the Snake and its surrounding territory from {\tt Gaia EDR3}. With the friends-of-friends algorithm, we identify 2694 and 9615 Snake member candidates from the two samples. Thirteen open clusters are embedded in these member candidates. By combining the spectroscopic data from multiple surveys, we investigate the comprehensive properties of the candidates and find that they \thj{are very likely to} belong to one sizable structure, since most of the components are well bridged in their spatial distributions, and follow a single stellar population with an age of $30-40$\,Myr and solar metallicity. This sizable structure is best explained as hierarchically primordial, and probably formed from a filamentary giant molecular cloud with unique formation history in localized regions. To analyze the dynamics of the Snake, we divide the structure into five groups according to their tangential velocities; we find that the groups are expanding at a coherent rate ($κ_X\sim3.0\,\times10^{-2}\,\rm km\,s^{-1}\,pc^{-1}$) along the length of the structure ($X$-direction). \thj{The corresponding expansion age ($τ\sim33$\,Myr) is highly consistent with the age of the Snake}. With over ten thousand member stars, the Snake is an ideal laboratory to study nearby coeval stellar formation, stellar physics, and environmental evolution over a large spatial extent.

preprint2021arXiv

The TMRT K Band Observations towards 26 Infrared Dark Clouds: NH$_{3}$, CCS, and HC$_{3}$N

We present one of the first Shanghai Tian Ma Radio Telescope (TMRT) K Band observations towards a sample of 26 infrared dark clouds (IRDCs). We observed the (1,1), (2,2), (3,3), and (4,4) transitions of NH$_{3}$ together with CCS (2$_{1}$-1$_{0}$) and HC$_{3}$N $J\,$=2-1, simultaneously. The survey dramatically increases the existing CCS-detected IRDC sample from 8 to 23, enabling a better statistical study of the ratios of carbon-chain molecules (CCM) to N-bearing molecules in IRDCs. With the newly developed hyperfine group ratio (HFGR) method of fitting NH$_{3}$ inversion lines, we found the gas temperature to be between 10 and 18 K. The column density ratios of CCS to NH$_{3}$ for most of the IRDCs are less than 10$^{-2}$, distinguishing IRDCs from low-mass star-forming regions. We carried out chemical evolution simulations based on a three-phase chemical model NAUTILUS. Our measurements of the column density ratios between CCM and NH$_{3}$ are consistent with chemical evolutionary ages of $\lesssim$10$^{5}$ yr in the models. Comparisons of the data and chemical models suggest that CCS, HC$_{3}$N, and NH$_{3}$ are sensitive to the chemical evolutionary stages of the sources.

preprint2020arXiv

A Fast Radio Burst discovered in FAST drift scan survey

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

preprint2020arXiv

ALMA ACA and Nobeyama observations of two Orion cores in deuterated molecular lines

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- I. Survey description and a first look at G9.62+0.19

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- II. Compact objects in ACA observations and star formation scaling relations

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

preprint2020arXiv

Constraints on individual supermassive binary black holes using observations of PSR J1909-3744

We perform a search for gravitational waves (GWs) from several supermassive binary black hole (SMBBH) candidates (NGC 5548, Mrk 231, OJ 287, PG 1302-102, NGC 4151, Ark 120 and 3C 66B) in long-term timing observations of the pulsar PSR J1909$-$3744 obtained using the Parkes radio telescope. No statistically significant signals were found. We constrain the chirp masses of those SMBBH candidates and find the chirp mass of NGC 5548 and 3C 66B to be less than $2.4 \times 10^9\,\rm M_{\odot}$ and $2.5 \times 10^9\,\rm M_{\odot}$ (with 95% confidence), respectively. Our upper limits remain a factor of 3 to 370 above the likely chirp masses for these candidates as estimated from other approaches. The observations processed here provide upper limits on the GW strain amplitude that improve upon the results from the first Parkes Pulsar Timing Array data release by a factor of 2 to 7. We investigate how information about the orbital parameters can help improve the search sensitivity for individual SMBBH systems. Finally, we show that these limits are insensitive to uncertainties in the Solar System ephemeris model.

preprint2020arXiv

Discovery and timing of pulsars in the globular cluster M13 with FAST

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

preprint2020arXiv

Electrokinetically enhanced cross-stream particle migration in viscoelastic flows

Advancements in understanding the lateral migration of particles have helped in enhanced separation in microfluidic devices. In this work, we investigate the effects of electrokinetics on the particle migration in a viscoelastic flow, where the electric field is applied parallel to the flow. Through experiments and theoretical analysis, we show that the interaction of electrokinetic and rheological effects can result in an enhancement in the migration by an order of magnitude. The theoretical analysis, in agreement with experiments, demonstrates that the particles can be focused at different equilibrium positions based on their intrinsic electrical properties.

preprint2020arXiv

Filament Intersections and Cold Dense Cores in Orion A North

We studied the filament structures and dense cores in OMC-2,3 region in Orion A North molecular cloud using the high-resolution N2H+ (1-0) spectral cube observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). The filament network over a total length of 2 pc is found to contain 170 intersections and 128 candidate dense cores. The dense cores are all displaced from the infrared point sources (possible young stars), and the major fraction of cores (103) are located around the intersections. Towards the intersections, there is also an increasing trend for the total column density Ntot as well as the the power-law index of the column-density Probability Distribution Function (N-PDF), suggesting that the intersections would in general have more significant gas assembly than the other part of the filament paths. The virial analysis shows that the dense cores mostly have virial mass ratio of alpha_vir=M_vir/M_gas<1.0, suggesting them to be bounded by the self gravity. In the mean time, only about 23 percent of the cores have critical mass ratio of alpha_crit=M_crit/M_gas<1.0, suggesting them to be unstable against core collapse. Combining these results, it shows that the major fraction of the cold starless and possible prestellar cores in OMC-2,3 are being assembled around the intersections, and currently in a gravitationally bound state. But more extensive core collapse and star formation may still require continuous core-mass growth or other perturbatio

preprint2020arXiv

First SETI Observations with China's Five-hundred-meter Aperture Spherical radio Telescope (FAST)

The Search for Extraterrestrial Intelligence (SETI) attempts to address the possibility of the presence of technological civilizations beyond the Earth. Benefiting from high sensitivity, large sky coverage, an innovative feed cabin for China's Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed the SETI first observations with FAST's newly commisioned 19-beam receiver; we report preliminary results in this paper. Using the data stream produced by the SERENDIP VI realtime multibeam SETI spectrometer installed at FAST, as well as its off-line data processing pipelines, we identify and remove four kinds of radio frequency interference(RFI): zone, broadband, multi-beam, and drifting, utilizing the Nebula SETI software pipeline combined with machine learning algorithms. After RFI mitigation, the Nebula pipeline identifies and ranks interesting narrow band candidate ET signals, scoring candidates by the number of times candidate signals have been seen at roughly the same sky position and same frequency, signal strength, proximity to a nearby star or object of interest, along with several other scoring criteria. We show four example candidates groups that demonstrate these RFI mitigation and candidate selection. This preliminary testing on FAST data helps to validate our SETI instrumentation techniques as well as our data processing pipeline.

preprint2020arXiv

From haloes to galaxies -- II. The fundamental relations in star formation and quenching

Star formation and quenching are two of the most important processes in galaxy formation and evolution. We explore in the local Universe the interrelationships among key integrated galaxy properties, including stellar mass $M_*$, star formation rate (SFR), specific SFR (sSFR), molecular gas mass $M_{\rm H_2}$, star formation efficiency (SFE) of the molecular gas and molecular gas to stellar mass ratio $μ$. We aim to identify the most fundamental scaling relations among these key galaxy properties and their interrelationships. We show the integrated $M_{\rm H_2}$-SFR, SFR-$M_*$ and $M_{\rm H_2}$-$M_*$ relation can be simply transformed from the $μ$-sSFR, SFE-$μ$ and SFE-sSFR relation, respectively. The transformation, in principle, can increase or decrease the scatter of each relation. Interestingly, we find the latter three relations all have significantly smaller scatter than the former three corresponding relations. We show the probability to achieve the observed small scatter by accident is extremely close to zero. This suggests that the smaller scatters of the latter three relations are driven by a more fundamental physical connection among these quantities. We then show the large scatters in the former relations are due to their systematic dependence on other galaxy properties, and on star formation and quenching process. We propose the sSFR-$μ$-SFE relation as the Fundamental Formation Relation (FFR), which governs the star formation and quenching process, and provides a simple framework to study galaxy evolution. Other scaling relations, including integrated Kennicutt-Schmidt law, star-forming main sequence and molecular gas main sequence, can all be derived from the FFR.

preprint2020arXiv

H I absorption towards radio Active Galactic Nuclei of different accretion modes

We present results of H I absorption experiment done using the Giant Metrewave Radio Telescope (GMRT) towards 27 low- and intermediate-luminosity ($P_{\rm 1.4 GHz}$ $\sim$10$^{23}$-10$^{26}$ W Hz$^{-1}$) radio active galactic nuclei (AGN), classified as either low excitation radio galaxies (LERGs) or high excitation radio galaxies (HERGs) and with WISE colour W2[4.6 $μ$m]$-$W3[12 $μ$m]$>$ 2. We report H I absorption detection towards seven radio AGNs, six of which are new. Combined with other sources from literature classified as LERGs or HERGs, we confirm our earlier result that compact radio AGNs with WISE colour W2$-$W3 $>$2 have higher detection rates compared to those with W2$-$W3 $<$2. We find that H I absorption detection rate is higher for HERGs (37.0$^{+15.8}_{-11.5}$ per cent) compared to LERGs (22.0$^{+3.9}_{-3.4}$ per cent), mainly due to a larger fraction of HERGs being gas and dust rich with a younger stellar population compared to LERGs. However, for similar compact radio structures and host galaxies with WISE colours W2$-$W3$>$2, we don't find any significant difference in detection rates of two types of AGNs implying detection of H I gas may not necessarily mean high excitation mode AGN. We further analysed the kinematics towards these sources. We find that while LERGs show a wide range in the shift of centroid velocities ($\sim$ $-$479 to $+$356 km s$^{-1}$) relative to the optical systemic velocity, most of the HERGs have centroid velocity shift less than 200 km s$^{-1}$, possibly due to differences in jet-interstellar medium interaction.

preprint2020arXiv

Molecular Cloud Cores with High Deuterium Fraction: Nobeyama Single-Pointing Survey

We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($λ$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the Nobeyama 45-m telescope. We find that early-type molecules (e.g., CCS) have low detection rates and that late-type molecules (e.g., N$_2$H$^+$, c-C$_3$H$_2$) and deuterated molecules (e.g., N$_2$D$^+$, DNC) have high detection rates, suggesting that most of the cores are chemically evolved. The deuterium fraction (D/H) is found to decrease with increasing distance, indicating that it suffers from differential beam dilution between the D/H pair of lines for distant cores ($>$1 kpc). For $λ$ Orionis, Orion A, and B located at similar distances, D/H is not significantly different, suggesting that there is no systematic difference in the observed chemical properties among these three regions. We identify at least eight high D/H cores in the Orion region and two at high latitudes, which are most likely to be close to the onset of star formation. There is no clear evidence of the evolutionary change in turbulence during the starless phase, suggesting that the dissipation of turbulence is not a major mechanism for the beginning of star formation as judged from observations with a beam size of 0.04 pc.

preprint2020arXiv

Molecular Oxygen in the nearest QSO Mrk 231

We report the detection of an emission feature at the 12 sigma level with FWHM line width of about 450 km/s toward the nearest quasi-stellar object, QSO Mrk 231. Based on observations with the IRAM 30 m telescope and the NOEMA Interferometer, the 11-10 transition of molecular oxygen is the likely origin of line with rest frequency close to 118.75 GHz. The velocity of the O2 emission in Mrk 231 coincides with the red wing seen in CO emission, suggesting that it is associated with the outflowing molecular gas, located mainly at about ten kpc away from the central AGN. This first detection of extragalactic molecular oxygen provides an ideal tool to study AGN-driven molecular outflows on dynamic time scales of tens of Myr. O2 may be a significant coolant for molecular gas in such regions affected by AGN-driven outflows. New astrochemical models are needed to explain the implied high molecular oxygen abundance in such regions several kpc away from the center of galaxies.

preprint2020arXiv

Observational features of exoplanetary synchrotron radio bursts

Magnetic fields of exoplanets are important in shielding the planets from cosmic rays and interplanetary plasma. Due to the interaction with the electrons from their host stars, the exoplanetary magnetospheres are predicted to have both cyclotron and synchrotron radio emissions, of which neither has been definitely identified in observations yet. As the coherent cyclotron emission has been extensively studied in literatures, here we focus on the planetary synchrotron radiation with bursty behaviors (i.e., radio flares) caused by the outbreaks of energetic electron ejections from the host star. Two key parameters of the bursty synchrotron emissions, namely the flux density and burst rate, and two key features namely the burst light curve and frequency shift, are predicted for star - hot Jupiter systems. The planetary orbital phase - burst rate relation is also considered as the signature of star-planet interactions (SPI). As examples, previous X-ray and radio observations of two well studied candidate systems, HD 189733 and V830 τ, are adopted to predict their specific burst rates and fluxes of bursty synchrotron emissions for further observational confirmations. The detectability of such emissions by current and upcoming radio telescopes shows that we are at the dawn of discoveries.

preprint2020arXiv

Opportunities to Search for Extra-Terrestrial Intelligence with the Five-hundred-meter Aperture Spherical radio Telescope

The discovery of ubiquitous habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, has ushered in a renaissance in the search for extra-terrestrial intelligence (SETI). Large scale SETI activities are now underway at numerous international facilities. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is the largest single-aperture radio telescope in the world, well positioned to conduct sensitive searches for radio emission indicative of exo-intelligence. SETI is one of the five key science goals specified in the original FAST project plan. A collaboration with the Breakthrough Listen Initiative has been initiated in 2016 with a joint statement signed both by Dr. Jun Yan, the then director of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), and Dr. Peter Worden, the Chairman of the Breakthrough Prize Foundation. In this paper, we highlight some of the unique features of FAST that will allow for novel SETI observations. We identify and describe three different signal types indicative of a technological source, namely, narrow-band, wide-band artificially dispersed, and modulated signals. We here propose observations with FAST to achieve sensitivities never before explored.

preprint2020arXiv

Readout for Kinetic-Inductance-Detector-Based Submillimeter Radio Astronomy

A substantial amount of important scientific information is contained within astronomical data at the submillimeter and far-infrared (FIR) wavelengths, including information regarding dusty galaxies, galaxy clusters, and star-forming regions; however, these wavelengths are among the least-explored fields in astronomy because of the technological difficulties involved in such research. Over the past 20 years, considerable efforts have been devoted to developing submillimeter- and millimeter-wavelength astronomical instruments and telescopes. The number of detectors is an important property of such instruments and is the subject of the current study. Future telescopes will require as many as hundreds of thousands of detectors to meet the necessary requirements in terms of the field of view, scan speed, and resolution. A large pixel count is one benefit of the development of multiplexable detectors that use kinetic inductance detector (KID) technology. This paper presents the development of all aspects of the readout electronics for a KID-based instrument, which enabled one of the largest detector counts achieved to date in submillimeter-/millimeter-wavelength imaging arrays: a total of 2304 detectors. The work presented in this paper had been implemented in the MUltiwavelength Submillimeter Inductance Camera (MUSIC), a instrument for the Caltech Submillimeter Observatory (CSO) between 2013 and 2015.

preprint2020arXiv

Resolution-dependent Subsonic Non-thermal Line Dispersion Revealed by ALMA

We report here Atacama Large Millimeter/submillimeter Array (ALMA) N$_2$H$^+$ (1-0) images of the Orion Molecular Cloud 2 and 3 (OMC-2/3) with high angular resolution (3'' or 1200 au) and high spatial dynamic range. Combining dataset from the ALMA main array, ALMA Compact Array (ACA), the Nobeyama 45m Telescope, and the JVLA (providing temperature measurement on matching scales), we find that most of the dense gas in OMC-2/3 is subsonic ($\rm σ_{NT}/c_{s}$ = 0.62) with a mean line width ($Δ\upsilon$) of 0.39 km s$^{-1}$ FWHM. This is markedly different from the majority of previous observations of massive star-forming regions. In contrast, line widths from the Nobeyama Telescope are transonic at 0.69 km s$^{-1}$ ($\rm σ_{NT}/c_{s}$ = 1.08). We demonstrated that the larger line widths obtained by the single-dish telescope arose from unresolved sub-structures within their respective beams. The dispersions from larger scales $σ_{ls}$ (as traced by the Nobeyama Telescope) can be decomposed into three components $\rm σ_{ls}^2 = σ_{ss}^2+ σ_{bm}^2+ σ_{rd}^2$, where small-scale $σ_{ss}$ is the line dispersion of each ALMA beam, bulk motion $σ_{bm}$ is dispersion between peak velocity of each ALMA beam, and $σ_{rd}$ is the residual dispersion. Such decomposition, though purely empirical, appears to be robust throughout our data cubes. Apparent supersonic line widths, commonly found in massive molecular clouds, are thus likely due to the effect of poor spatial resolution. The observed non-thermal line dispersion (sometimes referred to as 'turbulence') transits from supersonic to subsonic at $\sim 0.05$ pc scales in OMC-2/3 region. Such transition could be commonly found with sufficient spatial (not just angular) resolution, even in regions with massive young clusters, such as Orion molecular clouds studied here.

preprint2020arXiv

Rotation of Two Micron All Sky Survey Clumps in Molecular Clouds

We have analyzed the rotational properties of 12 clumps using $^{13}$CO (1--0) and C$^{18}$O (1--0) maps of the Five College Radio Astronomy Observatory 13.7 m radio telescope. The clumps, located within molecular clouds, have radii ($R$) in the range of 0.06 -- 0.27\,pc. The direction of clump elongation is not correlated with the direction of the velocity gradient. We measured the specific angular momentum (J/M) to be between 0.0022 and 0.025 pc\,km\,s$^{-1}$ based on $^{13}$CO images, and between 0.0025 and 0.021 pc\,km\,s$^{-1}$ based on C$^{18}$O images. The consistency of $J/M$ based on different tracers indicates the $^{13}$CO and C$^{18}$O in dense clumps trace essentially the same material despite significantly different opacities. We also found that $J/M$ increases monotonically as a function of $R$ in power--law form, $J/M~\propto~R^{1.58~\pm~0.11}$. The ratio between rotation energy and gravitational energy, $β$, ranges from 0.0012 to 0.018. The small values of $β$ imply that rotation alone is not sufficient to support the clump against gravitational collapse.

preprint2020arXiv

Supermassive Binary Black Hole Evolution can be traced by a small SKA Pulsar Timing Array

Supermassive black holes are commonly found in the center of galaxies and evolve with their hosts. The supermassive binary black holes (SMBBH) are thus expected to exist in close galaxy pairs, however, none has been unequivocally detected. The square kilometre array (SKA) is a multi-purpose radio telescope with a collecting area approaching 1 million square metres, with great potential for detecting nanohertz gravitational waves (GWs). In this paper, we quantify the GW detectability by SKA for a realistic SMBBH population using pulsar timing array (PTA) technique and quantify its impact on revealing SMBBH evolution with redshift for the first time. With only $\sim20$ pulsars, much smaller a requirement than in previous work, the SKA PTA is expected to obtain detection within about 5 years of operation and to achieve a detection rate of more than 100 SMBBHs/yr after about 10 years. Although beyond the scope of this paper, we must acknowledge that the presence of persistent red noise will reduce the number of expected detections here. It is thus imperative to understand and mitigate red noise in the PTA data. The GW signatures from a few well-known SMBBH candidates, such as OJ 287, 3C 66B, NGC 5548 and Ark 120, will be detected given the currently best-known parameters of each system. Within 30 years of operation, about 60 individual SMBBHs detection with $z<0.05$ and more than $10^4$ with $z<1$ are expected. The detection rate drops precipitately beyond $z=1$. The substantial number of expected detections and their discernible evolution with redshift by SKA PTA will make SKA a significant tool for studying SMBBHs.

preprint2020arXiv

The FAST discovery of an Eclipsing Binary Millisecond Pulsar in the Globular Cluster M92 (NGC 6341)

We report the discovery of an eclipsing binary millisecond pulsar in the globular cluster M92 (NGC6341) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J1717+4308A, or M92A, has a pulse frequency of 316.5~Hz (3.16~ms) and a dispersion measure of 35.45 pc cm$^{-3}$. The pulsar is a member of a binary system with an orbital period of 0.20~days around a low-mass companion which has a median mass of $\sim$0.18~\Ms. From observations so far, at least two eclipsing events have been observed in each orbit. The longer one lasted for ~5000~s in the orbital phase range 0.1--0.5. The other lasted for ~500~s and occurred between 1000--2000~s before or after the longer eclipsing event. The lengths of these two eclipsing events also change. These properties suggest that J1717+4308A is a ``red-back'' system with a low-mass main sequence or sub-giant companion. Timing observations of the pulsar and further searches of the data for additional pulsars are ongoing.

preprint2020arXiv

The Fundamental Performance of FAST with 19-beam Receiver at L Band

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) passed national acceptance and is taking pilot cycle of 'Shared-Risk' observations. The 19-beam receiver covering 1.05-1.45 GHz was used for most of these observations. The electronics gain fluctuation of the system is better than 1\% over 3.5 hours, enabling enough stability for observations. Pointing accuracy, aperture efficiency and system temperature are three key parameters of FAST. The measured standard deviation of pointing accuracy is 7.9$''$, which satisfies the initial design of FAST. When zenith angle is less than 26.4$^\circ$, the aperture efficiency and system temperature around 1.4 GHz are $\sim$ 0.63 and less than 24 K for central beam, respectively. The measured value of these two parameters are better than designed value of 0.6 and 25 K, respectively. The sensitivity and stability of the 19-beam backend are confirmed to satisfy expectation by spectral HI observations toward N672 and polarization observations toward 3C286. The performance allows FAST to take sensitive observations in various scientific goals, from studies of pulsar to galaxy evolution.

preprint2020arXiv

The M31/M33 tidal interaction: A hydrodynamic simulation of the extended gas distribution

We revisit the orbital history of the Triangulum galaxy (M33) around the Andromeda galaxy (M31) in view of the recent Gaia Data Release 2 proper motion measurements for both Local Group galaxies. Earlier studies consider highly idealised dynamical friction, but neglect the effects of dynamical mass loss. We show the latter process to be important using mutually consistent orbit integration and N-body simulations. Following this approach we find an orbital solution that brings these galaxies to within ~50 kpc of each other in the past, ~6.5 Gyr ago. We explore the implications of their interaction using an N-body/hydrodynamical simulation with a focus on the origin of two prominent features: 1) M31's Giant Stellar Stream; and 2) the M31-M33 HI filament. We find that the tidal interaction does not produce a structure reminiscent of the stellar stream that survives up to the present day. In contrast, the M31-M33 HI filament is likely a fossil structure dating back to the time of the ancient encounter between these galaxies. Similarly, the observed outer disc warp in M33 may well be a relic of this past event. Our model suggests the presence of a tidally induced gas envelope around these galaxies, and the existence of a diffuse gas stream, the 'Triangulum stream', stretching for tens of kpc from M33 away from M31. We anticipate upcoming observations with the recently commissioned, Five-hundred-meter Aperture Spherical radio Telescope (FAST) that will target the putative stream in its first years of operation.

preprint2020arXiv

Toward a direct measurement of the cosmic acceleration: roadmap and forecast on FAST

HI absorption systems are great targets for direct measurement of the Sandage-Loeb (SL) effect throughout a wide range of redshift for ground-based radio telescopes. We demonstrate the significance of improving the accuracy of SL effect measurement in cosmological model selection. With its wide sky coverage and high sensitivity, we forecast that for 1 year of the upcoming commensal survey (CRAFTS) the FAST telescope is capable of discovering about 800 HI absorption systems thereby improving the SL measurement accuracy. Aiming to measurement the cosmic redshift drift rate at the precision of $\dot{z} \sim 10^{-10} \mathrm{decade^{-1}}$, we propose combined observation mode with blind-searching and targeted observation. For a decade of consecutive targeted spectroscopic observation with the frequency resolution at a level of sub-$0.1\ \rm Hz$, we could detect the first-order derivative of the cosmological redshift with the expected precision.

preprint2020arXiv

Toward Eurasian SubMillimeter Telescopes: the concept of multicolor subTHz MKID-array demo camera MUSICAM and its instrumental testing

New challenges in submillimeter wave astronomy require instruments with a combination of high sensitivity and angular resolution, wide field of view and multiwave (multicolor) spectral range. New large single mm/submm telescopes are in high demand, as well as their inclusion in the global Event Horizon Telescope (EHT) VLBI network. At the same time, there are no large mm/submm telescopes in Asia at all while appropriate sites exist and their appearance in Asia or Eurasia is long overdue. Kinetic inductance detectors (KID) are ideal for large-format array implementation, which will be necessary for future telescope development. Concept of multicolor subTHz KID-array MUSICAM demo camera and its instrumental testing is given. It allows us to perform some necessary steps toward the creation of the Eurasian SubMillimeter Telescopes (ESMT), which concept and scientific tasks are presented as well.

preprint2019arXiv

ALMA observations reveal no preferred outflow--filament and outflow--magnetic field orientations

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $μm$ continuum emissions. The angles between outflow axes and corresponding filaments ($γ_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters.

preprint2019arXiv

Mapping Observations of complex organic molecules around Sagittarius B2 with ARO 12m telescope

We performed high-sensitivity mapping observations of several complex organic molecules around Sagittarius B2 with ARO 12m telescope at 3-mm wavelength. Based on their spatial distribution, molecules can be classified as either "extended" that detected not only in Sgr B2(N) and Sgr B2(M), or "compact" that only detected toward or near to Sgr B2(N) and Sgr B2(M). The "extended" molecules including glycolaldehyde (CH2OHCHO), methyl formate (CH3OCHO), formic acid (t-HCOOH), ethanol (C2H5OH) and methyl amine (CH3NH2), while the "compact" molecules including dimethyl ether (CH3OCH3), ethyl cyanide (C2H5CN), and amino acetonitrile (H2NCH2CN). These "compact" molecules are likely produced under strong UV radiation, while "extended" molecules are likely formed under low-temperature, via gas-phase or grain surface reactions. The spatial distribution of "warm" CH2OHCHO at 89 GHz differ from the spatial distribution of "cold" CH2OHCHO observed at 13 GHz. We found evidence for an overabundance of CH2OHCHO compared to that expected from the gas-phase model, which indicates that grain-surface reactions are necessary to explain the origin of CH2OHCHO in Sagittarius B2. Grain-surface reactions are also needed to explain the correlation between the abundances of "cold" CH2OHCHO and C2H5OH. These results demonstrate the importance of grain-surface chemistry in the production of complex organic molecules.

preprint2019arXiv

Pilot HI Survey of Planck Galactic Cold Clumps with FAST

We present a pilot HI survey of 17 Planck Galactic Cold Clumps (PGCCs) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). HI Narrow Self-Absorption (HINSA) is an effective method to detect cold HI being mixed with molecular hydrogen H$_2$ and improves our understanding of the atomic to molecular transition in the interstellar medium. HINSA was found in 58\% PGCCs that we observed. The column density of HINSA was found to have an intermediate correlation with that of $^{13}$CO, following $\rm log( N(HINSA)) = (0.52\pm 0.26) log(N_{^{13}CO}) + (10 \pm 4.1) $. HI abundance relative to total hydrogen [HI]/[H] has an average value of $4.4\times 10^{-3}$, which is about 2.8 times of the average value of previous HINSA surveys toward molecular clouds. For clouds with total column density N$\rm_H >5 \times 10^{20}$ cm$^{-2}$, an inverse correlation between HINSA abundance and total hydrogen column density is found, confirming the depletion of cold HI gas during molecular gas formation in more massive clouds. Nonthermal line width of $^{13}$CO is about 0-0.5 km s$^{-1}$ larger than that of HINSA. One possible explanation of narrower nonthermal width of HINSA is that HINSA region is smaller than that of $^{13}$CO. Based on an analytic model of H$_2$ formation and H$_2$ dissociation by cosmic ray, we found the cloud ages to be within 10$^{6.7}$-10$^{7.0}$ yr for five sources.

preprint2019arXiv

Probing the emission states of PSR J1107-5907

The emission from PSR J1107-5907 is erratic. Sometimes the radio pulse is undetectable, at other times the pulsed emission is weak, and for short durations the emission can be very bright. In order to improve our understanding of these state changes, we have identified archival data sets from the Parkes radio telescope in which the bright emission is present, and find that the emission never switches from the bright state to the weak state, but instead always transitions to the off state. Previous work had suggested the identification of the off state as an extreme manifestation of the weak state. However, the connection between the off and bright emission reported here suggests that the emission can be interpreted as undergoing only two emission states: a bursting state consisting of both bright pulses and nulls as well as the weak-emission state.

preprint2019arXiv

Revealing The CO X-factor In Dark Molecular Gas through Sensitive ALMA Absorption Observations

Carbon-bearing molecules, particularly CO, have been widely used as tracers of molecular gas in the interstellar medium (ISM). In this work, we aim to study the properties of molecules in diffuse, cold environments, where CO tends to be under-abundant and/or sub-thermally excited. We performed one of the most sensitive (down to $\mathrm{τ_{rms}^{CO} \sim 0.002}$ and $\mathrm{τ_{rms}^{HCO^+} \sim 0.0008}$) sub-millimeter molecular absorption line observations towards 13 continuum sources with the ALMA. CO absorption was detected in diffuse ISM down to $\mathrm{A_v< 0.32\,mag}$ and \hcop was down to $\mathrm{A_v < 0.2\,mag}$, where atomic gas and dark molecular gas (DMG) starts to dominate. Multiple transitions measured in absorption toward 3C454.3 allow for a direct determination of excitation temperatures $\mathrm{T_{ex}}$ of 4.1\,K and 2.7\,K, for CO and for \hcop, respectively, which are close to the cosmic microwave background (CMB) and provide explanation for their being undercounted in emission surveys. A stronger linear correlation was found between $\mathrm{N_{HCO^+}}$ and $\mathrm{N_{H_2}}$ (Pearson correlation coefficient P $\sim$ 0.93) than that of $\mathrm{N_{CO}}$ and $\mathrm{N_{H_2}}$ (P $\sim$ 0.33), suggesting \hcop\ being a better tracer of H$_2$ than CO in diffuse gas. The derived CO-to-\h2 conversion factor (the CO X-factor) of (14 $\pm$ 3) $\times$ 10$^{20}$ cm$^{-2}$ (K \kms)$^{-1}$ is approximately 6 times larger than the average value found in the Milky Way.

preprint2019arXiv

The design of China Reconfigurable Analog-digital backEnd for FAST

The Five-hundred-meter Aperture Spherical radio Telescope(FAST) was launched on September 25,2016.From early 2017,we began to use the FAST wideband receiver,which was designed,constructed and installed on the FAST in Guizhou,China.The front end of the receiver is composed an uncooled Quad Ridge Flared Horn feed(QRFH) with the frequency range of 270 to 1620 MHz,and a cryostat operating at 10 K.Stephen et al. 2016We have coop-erated with the Institute of Automation of the Chinese Academy of Sciences to developed the China Reconfigurable ANalog-digital backEnd.The system covers the 3 GHz operating band of FAST.The hardware part of the backend includes an Analog Front-end Board,a wideband high precision Analog Digital Converter,and a FAST Digital Back-end.Analog circuit boards, field programmable gate arrays, and control computers form a set of hardware, software, and firmware platforms to achieve flexible bandwidth requirements through parameter changes. It is also suitable for the versatility of different astronomical observations, and can meet specific requirements. This paper briefly introduces the hardware and software of CRANE, as well as some observations of the system.

preprint2019arXiv

Tracing the Formation of Molecular Clouds in a Low-Metallicity Galaxy: A HI Narrow Self-Absorption Survey of the Large Magellanic Cloud

Cold atomic hydrogen clouds are the precursors of molecular clouds. Due to self-absorption, the opacity of cold atomic hydrogen may be high, and this gas may constitute an important mass component of the interstellar medium (ISM). Atomic hydrogen gas can be cooled to temperatures much lower than found in the cold neutral medium (CNM) through collisions with molecular hydrogen. In this paper, we search for HI Narrow Self-Absorption (HINSA) features in the Large Magellanic Cloud (LMC) as an indicator of such cold HI clouds, and use the results to quantify atomic masses and atomic-to-molecular gas ratio. Our search for HINSA features was conducted towards molecular clouds in the LMC using the ATCA+Parkes HI survey and the MAGMA CO survey. HINSA features are prevalent in the surveyed sightlines. This is the first detection of HINSA in an external galaxy. The HINSA-HI/$\rm{H}_{2}$ ratio in the LMC varies from 0.5\e{-3} to 3.4\e{-3} (68\% interval), with a mean value of $(1.31 \pm 0.03)$\e{-3}, after correcting for the effect of foreground HI gas. This is similar to the Milky Way value and indicates that similar fractions of cold gas exist in the LMC and the Milky Way, despite their differing metallicities, dust content and radiation fields. The low ratio also confirms that, as with the Milky Way, the formation timescale of molecular clouds is short. The ratio shows no radial gradient, unlike the case for stellar metallicity. No correlation is found between our results and those from previous HI absorption studies of the LMC.

preprint2019arXiv

Wide Bandwidth Observations of Pulsars C, D and J in 47 Tucanae

We report the first wideband observations of pulsars C, D and J in the globular cluster 47Tucanae (NGC 104) using the Ultra-Wideband Low (UWL) receiver system recently installed on the Parkes 64 m radio telescope. The wide frequency range of the UWL receiver (704-4032 MHz), along with the well-calibrated system, allowed us to obtain flux density measurements and polarization pulse profiles. The mean pulse profiles have significant linear and circular polarization, allowing for determination of the Faraday rotation measure for each pulsar. Precise measurements of the dispersion measures show a significant deviation in the value for pulsar D compared to earlier results. Searches for new pulsars in the cluster are on-going and we have determined optimal bands for such searches using the Parkes UWL receiver system.

preprint2016arXiv

Cloud Structure of Galactic OB Cluster Forming Regions from Combining Ground and Space Based Bolometric Observations

We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, G10.3-0.1, we have performed single-component, modified black-body fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The $\sim$10$"$ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3-kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed $\sim$10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function, the two-point correlation function of column density, and the probability distribution function of clump/core separations).

preprint2016arXiv

Discovery of an extremely wide-angle bipolar outflow in AFGL 5142

Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle ($\sim$180$\arcdeg$) bipolar outflow ("EWBO") in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3-2) and HCO$^{+}$ (3-2) lines. This bipolar outflow is along a north-west to south-east direction with a line-of-sight flow velocity of about 3 km~s$^{-1}$ and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high velocity jet-like outflows and the "EWBO" are $\sim$24.5 M$_{\sun}$ and $\sim1.7\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, respectively. Global collapse of the clump is revealed by the "blue profile" in the HCO$^{+}$ (1-0) line. A hierarchical network of filaments was identified in NH$_{3}$ (1,1) emission. Clear velocity gradients of the order of 10 km~s$^{-1}$~pc$^{-1}$ are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is $\sim$3.1$\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

preprint2016arXiv

Discovery of Two New Pulsars in 47 Tucanae (NGC 104)

We report the discovery of two new millisecond pulsars (PSRs J0024$-$7204aa and J0024$-$7204ab) in the globular cluster 47\,Tucanae (NGC 104). Our results bring the total number of pulsars in 47\,Tucanae to 25. These pulsars were discovered by reprocessing archival observations from the Parkes radio telescope. We reprocessed the data using a standard search procedure based on the PRESTO software package as well as using a new method in which we incoherently added the power spectra corresponding to $\sim$1100\,hr of observations. The newly discovered PSR~J0024$-$7204aa, has a pulse frequency of $\rm \sim$541\,Hz (corresponding to a $\rm \sim$1.84 ms period), which is higher than any other pulsars currently known in the cluster and ranks 12$^{\rm{th}}$ amongst all the currently known pulsars. The dispersion measure of this pulsar, 24.941(7)\,cm$^{-3}$ pc, is the highest in the cluster. The second discovered pulsar, PSR~J0024$-$7204ab, is an isolated pulsar with a pulse frequency of $\rm \sim$270\,Hz (corresponding to a period of $\rm \sim$3.70 ms).

preprint2016arXiv

Evolution of OH and CO-dark Molecular Gas Fraction Across a Molecular Cloud Boundary In Taurus

We present observations of 12CO J=1-0, 13CO J=1-0, HI, and all four ground-state transitions of the hydroxyl (OH) radical toward a sharp boundary region of the Taurus molecular cloud. Based on a PDR model that reproduces CO and [CI] emission from the same region, we modeled the three OH transitions, 1612, 1665, 1667 MHz successfully through escape probability non-LTE radiative transfer model calculations. We could not reproduce the 1720 MHz observations, due to un-modeled pumping mechanisms, of which the most likely candidate is a C-shock. The abundance of OH and CO-dark molecular gas (DMG) are well constrained. The OH abundance [OH]/[H2] decreases from 8*10-7 to 1*10-7 as Av increases from 0.4 to 2.7 mag, following an empirical law [OH]/[H2]= 1.5 * 10^{-7} + 9.0 * 10^{-7} * exp(-Av/0.81), which is higher than PDR model predictions for low extinction regions by a factor of 80. The overabundance of OH at extinctions at or below 1 mag is likely the result of a C-shock. The dark gas fraction (DGF, defined as fraction of molecular gas without detectable CO emission) decreases from 80% to 20%, following a gaussian profile DGF= 0.90 * exp(-( Av -0.79 )/0.71)^2) This trend of the DGF is consistent with our understanding that the DGF drops at low visual extinction due to photodissociation of H2 and drops at high visual extinction due to CO formation. The DGF peaks in the extinction range where H2 has already formed and achieved self-shielding but 12CO has not. Two narrow velocity components with a peak-to-peak spacing of ~ 1 km s-1 were clearly identified. Their relative intensity and variation in space and frequency suggest colliding streams or gas flows at the boundary region.

preprint2016arXiv

Gas Kinematics and Star Formation in the Filamentary IRDC G34.43+0.24

We performed a multiwavelength study toward infrared dark cloud (IRDC) G34.43+0.24. New maps of 13CO $J$=1-0 and C18}O J=1-0 were obtained from the Purple Mountain Observatory (PMO) 13.7 m radio telescope. At 8 um (Spitzer - IRAC), IRDC G34.43+0.24 appears to be a dark filament extended by 18 arcmin along the north-south direction. Based on the association with the 870 um and C18O J=1-0 emission, we suggest that IRDC G34.43+0.24 should not be 18 arcmin in length, but extend by 34 arcmin. IRDC G34.43+0.24 contains some massive protostars, UC H II regions, and infrared bubbles. The spatial extend of IRDC G34.43+0.24 is about 37 pc assuming a distance of 3.7 kpc. IRDC G34.43+0.24 has a linear mass density of 1.6*10^{3} M_{sun} pc^{-1}, which is roughly consistent with its critical mass to length ratio. The turbulent motion may help stabilizing the filament against the radial collapse. Both infrared bubbles N61 and N62 show a ringlike structure at 8 um. Particularly, N61 has a double-shell structure, which has expanded into IRDC G34.43+0.24. The outer shell is traced by 8 um and 13}CO J=1-0 emission, while the inner shell is traced by 24 um and 20 cm emission. We suggest that the outer shell (9.9*10^{5} yr) is created by the expansion of H II region G34.172+0.175, while the inner shell (4.1-6.3*10^{5} yr) may be produced by the energetic stellar wind of its central massive star. From GLIMPSE I catalog, we selected some Class I sources with an age of 10^{5} yr. These Class I sources are clustered along the filamentary molecular cloud.

preprint2016arXiv

Massive Quiescent Cores in Orion: VI. The Internal Structures and a Candidate of Transiting Core in NGC 2024 Filament

We present a multi-wavelength observational study of the NGC 2024 filament using infrared to sub-millimeter continuum and the NH$_3$ $(1,1)$ and $(2,2)$ inversion transitions centered on FIR-3, the most massive core therein. FIR-3 is found to have no significant infrared point sources in the Spitzer/IRAC bands. But the NH$_3$ kinetic temperature map shows a peak value at the core center with $T_{\rm k}=25$ K which is significantly higher than the surrounding level ($T_{\rm k}=15-19$ K). Such internal heating signature without an infrared source suggests an ongoing core collapse possibly at a transition stage from first hydrostatic core (FHSC) to protostar. The eight dense cores in the filament have dust temperatures between 17.5 and 22 K. They are much cooler than the hot ridge ($T_{\rm d}=55$ K) around the central heating star IRS-2b. Comparison with a dust heating model suggests that the filament should have a distance of $3-5$ pc from IRS-2b. This value is much larger than the spatial extent of the hot ridge, suggesting that the filament is spatially separated from the hot region along the line of sight.

preprint2016arXiv

Physical properties of CO-dark molecular gas traced by C$^+$

Neither HI nor CO emission can reveal a significant quantity of so-called dark gas in the interstellar medium (ISM). It is considered that CO-dark molecular gas (DMG), the molecular gas with no or weak CO emission, dominates dark gas. We identified 36 DMG clouds with C$^+$ emission (data from Galactic Observations of Terahertz C+ (GOT C+) project) and HINSA features. Based on uncertainty analysis, optical depth of HI $τ\rm_{HI}$ of 1 is a reasonable value for most clouds. With the assumption of $τ\rm_{HI}=1$, these clouds were characterized by excitation temperatures in a range of 20 K to 92 K with a median value of 55 K and volume densities in the range of $6.2\times10^1$ cm$^{-3}$ to $1.2\times 10^3$ cm$^{-3}$ with a median value of $2.3\times 10^2$ cm$^{-3}$. The fraction of DMG column density in the cloud ($f\rm_{DMG}$) decreases with increasing excitation temperature following an empirical relation $f\rm_{DMG}=-2.1\times 10^{-3}T_(ex,τ_{HI}=1)$+1.0. The relation between $f\rm_{DMG}$ and total hydrogen column density $N_H$ is given by $f\rm_{DMG}$=$1.0-3.7\times 10^{20}/N_H$. The values of $f\rm_{DMG}$ in the clouds of low extinction group ($A\rm_V \le 2.7$ mag) are consistent with the results of the time-dependent, chemical evolutionary model at the age of ~ 10 Myr. Our empirical relation cannot be explained by the chemical evolutionary model for clouds in the high extinction group ($A\rm_V > 2.7$ mag). Compared to clouds in the low extinction group ($A\rm_V \le 2.7$ mag), clouds in the high extinction group ($A\rm_V > 2.7$ mag) have comparable volume densities but excitation temperatures that are 1.5 times lower. Moreover, CO abundances in clouds of the high extinction group ($A\rm_V > 2.7$ mag) are $6.6\times 10^2$ times smaller than the canonical value in the Milky Way. #[Full version of abstract is shown in the text.]#

preprint2015arXiv

A New Method for Constraining Molecular Cloud Thickness: A study of Taurus, Perseus and Ophiuchus

The core velocity dispersion (CVD) is a potentially useful tool for studying the turbulent velocity field of molecular clouds. CVD is based on centroid velocities of dense gas clumps, thus is less prone to density fluctuation and reflects more directly the cloud velocity field. Prior work demonstrated that the Taurus molecular cloud CVD resembles the well-known Larson's linewidth-size relation of molecular clouds. In this work, we studied the dependence of the CVD on the line-of-sight thickness of molecular clouds, a quantity which cannot be measured by direct means. We produced a simple statistical model of cores within clouds and analyzed the CVD of a variety of hydrodynamical simulations. We show that the relation between the CVD and the 2D projected separation of cores ($L_{2D}$) is sensitive to the cloud thickness. When the cloud is thin, the index of CVD-$L_{2D}$ relation ($γ$ in the relation CVD$\sim L_{2D}^γ$) reflects the underlying energy spectrum ($E(k)\sim k^{-β}$) in that $γ\sim(β-1)/2$. The CVD-$L_{2D}$ relation becomes flatter ($γ\to 0$) for thicker clouds. We used this result to constrain the thicknesses of Taurus, Perseus, and Ophiuchus. We conclude that Taurus has a ratio of cloud depth to cloud length smaller than about 1/10-1/8, i.e. it is a sheet. A simple geometric model fit to the linewidth-size relation indicates that the Taurus cloud has a $\sim 0.7$ pc line-of-sight dimension. In contrast, Perseus and Ophiuchus are thicker and have ratios of cloud depth to cloud length larger than about 1/10-1/8.

preprint2015arXiv

Distributed Bayesian Quickest Change Detection in Sensor Networks via Two-layer Large Deviation Analysis

We propose a distributed Bayesian quickest change detection algorithm for sensor networks, based on a random gossip inter-sensor communication structure. Without a control or fusion center, each sensor executes its local change detection procedure in a parallel and distributed fashion, interacting with its neighbor sensors via random inter-sensor communications to propagate information. By modeling the information propagation dynamics in the network as a Markov process, two-layer large deviation analysis is presented to analyze the performance of the proposed algorithm. The first-layer analysis shows that the relation between the probability of false alarm and the conditional averaged detection delay satisfies the large deviation principle, implying that the probability of false alarm according to a rare event decays to zero at an exponentially fast rate when the conditional averaged detection decay increases, where the Kullback-Leibler information number is established as a crucial factor. The second-layer analysis shows that the probability of the rare event that not all observations are available at a sensor decays to zero at an exponentially fast rate when the averaged number of communications increases, where the large deviation upper and lower bounds for this rate are also derived, based on which we show that the performance of the distributed algorithm converges exponentially fast to that of the centralized one, by proving that the defined distributed Kullback-Leibler information number converges to the centralized Kullback-Leibler information number.

preprint2015arXiv

Distributed Kalman Filtering over Massive Data Sets: Analysis Through Large Deviations of Random Riccati Equations

This paper studies the convergence of the estimation error process and the characterization of the corresponding invariant measure in distributed Kalman filtering for potentially unstable and large linear dynamic systems. A gossip network protocol termed Modified Gossip Interactive Kalman Filtering (M-GIKF) is proposed, where sensors exchange their filtered states (estimates and error covariances) and propagate their observations via inter-sensor communications of rate $\overlineγ$; $\overlineγ$ is defined as the averaged number of inter-sensor message passages per signal evolution epoch. The filtered states are interpreted as stochastic particles swapped through local interaction. The paper shows that the conditional estimation error covariance sequence at each sensor under M-GIKF evolves as a random Riccati equation (RRE) with Markov modulated switching. By formulating the RRE as a random dynamical system, it is shown that the network achieves weak consensus, i.e., the conditional estimation error covariance at a randomly selected sensor converges weakly (in distribution) to a unique invariant measure. Further, it is proved that as $\overlineγ \rightarrow \infty$ this invariant measure satisfies the Large Deviation (LD) upper and lower bounds, implying that this measure converges exponentially fast (in probability) to the Dirac measure $δ_{P^*}$, where $P^*$ is the stable error covariance of the centralized (Kalman) filtering setup. The LD results answer a fundamental question on how to quantify the rate at which the distributed scheme approaches the centralized performance as the inter-sensor communication rate increases.

preprint2015arXiv

Extremely Energetic Outflow and Decelerated Expansion in W49N

W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high-mass star formation studies. Due to its large distance (11.1$_{-0.7}^{+0.9}$ kpc), the kinematics inside and between the dense molecular clumps in W49N are far from well understood. The SMA observations resolved the continuum emission into two clumps. The molecular line observation of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) suggests that the two clumps have a velocity difference of $\sim$7 km~s$^{-1}$. The eastern clump is very close to two radio sources "G1" and "G2", and the western clump coincides with a radio source "B". The HCN (3-2) line reveals an extremely energetic outflow, which is among the most energetic molecular outflows in the Milky Way. This is the first report of high-velocity molecular outflow detection in W49N. The outflow jet might be in precession, which could account for the distribution, velocity and rotation of water maser spots. Three absorption systems are identified in HCO$^{+}$ (3-2) spectra. The absorption features are blueshifted with respect to the emission of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) lines, indicating that a cold layer is expanding in front of the warm gas. Further analysis indicates that the expansion is decelerated from the geometric expansion centers.

preprint2015arXiv

Galactic and Magellanic Evolution with the SKA

As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more detail than anywhere else in the Universe how galaxies acquire fresh gas to fuel their continuing star formation, how they exchange gas with the surrounding intergalactic medium, and turn warm, diffuse gas into molecular clouds and ultimately stars. The $λ$21-cm line of atomic hydrogen (HI) is an excellent tracer of these physical processes. With the SKA we will finally have the combination of surface brightness sensitivity, point source sensitivity and angular resolution to transform our understanding of the evolution of gas in the Milky Way, all the way from the halo down to the formation of individual molecular clouds.

preprint2015arXiv

Learning-Based Distributed Detection-Estimation in Sensor Networks with Unknown Sensor Defects

We consider the problem of distributed estimation of an unknown deterministic scalar parameter (the target signal) in a wireless sensor network (WSN), where each sensor receives a single snapshot of the field. We assume that the observation at each node randomly falls into one of two modes: a valid or an invalid observation mode. Specifically, mode one corresponds to the desired signal plus noise observation mode (\emph{valid}), and mode two corresponds to the pure noise mode (\emph{invalid}) due to node defect or damage. With no prior information on such local sensing modes, we introduce a learning-based distributed procedure, called the mixed detection-estimation (MDE) algorithm, based on iterative closed-loop interactions between mode learning (detection) and target estimation. The online learning step re-assesses the validity of the local observations at each iteration, thus refining the ongoing estimation update process. The convergence of the MDE algorithm is established analytically. Asymptotic analysis shows that, in the high signal-to-noise ratio (SNR) regime, the MDE estimation error converges to that of an ideal (centralized) estimator with perfect information about the node sensing modes. This is in contrast to the estimation performance of a naive average consensus based distributed estimator (without mode learning), whose estimation error blows up with an increasing SNR.

preprint2015arXiv

Outflows and Bubbles in Taurus: Star-formation Feedback Sufficient to Maintain Turbulence

We have identified outflows and bubbles in the Taurus molecular cloud based on the $\sim 100$ deg$^2$ Five College Radio Astronomy Observatory $^{12}$CO(1-0) and $^{13}$CO(1-0) maps and the Spitzer young stellar object catalogs. In the main 44 deg$^2$ area of Taurus we found 55 outflows, of which 31 were previously unknown. We also found 37 bubbles in the entire 100 deg$^2$ area of Taurus, all of which had not been found before. The total kinetic energy of the identified outflows is estimated to be $\bf \sim 3.9 \times 10^{45}$ erg, which is \textbf{1\%} of the cloud turbulent energy. The total kinetic energy of the detected bubbles is estimated to be $\sim 9.2 \times 10^{46}$ erg, which is 29\% of the turbulent energy of Taurus. The energy injection rate from outflows is $\bf \sim 1.3 \times 10^{33}~\rm erg\ s^{-1}$, \textbf{0.4 - 2 times} the dissipation rate of the cloud turbulence. The energy injection rate from bubbles is $\sim 6.4 \times 10^{33}$ erg s$^{-1}$, \textbf{2 - 10 times} the turbulent dissipation rate of the cloud. The gravitational binding energy of the cloud is $\bf \sim 1.5 \times 10^{48}$ {\bf erg}, \textbf{385} and 16 times the energy of outflows and bubbles, respectively. We conclude that neither outflows nor bubbles can \textbf{provide enough energy to balance the overall gravitational binding energy and the turbulent energy of Taurus. However,} in the current epoch, stellar feedback is sufficient to maintain the observed turbulence in Taurus.

preprint2015arXiv

Planck cold clumps in the $λ$ Orionis complex: I. Discovery of an extremely young Class 0 protostellar object and a proto-brown dwarf candidate in a bright rimmed clump PGCC G192.32-11.88

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $λ$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

preprint2015arXiv

Quantifying Dark Gas

A growing body of evidence has been supporting the existence of so-called "dark molecular gas" (DMG), which is invisible in the most common tracer of molecular gas, i.e., CO rotational emission. DMG is believed to be the main gas component of the intermediate extinction region between A$\rm_v$$\sim$0.05-2, roughly corresponding to the self-shielding threshold of H$_2$ and $^{13}$CO. To quantify DMG relative to HI and CO, we are pursuing three observational techniques, namely, HI self-absorption, OH absorption, and TeraHz C$^+$ emission. In this paper, we focus on preliminary results from a CO and OH absorption survey of DMG candidates. Our analysis show that the OH excitation temperature is close to that of the Galactic continuum background and that OH is a good DMG tracer co-existing with molecular hydrogen in regions without CO. Through systematic "absorption mapping" by Square Kilometer Array (SKA) and ALMA, we will have unprecedented, comprehensive knowledge of the ISM components including DMG in terms of their temperature and density, which will impact our understanding of galaxy evolution and star formation profoundly.

preprint2015arXiv

Radio Afterglows and Host Galaxies of Gamma-Ray Bursts

Considering the contribution of the emission from the host galaxies of gamma-ray bursts (GRBs) to the radio afterglows, we investigate the effect of host galaxies on observations statistically. For the three types of events, e.g. low-luminosity, standard and high-luminosity GRBs, it is found that a tight correlation exists between the ratio of the radio flux (RRF) of host galaxy to the total radio peak emission and the observational frequency. Especially, toward lower frequencies, the contribution from the host increases significantly. The correlation can be used to get a useful estimate for the radio brightness of those host galaxies which only have very limited radio afterglow data. Using this prediction, we re-considered the theoretical radio afterglow light curves for four kinds of events, i.e. high-luminosity, low-luminosity, standard and failed GRBs, taking into account the contribution from the host galaxies and aiming at exploring the detectability of these events by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Lying at a typical redshift of $z=1$, most of the events can be detected by FAST easily. For the less fierce low-luminosity GRBs, their radio afterglows are not strong enough to exceed the sensitivity limit of FAST at such distances. However, since a large number of low luminosity bursts actually happen very near to us, it is expected that FAST will still be able to detect many of them.

preprint2015arXiv

Sulfur-Bearing Molecules In Massive Star-Forming Regions: Observations Of OCS, CS, H2S and SO

We studied sulfur chemistry of massive star-forming regions through single-dish submillimeter spectroscopy. OCS, O13CS, 13CS, H2S and SO transitions were observed toward a sample of massive star-forming regions with embedded UCH II or CH II regions. These sources could be divided into H II-hot core and H II-only sources based on their CH3CN emission. Our results show that the OCS line of thirteen sources is optically thick, with optical depth ranging from 5 to 16. Column densities of these molecules were computed under LTE conditions. CS column densities were also derived using its optically thin isotopologue 13CS. H2S is likely to be the most abundant gas-phase sulfuretted molecules in hot assive cores. Both the column density and abundance of sulfur-bearing molecules decrease significantly from H II-hot core to H II-only sources. Ages derived from hot core models appear to be consistent with star-formation theories, suggesting that abundance ratios of [CS]/[SO], [SO]/[OCS] and [OCS]/[CS] could be used as chemical clocks in massive star-forming regions.

preprint2015arXiv

Superconductivity in efficient thermoelectric Cu$_3$Sb$_{0.98}$Al$_{0.02}$Se$_4$

Both superconductivity and thermoelectricity offer promising prospects for daily energy efficiency applications. The advancements of thermoelectric materials have led to the huge improvement of the thermoelectric figure of merit in the past decade. By applying pressure on a highly efficient thermoelectric material Cu$_{3}$Sb$_{0.98}$Al$_{0.02}$Se$_4$, we achieve dome-shape superconductivity developing at around 8.5 GPa but having a maximum critical temperature of 3.2 K at pressure of 12.7 GPa. The novel superconductor is realized through the first-order structural transformation from its initial phase to an orthorhombic one. The superconducting phase is determined in the ultimate formation of the Cu-Al-Sb-Se alloy.

preprint2014arXiv

350 $μ$m map of the Ophiuchus molecular cloud: core mass function

Stars are born in dense cores of molecular clouds. The core mass function (CMF), which is the mass distribution of dense cores, is important for understanding the stellar initial mass function (IMF). We obtained 350 $μ$m dust continuum data using the SHARC-II camera at the Caltech Submillimeter Observatory (CSO) telescope. A 350 $μ$m map covering 0.25 ${deg}^{2}$ of the Ophiuchus molecular cloud was created by mosaicing 56 separate scans. The CSO telescope had an angular resolution of 9 $^{\prime\prime}$, corresponding to $1.2\times {10}^{3}\ $AU at the distance of the Ophiuchus molecular cloud (131 pc). The data was reduced using the Comprehensive Reduction Utility for SHARC-II (CRUSH). The flux density map was analyzed using the GaussClumps algorithm, within which 75 cores has been identified. We used the Spitzer c2d catalogs to separate the cores into 63 starless cores and 12 protostellar cores. By locating Jeans instabilities, 55 prestellar cores (a subcategory of starless cores) were also identified. The excitation temperatures, which were derived from FCRAO ${}^{12}$CO data, help to improve the accuracy of the masses of the cores. We adopted a Monte Carlo approach to analyze the CMF with two types of functional forms; power law and log-normal. The whole and prestellar CMF are both well fitted by a log-normal distribution, with $μ=-1.18\pm0.10,\ σ=0.58\pm0.05$ and $μ=1.40\pm0.10,\ σ=0.50\pm0.05$ respectively. This finding suggests that turbulence influences the evolution of the Ophiuchus molecular cloud.

preprint2014arXiv

A CO observation of the galactic methanol masers

Context: We investigated the molecular gas associated with 6.7 GHz methanol masers throughout the Galaxy using a J=1-0 transition of the CO isotopologues. Methods:Using the 13.7-meter telescope at the Purple Mountain Observatory (PMO), we have obtained ^{12}CO and ^{13}CO (1-0) lines for 160 methanol masers sources from the first to the third Galactic quadrants. We made efforts to resolve the distance ambiguity by careful comparison with the radio continuum and HI 21 cm observations. Results: First, the maser sources show increased ^{13}CO line widths toward the Galactic center, suggesting that the molecular gas are more turbulent toward the Galactic center. This trend can be noticeably traced by the ^{13}CO line width. Second, the ^{12}CO excitation temperature shows a noticeable correlation with the H_2 column density. A possible explanation consistent with the collapse model is that the higher surface-density gas is more efficient to the stellar heating and/or has a higher formation rate of high-mass stars. Third, comparing the IRDCs, the maser sources on average have significantly lower H_2 column densities, moderately higher temperatures, and similar line widths. Fourth, in the mapped regions, 51 ^{13}CO cores have been revealed. Only 17 coincide with the radio continuum emission (F_{cm}>6 mJy), while a larger fraction (30 cores) coincide with the infrared emissions. The IR-bright and radio-bright sources exhibit significantly higher $^{12}$CO excitation temperatures than the IR-faint and radio-faint sources, respectively. Conclusions: The 6.7 GHz masers show a moderately low ionization rate but have a common-existing stellar heating that generates the IR emissions. The relevant properties can be characterized by the ^{12}CO and ^{13}CO (1-0) emissions in several aspects as described above.

preprint2014arXiv

Characterization of Molecular Outflows in The Substellar Domain

We report here our latest search for molecular outflows from young brown dwarfs and very low-mass stars in nearby star-forming regions. We have observed three sources in Taurus with the Submillimeter Array and the Combined Array for Research in Millimeter-wave Astronomy at 230 GHz frequency to search for CO J=2-1 outflows. We obtain a tentative detection of a redshifted and extended gas lobe at about 10 arcsec from the source GM Tau, a young brown dwarf in Taurus with an estimated mass of 73 M_J, which is right below the hydrogen-burning limit. No blueshifted emission around the brown dwarf position is detected. The redshifted gas lobe that is elongated in the northeast direction suggests a possible bipolar outflow from the source with a position angle of about 36 degrees. Assuming that the redshifted emission is outflow emission from GM Tau, we then estimate a molecular outflow mass in the range from 1.9x10^-6 M_Sun to 2.9x10^-5 M_Sun and an outflow mass-loss rate from 2.7x10^-9 M_Sun yr^-1 to 4.1x10^-8 M_Sun yr^-1. These values are comparable to those we have observed in the young brown dwarf ISO-Oph 102 of 60 M_J in rho Ophiuchi and the very low-mass star MHO 5 of 90 M_J in Taurus. Our results suggest that the outflow process in very low-mass objects is episodic with duration of a few thousand years and the outflow rate of active episodes does not significantly change for different stages of the formation process of very low-mass objects. This may provide us with important implications that clarify the formation process of brown dwarfs.

preprint2014arXiv

Detecting Radio Afterglows of Gamma-Ray Bursts with FAST

Using the generic hydrodynamic model of Gamma-Ray Burst (GRB) afterglows, we calculate the radio afterglow light curves of low luminosity, high luminosity, failed and standard GRBs in different observational bands of FAST's energy window. The GRBs are assumed to be located at different distances from us. Our results rank the detectability of GRBs in descending order as high luminosity, standard, failed and low luminosity GRBs. We predict that almost all types of radio afterglows except that of low luminosity GRBs could be observed by the large radio telescope as long as the domains of time and frequency are appropriate. It is important to note that FAST can detect relatively weak radio afterglows at a higher frequency of 2.5 GHz for very high redshift up to z=15 or even more. Radio afterglows of low luminosity GRBs can be detected only after the completion of the 2nd phase of FAST. FAST is expected to significantly expand the sample of GRB radio afterglows in the near future.

preprint2014arXiv

Evidence for Large Grains in the Star-forming Filament OMC-2/3

We present a new 3.3 mm continuum map of the OMC-2/3 region. When paired with previously published maps of 1.2mm continuum and NH3-derived temperature, we derive the emissivity spectral index of dust emission in this region, tracking its changes across the filament and cores. We find that the median value of the emissivity spectral index is 0.9, much shallower than previous estimates in other nearby molecular clouds. We find no significant difference between the emissivity spectral index of dust in the OMC-2/3 filament and the starless or protostellar cores. Furthermore, the temperature and emissivity spectral index, beta, are anti-correlated at the 4 sigma level. The low values of the emissivity spectral index found in OMC-2/3 can be explained by the presence of millimeter-sized dust grains in the dense regions of the filaments to which these maps are most sensitive. Alternatively, a shallow dust emissivity spectral index may indicate non-powerlaw spectral energy distributions, significant free-free emission, or anomalous microwave emission. We discuss the possible implications of millimeter-sized dust grains compared to the alternatives.

preprint2014arXiv

Follow-up observations toward Planck cold clumps with ground-based radio telescopes

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

preprint2014arXiv

IC348-SMM2E: a Class 0 proto-brown dwarf candidate forming as a scaled-down version of low-mass stars

We report on Submillimeter Array observations of the 870 micron continuum and CO(3-2), 13CO(2-1) and C18O(2-1) line emission of a faint object, SMM2E, near the driving source of the HH797 outflow in the IC348 cluster. The continuum emission shows an unresolved source for which we estimate a mass of gas and dust of 30 Mjup, and the CO(3-2) line reveals a compact bipolar outflow centred on SMM2E, and barely seen also in 13CO(2-1). In addition, C18O(2-1) emission reveals hints of a possible rotating envelope/disk perpendicular to the outflow, for which we infer a dynamical mass of ~16 Mjup. In order to further constrain the accreted mass of the object, we gathered data from Spitzer, Herschel, and new and archive submillimetre observations, and built the Spectral Energy Distribution (SED). The SED can be fitted with one single modified black-body from 70 micron down to 2.1 cm, using a dust temperature of ~24 K, a dust emissivity index of 0.8, and an envelope mass of ~35 Mjup. The bolometric luminosity is 0.10 Lsun, and the bolometric temperature is 35 K. Thus, SMM2E is comparable to the known Class 0 objects in the stellar domain. An estimate of the final mass indicates that SMM2E will most likely remain substellar, and the SMM2E outflow force matches the trend with luminosity known for young stellar objects. Thus, SMM2E constitutes an excellent example of a Class 0 proto-brown dwarf candidate which forms as a scaled-down version of low-mass stars. Finally, SMM2E seems to be part of a wide (~2400 AU) multiple system of Class 0 sources.

preprint2014arXiv

Massive Quiescent Cores in Orion: V. The Physical Structures on Sub-Jeans Scale and the Chemical Properties in Two Extremely Dense Cores

We present a high-resolution $1.5"$ observational study towards two massive dust-and-gas cores, ORI8nw\_2 and ORI2\_6 in Orion Molecular Cloud using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). In both regions the 3.2 mm continuum emission reveals a dense and compact dust core of 1 to 3 solar masses. The cores are estimated to have number densities more than $10^9$ cm$^{-3}$, which are among the highest volume density ever published for the star-forming cores. In both regions the ${\rm N_2H^+}$ shows multiple gas clumps which are spatially displaced from the dense gas traced by the 3.2 mm continuum as well as the ${\rm HCO^+}$. The ${\rm N_2H^+}$ clumps in ORI8nw\_2 are gravitationally unbound. The clumps have masses comparable to the thermal Jeans mass, and are separated with intervals comparable with the Jeans length, thus likely result from thermal Jeans fragmentation. The fragmentation may have prohibited the high-mass star formation in the two cores which originally have larger total masses. In both regions the abundance ratio of ${\rm [N_2H^+]/[HCO^+]}$ appears to be higher than those in infrared dark clouds as a result of CO depletion. The extreme high central density in these Orion cores implies a very small Jeans scale ($<500$ AU). To further reveal the evolution of such dense cores requires higher resolution and sensitivity.

preprint2014arXiv

Mobile Adhoc Offloading

This problem is a series of biddings and auctions. Each round of bidding and auction are different from previous ones because of the change of network topology, variance of budget set by the sender, and possible evolution of strategies of other nodes. The huge strategy space of relay nodes makes the formulation to a game very difficult. We present a brief qualitative analysis in this paper, and propose a bidding strategy based on learning algorithms.

preprint2013arXiv

Quickest Change Point Detection and Identification Across a Generic Sensor Array

In this paper, we consider the problem of quickest change point detection and identification over a linear array of $N$ sensors, where the change pattern could first reach any of these sensors, and then propagate to the other sensors. Our goal is not only to detect the presence of such a change as quickly as possible, but also to identify which sensor that the change pattern first reaches. We jointly design two decision rules: a stopping rule, which determines when we should stop sampling and claim a change occurred, and a terminal decision rule, which decides which sensor that the change pattern reaches first, with the objective to strike a balance among the detection delay, the false alarm probability, and the false identification probability. We show that this problem can be converted to a Markov optimal stopping time problem, from which some technical tools could be borrowed. Furthermore, to avoid the high implementation complexity issue of the optimal rules, we develop a scheme with a much simpler structure and certain performance guarantee.

preprint2012arXiv

13CO Cores in Taurus Molecular Cloud

Young stars form in molecular cores, which are dense condensations within molecular clouds. We have searched for molecular cores traced by $^{13}$CO $J=1\to 0$ emission in the Taurus molecular cloud and studied their properties. Our data set has a spatial dynamic range (the ratio of linear map size to the pixel size) of about 1000 and spectrally resolved velocity information, which together allow a systematic examination of the distribution and dynamic state of $^{13}$CO cores in a large contiguous region. We use empirical fit to the CO and CO$_2$ ice to correct for depletion of gas-phase CO. The $^{13}$CO core mass function ($^{13}$CO CMF) can be fitted better with a log-normal function than with a power law function. We also extract cores and calculate the $^{13}$CO CMF based on the integrated intensity of $^{13}$CO and the CMF from 2MASS. We demonstrate that there exists core blending, i.e.\ combined structures that are incoherent in velocity but continuous in column density. The core velocity dispersion (CVD), which is the variance of the core velocity difference $δv$, exhibits a power-law behavior as a function of the apparent separation $L$:\ CVD (km/s) $\propto L ({\rm pc})^{0.7}$. This is similar to Larson's law for the velocity dispersion of the gas. The peak velocities of $^{13}$CO cores do not deviate from the centroid velocities of the ambient $^{12}$CO gas by more than half of the line width. The low velocity dispersion among cores, the close similarity between CVD and Larson's law, and the small separation between core centroid velocities and the ambient gas all suggest that molecular cores condense out of the diffuse gas without additional energy from star formation or significant impact from converging flows.

preprint2012arXiv

Gas emissions in Planck cold dust clumps---A Survey of the J=1-0 Transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O

A survey toward 674 Planck cold clumps of the Early Cold Core Catalogue (ECC) in the J=1-0 transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O has been carried out using the PMO 13.7 m telescope. 673 clumps were detected with the $^{12}$CO and $^{13}$CO, and 68% of the samples have C$^{18}$O emission. Additional velocity components were also identified.A close consistency of the three line peak velocities was revealed for the first time. Kinematic distances are given out for all the velocity components and half of the clumps are located within 0.5 and 1.5 kpc. Excitation temperatures range from 4 to 27 K, slightly larger than those of $T_d$. Line width analysis shows that the majority of ECC clumps are low mass clumps. Column densities N$_{H_{2}}$ span from 10$^{20}$ to 4.5$\times10^{22}$ cm$^{-2}$ with an average value of (4.4$\pm$3.6)$\times10^{21}$ cm$^{-2}$. N$_{H_{2}}$ cumulative fraction distribution deviates from the lognormal distribution, which is attributed to optical depth. The average abundance ratio of the $^{13}$CO to C$^{18}$O in these clumps is 7.0$\pm$3.8, higher than the terrestrial value. Dust and gas are well coupled in 95% of the clumps. Blue profile, red profile and line asymmetry in total was found in less than 10% of the clumps, generally indicating star formation is not developed yet. Ten clumps were mapped. Twelve velocity components and 22 cores were obtained. Their morphologies include extended diffuse, dense isolated, cometary and filament, of which the last is the majority. 20 cores are starless.Only 7 cores seem to be in gravitationally bound state. Planck cold clumps are the most quiescent among the samples of weak-red IRAS, infrared dark clouds, UC H{\sc ii} region candidates, EGOs and methanol maser sources, suggesting that Planck cold clumps have expanded the horizon of cold Astronomy.

preprint2012arXiv

Herschel Search for O2 Toward the Orion Bar

We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 - 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams.

preprint2012arXiv

Is the Taurus B213 Region a True Filament?: Observations of Multiple Cyanoacetylene Transitions

We have obtained spectra of the J=2-1 and J=10-9 transitions of cyanoacetylene (\hc3n) toward a collection of positions in the most prominent filament, B213, in the Taurus molecular cloud. The analysis of the excitation conditions of these transitions reveals an average gas H$_2$ volume density of $(1.8\pm 0.7) \times10^{4} $ \cc. Based on column density derived from 2MASS and this volume density, the line of sight dimension of the high density portion of B213 is found to be $\simeq$ 0.12 pc, which is comparable to the smaller projected dimension and much smaller than the elongated dimension of B213 ($\sim$2.4 pc). B213 is thus likely a true cylinder--like filament rather than a sheet seen edge-on. The line width and velocity gradient seen in \hc3n are also consistent with Taurus B213 being a self-gravitating filament in the early stage of either fragmentation and/or collape.

preprint2012arXiv

The Five-hundred-meter Aperture Spherical Radio Telescope Project and its Early Science Opportunities

The National Astronomical Observatories, Chinese Academy of Science (NAOC), has started building the largest antenna in the world. Known as FAST, the Five-hundred-meter Aperture Spherical radio Telescope is a Chinese mega-science project funded by the National Development and Reform Commission (NDRC). FAST also represents part of Chinese contribution to the international efforts to build the square kilometer array (SKA). Upon its finishing around September of 2016, FAST will be the most sensitive single-dish radio telescope in the low frequency radio bands between 70 MHz and 3 GHz. The design specifications of FAST, its expected capabilities, and its main scientific aspirations were described in an overview paper by Nan et al. (2011). In this paper, we briefly review the design and the key science goals of FAST, speculate the likely limitations at the initial stages of FAST operation, and discuss the opportunities for astronomical discoveries in the so-called early science phase.

preprint2011arXiv

Herschel Measurements of Molecular Oxygen in Orion

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

preprint2011arXiv

The Five-Hundred-Meter Aperture Spherical Radio Telescope (FAST) Project

Five-hundred-meter Aperture Spherical radio Telescope (FAST) is a Chinese mega-science project to build the largest single dish radio telescope in the world. Its innovative engineering concept and design pave a new road to realize a huge single dish in the most effective way. FAST also represents Chinese contribution in the international efforts to build the square kilometer array (SKA). Being the most sensitive single dish radio telescope, FAST will enable astronomers to jump-start many science goals, for example, surveying the neutral hydrogen in the Milky Way and other galaxies, detecting faint pulsars, looking for the first shining stars, hearing the possible signals from other civilizations, etc. The idea of sitting a large spherical dish in a karst depression is rooted in Arecibo telescope. FAST is an Arecibo-type antenna with three outstanding aspects: the karst depression used as the site, which is large to host the 500-meter telescope and deep to allow a zenith angle of 40 degrees; the active main reflector correcting for spherical aberration on the ground to achieve a full polarization and a wide band without involving complex feed systems; and the light-weight feed cabin driven by cables and servomechanism plus a parallel robot as a secondary adjustable system to move with high precision. The feasibility studies for FAST have been carried out for 14 years, supported by Chinese and world astronomical communities. The project time is 5.5 years from the commencement of work in March of 2011 and the first light is expected to be in 2016. This review intends to introduce FAST project with emphasis on the recent progress since 2006. In this paper, the subsystems of FAST are described in modest details followed by discussions of the fundamental science goals and examples of early science projects.

preprint2011arXiv

The Magnetic Field in Taurus Probed by Infrared Polarization

We present maps of the plane-of-sky magnetic field within two regions of the Taurus molecular cloud: one in the dense core L1495/B213 filament, the other in a diffuse region to the west. The field is measured from the polarization of background starlight seen through the cloud. In total, we measured 287 high-quality near-infrared polarization vectors in these regions. In L1495/B213, the percent polarization increases with column density up to Av ~ 9 mag, the limits of our data. The Radiative Torques model for grain alignment can explain this behavior, but models that invoke turbulence are inconsistent with the data. We also combine our data with published optical and near-infrared polarization measurements in Taurus. Using this large sample, we estimate the strength of the plane-of-sky component of the magnetic field in nine subregions. This estimation is done with two different techniques that use the observed dispersion in polarization angles. Our values range from 5-82 microgauss and tend to be higher in denser regions. In all subregions, the critical index of the mass-to-magnetic flux ratio is sub-unity, implying that Taurus is magnetically supported on large scales (~2 pc). Within the region observed, the B213 filament makes a sharp turn to the north and the direction of the magnetic field also takes a sharp turn, switching from being perpendicular to the filament to becoming parallel. This behavior can be understood if we are observing the rim of a bubble. We argue that it has resulted from a supernova remnant associated with a recently discovered nearby gamma-ray pulsar.

preprint2010arXiv

A Selection Region Based Routing Protocol for Random Mobile ad hoc Networks

We propose a selection region based multi-hop routing protocol for random mobile ad hoc networks, where the selection region is defined by two parameters: a reference distance and a selection angle. At each hop, a relay is chosen as the nearest node to the transmitter that is located within the selection region. By assuming that the relay nodes are randomly placed, we derive an upper bound for the optimum reference distance to maximize the expected density of progress and investigate the relationship between the optimum selection angle and the optimum reference distance. We also note that the optimized expected density of progress scales as $Θ(\sqrtλ)$, which matches the prior results in the literature. Compared with the spatial-reuse multi-hop protocol in \cite{Baccelli:Aloha} recently proposed by Baccelli \emph{et al.}, in our new protocol the amount of nodes involved and the calculation complexity for each relay selection are reduced significantly, which is attractive for energy-limited wireless ad hoc networks (e.g., wireless sensor networks).

preprint2010arXiv

A Selection Region Based Routing Protocol for Random Mobile ad hoc Networks with Directional Antennas

In this paper, we propose a selection region based multihop routing protocol with directional antennas for wireless mobile ad hoc networks, where the selection region is defined by two parameters: a reference distance and the beamwidth of the directional antenna. At each hop, we choose the nearest node to the transmitter within the selection region as the next hop relay. By maximizing the expected density of progress, we present an upper bound for the optimum reference distance and derive the relationship between the optimum reference distance and the optimum transmission probability. Compared with the results with routing strategy using omnidirectional antennas in \cite{Di:Relay-Region}, we find interestingly that the optimum transmission probability is a constant independent of the beamwidth, the expected density of progress with the new routing strategy is increased significantly, and the computational complexity involved in the relay selection is also greatly reduced.

preprint2010arXiv

Distance to the SNR CTB 109/AXP 1E 2259+586 by HI absorption and self-absorption

We suggest a revised distance to the supernova remnant (SNR) G109.1-1.0 (CTB 109) and its associated anomalous X-ray pulsar (AXP) 1E 2259+586 by analyzing 21cm HI-line and 12CO-line spectra of CTB 109, HII region Sh 152, and the adjacent molecular cloud complex. CTB 109 has been established to be interacting with a large molecular cloud (recession velocity at v=-55 km s^-1). The highest radial velocities of absorption features towards CTB 109 (-56 km s^-1) and Sh 152 (-65 km s^-1) are larger than the recombination line velocity (-50 km s^-1) of Sh 152 demonstrating the velocity reversal within the Perseus arm. The molecular cloud has cold HI column density large enough to produce HI self-absorption (HISA) and HI narrow self-absorption (HINSA) if it was at the near side of the velocity reversal. Absence of both HISA and HINSA indicates that the cloud is at the far side of the velocity reversal within the Perseus Arm, so we obtain a distance for CTB 109 of 4+/-0.8 kpc. The new distance still leads to a normal explosion energy for CTB 109/AXP 1E 2259+586.

preprint2010arXiv

Distribution of Water Vapor in Molecular Clouds

We report the results of a large-area study of water vapor along the Orion Molecular Cloud ridge, the purpose of which was to determine the depth-dependent distribution of gas-phase water in dense molecular clouds. We find that the water vapor measured toward 77 spatial positions along the face-on Orion ridge, excluding positions surrounding the outflow associated with BN/KL and IRc2, display integrated intensities that correlate strongly with known cloud surface tracers such as CN, C2H, 13CO J =5-4, and HCN, and less well with the volume tracer N2H+. Moreover, at total column densities corresponding to Av < 15 mag., the ratio of H2O to C18O integrated intensities shows a clear rise approaching the cloud surface. We show that this behavior cannot be accounted for by either optical depth or excitation effects, but suggests that gas-phase water abundances fall at large Av. These results are important as they affect measures of the true water-vapor abundance in molecular clouds by highlighting the limitations of comparing measured water vapor column densities with such traditional cloud tracers as 13CO or C18O. These results also support cloud models that incorporate freeze-out of molecules as a critical component in determining the depth-dependent abundance of water vapor.

preprint2010arXiv

The relation between gas and dust in the Taurus Molecular Cloud

(abridged) We report a study of the relation between dust and gas over a 100deg^2 area in the Taurus molecular cloud. We compare the H2 column density derived from dust extinction with the CO column density derived from the 12CO and 13CO J= 1-0 lines. We derive the visual extinction from reddening determined from 2MASS data. The comparison is done at an angular size of 200", corresponding to 0.14pc at a distance of 140pc. We find that the relation between visual extinction Av and N(CO) is linear between Av~3 and 10 mag in the region associated with the B213--L1495 filament. In other regions the linear relation is flattened for Av > 4 mag. We find that the presence of temperature gradients in the molecular gas affects the determination of N(CO) by ~30--70% with the largest difference occurring at large column densities. Adding a correction for this effect and accounting for the observed relation between the column density of CO and CO2 ices and Av, we find a linear relationship between the column of carbon monoxide and dust for observed visual extinctions up to the maximum value in our data 23mag. We have used these data to study a sample of dense cores in Taurus. Fitting an analytical column density profile to these cores we derive an average volume density of about 1.4e4 cm^-3 and a CO depletion age of about 4.2e5 years. We estimate the H2 mass of Taurus to be about 1.5e4 M_sun, independently derived from the Av and N(CO) maps. We derive a CO integrated intensity to H2 conversion factor of about 2.1e20 cm^-2 (K km/s)^-1, which applies even in the region where the [CO]/[H_2] ratio is reduced by up to two orders of magnitude. The distribution of column densities in our Taurus maps resembles a log--normal function but shows tails at large and low column densities.