Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

We performed a multi-wavelength observation toward LkHa 101 embedded cluster and its adjacent 85arcmin*60arcmin region. The LkHa 101 embedded cluster is the first and only one significant cluster in California molecular cloud (CMC). These observations have revealed that the LkHa 101 embedded cluster is just located at the projected intersectional region of two filaments. One filament is the highest-density section of the CMC, the other is a new identified filament with a low-density gas emission. Toward the projected intersection, we find the bridging features connecting the two filaments in velocity, and identify a V-shape gas structure. These agree with the scenario that the two filaments are colliding with each other. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we measured that the RRL velocity of the LkHa 101 H II region is 0.5 km/s, which is related to the velocity component of the CMC filament. Moreover, there are some YSOs distributed outside the intersectional region. We suggest that the cloud-cloud collision together with the fragmentation of the main filament may play an important role in the YSOs formation of the cluster.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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.

Relativistic, magnetized jets are observed to propagate to very large distances in many Active Galactic Nuclei (AGN). We use 3D relativistic MHD (RMHD) simulations to study the propagation of Poynting flux-driven jets in AGN. These jets are assumed already being launched from the vicinity ($\sim 10^3$ gravitational radii) of supermassive black holes. Jet injections are characterized by a model described in Li et al. (2006) and we follow the propagation of these jets to ~ parsec scales. We find that these current-carrying jets are always collimated and mildly relativistic. When $\alpha$, the ratio of toroidal-to-poloidal magnetic flux injection, is large the jet is subject to non-axisymmetric current-driven instabilities (CDI) which lead to substantial dissipation and reduced jet speed. However, even with the presence of instabilities, the jet is not disrupted and will continue to propagate to large distances. We suggest that the relatively weak impact by the instability is due to the nature of the instability being convective and the fact that the jet magnetic fields are rapidly evolving on Alfv\'enic timescale. We present the detailed jet properties and show that far from the jet launching region, a substantial amount of magnetic energy has been transformed into kinetic energy and thermal energy, producing a jet magnetization number $\sigma < 1$. In addition, we have also studied the effects of a gas pressure supported "disk" surrounding the injection region and qualitatively similar global jet behaviors were observed. We stress that jet collimation, CDIs, and the subsequent energy transitions are intrinsic features of current-carrying jets.

Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by the magneto-rotational instability (MRI). Using them to predict accurately quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: Qz and Qphi, the ratio of the grid zone size to the characteristic MRI wavelength, alpha_mag, the ratio of the Maxwell stress to the magnetic pressure, and the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by Qz > 15, Qphi > 20, alpha_mag = 0.45, and a field energy ratio of 0.2. These values are associated with > 35 zones per scaleheight, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the HLLE flux solver or second-order interpolation can significantly degrade the effective resolution compared to the HLLD flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation.

Methods: We used the dual-band receiver GREAT on board the SOFIA airborne telescope to perform observations of the [C II] 158 {\mu}m fine-structure line at the postitions of two giant molecular clouds (GMC) in the center of IC 342 (GMCs C and E) and compared the spectra with corresponding ground-based data for low- and mid-J CO and [C I]. We performed model calculations assuming a clumpy photo-dissociation region (PDR) environment using the KOSMA-tau PDR model code to derive physical parameters of the local medium. Results: The [C II] 158 {\mu}m emission resembles the spectral signature of ground-based atomic and molecular lines, which indicates a common origin. The emission from GMC E can be decomposed into a cool, molecular component with weak far-ultraviolet (FUV) fields and low, mean densities of 103 cm^-3 and a strongly excited starburst/PDR region with higher densities of 104 cm^-3 and FUV intensities of 250-300 Draine fields. The emission from GMC C is consistent with gas densities of 5000 cm^-3, FUV intensities of a few Draine fields and total gas masses of 20\times10^6 M$_\odot$. Conclusions: The high spectral resolution of the GREAT receiver allowed us to decompose the [C II] emission of the GMC E into a strongly excited gas component resembling a PDR/starburst environment and a quieter, less excited gas component and to analyze the different components within a single beam individually.

Context. The HII region/PDR/molecular cloud complex S106 is excited by a single O-star. The full extent of the warm and dense gas close to the star has not been mapped in spectrally resolved high-J CO or [CII] lines, so the kinematics of the warm, partially ionized gas, are unknown. Whether the prominent dark lane bisecting the hourglass-shaped nebula is due solely to the shadow cast by a small disk around the exciting star or also to extinction in high column foreground gas was an open question until now. Aims. To disentangle the morphology and kinematics of warm neutral and ionized gas close to the star, study their relation to the bulk of the molecular gas, and to investigate the nature of the dark lane. Methods. We use the heterodyne receiver GREAT on board SOFIA to observe velocity resolved spectral lines of [C II] and CO 11-10 in comparison with so far unpublished submm continuum data at 350 micron (SHARC-II) and complementary molecular line data. Results. The high angular and spectral resolution observations show a very complex morphology and kinematics of the inner S106 region, with many different components at different excitation conditions contributing to the observed emission. The [C II] lines are found to be bright and very broad, tracing high velocity gas close to the interface of molecular cloud and HII region. CO 11-10 emission is more confined, both spatially and in velocity, to the immediate surroundings of S106 IR showing the presence of warm, high density (clumpy) gas. Our high angular resolution submm continuum observations rule out the scenario where the dark lane separating the two lobes is due solely to the shadow cast by a small disk close to the star. The lane is clearly seen also as warm, high column density gas at the boundary of the molecular cloud and HII region.

Context: We analyze the NGC 2024 HII region and molecular cloud interface using [12CII] and [13CII] observations. Aims: We attempt to gain insight into the physical structure of the interface layer between the molecular cloud and the HII region. Methods. Observations of [12CII] and [13CII] emission at 158 {\mu}m with high spatial and spectral resolution allow us to study the detailed structure of the ionization front and estimate the column densities and temperatures of the ionized carbon layer in the PDR. Results: The [12CII] emission closely follows the distribution of the 8 mum continuum. Across most of the source, the spectral lines have two velocity peaks similar to lines of rare CO isotopes. The [13CII] emission is detected near the edge-on ionization front. It has only a single velocity component, which implies that the [12CII] line shape is caused by self-absorption. An anomalous hyperfine line-intensity ratio observed in [13CII] cannot yet be explained. Conclusions: Our analysis of the two isotopes results in a total column density of N(H)~1.6\times10^23 cm^-2 in the gas emitting the [CII] line. A large fraction of this gas has to be at a temperature of several hundred K. The self-absorption is caused by a cooler (T<=100 K) foreground component containing a column density of N(H)~10^22 cm^-2.

The GREAT observations need frequency-selective calibration across the passband for the residual atmospheric opacity at flight altitude. At these altitudes the atmospheric opacity has both narrow and broad spectral features. To determine the atmospheric transmission at high spectral resolution, GREAT compares the observed atmospheric emission with atmospheric model predictions, and therefore depends on the validity of the atmospheric models. We discusse the problems identified in this comparison with respect to the observed data and the models, and describe the strategy used to calibrate the science data from GREAT/SOFIA during the first observing periods.

Discretized numerical simulations are a powerful tool for investigation of nonlinear MHD turbulence in accretion disks. However, confidence in their quantitative predictions requires a demonstration that further refinement of the spatial gridscale would not result in any significant change. This has yet to be accomplished, particularly for global disk simulations. In this paper, we combine data from previously published stratified shearing box simulations and new global disk simulations to calibrate several quantitative diagnostics by which one can estimate progress toward numerical convergence. Using these diagnostics, we find that the established criterion for an adequate numerical description of linear growth of the magneto-rotational instability (the number of cells across a wavelength of the fastest-growing vertical wavenumber mode) can be extended to a criterion for adequate description of nonlinear MHD disk turbulence, but the standard required is more stringent. We also find that azimuthal resolution, which has not often been extensively examined in previous studies, can significantly affect the evolution of the poloidal magnetic field. We further analyze the comparative resolution requirements of a small sample of initial magnetic field geometries; not surprisingly, more complicated initial field geometries require higher spatial resolution. Otherwise, they tend to evolve to qualitatively similar states if evolved for sufficient time. Applying our quantitative resolution criteria to a sample of previously published global simulations, we find that, with perhaps a single exception, they are significantly under-resolved, and therefore underestimate the magnetic turbulence and resulting stress levels throughout the accretion flow.

We present a high-resolution study of a massive dense core JCMT 18354-0649S with the Submillimeter Array. The core is mapped with continuum emission at 1.3 mm, and molecular lines including CH$_{3}$OH ($5_{23}$-$4_{13}$) and HCN (3-2). The dust core detected in the compact configuration has a mass of $47 M_{\odot}$ and a diameter of $2\arcsec$ (0.06 pc), which is further resolved into three condensations with a total mass of $42 M_{\odot}$ under higher spatial resolution. The HCN (3-2) line exhibits asymmetric profile consistent with infall signature. The infall rate is estimated to be $2.0\times10^{-3} M_{\odot}\cdot$yr$^{-1}$. The high velocity HCN (3-2) line wings present an outflow with three lobes. Their total mass is $12 M_{\odot}$ and total momentum is $121 M_{\odot}\cdot$km s$^{-1}$, respectively. Analysis shows that the N-bearing molecules especially HCN can trace both inflow and outflow.

We consider local, stratified, numerical models of isothermal accretion disks. The novel feature of our treatment is that radial extent L_x and azimuthal extent L_y satisfy H << L_x, L_y << R, where H is the scale height and R is the local radius. This enables us to probe mesoscale structure in stratified thin disks. We evolve the model at several resolutions, sizes, and initial magnetic field strengths. Consistent with earlier work, we find that the saturated, turbulent state consists of a weakly magnetized disk midplane coupled to a strongly magnetized corona, with a transition at |z| ~ 2H. The saturated \alpha ~ 0.01 - 0.02. A two-point correlation function analysis reveals that the central 4H of the disk is dominated by small scale turbulence that is statistically similar to unstratified disk models, while the coronal magnetic fields are correlated on scales ~ 10 H. Nevertheless angular momentum transport through the corona is small. A study of magnetic field loops in the corona reveals few open field lines and predominantly toroidal loops with a characteristic distance between footpoints that is ~ H. Finally we find quasi-periodic oscillations with characteristic timescale ~ 30 \Omega^{-1} in the magnetic field energy density. These oscillations are correlated with oscillations in the mean azimuthal field; we present a phenomenological, alpha-dynamo model that captures most aspects of the oscillations.

The magnetic Prandtl number Pr_M is the ratio of viscosity to resistivity. In astrophysical disks the diffusion of angular momentum (viscosity) and magnetic fields (resistivity) are controlled by turbulence. Phenomenological models of the evolution of large scale poloidal magnetic fields in disks suggest that the turbulent magnetic Prandtl number Pr_{M,T} controls the rate of escape of vertical field from the disk; for Pr_{M,T} leq R/H vertical field diffuses outward before it can be advected inward by accretion. Here we measure field diffusion and angular momentum transport due to MHD turbulence in a shearing box, and thus Pr_{M,T}, by studying the evolution of a sinusoidal perturbation in the magnetic field that is injected into a turbulent background. We show that the perturbation is always stable, decays approximately exponentially, has decay rate proportional to k^2, and that the implied Pr_{M,T} ~ 1.

We numerically evolve turbulence driven by the magnetorotational instability (MRI) in a 3D, unstratified shearing box and study its structure using two-point correlation functions. We confirm Fromang and Papaloizou's result that shearing box models with zero net magnetic flux are not converged; the dimensionless shear stress $\alpha$ is proportional to the grid scale. We find that the two-point correlation of the magnetic field shows that it is composed of narrow filaments that are swept back by differential rotation into a trailing spiral. The correlation lengths along each of the correlation function principal axes decrease monotonically with the grid scale. For mean azimuthal field models, which we argue are more relevant to astrophysical disks than the zero net field models, we find that: $\alpha$ increases weakly with increasing resolution at fixed box size; $\alpha$ increases slightly as the box size is increased; $\alpha$ increases linearly with net field strength, confirming earlier results; the two-point correlation function of the magnetic field is resolved and converged, and is composed of narrow filaments swept back by the shear; the major axis of the two-point increases slightly as the box size is increased; these results are code independent, based on a comparison of ATHENA and ZEUS runs. The velocity, density, and magnetic fields decorrelate over scales larger than $\sim H$, as do the dynamical terms in the magnetic energy evolution equations. We conclude that MHD turbulence in disks is localized, subject to the limitations imposed by the absence of vertical stratification, the use of an isothermal equation of state, finite box size, finite run time, and finite resolution

In numerical models of thin astrophysical disks that use an Eulerian scheme, gas orbits supersonically through a fixed grid. As a result the time step is sharply limited by the Courant condition. Also, because the mean flow speed with respect to the grid varies with position, the truncation error varies systematically with position. For hydrodynamic (unmagnetized) disks an algorithm called FARGO has been developed that advects the gas along its mean orbit using a separate interpolation substep. This relaxes the constraint imposed by the Courant condition, which now depends only on the peculiar velocity of the gas, and results in a truncation error that is more nearly independent of position. This paper describes a FARGO-like algorithm suitable for evolving magnetized disks. Our method is second order accurate on a smooth flow and preserves the divergence-free constraint to machine precision. The main restriction is that the magnetic field must be discretized on a staggered mesh. We give a detailed description of an implementation of the code and demonstrate that it produces the expected results on linear and nonlinear problems. We also point out how the scheme might be generalized to make the integration of other supersonic/super-fast flows more efficient. Although our scheme reduces the variation of truncation error with position, it does not eliminate it. We show that the residual position dependence leads to characteristic radial variations in the density over long integrations.

Using the IRAM 30 m telescope, a mapping survey in optically thick and thin lines was performed towards 46 high mass star-forming regions. The sample includes UC H{\sc ii} precursors and UC H{\sc ii} regions. Seventeen sources are found to show "blue profiles", the expected signature of collapsing cores. The excess of sources with blue over red profiles ([$N_{\rm blue}$ -- $N_{\rm red}$]/$N_{\rm total}$) is 29% in the HCO$^+$ $J$=1--0 line, with a probability of 0.6% that this is caused by random fluctuations. UC H{\sc ii} regions show a higher excess (58%) than UC H{\sc ii} precursors (17%), indicating that material is still accreted after the onset of the UC H{\sc ii} phase. Similar differences in the excess of blue profiles as a function of evolutionary state are not observed in low mass star-forming regions. Thus, if confirmed for high mass star-forming sites, this would point at a fundamental difference between low- and high-mass star formation. Possible explanations are inadequate thermalization, stronger influence of outflows in massive early cores, larger gas reserves around massive stellar objects or different trigger mechanisms between low- and high- mass star formation.