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.

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.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

Stellar kinematics provides the key to understanding the formation process and dynamical evolution of stellar systems. Here, we present a kinematic study of the massive star-forming region W4 in the Cassiopeia OB6 association using the Gaia Data Release 2 and high-resolution optical spectra. This star-forming region is composed of a core cluster (IC 1805) and a stellar population distributed over 20 pc, which is a typical structural feature found in many OB associations. According to a classical model, this structural feature can be understood in the context of the dynamical evolution of a star cluster. The core-extended structure exhibits internally different kinematic properties. Stars in the core have an almost isotropic motion, and they appear to reach virial equilibrium given their velocity dispersion (0.9 +/- 0.3 km/s) comparable to that in a virial state (~0.8 km/s). On the other hand, the distributed population shows a clear pattern of radial expansion. From the N-body simulation for the dynamical evolution of a model cluster in subvirial state, we reproduce the observed structure and kinematics of stars. This model cluster experiences collapse for the first 2 Myr. Some members begin to radially escape from the cluster after the initial collapse, eventually forming a distributed population. The internal structure and kinematics of the model cluster appear similar to those of W4. Our results support the idea that the stellar population distributed over 20 pc in W4 originate from the dynamical evolution of IC 1805.

Gravitational interactions between a protoplanetary disk and its embedded planet is one of the formation mechanisms of gaps and rings found in recent ALMA observations. To quantify the gap properties measured in not only surface density but also rotational velocity profiles, we run two-dimensional hydrodynamic simulations of protoplanetary disks by varying three parameters: the mass ratio $q$ of a planet to a central star, the ratio of the disk scale height $h_p$ to the orbital radius $r_p$ of the planet, and the viscosity parameter $\alpha$. We find the gap depth $\delta_\Sigma$ in the gas surface density depends on a single dimensionless parameter $K\equiv q^2(h_p/r_p)^{-5}\alpha^{-1}$ as $\delta_\Sigma = (1 + 0.046K)^{-1}$, consistent with the previous results of Kanagawa et al. (2015a). The gap depth $\delta_V$ in the rotational velocity is given by $\delta_V= 0.007 (h_p/r_p) K^{1.38}/(1 +0.06K^{1.03})$. The gap width, in both surface density and rotational velocity, has a minimum of about $4.7 h_p$ when the planet mass $M_p$ is around the disk thermal mass $M_{\text{th}}$, while it increases in a power-law fashion as $M_p/M_{\text{th}}$ increases or decrease from unity. Such a minimum in the gap width arises because spirals from sub-thermal planets have to propagate before they shock the disk gas and open a gap. We compare our relations for the gap depth and width with the previous results, and discuss their applicability to observations.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

In this study we present the study of solar active regions based on BOAO vector magnetograms and $H\alpha$ filtergrams. With the new calibration method we analyzed BOAO vector magnetograms taken from the SOFT observational system to compare with those of other observing systems. In this study it has been demonstrated that (1) our longitudinal magnetogram matches very well the corresponding Mitaka's magnetogram to the extent that the maximum correlation yields r=0.962 between our re-scaled longitudinal magnetogram and the Mitaka's magnetogram; (2) according to a comparison of our magnetograms of AR 8422 with those taken at Mitaka solar observatory their longitudinal fields are very similar to each other while transverse fields are a little different possibly due to large noise level; (3) main features seen by our longitudinal magnetograms of AR 8422 and AR 8419 and the corresponding Kitt Peak magnetograms are very similar to each other; (4) time series of our vector magnetograms and H-alpha observations of AR 8419 during its flaring (M3.1/1B) activity show that the filament eruption followed the sheared inversion line of the quadrupolar configuration of sunspots, indicating that the flare should be associated with the quadrupolar field configuration and its interaction with new filament eruption. Finally, it may be concluded that the Solar Flare Telescope at BOAO works normally and it is ready to do numerous observational and theoretical works associated with solar activities such as flares.

The evolution of nonpotential characteristics of magnetic fields in AR 5747 is presented using Mees Solar Observatory magnetograms taken on Oct. 20, 1989 to Oct. 22, 1989. The active region showed such violent flaring activities during the observational span that strong X-ray flares took place including a 2B/X3 flare. The magnetogram data were obtained by the Haleakala Stokes Polarimeter which provides simultaneous Stokes profiles of the Fe I doublet 6301.5 and 6302.5. A nonlinear least square method was adopted to derive the magnetic field vectors from the observed Stokes profiles and a multi-step ambiguity solution method was employed to resolve the 180 degree ambiguity. From the ambiguity-resolved vector magnetograms, we have derived a set of physical quantities characterizing the field configuration, which are magnetic flux, vertical current density, magnetic shear angle, angular shear, magnetic free energy density and a measure of magnetic field discontinuity MAD (Maximum Angular Difference between two adjacent field vectors). In our results, all the physical parameters decreased with time, which implies that the active region was in a relaxation stage of its evolution. To examine the force-free characteristics of the field, we calculated the integrated Lorentz force and and also compared the longitudinal field component with the corresponding vertical current density. In this investigation, we found that the magnetic field in this active region was approximately linearly force-free throughout the observing period. The time variation of the linear force-free coefficient is consistent with the evolutionary trend of other nonpotentiality parameters. This suggests that the linear force-free coefficient could be a good indicator of the evolutionary status of active regions.

We present the evolution of magnetic nonpotentiality associated with an X-class flare in AR 6919 using a set of MSO (Mees Solar Observatory) magnetograms. The magnetogram data were obtained before and after the flare, using the Haleakala Stokes Polarimeter which provides simultaneous Stokes profiles of the Fe I doublet 6301.5 and 6302.5. A nonlinear least square method was adopted to derive the magnetic field configuration from the observed Stokes profiles and a multi-step ambiguity solution method was employed to resolve the 180 degree ambiguity. From the ambiguity-resolved vector magnetograms, we have derived a set of physical quantities characterizing the field configuration such as vertical current density, magnetic shear angle, angular shear and magnetic free energy density. We have examined their changes before and after the flare occurrence and obtained the following results. 1) There was a flux decrease in both polarities around a delta-type sunspot region, where an X-class flare occurred. 2) The vertical current near the delta sunspot region was strongly enhanced before the flare. 3) The magnetic shear near the delta sunspot region increased before the flare and then decreased after it. 4) The sum of magnetic energy density significantly decreased before the flare onset, implying that magnetic free energy was released through the flaring processes.