Search

The major satellites of Jupiter and Saturn are believed to have formed in circumplanetary discs, which orbit forming giant protoplanets. Gas and dust in CPDs have different distributions and affect each other by drag, which varies with grain size. Yet simulations of multiple dust grain sizes with separate dynamics have not been done before. We seek to assess how much dust of each grain size there is in circumplanetary discs. We run multifluid 3D hydrodynamical simulations including gas and four discrete grain sizes of dust from 1$\mu$m to 1mm, representing a continuous distribution. We consider a 1 $M_\mathrm{Jup}$ protoplanet embedded in a protoplanetary disc around a 1 $M_{\odot}$ star. Our results show a truncated MRN distribution at smaller grain sizes, which starts to tail off by $a=100\mu$m and is near zero at 1mm. Large dust grains, which hold most of the dust mass, have very inefficient accretion to the CPD, due to dust filtration. Therefore CPDs' dust masses must be small, with mass ratio ~ a few $\times 10^{-6}$ to the protoplanet. These masses and the corresponding millimetre opacities are in line with CPD fluxes observed to date., Comment: 14 pages, 15 figures; accepted for publication in MNRAS

Dust production is one of the more curious phenomena observed in massive binary systems with interacting winds. The high temperatures, UV photon flux and violent shocks should destroy any dust grains that condense. However, in some extreme cases dust production yields of approximately 30% of the total mass loss rate of the stellar winds have been observed. In order to better understand this phenomenon a parameter space exploration was performed using a series of numerical models of dust producing carbon phase Wolf-Rayet (WCd) systems. These models incorporated a passive scalar dust model simulating dust growth, destruction and radiative cooling. We find that reasonable dust yields were produced by these simulations. Significant changes in the dust yield were caused by changing the mass loss rates of the stars, with a greater mass loss rate contributing to increased dust yields. Similarly, a close orbit between the stars also resulted in higher dust yields. Finally, a high velocity wind shear, which induces Kelvin-Helmholtz (KH) instabilities and wind mixing, drastically increases the dust yields., 13 pages, 18 figures, 10 tables

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., 18 pages, accepted for publication in ApJ

Tao-Chung Ching et al., We present 850 μm dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the subfilaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1–10 pc scales. The magnetic fields revealed in the Planck data are well-aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6–1.0 mG in the DR21 filament and ∼0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by James Clerk Maxwell Telescope. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and subfilaments., This work is supported by National Natural Science Foundation of China (NSFC) grant Nos. 11988101, U1931117, 11725313, and 12073061 and the CAS International Partnership Program of Chinese Academy of Sciences grant No. 114A11KYSB20160008. T.-C.C. is funded by Chinese Academy of Sciences Taiwan Young Talent Program grant No. 2018TW2JB0002. T.-C.C. and C.E. were supported by Special Funding for Advanced Users, budgeted and administrated by Center for Astronomical Mega-Science (CAMS), Chinese Academy of Sciences. K.P. is a Royal Society University Research Fellow, supported by grant No. URF\R1\211322. D.J. is supported by the National Research Council of Canada and by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant. P.M.K. is supported by the Ministry of Science and Technology (MoST) in Taiwan through grants 109-2112-M-001-022 and 110-2112-M-001-057. C.W.L. is supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF-2019R1A2C1010851), and by the Korea Astronomy and Space Science Institute grant funded by the Korea government (MSIT; Project No. 2022-1-840-05). T.H. is supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) through the Mid-career Research Program (2019R1A2C1087045). W.K. was supported by the NRF grant funded by the MSIT (2021R1F1A1061794). C.E. acknowledges the financial support from grant RJF/2020/000071 as a part of Ramanujan Fellowship awarded by the Science and Engineering Research Board (SERB), Department of Science and Technology (DST), Govt. of India. F.P. acknowledges support from the Spanish State Research Agency (AEI) under grant No. PID2019-105552RB-C43. M.T. is supported by JSPS KAKENHI grant Nos. 18H05442, 15H02063, and 22000005. J.K. is supported by JSPS KAKENHI grant No. 19K14775. L.F. and F.K. acknowledge the support by the MoST in Taiwan through grant 107-2119-M-001-031-MY3 and Academia Sinica through grant AS-IA-106-M03. L.F. acknowledges the support by the MoST in Taiwan through grants 111-2811-M-005-007 and 109-2112-M-005-003-MY3. C.L.H.H. acknowledges the support of the NAOJ Fellowship and JSPS KAKENHI grants 18K13586 and 20K14527. F.K. is supported by the Spanish program Unidad de Excelencia María de Maeztu CEX2020-001058-M, financed by MCIN/AEI/10.13039/501100011033. K.Q. is partially supported by National Key R&D Program of China No. 2022YFA1603100, and acknowledges the National Natural Science Foundation of China (NSFC) grant U1731237. S.P.L. acknowledges grants from the Ministry of Science and Technology of Taiwan 106-2119-M-007-021-MY3 and 109-2112-M-007-010-MY3. Y.D. acknowledges the support of JSPS KAKENHI grants 25247016 and 18H01250. Y.S.D. is supported by the National Key R&D Program of China for grant No. 2022YFA1605300, and NSFC grants Nos. 12273051, 11933003.

Hwang et al., We present and analyze observations of polarized dust emission at 850 μm toward the central 1 × 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields in Star-forming Region Observations survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis–Chandrasekhar–Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from Herschel data and the C18O (J = 3 − 2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 ± 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 ± 0.02. Additionally, the mean Alfvén Mach number is 0.35 ± 0.01. This suggests that, in Mon R2, the magnetic fields provide resistance against large-scale gravitational collapse, and the magnetic pressure exceeds the turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically subcritical., J.H. is supported by the University of Science and Technology (UST) Overseas Training Program 2022, funded by the University of Science and Technology, Korea (No. 2022-017). K.P. is a Royal Society University Research Fellow, supported by grant number URF\R1\211322. C.W.L. was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF- 2019R1A2C1010851), and by the Korea Astronomy and Space Science Institute grant funded by the Korea government (MSIT; Project No. 2022-1-840-05). D.J. and J.d.F. are supported by the National Research Council of Canada and by Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grants. G.P. is supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF; 2020R1A6A3A01100208) W.K. was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT; NRF-2021R1F1A1061794). M.T. is supported by JSPS KAKENHI grant Nos. 18H05442, 15H02063, and 22000005. J.K. is supported by JSPS KAKENHI grant No. 19K14775. M.M. is supported by JSPS KAKENHI grant No. 20K03276. F.P. acknowledges support from the Spanish State Research Agency Agencia Estatal de Investigación (AEI) under grant number PID2019-105552RB-C43. S.P.L., H.Y.D., S.J.L., and M.Z.Y. acknowledge grants from the Ministry of Science and Technology of Taiwan 106-2119-M-007-021-MY3 and 109-2112-M-007-010-MY3. C.E. acknowledges the financial support from grant RJF/2020/000071 as a part of a Ramanujan Fellowship awarded by Science and Engineering Research Board (SERB), Department of Science and Technology (DST), Govt. of India. L.F. and F.K. acknowledge support from the Ministry of Science and Technology of Taiwan, under grant No. MoST107-2119-M-001-031- MY3 and from Academia Sinica under grant No. AS-IA-106-M03. L.F. acknowledges the support by the MoST in Taiwan through grants 111-2811-M-005-007 and 109-2112-M-005-003-MY3. F.K. acknowledges support from the Spanish program Unidad de Excelencia María de Maeztu CEX2020-001058-M, financed by MCIN/AEI/10.13039/501100011033. Y.S.D. would like to acknowledge the support from NSFC grants 12273051 and 10878003. M.R. is supported by the international Gemini Observatory, a program of NSFs NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation, on behalf of the Gemini partnership of Argentina, Brazil, Canada, Chile, the Republic of Korea, and the United States of America. The JCMT is operated by the East Asian Observatory on behalf of The National Astronomical Observatory of Japan; Academia Sinica Institute of Astronomy and Astrophysics; the Korea Astronomy and Space Science Institute; the Operation, Maintenance and Upgrading Fund for Astronomical Telescopes and Facility Instruments, budgeted from the Ministry of Finance of China. Additional funding support is provided by the Science and Technology Facilities Council of the United Kingdom and participating universities and organizations in the United Kingdom, Canada, and Ireland. Additional funds for the construction of SCUBA-2 were provided by the Canada Foundation for Innovation.

The wind collision region (WCR) in a colliding wind binary (CWB) is a particularly violent place, as such, it is surprising that it is also a region where significant quantities of interstellar dust can form. In extreme cases, approximately 30% of the total mass loss rate of a system can be converted into dust. These regions are poorly understood, as observation and simulation of these systems are difficult. In our previous paper we simulated dust growth in CWB systems using an advected scalar model and found our model to be suitable for qualitative study. For this paper we simulated the periodic dust forming CWB (WCd) system WR140 with our dust model, to determine how dust growth changes over the systems periastron passage. We found that dust production increases significantly at periastron passage, which is consistent with IR emission of the surrounding dusty shell. We also find that the dust production rate of the system decreases rapidly as the stars recede from each other, though the rate of decrease is significantly lower than the rate of increase during periastron passage. This was found to be due to strong cooling and its associated thermal instabilities, resulting in cool, high-density pockets of gas in the WCR where dust forms. The WCR also shows a degree of hysteresis, resulting in a radiative post-shock flow even when the stars are separated enough for the region to behave adiabatically., Comment: 9 pages, 9 figures, 5 tables

We present the results of simultaneous 450 $\mu$m and 850 $\mu$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could be due to a rotating toroidal disk-like structure and a bipolar outflow originating from the central young stellar object, IRS 3. Using the modified Davis-Chandrasekhar-Fermi method, we obtain a plane-of-sky magnetic field strength of 563$\pm$421 $\mu$G in the central $\sim$0.12 pc region from 850 $\mu$m polarization data. The corresponding magnetic energy density of 2.04$\times$10$^{-8}$ erg cm$^{-3}$ is comparable to the turbulent and gravitational energy densities in the region. We find that the magnetic field direction is very well aligned with the whole of the IRS 3 bipolar outflow structure. We find that the median value of polarization fractions, 3.0 \%, at 450 $\mu$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $\mu$m. The trend could be due to the better alignment of warmer dust in the strong radiation environment. We also find that polarization fractions decrease with intensity at both wavelengths, with slopes, determined by fitting a Rician noise model, of $0.59 \pm 0.03$ at 450 $\mu$m and $0.36 \pm 0.04$ at 850 $\mu$m, respectively. We think that the shallow slope at 850 $\mu$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects., Comment: 17 pages, 11 figures, 4 tables

We have obtained sensitive dust continuum polarization observations at 850 μm in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope as part of the B-fields in STar-forming Region Observations (BISTRO) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution (∼2000 au or ∼0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis–Chandrasekhar–Fermi method, we estimate the B-field strengths in K04166, K04169, and Miz-8b to be 38 ± 14, 44 ± 16, and 12 ± 5 μG, respectively. These cores show distinct mean B-field orientations. The B-field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B-field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. The B-field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B-field and not well correlated with other axes. In contrast, Miz-8b exhibits a disordered B-field that shows no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B-field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B-field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux.

We report the first high spatial resolution measurement of magnetic fields surrounding LkH$\alpha$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $\mu$m is found to be mostly associated with the red-shifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis-Chandrasekhar-Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of $B_{\rm POS}\sim 115$ $\mathrm{\mu}$G. The mass-to-magnetic-flux ratio in critical value units, $\lambda\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$\alpha$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $\delta B/B_0\sim 0.3$ implies that the large scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power-law with an index of $\alpha=0.82\pm0.03$, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH$\alpha$ 101) in the region. The magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques are potential of explaining such a decreasing trend., Comment: 25 pages, 19 figures, Accepted for publication in The Astrophysical Journal

We present new observations of the active star formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope 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 × 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 multiscale 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 toward the filament network., Erratum published in volume 899, issue 1, 5 p., April 10 2021. DOI: 10.3847/1538-4357/abf2b3

In the published article, we presented high-resolution polarimetry data obtained by using JCMT SCUBA-2/POL-2, and compared them with the larger-scale magnetic structure observed by Planck (Planck Collaboration et al. 2020). There was a miscalculation in the analysis of the Planck data for comparison, and the mean position angle of the Planck magnetic field should be corrected from -40° ± 7° to -48° ± 6°. Thus, we replace the descriptions in the published article as follows. In Section 4.2, paragraph 3, the second sentence should read: "The differences between the orientations of IRAS 4A, IRAS 4B, and IRAS 2A are not statistically significant, but they do show significantly different orientations from those of the global B-field observed by Planck (-48° ± 6°; Section 4.1)." In Section 4.1, paragraph 3, the second line should read: "The Planck B-field orientation shows a smoothly and slowly varying field distribution with a position angle of -48° ± 6° in our observed NGC 1333 area." Accordingly, we replace Figures 4, 5, 9, 10, and 17 to reflect the correct Planck data. As described above, this error of modest magnitude is related only to our derivation of the Planck polarization angle, and the JCMT observation results are unaffected. Therefore, all conclusions drawn in the published article are unchanged even after the above correction is applied. The 1 pc scale magnetic field observed by Planck shows a smooth distribution, and the interstellar magnetic field in molecular clouds increases the complexity significantly on the scale of less than 1 pc.

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 $, Comment: 24 pages, 15 figures, accepted for publication in ApJ

著者人数: 122名 (所属. 宇宙航空研究開発機構宇宙科学研究所(JAXA)(ISAS): 權, 靜美; 中西, 裕之; 中川, 貴雄）, Accepted: 2018-03-26, 資料番号: SA1170332000

著者人数: 124名（所属. 宇宙航空研究開発機構宇宙科学研究所(JAXA)(ISAS): 權, 靜美), Accepted: 2018-05-10, 資料番号: SA1180101000

We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation sub-grid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold; (2) "Magnetically-Regulated," i.e., fixed efficiency per free-fall time in regions that are magnetically supercritical. Variations of parameters associated with these models are also explored. In the non-colliding simulations, the overall level of star formation is sensitive to model parameter choices that relate to effective density thresholds. In the GMC collision simulations, the final star formation rates and efficiencies are relatively independent of these parameters. Between non-colliding and colliding cases, we compare the morphologies of the resulting star clusters, properties of star-forming gas, time evolution of the star formation rate (SFR), spatial clustering of the stars, and resulting kinematics of the stars in comparison to the natal gas. We find that typical collisions, by creating larger amounts of dense gas, trigger earlier and enhanced star formation, resulting in 10 times higher SFRs and efficiencies. The star clusters formed from GMC collisions show greater spatial sub-structure and more disturbed kinematics., 22 pages, 10 figures, Accepted for publication in ApJ

We have used the AMR hydrodynamic code, MG, to perform idealised 3D MHD simulations of the formation of clumpy and filamentary structure in a thermally unstable medium without turbulence. A stationary thermally unstable spherical diffuse atomic cloud with uniform density in pressure equilibrium with low density surroundings was seeded with random density variations and allowed to evolve. A range of magnetic field strengths threading the cloud have been explored, from beta=0.1 to beta=1.0 to the zero magnetic field case (beta=infinity), where beta is the ratio of thermal pressure to magnetic pressure. Once the density inhomogeneities had developed to the point where gravity started to become important, self-gravity was introduced to the simulation. With no magnetic field, clouds and clumps form within the cloud with aspect ratios of around unity, whereas in the presence of a relatively strong field (beta=0.1) these become filaments, then evolve into interconnected corrugated sheets that are predominantly perpendicular to the magnetic field. With magnetic and thermal pressure equality (beta=1.0), filaments, clouds and clumps are formed. At any particular instant, the projection of the 3D structure onto a plane parallel to the magnetic field, i.e. a line of sight perpendicular to the magnetic field, resembles the appearance of filamentary molecular clouds. The filament densities, widths, velocity dispersions and temperatures resemble those observed in molecular clouds. In contrast, in the strong field case beta=0.1, projection of the 3D structure along a line of sight parallel to the magnetic field reveals a remarkably uniform structure., 17 pages, 10 figures. Accepted to MNRAS, minor changes and additions to text, primarily in Section 6

We investigate giant molecular cloud (GMCs) collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of three dimensional, magnetohydrodynamics (MHD), adaptive mesh refinement (AMR) simulations to study how cloud collisions trigger formation of dense filaments and clumps. Heating and cooling functions are implemented based on photo-dissociation region (PDR) models that span the atomic to molecular transition and can return detailed diagnostic information. The clouds are initialized with supersonic turbulence and a range of magnetic field strengths and orientations. Collisions at various velocities and impact parameters are investigated. Comparing and contrasting colliding and non-colliding cases, we characterize morphologies of dense gas, magnetic field structure, cloud kinematic signatures, and cloud dynamics. We present key observational diagnostics of cloud collisions, especially: relative orientations between magnetic fields and density structures, like filaments; $^{13}$CO($J$=2-1), $^{13}$CO($J$=3-2), and $^{12}$CO($J$=8-7) integrated intensity maps and spectra; and cloud virial parameters. We compare these results to observed Galactic clouds., 27 pages, 17 figures, Accepted for publication in ApJ

We utilize magnetohydrodynamic (MHD) simulations to develop a numerical model for GMC-GMC collisions between nearly magnetically critical clouds. The goal is to determine if, and under what circumstances, cloud collisions can cause pre-existing magnetically subcritical clumps to become supercritical and undergo gravitational collapse. We first develop and implement new photodissociation region (PDR) based heating and cooling functions that span the atomic to molecular transition, creating a multiphase ISM and allowing modeling of non-equilibrium temperature structures. Then in 2D and with ideal MHD, we explore a wide parameter space of magnetic field strength, magnetic field geometry, collision velocity, and impact parameter, and compare isolated versus colliding clouds. We find factors of ~2-3 increase in mean clump density from typical collisions, with strong dependence on collision velocity and magnetic field strength, but ultimately limited by flux-freezing in 2D geometries. For geometries enabling flow along magnetic field lines, greater degrees of collapse are seen. We discuss observational diagnostics of cloud collisions, focussing on 13CO(J=2-1), 13CO(J=3-2), and 12CO(J=8-7) integrated intensity maps and spectra, which we synthesize from our simulation outputs. We find the ratio of J=8-7 to lower-J emission is a powerful diagnostic probe of GMC collisions., Comment: 22 pages, 23 figures. Accepted for publication in ApJ

The regulation of galactic-scale star formation rates (SFRs) is a basic problem for theories of galaxy formation and evolution: which processes are responsible for making observed star formation rates so inefficient compared to maximal rates of gas content divided by dynamical timescale? Here we study the effect of magnetic fields of different strengths on the evolution of giant molecular clouds (GMCs) within a kiloparsec patch of a disk galaxy and resolving scales down to $\simeq0.5\:{\rm{pc}}$. Including an empirically motivated prescription for star formation from dense gas ($n_{\rm{H}}>10^5\:{\rm{cm}^{-3}}$) at an efficiency of 2\% per local free-fall time, we derive the amount of suppression of star formation by magnetic fields compared to the nonmagnetized case. We find GMC fragmentation, dense clump formation and SFR can be significantly affected by the inclusion of magnetic fields, especially in our strongest investigated $B$-field case of $80\:{\rm��}$G. However, our chosen kpc-scale region, extracted from a global galaxy simulation, happens to contain a starbursting cloud complex that is only modestly affected by these magnetic fields and likely requires internal star formation feedback to regulate its SFR., 7 pages, 5 figures. Accepted by ApJ Letters

We examine the role of the gravitational instability in an isothermal, self-gravitating layer threaded by magnetic fields on the formation of filaments and dense cores. Using numerical simulation we follow the non-linear evolution of a perturbed equilibrium layer. The linear evolution of such a layer is described in the analytic work of Nagai et al (1998). We find that filaments and dense cores form simultaneously. Depending on the initial magnetic field, the resulting filaments form either a spiderweb-like network (for weak magnetic fields) or a network of parallel filaments aligned perpendicular to the magnetic field lines (for strong magnetic fields). Although the filaments are radially collapsing, the density profile of their central region (up to the thermal scale height) can be approximated by a hydrodynamical equilibrium density structure. Thus, the magnetic field does not play a significant role in setting the density distribution of the filaments. The density distribution outside of the central region deviates from the equilibrium. The radial column density distribution is then flatter than the expected power law of $r^{-4}$ and similar to filament profiles observed with Herschel. Our results does not explain the near constant filament width of $\sim 0.1$pc. However, our model does not include turbulent motions. It is expected that accretion-driven amplification of these turbulent motions provides additional support within the filaments against gravitational collapse. Finally, we interpret the filamentary network of the massive star forming complex G14.225-0.506 in terms of the gravitational instability model and find that the properties of the complex are consistent with being formed out of an unstable layer threaded by a strong, parallel magnetic field., Comment: Accepted by ApJ, 14 pages, 22 figures

We present hydrodynamic simulations of self-gravitating dense gas in a galactic disk, exploring scales ranging from 1 kpc down to $\sim 0.1$~pc. Our primary goal is to understand how dense filaments form in Giant Molecular Clouds (GMCs). These structures, often observed as Infrared Dark Clouds (IRDCs) in the Galactic plane, are thought to be the precursors to massive stars and star clusters, so their formation may be the rate limiting step controlling global star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. Our study follows on from Van Loo et al. (2013, Paper I), which carried out simulations to 0.5~pc resolution and examined global aspects of the formation of dense gas clumps and the resulting star formation rate. Here, using our higher resolution, we examine the detailed structural, kinematic and dynamical properties of dense filaments and clumps, including mass surface density ($\Sigma$) probability distribution functions, filament mass per unit length and its dispersion, lateral $\Sigma$ profiles, filament fragmentation, filament velocity gradients and infall, and degree of filament and clump virialization. Where possible, these properties are compared to observations of IRDCs. By many metrics, especially too large mass fractions of high $\Sigma>1\:{\rm g\:cm^{-2}}$ material, too high mass per unit length dispersion due to dense clump formation, too high velocity gradients and too high velocity dispersion for a given mass per unit length, the simulated filaments differ from observed IRDCs. We thus conclude that IRDCs do not form from global fast collapse of GMCs. Rather, we expect IRDC formation and collapse is slowed significantly by the influence of dynamically important magnetic fields, which may thus play a crucial role in regulating galactic star formation rates., Comment: 26 pages, 15 figures, accepted to ApJ

We present hydrodynamic simulations of the evolution of self-gravitating dense gas on scales of 1 kiloparsec down to < parsec in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These structures are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. We follow the thermal evolution of the gas down to ~5K using extinction-dependent heating and cooling functions. We do not yet include magnetic fields or localized stellar feedback, so the evolution of the GMCs and clumps is determined solely by self-gravity balanced by thermal and turbulent pressure support and the large scale galactic shear. While cloud structures and densities change significantly during the simulation, GMC virial parameters remain mostly above unity for time scales exceeding the free-fall time of GMCs indicating that energy from galactic shear and large-scale cloud motions continuously cascades down to and within the GMCs. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but even this yields global star formation rates that are about two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to over-production of dense gas clumps. We expect a combination of magnetic support and localized stellar feedback is required to inhibit dense clump formation to ~1% of the rate that results from the nonmagnetic, zero-feedback limit., Accepted by ApJ, 14 pages and 16 figures

All stars are born in molecular clouds, and most in giant molecular clouds (GMCs), which thus set the star formation activity of galaxies. We first review their observed properties, including measures of mass surface density, Sigma, and thus mass, M. We discuss cloud dynamics, concluding most GMCs are gravitationally bound. Star formation is highly clustered within GMCs, but overall is very inefficient. We compare properties of star-forming clumps with those of young stellar clusters (YSCs). The high central densities of YSCs may result via dynamical evolution of already-formed stars during and after star cluster formation. We discuss theoretical models of GMC evolution, especially addressing how turbulence is maintained, and emphasizing the importance of GMC collisions. We describe how feedback limits total star formation efficiency, epsilon, in clumps. A turbulent and clumpy medium allows higher epsilon, permitting formation of bound clusters even when escape speeds are less than the ionized gas sound speed., Comment: Invited review, IAU Symp. 292 Molecular Gas, Dust, and Star Formation in Galaxies

We report the results of a multiwavelength study of the non-thermal radio emitter HD 168112 (O5.5III(f+)). The detailed analysis of two quasi-simultaneous XMM-Newton and VLA observations reveals strong variability of this star both in the X-ray and radio ranges. The X-ray observations separated by five months reveal a decrease of the X-ray flux of ~30%. The radio emission on the other hand increases by a factor 5-7 between the two observations obtained roughly simultaneously with the XMM-Newton pointings. The X-ray data reveal a hard emission that is most likely produced by a thermal plasma at kT ~2-3 keV while the VLA data confirm the non-thermal status of this star in the radio waveband. Comparison with archive X-ray and radio data confirms the variability of this source in both wavelength ranges over a yet ill defined time scale. The properties of HD 168112 in the X-ray and radio domain point towards a binary system with a significant eccentricity and an orbital period of a few years. However, our optical spectra reveal no significant changes of the star's radial velocity suggesting that if HD 168112 is indeed a binary, it must be seen under a fairly low inclination. Based on observations with XMM-Newton, an ESA Science Mission with instruments and contributions directly funded by ESA Member states and the USA (NASA). Also based on observations collected with the VLA, an instrument of the National Radio Astronomy Observatory, which is a facility of the National Science Foundation operated by Associated Universities, Inc. Optical data were collected at the European Southern Observatory (La Silla, Chile), and at the Observatorio Astronómico Nacional of San Pedro Mártir (Mexico).