Showing 1–50 of 52 results for author: Rawlings, M G

First Very Long Baseline Interferometry Detections at 870μm

Authors: Alexander W. Raymond, Sheperd S. Doeleman, Keiichi Asada, Lindy Blackburn, Geoffrey C. Bower, Michael Bremer, Dominique Broguiere, Ming-Tang Chen, Geoffrey B. Crew, Sven Dornbusch, Vincent L. Fish, Roberto Garc�a, Olivier Gentaz, Ciriaco Goddi, Chih-Chiang Han, Michael H. Hecht, Yau-De Huang, Michael Janssen, Garrett K. Keating, Jun Yi Koay, Thomas P. Krichbaum, Wen-Ping Lo, Satoki Matsushita, Lynn D. Matthews, James M. Moran , et al. (254 additional authors not shown)

Abstract: The first very long baseline interferometry (VLBI) detections at 870$μ$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescop… ▽ More The first very long baseline interferometry (VLBI) detections at 870$μ$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$λ$ corresponding to an angular resolution, or fringe spacing, of 19$μ$as. The Allan deviation of the visibility phase at 870$μ$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$μ$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time. △ Less

Submitted 9 October, 2024; originally announced October 2024.

Comments: Corresponding author: S. Doeleman

Journal ref: The Astronomical Journal, Volume 168, Issue 3, id.130, 19 pp. 2024

Relative Alignments Between Magnetic Fields, Velocity Gradients, and Dust Emission Gradients in NGC 1333

Authors: Michael Chun-Yuan Chen, Laura M. Fissel, Sarah I. Sadavoy, Erik Rosolowsky, Yasuo Doi, Doris Arzoumanian, Pierre Bastien, Simon Coud�, James Di Francesco, Rachel Friesen, Ray S. Furuya, Jihye Hwang, Shu-ichiro Inutsuka, Doug Johnstone, Janik Karoly, Jungmi Kwon, Woojin Kwon, Valentin J. M. Le Gouellec, Hong-Li Liu, Steve Mairs, Takashi Onaka, Kate Pattle, Mark G. Rawlings, Mehrnoosh Tahani, Motohide Tamura , et al. (1 additional authors not shown)

Abstract: Magnetic fields play an important role in shaping and regulating star formation in molecular clouds. Here, we present one of the first studies examining the relative orientations between magnetic ($B$) fields and the dust emission, gas column density, and velocity centroid gradients on the 0.02 pc (core) scales, using the BISTRO and VLA+GBT observations of the NGC 1333 star-forming clump. We quant… ▽ More Magnetic fields play an important role in shaping and regulating star formation in molecular clouds. Here, we present one of the first studies examining the relative orientations between magnetic ($B$) fields and the dust emission, gas column density, and velocity centroid gradients on the 0.02 pc (core) scales, using the BISTRO and VLA+GBT observations of the NGC 1333 star-forming clump. We quantified these relative orientations using the Project Rayleigh Statistic (PRS) and found preferential global parallel alignment between the $B$ field and dust emission gradients, consistent with large-scale studies with Planck. No preferential global alignments, however, are found between the $B$ field and velocity gradients. Local PRS calculated for subregions defined by either dust emission or velocity coherence further revealed that the $B$ field does not preferentially align with dust emission gradients in most emission-defined subregions, except in the warmest ones. The velocity-coherent structures, on the other hand, also showed no preferred $B$ field alignments with velocity gradients, except for one potentially bubble-compressed region. Interestingly, the velocity gradient magnitude in NGC 1333 ubiquitously features prominent ripple-like structures that are indicative of magnetohydrodynamic (MHD) waves. Finally, we found $B$ field alignments with the emission gradients to correlate with dust temperature and anticorrelate with column density, velocity dispersion, and velocity gradient magnitude. The latter two anticorrelations suggest that alignments between gas structures and $B$ fields can be perturbed by physical processes that elevate velocity dispersion and velocity gradients, such as infall, accretions, and MHD waves. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: Accepted to the Monthly Notices of the Royal Astronomical Society Main Journal

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. R�der, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

On the Scarcity of Dense Cores ($n>10^{5}$ cm$^{-3}$) in High Latitude Planck Galactic Cold Clumps

Authors: Fengwei Xu, Ke Wang, Tie Liu, David Eden, Xunchuan Liu, Mika Juvela, Jinhua He, Doug Johnstone, Paul Goldsmith, Guido Garay, Yuefang Wu, Archana Soam, Alessio Traficante, Isabelle Ristorcelli, Edith Falgarone, Huei-Ru Vivien Chen, Naomi Hirano, Yasuo Doi, Woojin Kwon, Glenn J. White, Anthony Whitworth, Patricio Sanhueza, Mark G. Rawlings, Dana Alina, Zhiyuan Ren , et al. (12 additional authors not shown)

Abstract: High-latitude ($|b|>30^{\circ}$) molecular clouds have virial parameters that exceed 1, but whether these clouds can form stars has not been studied systematically. Using JCMT SCUBA-2 archival data, we surveyed 70 fields that target high-latitude Planck galactic cold clumps (HLPCs) to find dense cores with density of $10^{5}$-$10^{6}$ cm$^{-3}$ and size of $<0.1$ pc. The sample benefits from both… ▽ More High-latitude ($|b|>30^{\circ}$) molecular clouds have virial parameters that exceed 1, but whether these clouds can form stars has not been studied systematically. Using JCMT SCUBA-2 archival data, we surveyed 70 fields that target high-latitude Planck galactic cold clumps (HLPCs) to find dense cores with density of $10^{5}$-$10^{6}$ cm$^{-3}$ and size of $<0.1$ pc. The sample benefits from both the representativeness of the parent sample and covering densest clumps at the high column density end ($>1\times10^{21}$ cm$^{-2}$). At an average noise rms of 15 mJy/beam, we detected Galactic dense cores in only one field, G6.04+36.77 (L183), while also identifying 12 extragalactic objects and two young stellar objects. Compared to the low-latitude clumps, dense cores are scarce in HLPCs. With synthetic observations, the densities of cores are constrained to be $n_c\lesssim10^5$ cm$^{-3}$, should they exist in HLPCs. Low-latitude clumps, Taurus clumps, and HLPCs form a sequence where a higher virial parameter corresponds to a lower dense core detection rate. If HLPCs were affected by the Local Bubble, the scarcity should favor turbulence-inhibited rather than supernova-driven star formation. Studies of the formation mechanism of the L183 molecular cloud are warranted. △ Less

Submitted 22 February, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Comments: 9 pages for the main text. 4 figures, 1 table. Published in Astrophysical Journal Letter

Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations

Authors: Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, Fr�d�rick Poidevin, Mehrnoosh Tahani, Takashi Onaka, Mark G. Rawlings, Eun Jung Chung , et al. (132 additional authors not shown)

Abstract: We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from… ▽ More We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

The magnetic field in colliding filaments G202.3+2.5

Authors: Qi-Lao Gu, Tie Liu, Pak Shing Li, Zhi-Qiang Shen, Xunchuan Liu, Junhao Liu, Xing Lu, Julien Montillaud, Sihan Jiao, Mika Juvela, Mark G. Rawlings, Qizhou Zhang, Patrick Koch, Isabelle Ristorcelli, Jean-Sébastien Carriere, David Eden, Zhiyuan Ren, Ken'ichi Tatematsu, Naomi Hirano, Qiu-yi Luo, Xiaofeng Mai, Namitha Issac

Abstract: We observe the magnetic field morphology towards a nearby star-forming filamentary cloud, G202.3+2.5, by the JCMT/POL-2 850 μm thermal dust polarization observation with an angular resolution of 14.4" (~0.053 pc). The average magnetic field orientation is found to be perpendicular to the filaments while showing different behaviors in the four subregions, suggesting various effects from filaments'… ▽ More We observe the magnetic field morphology towards a nearby star-forming filamentary cloud, G202.3+2.5, by the JCMT/POL-2 850 μm thermal dust polarization observation with an angular resolution of 14.4" (~0.053 pc). The average magnetic field orientation is found to be perpendicular to the filaments while showing different behaviors in the four subregions, suggesting various effects from filaments' collision in these subregions. With the kinematics obtained by N2H+ observation by IRAM, we estimate the plane-of-sky (POS) magnetic field strength by two methods, the classical Davis-Chandrasekhar-Fermi (DCF) method and the angular dispersion function (ADF) method, B_{pos,dcf} and B_{pos,adf} are ~90 μG and ~53 μG. We study the relative importance between the gravity (G), magnetic field (B) and turbulence (T) in the four subregions, find G > T > B, G >= T > B, G ~ T > B and T > G > B in the north tail, west trunk, south root and east wing, respectively. In addition, we investigate the projection effect on the DCF and ADF methods based on a similar simulation case and find the 3D magnetic field strength may be underestimated by a factor of ~3 if applying the widely-used statistical B_{pos}-to-B_{3D} factor when using DCF or ADF method, which may further underestimate/overestimate related parameters. △ Less

Submitted 12 January, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: Accepted by ApJ. 20 pages, 9 figures

Charting Circumstellar Chemistry of Carbon-rich AGB Stars: I. ALMA 3 mm spectral surveys

Authors: R. Unnikrishnan, E. De Beck, L. A. Nyman, H. Olofsson, W. H. T. Vlemmings, D. Tafoya, M. Maercker, S. B. Charnley, M. A. Cordiner, I. de Gregorio, E. Humphreys, T. J. Millar, M. G. Rawlings

Abstract: AGB stars are major contributors to the chemical enrichment of the ISM through nucleosynthesis and extensive mass loss. Most of our current knowledge of AGB atmospheric and circumstellar chemistry, in particular in a C-rich environment, is based on observations of the carbon star IRC+10216. We aim to obtain a more generalised understanding of the chemistry in C-rich AGB CSEs by studying a sample o… ▽ More AGB stars are major contributors to the chemical enrichment of the ISM through nucleosynthesis and extensive mass loss. Most of our current knowledge of AGB atmospheric and circumstellar chemistry, in particular in a C-rich environment, is based on observations of the carbon star IRC+10216. We aim to obtain a more generalised understanding of the chemistry in C-rich AGB CSEs by studying a sample of three carbon stars, IRAS15194-5115, IRAS15082-4808, and IRAS07454-7112, and test the archetypal status often attributed to IRC+10216. We performed spatially resolved, unbiased spectral surveys in ALMA Band 3. We identify a total of 132 rotational transitions from 49 molecular species. There are two main morphologies of the brightness distributions: centrally-peaked (e.g. HCN) and shell-like (e.g. C$_2$H). We estimated the sizes of the molecular emitting regions using azimuthally-averaged radial profiles of the line brightness distributions, and derived abundance estimates. Of the shell distributions, the cyanopolyynes peak at slightly smaller radii than the hydrocarbons, and CN and HNC show the most extended emission. The emitting regions for each species are the smallest for IRAS07454-7112. We find that, within the uncertainties of the analysis, the three stars present similar abundances for most species, also compared to IRC+10216. We find that SiO is more abundant in our three stars compared to IRC+10216. Our estimated isotopic ratios match well the literature values for the sources. The observed circumstellar chemistry appears very similar across our sample and compared to that of IRC+10216, both in terms of the relative location of the emitting regions and molecular abundances. This implies that, to a first approximation, the chemical models tailored to IRC+10216 are able to reproduce the observed chemistry in C-rich envelopes across roughly an order of magnitude in wind density. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 21 pages. 13 figures and 8 tables in the main text. 5 appendices contain additional figures and tables. Appendix tables are available in electronic form at the CDS, along with the reprocessed ALMA cubes and spectra, at http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ or via anonymous ftp to cdsarc.u-strasbg.fr(130.79.128.5). Accepted for publication in Astronomy & Astrophysics

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Discovery of an extremely dense and compact object embedded in the prestellar core G208.68-19.92-N2

Authors: Naomi Hirano, Dipen Sahu, Sheng-Yaun Liu, Tie Liu, Ken'ichi Tatematsu, Somnath Dutta, Shanghuo Li, Chin-Fei Lee, Pak Shing Li, Shih-Ying Hsu, Sheng-Jun Lin, Doug Johnstone, Leonardo Bronfman, Huei-Ru Vivien Chen, David J. Eden, Yi-Jehng Kuan, Woojin Kwon, Chang Won Lee, Hong-Li Liu, Mark G. Rawlings, Isabelle Ristorcelli, Alessio Traficante

Abstract: The internal structure of the prestellar core G208.68-19.02-N2 (G208-N2) in the Orion Molecular Cloud 3 (OMC-3) region has been studied with the Atacama Large Millimeter/submillimeter Array (ALMA). The dust continuum emission revealed a filamentary structure with a length of $\sim$5000 au and an average H$_2$ volume density of $\sim$6 $\times$ 10$^7$ cm$^{-3}$. At the tip of this filamentary struc… ▽ More The internal structure of the prestellar core G208.68-19.02-N2 (G208-N2) in the Orion Molecular Cloud 3 (OMC-3) region has been studied with the Atacama Large Millimeter/submillimeter Array (ALMA). The dust continuum emission revealed a filamentary structure with a length of $\sim$5000 au and an average H$_2$ volume density of $\sim$6 $\times$ 10$^7$ cm$^{-3}$. At the tip of this filamentary structure, there is a compact object, which we call a ``nucleus", with a radius of $\sim$150--200 au and a mass of $\sim$0.1 M$_{\odot}$. The nucleus has a central density of $\sim$2 $\times$ 10$^9$ cm$^{-3}$ with a radial density profile of $r^{-1.87{\pm}0.11}$. The density scaling of the nucleus is $\sim$3.7 times higher than that of the singular isothermal sphere. This as well as the very low virial parameter of 0.39 suggest that the gravity is dominant over the pressure everywhere in the nucleus. However, there is no sign of CO outflow localized to this nucleus. The filamentary structure is traced by the N$_2$D$^+$ 3--2 emission, but not by the C$^{18}$O 2--1 emission, implying the significant CO depletion due to high density and cold temperature. Toward the nucleus, the N$_2$D$^+$ also shows the signature of depletion. This could imply either the depletion of the parent molecule, N$_2$, or the presence of the embedded very-low luminosity central source that could sublimate the CO in the very small area. The nucleus in G208-N2 is considered to be a prestellar core on the verge of first hydrostatic core (FHSC) formation or a candidate for the FHSC. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: 27 pages, 16 figures

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): The Warm-Envelope Origin of Hot Corinos

Authors: Shih-Ying Hsu, Sheng-Yuan Liu, Doug Johnstone, Tie Liu, Leonardo Bronfman, Huei-Ru Vivien Chen, Somnath Dutta, David J. Eden, Neal J. Evans II, Naomi Hirano, Mika Juvela, Yi-Jehng Kuan, Woojin Kwon, Chin-Fei Lee, Chang Won Lee, Jeong-Eun Lee, Shanghuo Li, Chun-Fan Liu, Xunchuan Liu, Qiuyi Luo, Sheng-Li Qin, Mark G. Rawlings, Dipen Sahu, Patricio Sanhueza, Hsien Shang , et al. (2 additional authors not shown)

Abstract: Hot corinos are of great interest due to their richness in interstellar complex organic molecules (COMs) and the consequent potential prebiotic connection to solar-like planetary systems. Recent surveys have reported an increasing number of hot corino detections in Class 0/I protostars; however, the relationships between their physical properties and the hot-corino signatures remain elusive. In th… ▽ More Hot corinos are of great interest due to their richness in interstellar complex organic molecules (COMs) and the consequent potential prebiotic connection to solar-like planetary systems. Recent surveys have reported an increasing number of hot corino detections in Class 0/I protostars; however, the relationships between their physical properties and the hot-corino signatures remain elusive. In this study, our objective is to establish a general picture of the detectability of the hot corinos by identifying the origin of the hot-corino signatures in the sample of young stellar objects (YSOs) obtained from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We apply spectral energy distribution (SED) modeling to our sample and identify the physical parameters of the modeled YSOs directly, linking the detection of hot-corino signatures to the envelope properties of the YSOs. Imaging simulations of the methanol emission further support this scenario. We, therefore, posit that the observed COM emission originates from the warm inner envelopes of the sample YSOs, based on both the warm region size and the envelope density profile. The former is governed by the source luminosity and is additionally affected by the disk and cavity properties, while the latter is related to the evolutionary stages. This scenario provides a framework for detecting hot-corino signatures toward luminous Class 0 YSOs, with fewer detections observed toward similarly luminous Class I sources. △ Less

Submitted 11 August, 2023; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: 28 pages, 11 figures

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): A forming quadruple system with continuum `ribbons' and intricate outflows

Authors: Qiu-yi Luo, Tie Liu, Aaron T. Lee, Stella S. R. Offner, James di Francesco, Doug Johnstone, Mika Juvela, Paul F. Goldsmith, Sheng-Li Qin, Xiaofeng Mai, Xun-chuan Liu, Patricio Sanhueza, Feng-Wei Xu, Ken'ichi Tatematsu, Somnath Dutta, Huei-Ru Vivien Chen, Shanghuo Li, Aiyuan Yang, Sheng-Yuan Liu, Chin-Fei Lee, Naomi Hirano, Chang Won Lee, Dipen Sahu, Hsien Shang, Shih-Ying Hsu , et al. (9 additional authors not shown)

Abstract: One of the most poorly understood aspects of low-mass star formation is how multiple-star systems are formed. Here we present the results of Atacama Large Millimeter/submillimeter Array (ALMA) Band-6 observations towards a forming quadruple protostellar system, G206.93-16.61E2, in the Orion B molecular cloud. ALMA 1.3 mm continuum emission reveals four compact objects, of which two are Class I you… ▽ More One of the most poorly understood aspects of low-mass star formation is how multiple-star systems are formed. Here we present the results of Atacama Large Millimeter/submillimeter Array (ALMA) Band-6 observations towards a forming quadruple protostellar system, G206.93-16.61E2, in the Orion B molecular cloud. ALMA 1.3 mm continuum emission reveals four compact objects, of which two are Class I young stellar objects (YSOs), and the other two are likely in prestellar phase. The 1.3 mm continuum emission also shows three asymmetric ribbon-like structures that are connected to the four objects, with lengths ranging from $\sim$500 au to $\sim$2200 au. By comparing our data with magneto-hydrodynamic (MHD) simulations, we suggest that these ribbons trace accretion flows and also function as gas bridges connecting the member protostars. Additionally, ALMA CO J=2-1 line emission reveals a complicated molecular outflow associated with G206.93-16.61E2 with arc-like structures suggestive of an outflow cavity viewed pole-on. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: The paper was accepted by APJL

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

Abstract: We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles… ▽ More We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

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

Authors: Woojin Kwon, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, Derek Ward-Thompson, James Di Francesco, Patrick M. Koch, Ray Furuya, Yasuo Doi, Valentin J. M. Le Gouellec, Jihye Hwang, A-Ran Lyo, Archana Soam, Xindi Tang, Thiem Hoang, Florian Kirchschlager, Chakali Eswaraiah, Lapo Fanciullo, Kyoung Hee Kim, Takashi Onaka, Vera Könyves, Ji-hyun Kang, Chang Won Lee, Motohide Tamura , et al. (127 additional authors not shown)

Abstract: 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 filament… ▽ More 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. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 18 pages, accepted for publication in ApJ

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The HASHTAG project: The First Submillimeter Images of the Andromeda Galaxy from the Ground

Authors: Matthew W. L. Smith, Stephen A. Eales, Thomas G. Williams, Bumhyun Lee, Zongnan Li, Pauline Barmby, Martin Bureau, Scott Chapman, Brian S. Cho, Aeree Chung, Eun Jung Chung, Hui-Hsuan Chung, Christopher J. R. Clark, David L. Clements, Timothy A. Davis, Ilse De Looze, David J. Eden, Gayathri Athikkat-Eknath, George P. Ford, Yu Gao, Walter Gear, Haley L. Gomez, Richard de Grijs, Jinhua He, Luis C. Ho , et al. (24 additional authors not shown)

Abstract: Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing… ▽ More Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

The JCMT BISTRO Survey: An 850/450$μ$m Polarization Study of NGC 2071IR in OrionB

Authors: A-Ran Lyo, Jongsoo Kim, Sarah Sadavoy, Doug Johnstone, David Berry, Kate Pattle, Woojin Kwon, Pierre Bastien, Takashi Onaka, James Di Francesco, Ji-Hyun Kang, Ray Furuya, Charles L. H. Hull, Motohide Tamura, Patrick M. Koch, Derek Ward-Thompson, Tetsuo Hasegawa, Thiem Hoang, Doris Arzoumanian, Chang Won Lee, Chin-Fei Lee, Do-Young Byun, Florian Kirchschlager, Yasuo Doi, Kee-Tae Kim , et al. (121 additional authors not shown)

Abstract: We present the results of simultaneous 450 $μ$m and 850 $μ$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 b… ▽ More We present the results of simultaneous 450 $μ$m and 850 $μ$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 $μ$G in the central $\sim$0.12 pc region from 850 $μ$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 $μ$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $μ$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 $μ$m and $0.36 \pm 0.04$ at 850 $μ$m, respectively. We think that the shallow slope at 850 $μ$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: 17 pages, 11 figures, 4 tables

Journal ref: The Astrophysical Journal, 2021, 918, 85

A Decade of SCUBA-2: A Comprehensive Guide to Calibrating 450 $μ$m and 850 $μ$m Continuum Data at the JCMT

Authors: Steve Mairs, Jessica T. Dempsey, Graham S. Bell, Harriet Parsons, Malcolm J. Currie, Per Friberg, Xue-Jian Jiang, Alexandra J. Tetarenko, Dan Bintley, Jamie Cookson, Shaoliang Li, Mark G. Rawlings, Jan Wouterloot, David Berry, Sarah Graves, Izumi Mizuno, Alexis Ann Acohido, Alyssa Clark, Jeff Cox, Miriam Fuchs, James Hoge, Johnathon Kemp, E'lisa Lee, Callie Matulonis, William Montgomerie , et al. (2 additional authors not shown)

Abstract: The Submillimetre Common User Bolometer Array 2 (SCUBA-2) is the James Clerk Maxwell Telescope's continuum imager, operating simultaneously at 450 and 850~$μ$m. SCUBA-2 was commissioned in 2009--2011 and since that time, regular observations of point-like standard sources have been performed whenever the instrument is in use. Expanding the calibrator observation sample by an order of magnitude com… ▽ More The Submillimetre Common User Bolometer Array 2 (SCUBA-2) is the James Clerk Maxwell Telescope's continuum imager, operating simultaneously at 450 and 850~$μ$m. SCUBA-2 was commissioned in 2009--2011 and since that time, regular observations of point-like standard sources have been performed whenever the instrument is in use. Expanding the calibrator observation sample by an order of magnitude compared to previous work, in this paper we derive updated opacity relations at each wavelength for a new atmospheric-extinction correction, analyze the Flux-Conversion Factors (FCFs) used to convert instrumental units to physical flux units as a function of date and observation time, present information on the beam profiles for each wavelength, and update secondary-calibrator source fluxes. Between 07:00 and 17:00 UTC, the portion of the night that is most stable to temperature gradients that cause dish deformation, the total-flux uncertainty and the peak-flux uncertainty measured at 450~$μ$m are found to be 14\% and 17\%, respectively. Measured at 850~$μ$m, the total-flux and peak-flux uncertainties are 6\%, and 7\%, respectively. The analysis presented in this work is applicable to all SCUBA-2 projects observed since 2011. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 31 pages, 15 figures, 8 tables. This is the Accepted Manuscript version of an article accepted for publication in the Astronomical Journal. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. This Accepted Manuscript is published under a CC BY licence

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

The JCMT BISTRO Survey: Magnetic Fields Associated with a Network of Filaments in NGC 1333

Authors: Yasuo Doi, Tetsuo Hasegawa, Ray S. Furuya, Simon Coudé, Charles L. H. Hull, Doris Arzoumanian, Pierre Bastien, Michael Chun-Yuan Chen, James di Francesco, Rachel Friesen, Martin Houde, Shu-ichiro Inutsuka, Steve Mairs, Masafumi Matsumura, Takashi Onaka, Sarah Sadavoy, Yoshito Shimajiri, Mehrnoosh Tahani, Kohji Tomisaka, Chakali Eswaraiah, Patrick M. Koch, Kate Pattle, Chang Won Lee, Motohide Tamura, David Berry , et al. (113 additional authors not shown)

Abstract: 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 x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary… ▽ More 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 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 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. △ Less

Submitted 21 July, 2020; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 33 pages, 21 figures, accepted for publication in ApJ

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

Abstract: Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabili… ▽ More Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 15 figures, accepted for publication in A&A

The MALATANG Survey: Dense Gas and Star Formation from High Transition HCN and HCO+ maps of NGC253

Authors: Xue-Jian Jiang, Thomas R. Greve, Yu Gao, Zhi-Yu Zhang, Qinghua Tan, Richard de Grijs, Luis C. Ho, Michal J. Michalowski, Malcolm J. Currie, Christine D. Wilson, Elias Brinks, Yiping Ao, Yinghe Zhao, Jinhua He, Nanase Harada, Chentao Yang, Qian Jiao, Aeree Chung, Bumhyun Lee, Matthew W. L. Smith, Daizhong Liu, Satoki Matsushita, Yong Shi, Masatoshi Imanishi, Mark G. Rawlings , et al. (4 additional authors not shown)

Abstract: To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission feat… ▽ More To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $Σ_\text{stellar}$, and SFE$_\text{dense}$ versus $Σ_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $Σ_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars. △ Less

Submitted 14 March, 2020; originally announced March 2020.

Comments: accepted to MNRAS

Multi-scale analysis of the Monoceros OB 1 star-forming region: II. Colliding filaments in the Monoceros OB1 molecular cloud

Authors: Julien Montillaud, Mika Juvela, Charlotte Vastel, J. H. He, Tie Liu, Isabelle Ristorcelli, David J. Eden, Sung-ju Kang, Kee-Tae Kim, Patrick M. Koch, Chang Won Lee, Mark G. Rawlings, Mika Saajasto, Patricio Sanhueza, Archana Soam, Sarolta Zahorecz, Dana Alina, Rebeka Bögner, David Cornu, Yasuo Doi, Johanna Malinen, Douglas Marshall, E. R. Micelotta, V. M. Pelkonen, L. V. Tóth , et al. (2 additional authors not shown)

Abstract: We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the… ▽ More We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: A&A, in press; IRAM and TRAO data cubes will be available at CDS

Journal ref: A&A 631, A3 (2019)

Multi-scale analysis of the Monoceros OB 1 star-forming region: I. The dense core population

Authors: Julien Montillaud, Mika Juvela, Charlotte Vastel, J. H. He, Tie Liu, Isabelle Ristorcelli, David J. Eden, Sung-ju Kang, Kee-Tae Kim, Patrick M. Koch, Chang Won Lee, Mark G. Rawlings, Mika Saajasto, Patricio Sanhueza, Archana Soam, Sarolta Zahorecz, Dana Alina, Rebeka Bögner, David Cornu, Yasuo Doi, Johanna Malinen, Douglas Marshall, E. R. Micelotta, V. M. Pelkonen, L. V. Tóth , et al. (2 additional authors not shown)

Abstract: Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures.… ▽ More Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Journal ref: A&A 631, L1 (2019)

The JCMT BISTRO Survey: The Magnetic Field of the Barnard 1 Star-Forming Region

Authors: Simon Coudé, Pierre Bastien, Martin Houde, Sarah Sadavoy, Rachel Friesen, James Di Francesco, Doug Johnstone, Steve Mairs, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Keping Qiu, Derek Ward-Thompson, David Berry, Michael Chun-Yuan Chen, Jason Fiege, Erica Franzmann, Jennifer Hatchell, Kevin Lacaille, Brenda C. Matthews, Gerald H. Moriarty-Schieven, Andy Pon, Philippe André, Doris Arzoumanian, Yusuke Aso , et al. (96 additional authors not shown)

Abstract: We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We th… ▽ More We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $μ$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $μ$G and a criticality criterion $λ_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component. △ Less

Submitted 23 April, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

Comments: 22 pages, 7 figures, accepted in the Astrophysical Journal

Near-Infrared Diffuse Interstellar Bands towards Her 36

Authors: M. G. Rawlings, A. J. Adamson, C. C. M. Marshall, P. J. Sarre

Abstract: Discovered almost a century ago, the Diffuse Interstellar Bands (DIBs) still lack convincing and comprehensive identification. Hundreds of DIBs have now been observed in the near-ultraviolet (NUV), visible and near-infrared (NIR). They are widely held to be molecular in origin, and modelling of their band profiles offers powerful constraints on molecular constants. Herschel 36, the illuminating st… ▽ More Discovered almost a century ago, the Diffuse Interstellar Bands (DIBs) still lack convincing and comprehensive identification. Hundreds of DIBs have now been observed in the near-ultraviolet (NUV), visible and near-infrared (NIR). They are widely held to be molecular in origin, and modelling of their band profiles offers powerful constraints on molecular constants. Herschel 36, the illuminating star of the Lagoon Nebula, has been shown to possess unusually broad and asymmetric DIB profiles in the visible, and is also bright enough for NIR observation. We present here high-resolution spectroscopic observations targeting the two best-known NIR DIBs at 11797.5 and 13175 A toward this object and a nearby comparison O-star, 9 Sgr, using the GNIRS instrument on Gemini North. We show a clear detection of the 13175 A DIB in both stars, and find (i) that it does not exhibit the unusual wing structure of some of the visual DIBs in Her 36 and (ii) that the depth of the band in the two objects is very similar, also contrary to the behaviour of the visual DIBs. We discuss the implications of these results for multiple DIB carrier candidates, and the location of their carriers along the observed lines of sight. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: 5 pages, 2 figures, Accepted by MNRAS

SCOPE: SCUBA-2 Continuum Observations of Pre-protostellar Evolution - Survey Description and Compact Source Catalogue

Authors: D. J. Eden, Tie Liu, Kee-Tae Kim, S. -Y. Liu, K. Tatematsu, J. Di Francesco, K. Wang, Y. Wu, M. A. Thompson, G. A. Fuller, Di Li, I. Ristorcelli, Sung-ju Kang, N. Hirano, D. Johnstone, Y. Lin, J. H. He, P. M. Koch, Patricio Sanhueza, S. -L. Qin, Q. Zhang, P. F. Goldsmith, N. J. Evans II, J. Yuan, C. -P. Zhang , et al. (136 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 16 pages, 13 figures, 1 table. Accepted for publication in MNRAS

The JCMT BISTRO Survey: The Magnetic Field In The Starless Core $ρ$ Ophiuchus C

Authors: Junhao Liu, Keping Qiu, David Berry, James Di Francesco, Pierre Bastien, Patrick M. Koch, Ray S. Furuya, Kee-Tae Kim, Simon Coudé, Chang Won Lee, Archana Soam, Chakali Eswaraiah, Di Li, Jihye Hwang, A-Ran Lyo, Kate Pattle, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Derek Ward-Thompson, Tao-Chung Ching, Zhiwei Chen, Qilao Gu, Dalei Li, Hua-bai Li , et al. (106 additional authors not shown)

Abstract: We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic… ▽ More We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $μ$G, 136 $\pm$ 69 $μ$G, and 213 $\pm$ 115 $μ$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: 22 pages, 12 figures

JCMT BISTRO survey: Magnetic Fields within the Hub-Filament Structure in IC 5146

Authors: Jia-Wei Wang, Shih-Ping Lai, Chakali Eswaraiah, Kate Pattle, James Di Francesco, Doug Johnstone, Patrick M. Koch, Tie Liu, Motohide Tamura, Ray S. Furuya, Takashi Onaka, Derek Ward-Thompson, Archana Soam, Kee-Tae Kim, Chang Won Lee, Chin-Fei Lee, Steve Mairs, Doris Arzoumanian, Gwanjeong Kim, Thiem Hoang, Jihye Hwang, Sheng-Yuan Liu, David Berry, Pierre Bastien, Tetsuo Hasegawa , et al. (108 additional authors not shown)

Abstract: We present the 850 $μ$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-scal… ▽ More We present the 850 $μ$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é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. △ Less

Submitted 27 March, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

Comments: 24 pages, 15 figures, accepted for publication in ApJ

Compressed magnetic field in the magnetically-regulated global collapsing clump of G9.62+0.19

Authors: Tie Liu, Kee-Tae Kim, Sheng-Yuan Liu, Mika Juvela, Qizhou Zhang, Yuefang Wu, Pak Shing Li, Harriet Parsons, Archana Soam, Paul F. Goldsmith, Yu-Nung Su, Kenichi Tatematsu, Sheng-Li Qin, Guido Garay, Tomoya Hirota, Jan Wouterloot, Huei-Ru Chen, Neal J. Evans II, Sarah Graves, Sung-ju Kang, Di Li, Diego Mardones, Mark G. Rawlings, Zhiyuan Ren, Ke Wang

Abstract: How stellar feedback from high-mass stars (e.g., H{\sc ii} regions) influences the surrounding interstellar medium and regulates new star formation is still unclear. To address this question, we observed the G9.62+0.19 complex in 850 $μ$m continuum with the JCMT/POL-2 polarimeter. An ordered magnetic field has been discovered in its youngest clump, the G9.62 clump. The magnetic field strength is d… ▽ More How stellar feedback from high-mass stars (e.g., H{\sc ii} regions) influences the surrounding interstellar medium and regulates new star formation is still unclear. To address this question, we observed the G9.62+0.19 complex in 850 $μ$m continuum with the JCMT/POL-2 polarimeter. An ordered magnetic field has been discovered in its youngest clump, the G9.62 clump. The magnetic field strength is determined to be $\sim$1 mG. Magnetic field plays a larger role than turbulence in supporting the clump. However, the G9.62 clump is still unstable against gravitational collapse even if thermal, turbulent, and magnetic field support are taken into account all together. The magnetic field segments in the outskirts of the G9.62 clump seem to point toward the clump center, resembling a dragged-in morphology, indicating that the clump is likely undergoing magnetically-regulated global collapse. However, The magnetic field in its central region is aligned with the shells of the photodissociation regions (PDRs) and is approximately parallel to the ionization (or shock) front, indicating that the magnetic field therein is likely compressed by the expanding H{\sc ii} regions that formed in the same complex. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: 8 pages, 3 figures, accepted to be published on ApJ letter

Dust spectrum and polarisation at 850 um in the massive IRDC G035.39-00.33

Authors: Mika Juvela, Vincent Guillet, Tie Liu, Isabelle Ristorcelli, Veli-Matti Pelkonen, Dana Alina, Leonardo Bronfman, David J. Eden, Kee Tae Kim, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Johanna Malinen, Elisabetta Micelotta, Julien Montillaud, Mark G. Rawlings, Patricio Sanhueza, Archana Soam, Alessio Traficante, Nathalie Ysard, Chuan-Peng Zhang

Abstract: The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation… ▽ More The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts. △ Less

Submitted 6 September, 2018; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 620, A26 (2018)

Planck Cold Clumps in the $λ$ Orionis complex. II. Environmental effects on core formation

Authors: Hee-Weon Yi, Jeong-Eun Lee, Tie Liu, Kee-Tae Kim, Minho Choi, David Eden, Neal J. Evans II, James Di Francesco, Gary Fuller, N. Hirano, Mika Juvela, Sung-ju Kang, Gwanjeong Kim, Patrick M. Koch, Chang Won Lee, Di Li, H. -Y. B. Liu, Hong-Li Liu, Sheng-Yuan Liu, Mark G. Rawlings, I. Ristorcelli, Patrico Sanhueza, Archana Soam, Ken'ichi Tatematsu, Mark Thompson , et al. (5 additional authors not shown)

Abstract: Based on the 850 $μ$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $λ$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $λ$ Orionis cloud has a different environment as a consequence of the interaction wi… ▽ More Based on the 850 $μ$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $λ$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $λ$ Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant Hii region. PGCCs in the $λ$ Orionis cloud have higher dust temperatures ($Td=16.13\pm0.15$ K) and lower values of dust emissivity spectral index ($ β=1.65\pm0.02$) than PGCCs in the Orion A (Td=13.79$\pm 0.21$K, $β=2.07\pm0.03$) and Orion B ($Td=13.82\pm0.19$K, $β=1.96\pm0.02$) clouds. We find 119 sub-structures within the 40 detected PGCCs and identify them as cores. Of total 119 cores, 15 cores are discovered in the $λ$ Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the $λ$ Orionis cloud show much lower mean values of size R=0.08 pc, column density N(H2)=$(9.5\pm1.2) \times 10^{22}$ cm$^{-2}$, number density n(H2)=$(2.9 \pm 0.4)\times10^{5}$ cm$^{-3}$, and mass $M_{core}$=$1.0\pm0.3$ M$_{\odot}$ compared to the cores in the Orion A (R=0.11pc, $N(H2)=(2.3\pm0.3) \times 10^{23}$ cm$^{-2}$, n(H2)=$(3.8\pm0.5) \times 10^{5}$cm$^{-3}$, and $M_{core}$=$2.4 \pm 0.3$ M$_{\odot}$) and Orion B (R=0.16pc, N(H2)=$(3.8 \pm 0.4) \times 10^{23}$cm$^{-2}$, n(H2)=$(15.6\pm1.8)\times10^{5}$ cm$^{-3}$, and $M_{core}$= $2.7\pm0.3$ M$_{\odot}$) clouds. These core properties in the $λ$ Orionis cloud can be attributed to the photodissociation and external heating by the nearby Hii region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation. △ Less

Submitted 15 May, 2018; originally announced May 2018.

The MALATANG Survey: the L_gas-L_IR correlation on sub-kiloparsec scale in six nearby star-forming galaxies as traced by HCN J=4-3 and HCO^+ J=4-3

Authors: Qing-Hua Tan, Yu Gao, Zhi-Yu Zhang, Thomas R. Greve, Xue-Jian Jiang, Christine D. Wilson, Chen-Tao Yang, Ashley Bemis, Aeree Chung, Satoki Matsushita, Yong Shi, Yi-Ping Ao, Elias Brinks, Malcolm J. Currie, Timothy A. Davis, Richard de Grijs, Luis C. Ho, Masatoshi Imanishi, Kotaro Kohno, Bumhyun Lee, Harriet Parsons, Mark G. Rawlings, Dimitra Rigopoulou, Erik Rosolowsky, Joanna Bulger , et al. (21 additional authors not shown)

Abstract: We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where… ▽ More We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 μm/100 μm flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: 29 pages, 12 figures, accepted for publication in ApJ

A holistic perspective on the dynamics of G035.39-00.33: the interplay between gas and magnetic fields

Authors: Tie Liu, Pak Shing Li, Mika Juvela, Kee-Tae Kim, Neal J. Evans II, James Di Francesco, Sheng-Yuan Liu, Jinghua Yuan, Ken'ichi Tatematsu, Qizhou Zhang, Derek Ward-Thompson, Gary Fuller, Paul F. Goldsmith, P. M. Koch, Patricio Sanhueza, I. Ristorcelli, Sung-ju Kang, Huei-Ru Chen, N. Hirano, Yuefang Wu, Vlas Sokolov, Chang Won Lee, Glenn J. White, Ke Wang, David Eden , et al. (32 additional authors not shown)

Abstract: Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at… ▽ More Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $μ$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament. △ Less

Submitted 18 June, 2018; v1 submitted 26 March, 2018; originally announced March 2018.

Comments: Published on ApJ, 27 pages

Journal ref: Liu, T., Li, P. S., Juvela M., et al., 2018, ApJ, 859 ,151

The properties of Planck Galactic cold clumps in the L1495 dark cloud

Authors: Mengyao Tang, Tie Liu, Sheng-Li Qin, Kee-Tae Kim, Yuefang Wu, Ken'ichi Tatematsu, Jinghua Yuan, Ke Wang, Harriet Parsons, Patrick M. Koch, Patricio Sanhueza, D. Ward-Thompson, L. Viktor Tóth, Archana Soam, Chang Won Lee, David Eden, James Di Francesco, Jonathan Rawlings, Mark G. Rawlings, Julien Montillaud, Chuan-Peng Zhang, M. R. Cunningham

Abstract: Planck Galactic Cold Clumps (PGCCs) possibly represent the early stages of star formation. To understand better the properties of PGCCs, we studied 16 PGCCs in the L1495 cloud with molecular lines and continuum data from Herschel, JCMT/SCUBA-2 and the PMO 13.7 m telescope. Thirty dense cores were identified in 16 PGCCs from 2-D Gaussian fitting. The dense cores have dust temperatures of… ▽ More Planck Galactic Cold Clumps (PGCCs) possibly represent the early stages of star formation. To understand better the properties of PGCCs, we studied 16 PGCCs in the L1495 cloud with molecular lines and continuum data from Herschel, JCMT/SCUBA-2 and the PMO 13.7 m telescope. Thirty dense cores were identified in 16 PGCCs from 2-D Gaussian fitting. The dense cores have dust temperatures of $T_{\rm d}$ = 11-14 K, and H$_{2}$ column densities of $N_{\rm H_{2}}$ = 0.36-2.5$\times10^{22}$ cm$^{-2}$. We found that not all PGCCs contain prestellar objects. In general, the dense cores in PGCCs are usually at their earliest evolutionary stages. All the dense cores have non-thermal velocity dispersions larger than the thermal velocity dispersions from molecular line data, suggesting that the dense cores may be turbulence-dominated. We have calculated the virial parameter $α$ and found that 14 of the dense cores have $α$ $<$ 2, while 16 of the dense cores have $α$ $>$ 2. This suggests that some of the dense cores are not bound in the absence of external pressure and magnetic fields. The column density profiles of dense cores were fitted. The sizes of the flat regions and core radii decrease with the evolution of dense cores. CO depletion was found to occur in all the dense cores, but is more significant in prestellar core candidates than in protostellar or starless cores. The protostellar cores inside the PGCCs are still at a very early evolutionary stage, sharing similar physical and chemical properties with the prestellar core candidates. △ Less

Submitted 14 February, 2018; originally announced February 2018.

Herschel and SCUBA-2 observations of dust emission in a sample of Planck cold clumps

Authors: Mika Juvela, Jinhua He, Katherine Pattle, Tie Liu, George Bendo, David J. Eden, Orsolya Feher, Michel Fich, Gary Fuller, Naomi Hirano, Kee-Tae Kim, Di Li, Sheng-Yuan Liu, Johanna Malinen, Douglas J. Marshall, Deborah Paradis, Harriet Parsons, Veli-Matti Pelkonen, Mark G. Rawlings, Isabelle Ristorcelli, Manash R. Samal, Ken'ichi Tatematsu, Mark Thompson, Alessio Traficante, Ke Wang , et al. (4 additional authors not shown)

Abstract: Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs… ▽ More Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths. △ Less

Submitted 27 December, 2017; v1 submitted 26 November, 2017; originally announced November 2017.

Comments: 37 pages, accepted to A&A, final version with language corrections

Journal ref: A&A 612, A71 (2018)

The ALMA detection of CO rotational line emission in AGB stars in the Large Magellanic Cloud

Authors: M. A. T. Groenewegen, W. H. T. Vlemmings, P. Marigo, G. C. Sloan, L. Decin, M. W. Feast, S. R. Goldman, K. Justtanont, F. Kerschbaum, M. Matsuura, I. McDonald, H. Olofsson, R. Sahai, J. Th. van Loon, P. R. Wood, A. A. Zijlstra, J. Bernard-Salas, M. L. Boyer, L. Guzman-Ramirez, O. C. Jones, E. Lagadec, M. Meixner, M. G. Rawlings, S. Srinivasan

Abstract: Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB… ▽ More Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models. △ Less

Submitted 30 September, 2016; originally announced September 2016.

Comments: Accepted for publication by A&A

Journal ref: A&A 596, A50 (2016)

An Overview of the 2014 ALMA Long Baseline Campaign

Authors: ALMA Partnership, E. B. Fomalont, C. Vlahakis, S. Corder, A. Remijan, D. Barkats, R. Lucas, T. R. Hunter, C. L. Brogan, Y. Asaki, S. Matsushita, W. R. F. Dent, R. E. Hills, N. Phillips, A. M. S. Richards, P. Cox, R. Amestica, D. Broguiere, W. Cotton, A. S. Hales, R. Hiriart, A. Hirota, J. A. Hodge, C. M. V. Impellizzeri, J. Kern , et al. (224 additional authors not shown)

Abstract: A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and… ▽ More A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy. △ Less

Submitted 24 April, 2015; v1 submitted 19 April, 2015; originally announced April 2015.

Comments: 11 pages, 7 figures, 2 tables; accepted for publication in the Astrophysical Journal Letters; this version with small changes to affiliations

Journal ref: The Astrophysical Journal Letters, volume 808, issue 1, article id. L1, 11 pp. (2015)

A High-Resolution Study of Near-Infrared Diffuse Interstellar Bands

Authors: M. G. Rawlings, A. J. Adamson, T. H. Kerr

Abstract: We present high-resolution echelle spectroscopic observations of the two near- infrared (NIR) Diffuse Interstellar Bands (DIBs) at 13175 A and 11797.5 A. The DIBs have been observed in a number of diffuse interstellar medium sightlines that exhibit a wide range of visual extinctions. Band profiles are similar to those seen in narrow DIBs, clearly asymmetric and can be closely fitted in most cases… ▽ More We present high-resolution echelle spectroscopic observations of the two near- infrared (NIR) Diffuse Interstellar Bands (DIBs) at 13175 A and 11797.5 A. The DIBs have been observed in a number of diffuse interstellar medium sightlines that exhibit a wide range of visual extinctions. Band profiles are similar to those seen in narrow DIBs, clearly asymmetric and can be closely fitted in most cases using two simple Gaussian components. Gaussian fits were generally found to be more successful than fits based on a multiple-cloud model using a template DIB profile. For a sample of 9 objects in which both bands are observed, the strength of both NIR DIBs generally increases with A(V), and we report a correlation between the two observed bands over a large A(V) range and widely-separated lines of sight. The strength of the two bands is also compared against those of two visual DIBs and the diffuse ISM aliphatic dust absorption feature at 3.4 microns previously detected in the same sightlines. We find that the NIR DIBs do not exhibit notable (anti)correlations with either. Implications of these observations on possible DIB carrier species are discussed. △ Less

Submitted 25 September, 2014; originally announced September 2014.

Comments: 35 pages; contains 8 figures (including a 3-page Figure 1) and 3 tables. Accepted for publication in the Astrophysical Journal

Journal ref: 2014, ApJ, 796, 58

Mapping of interstellar clouds with infrared light scattered from dust: TMC-1N

Authors: J. Malinen, M. Juvela, V. -M. Pelkonen, M. G. Rawlings

Abstract: Mapping of near-infrared (NIR) scattered light is a recent method for the study of interstellar clouds, complementing other, more commonly used methods, like dust emission and extinction. Our goal is to study the usability of this method on larger scale, and compare the properties of a filament using NIR scattering and other methods. We also study the radiation field and differences in grain emiss… ▽ More Mapping of near-infrared (NIR) scattered light is a recent method for the study of interstellar clouds, complementing other, more commonly used methods, like dust emission and extinction. Our goal is to study the usability of this method on larger scale, and compare the properties of a filament using NIR scattering and other methods. We also study the radiation field and differences in grain emissivity between diffuse and dense areas. We have used scattered J, H, and K band surface brightness WFCAM-observations to map filament TMC-1N in Taurus, covering an area of 1dx1d corresponding to ~(2.44 pc)^2. We have converted the data into optical depth and compared the results with NIR extinction and Herschel observations of submm dust emission. We see the filament in scattered light in all three NIR bands. We note that our WFCAM observations in TMC-1N show notably lower intensity than previous results in Corona Australis using the same method. We show that 3D radiative transfer simulations predict similar scattered surface brightness levels as seen in the observations. However, changing the assumptions about the background can change the results of simulations notably. We derive emissivity by using optical depth in the J band as an independent tracer of column density. We obtain opacity sigma(250um) values 1.7-2.4x10^-25 cm^2/H, depending on assumptions of the extinction curve, which can change the results by over 40%. These values are twice as high as obtained for diffuse areas, at the lower limit of earlier results for denser areas. We show that NIR scattering can be a valuable tool in making high resolution maps. We conclude, however, that NIR scattering observations can be complicated, as the data can show relatively low-level artefacts. This suggests caution when planning and interpreting the observations. △ Less

Submitted 9 April, 2014; originally announced April 2014.

Comments: abstract shortened and figures reduced for astroph

Journal ref: A&A 558, A44 (2013)

Silicon isotopic abundance toward evolved stars and its application for presolar grains

Authors: T. -C. Peng, E. M. L. Humphreys, L. Testi, A. Baudry, M. Wittkowski, M. G. Rawlings, I. de Gregorio-Monsalvo, W. Vlemmings, L. -A. Nyman, M. D. Gray, C. de Breuck

Abstract: Galactic chemical evolution (GCE) is important for understanding the composition of the present-day interstellar medium (ISM) and of our solar system. In this paper, we aim to track the GCE by using the 29Si/30Si ratios in evolved stars and tentatively relate this to presolar grain composition. We used the APEX telescope to detect thermal SiO isotopologue emission toward four oxygen-rich M-type st… ▽ More Galactic chemical evolution (GCE) is important for understanding the composition of the present-day interstellar medium (ISM) and of our solar system. In this paper, we aim to track the GCE by using the 29Si/30Si ratios in evolved stars and tentatively relate this to presolar grain composition. We used the APEX telescope to detect thermal SiO isotopologue emission toward four oxygen-rich M-type stars. Together with the data retrieved from the Herschel science archive and from the literature, we were able to obtain the 29Si/30Si ratios for a total of 15 evolved stars inferred from their optically thin 29SiO and 30SiO emission. These stars cover a range of masses and ages, and because they do not significantly alter 29Si/30Si during their lifetimes, they provide excellent probes of the ISM metallicity (or 29Si/30Si ratio) as a function of time. The 29Si/30Si ratios inferred from the thermal SiO emission tend to be lower toward low-mass oxygen-rich stars (e.g., down to about unity for W Hya), and close to an interstellar or solar value of 1.5 for the higher-mass carbon star IRC+10216 and two red supergiants. There is a tentative correlation between the 29Si/30Si ratios and the mass-loss rates of evolved stars, where we take the mass-loss rate as a proxy for the initial stellar mass or current stellar age. This is consistent with the different abundance ratios found in presolar grains. We found that older objects (up to possibly 10 Gyr old) in our sample trace a previous, lower 29Si/30Si value of about 1. Material with this isotopic ratio is present in two subclasses of presolar grains, providing independent evidence of the lower ratio. Therefore, the 29Si/30Si ratio derived from the SiO emission of evolved stars is a useful diagnostic tool for the study of the GCE and presolar grains. △ Less

Submitted 14 November, 2013; originally announced November 2013.

Comments: 7 pages, 4 figures

Journal ref: A&A, 2013, 559, L8

ALMA imaging of the CO snowline of the HD 163296 disk with DCO+

Authors: G. S. Mathews, P. D. Klaassen, A. Juhasz, D. Harsono, E. Chapillon, E. F. van Dishoeck, D. Espada, I. de Gregorio-Monsalvo, A. Hales, M. R. Hogerheijde, J. C. Mottram, M. G. Rawlings, S. Takahashi, L. Testi

Abstract: The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 1… ▽ More The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 19 K, making it a useful probe of the cold disk midplane. We compare ALMA line observations of HD 163296 to a grid of models. We vary the upper- and lower-limit temperatures of the region in which DCO+ is present as well as the abundance of DCO+ in order to fit channel maps of the DCO+ J=5-4 line. To determine the abundance enhancement compared to the general interstellar medium, we carry out similar fitting to HCO+ J=4-3 and H13CO+ J=4-3 observations. ALMA images show centrally peaked extended emission from HCO+ and H13CO+. DCO+ emission lies in a resolved ring from ~110 to 160 AU. The outer radius approximately corresponds to the size of the CO snowline as measured by previous lower resolution observations of CO lines in this disk. The ALMA DCO+ data now resolve and image the CO snowline directly. In the best fitting models, HCO+ exists in a region extending from the 19 K isotherm to the photodissociation layer with an abundance of 3x10^-10 relative to H2. DCO+ exists within the 19-21 K region of the disk with an abundance ratio [DCO+] / [HCO+] = 0.3. This represents a factor of 10^4 enhancement of the DCO+ abundance within this narrow region of the HD 163296 disk. Such a high enhancement has only previously been seen in prestellar cores. The inferred abundances provide a lower limit to the ionization fraction in the midplane of the cold outer disk (approximately greater than 4x10^-10), and suggest the utility of DCO+ as a tracer of its parent molecule H2D+. Abridged △ Less

Submitted 12 July, 2013; originally announced July 2013.

Comments: 10 pages, accepted to A&A

Unveiling the gas and dust disk structure in HD 163296 using ALMA observations

Authors: I. de Gregorio-Monsalvo, F. Ménard, W. Dent, C. Pinte, C. López, P. Klaassen, A. Hales, P. Cortés, M. G. Rawlings, K. Tachihara, L. Testi, S. Takahashi, E. Chapillon, G. Mathews, A. Juhasz, E. Akiyama, A. E. Higuchi, M. Saito, L. -Å. Nyman, N. Phillips, J. Rodń, S. Corder, T. Van Kempen

Abstract: Aims: The aim of this work is to study the structure of the protoplanetary disk surrounding the Herbig Ae star HD 163296. Methods: We have used high-resolution and high-sensitivity ALMA observations of the CO(3-2) emission line and the continuum at 850 microns, as well as the 3- dimensional radiative transfer code MCFOST to model the data presented in this work. Results: The CO(3-2) emission unvei… ▽ More Aims: The aim of this work is to study the structure of the protoplanetary disk surrounding the Herbig Ae star HD 163296. Methods: We have used high-resolution and high-sensitivity ALMA observations of the CO(3-2) emission line and the continuum at 850 microns, as well as the 3- dimensional radiative transfer code MCFOST to model the data presented in this work. Results: The CO(3-2) emission unveils for the first time at sub-millimeter frequencies the vertical structure details of a gaseous disk in Keplerian rotation, showing the back- and the front-side of a flared disk. Continuum emission at 850 microns reveals a compact dust disk with a 240 AU outer radius and a surface brightness profile that shows a very steep decline at radius larger than 125 AU. The gaseous disk is more than two times larger than the dust disk, with a similar critical radius but with a shallower radial profile. Radiative transfer models of the continuum data confirms the need for a sharp outer edge to the dust disk. The models for the CO(3-2) channel map require the disk to be slightly more geometrically thick than previous models suggested, and that the temperature at which CO gas becomes depleted (frozen-out) from the outer regions of the disk midplane is T < 20 K, in agreement with previous studies. △ Less

Submitted 4 July, 2013; originally announced July 2013.

Comments: Accepted for publication in Astronomy and Astrophysics

Observations of 6 - 200 μm emission of the Ophiuchus cloud LDN 1688

Authors: M. G. Rawlings, M. Juvela, K. Lehtinen, K. Mattila, D. Lemke

Abstract: We examine two positions, ON1 and ON2, within the Ophiuchus cloud LDN 1688 using observations made with the ISOPHOT instrument aboard the ISO satellite. The data include mid-IR spectra (~6-12μm) and several photometric bands up to 200μm. The data probe the emission from molecular PAH-type species, transiently-heated Very Small Grains (VSGs), and large classical dust grains. We compare the observat… ▽ More We examine two positions, ON1 and ON2, within the Ophiuchus cloud LDN 1688 using observations made with the ISOPHOT instrument aboard the ISO satellite. The data include mid-IR spectra (~6-12μm) and several photometric bands up to 200μm. The data probe the emission from molecular PAH-type species, transiently-heated Very Small Grains (VSGs), and large classical dust grains. We compare the observations to earlier studies, especially those carried out towards an isolated translucent cloud in Chamaeleon (Paper I). The spectra towards the two LDN 1688 positions are very similar to each other, in spite of position ON1 having a larger column density and probably being subjected to a stronger radiation field. The ratios of the mid-IR features are similar to those found in other diffuse and translucent clouds. Compared to paper I, the 7.7/11.3μm band ratios are lower, ~2.0, at both LDN 1688 positions. A continuum is detected in the ~10μm region. This is stronger towards the position ON1 but still lower than on any of the sightlines in Paper I. The far-infrared opacities are higher than for diffuse medium. The value of the position ON2, τ200/N(H) = 3.9 x 10^{-25} cm^2/H, is twice the value found for ON1. The radiation field of LDN 1688 is dominated by the two embedded B type double stars, ρ Oph AB and HD 147889, with an additional contribution from the Upper Sco OB association. The strong heating is reflected in the high colour temperature, ~24 K, of the large grain emission. Radiative transfer modelling confirms a high level of the radiation field and points to an increased abundance of PAH grains. However, when the hardening of the radiation field caused by the local B-stars is taken into account, the observations can be fitted with almost no change to the standard dust models. However, all the examined models underestimate the level of the mid-IR continuum. △ Less

Submitted 17 October, 2012; originally announced October 2012.

Comments: 12 pages, 7 figures, 3 tables. Accepted for publication in MNRAS

Journal ref: 2013, MNRAS, 428, 2617

Profiling filaments: comparing near-infrared extinction and submillimetre data in TMC-1

Authors: J. Malinen, M. Juvela, M. G. Rawlings, D. Ward-Thompson, P. Palmeirim, Ph. Andre

Abstract: Interstellar filaments are an important part of star formation. To understand the structure of filaments, cross-section profiles are often fitted with Plummer profiles. This profiling is often done with submm studies, such as Herschel. It would be convenient if filament properties could also be studied using groundbased NIR data. We compare the filament profiles obtained by NIR extinction and subm… ▽ More Interstellar filaments are an important part of star formation. To understand the structure of filaments, cross-section profiles are often fitted with Plummer profiles. This profiling is often done with submm studies, such as Herschel. It would be convenient if filament properties could also be studied using groundbased NIR data. We compare the filament profiles obtained by NIR extinction and submm observations to find out if reliable profiles can be derived using NIR data. We use J-, H-, and K-band data of a filament north of TMC-1 to derive an extinction map from colour excesses of background stars. We compare the Plummer profiles obtained from extinction maps with Herschel dust emission maps. We present 2 methods to estimate profiles from NIR: Plummer profile fits to median Av of stars or directly to the Av of individual stars. We compare the methods by simulations. In simulations extinction maps and the new methods give correct results to within ~10-20 for modest densities. Direct fit to data on individual stars gives more accurate results than extinction map, and can work in higher density. In profile fits to real observations, values of Plummer parameters are generally similar to within a factor of ~2. Although parameter values can vary significantly, estimates of filament mass usually remain accurate to within some tens of per cent. Our results for TMC-1 are in agreement with earlier results. High resolution NIR data give more details, but 2MASS data can be used to estimate profiles. NIR extinction can be used as an alternative to submm observations to profile filaments. Direct fits of stars can also be a valuable tool. Plummer profile parameters are not always well constrained, and caution should be taken when making fits. In the evaluation of Plummer parameters, one can use the independence of dust emission and NIR data and the difference in the shapes of the confidence regions. △ Less

Submitted 21 June, 2012; originally announced June 2012.

Comments: accepted to Astronomy & Astrophysics; abstract has been shortened for astroph

[CII] line emission in massive star-forming galaxies at z=4.7

Authors: J. Wagg, T. Wiklind, C. L. Carilli, D. Espada, A. Peck, D. Riechers, F. Walter, A. Wootten, M. Aravena, D. Barkats, J. R. Cortes, R. Hills, J. Hodge, C. M V. Impellizzeri, D. Iono, A. Leroy, S. Martin, M. G. Rawlings, R. Maiolino, R. G. McMahon, K. S. Scott, E. Villard, C. Vlahakis

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [CII] 157.7micron fine structure line and thermal dust continuum emission from a pair of gas-rich galaxies at z=4.7, BR1202-0725. This system consists of a luminous quasar host galaxy and a bright submm galaxy (SMG), while a fainter star-forming galaxy is also spatially coincident within a 4" (25 kpc) region. All th… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [CII] 157.7micron fine structure line and thermal dust continuum emission from a pair of gas-rich galaxies at z=4.7, BR1202-0725. This system consists of a luminous quasar host galaxy and a bright submm galaxy (SMG), while a fainter star-forming galaxy is also spatially coincident within a 4" (25 kpc) region. All three galaxies are detected in the submm continuum, indicating FIR luminosities in excess of 10^13 Lsun for the two most luminous objects. The SMG and the quasar host galaxy are both detected in [CII] line emission with luminosities, L([CII]) = (10.0 +/- 1.5)x10^9 Lsun and L([CII]) = (6.5+/-1.0)x10^9 Lsun, respectively. We estimate a luminosity ratio, L([CII])/L(FIR) = (8.3+/-1.2)x10^-4 for the starburst SMG to the North, and L([CII])/L(FIR) = (2.5+/-0.4)x10^-4 for the quasar host galaxy, in agreement with previous high-redshift studies that suggest lower [CII]-to-FIR luminosity ratios in quasars than in starburst galaxies. The third fainter object with a flux density, S(340GHz) = 1.9+/-0.3 mJy, is coincident with a Ly-Alpha emitter and is detected in HST ACS F775W and F814W images but has no clear counterpart in the H-band. Even if this third companion does not lie at a similar redshift to BR1202-0725, the quasar and the SMG represent an overdensity of massive, infrared luminous star-forming galaxies within 1.3 Gyr of the Big Bang. △ Less

Submitted 18 May, 2012; v1 submitted 15 May, 2012; originally announced May 2012.

Comments: 14 pages, accepted for publication in ApJ Letters

The nature of the infrared counterpart of IGR J19140+0951

Authors: D. C. Hannikainen, M. G. Rawlings, P. Muhli, O. Vilhu, J. Schultz, J. Rodriguez

Abstract: The INTEGRAL observatory has been (re-)discovering new X-ray sources since the beginning of nominal operations in early 2003. These sources include X-ray binaries, Active Galactic Nuclei, cataclysmic variables, etc. Amongst the X-ray binaries, the true nature of many of these sources has remained largely elusive, though they seem to make up a population of highly absorbed high-mass X-ray binarie… ▽ More The INTEGRAL observatory has been (re-)discovering new X-ray sources since the beginning of nominal operations in early 2003. These sources include X-ray binaries, Active Galactic Nuclei, cataclysmic variables, etc. Amongst the X-ray binaries, the true nature of many of these sources has remained largely elusive, though they seem to make up a population of highly absorbed high-mass X-ray binaries. One of these new sources, IGR J19140+0951, was serendipitously discovered on 2003 Mar 6 during an observation of the galactic microquasar GRS 1915+105. We observed IGR J19140+0951 with UKIRT in order to identify the infrared counterpart. Here we present the H- and K-band spectra. We determined that the companion is a B0.5-type bright supergiant in a wind-fed system, at a distance $\la$ 5 kpc. △ Less

Submitted 8 June, 2007; originally announced June 2007.

Comments: 5 pages, 2 figures, accepted for publication in MNRAS

The UKIRT Infrared Deep Sky Survey Second Data Release

Authors: S. J. Warren, N. J. G. Cross, S. Dye, N. C. Hambly, O. Almaini, A. C. Edge, P. C. Hewett, S. T. Hodgkin, M. J. Irwin, R. F. Jameson, A. Lawrence, P. W. Lucas, D. J. Mortlock, A. J. Adamson, J. Bryant, R. S. Collins, C. J. Davis, J. P. Emerson, D. W. Evans, E. A. Gonzales-Solares, P. Hirst, T. H. Kerr, J. R. Lewis, R. G. Mann, M. G. Rawlings , et al. (4 additional authors not shown)

Abstract: The UKIRT Infrared Deep Sky Survey (UKIDSS) is a set of five large near-infrared surveys, covering a complementary range of areas, depths, and Galactic latitudes. The UKIDSS Second Data Release (DR2) includes the First Data Release (DR1), with minor improvements, plus new data for the LAS, GPS, GCS, and DXS, from observations made over 2006 May through July (when the UDS was unobservable). DR2 w… ▽ More The UKIRT Infrared Deep Sky Survey (UKIDSS) is a set of five large near-infrared surveys, covering a complementary range of areas, depths, and Galactic latitudes. The UKIDSS Second Data Release (DR2) includes the First Data Release (DR1), with minor improvements, plus new data for the LAS, GPS, GCS, and DXS, from observations made over 2006 May through July (when the UDS was unobservable). DR2 was staged in two parts. The first part excluded the GPS, and took place on 2007 March 1. The GPS was released on 2007 April 12. DR2 includes 282 sq. degs of multicolour data to (Vega) K=18, complete in the full YJHK set for the LAS, 57 sq. degs in the ZYJHK set for the GCS, and 236 sq. degs in the JHK set for the GPS. DR2 includes nearly 7 sq. degs of deep JK data (DXS, UDS) to an average depth K=21. In addition the release includes a comparable quantity of data where coverage of the filter set for any survey is incomplete. We document changes that have occurred since DR1 to the pipeline, calibration, and archive procedures. The two most noteworthy changes are presentation of the data in a single database (compared to two previously), and provision of additional error flags for detected sources, flagging potentially spurious artifacts, corrupted data and suspected cross-talk sources. We summarise the contents of each of the surveys in terms of filters, areas, and depths. △ Less

Submitted 26 July, 2007; v1 submitted 1 March, 2007; originally announced March 2007.

Comments: submitted to astro-ph only (not a journal publication) revised in v2 to include GPS details

The UKIRT Infrared Deep Sky Survey First Data Release

Authors: S. J. Warren, N. C. Hambly, S. Dye, O. Almaini, N. J. G. Cross, A. C. Edge, S. Foucaud, P. C. Hewett, S. T. Hodgkin, M. J. Irwin, R. F. Jameson, A. Lawrence, P. W. Lucas, A. J. Adamson, R. M. Bandyopadhyay, J. Bryant, R. S. Collins, C. J. Davis, J. S. Dunlop, J. P. Emerson, D. W. Evans, E. A. Gonzales-Solares, P. Hirst, M. J. Jarvis, T. R. Kendall , et al. (15 additional authors not shown)

Abstract: The First Data Release (DR1) of the UKIRT Infrared Deep Sky Survey (UKIDSS) took place on 2006 July 21. UKIDSS is a set of five large near-infrared surveys, covering a complementary range of areas, depths, and Galactic latitudes. DR1 is the first large release of survey-quality data from UKIDSS and includes 320 sq degs of multicolour data to (Vega) K=18, complete (depending on the survey) in thr… ▽ More The First Data Release (DR1) of the UKIRT Infrared Deep Sky Survey (UKIDSS) took place on 2006 July 21. UKIDSS is a set of five large near-infrared surveys, covering a complementary range of areas, depths, and Galactic latitudes. DR1 is the first large release of survey-quality data from UKIDSS and includes 320 sq degs of multicolour data to (Vega) K=18, complete (depending on the survey) in three to five bands from the set ZYJHK, together with 4 sq degs of deep JK data to an average depth K=21. In addition the release includes a similar quantity of data with incomplete filter coverage. In JHK, in regions of low extinction, the photometric uniformity of the calibration is better than 0.02 mag. in each band. The accuracy of the calibration in ZY remains to be quantified, and the same is true of JHK in regions of high extinction. The median image FWHM across the dataset is 0.82 arcsec. We describe changes since the Early Data Release in the implementation, pipeline and calibration, quality control, and archive procedures. We provide maps of the areas surveyed, and summarise the contents of each of the five surveys in terms of filters, areas, and depths. DR1 marks completion of 7 per cent of the UKIDSS 7-year goals. △ Less

Submitted 30 July, 2007; v1 submitted 6 October, 2006; originally announced October 2006.

Comments: Replaced with accepted version. 12 figures. 6 figures are low resolution. Version of paper with full-resolution figures available at http://astro.ic.ac.uk/~sjw4/ukidss_dr1_accepted.ps.gz

Journal ref: Mon.Not.Roy.Astron.Soc.372:1227-1252,2006; Mon.Not.Roy.Astron.Soc.375:213-226,2007