Showing 101–150 of 372 results for author: Goodman, A

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Adam Ginsburg, Stella S. R. Offner, Erik Rosolowsky, Rachel K. Friesen, Felipe O. Alves, Ana Chacón-Tanarro, Anna Punanova, Elena Redaelli, Helen Kirk, Philip C. Myers, Peter G. Martin, Yancy Shirley, Michael Chun-Yuan Chen, Alyssa A. Goodman, James Di Francesco

Abstract: Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity tow… ▽ More Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations. △ Less

Submitted 20 July, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: Accepted for publication in Astronomy & Astrophysics on 06/07/2020. 15 pages, 16 figures, 1 table. Language edits from previous version

Journal ref: A&A 640, L6 (2020)

Turbulence, Coherence and Collapse: Three Phases for Core Evolution

Authors: Stella S. R. Offner, Josh Taylor, Carleen Markey, Hope How-Huan Chen, Jaime E. Pineda, Alyssa A. Goodman, Andreas Burkert, Adam Ginsburg, Spandan Choudhury

Abstract: We study the formation, evolution and collapse of dense cores by tracking structures in a magnetohydrodynamic simulation of a star-forming cloud. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including Neural Gas prototype learning and Fuzzy $c$-means clustering, to analyze the density and velocity dispersion profiles of cores together with six bulk prop… ▽ More We study the formation, evolution and collapse of dense cores by tracking structures in a magnetohydrodynamic simulation of a star-forming cloud. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including Neural Gas prototype learning and Fuzzy $c$-means clustering, to analyze the density and velocity dispersion profiles of cores together with six bulk properties. We produce a 2-d visualization using a Uniform Manifold Approximation and Projection (UMAP), which facilitates the connection between physical properties and three partially-overlapping phases: i) unbound turbulent structures (Phase I), ii) coherent cores that have low turbulence (Phase II), and iii) bound cores, many of which become protostellar (Phase III). Within Phase II we identify a population of long-lived coherent cores that reach a quasi-equilibrium state. Most prestellar cores form in Phase II and become protostellar after evolving into Phase III. Due to the turbulent cloud environment, the initial core properties do not uniquely predict the eventual evolution, i.e., core evolution is stochastic, and cores follow no one evolutionary path. The phase lifetimes are 1.0$\pm$0.1$\times$10$^5$ yr, 1.3$\pm$0.2$\times$10$^5$ yr, and 1.8$\pm$0.3$\times$10$^5$ yr for Phase I, II, and III, respectively. We compare our results to NH$_3$ observations of dense cores. Known coherent cores predominantly map into Phase II, while most turbulent pressure-confined cores map to Phase I or III. We predict that a significant fraction of observed starless cores have unresolved coherent regions and that $\gtrsim 20$% of observed starless cores will not form stars. Measurements of core radial profiles, in addition to the usual bulk properties, will enable more accurate predictions of core evolution. △ Less

Submitted 26 September, 2022; v1 submitted 12 June, 2020; originally announced June 2020.

Comments: 26 pages, 16 figures. Accepted to MNRAS

HAWC J2227+610 and its association with G106.3+2.7, a new potential Galactic PeVatron

Authors: A. Albert, R. Alfaro, C. Alvarez, J. R. Angeles Camacho, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, V. Baghmanyan, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti, J. Cotzomi, S. Coutiño de León, E. De la Fuente, L. Diaz-Cruz, B. L. Dingus, M. A. DuVernois, J. C. Díaz-Vélez, R. W. Ellsworth , et al. (75 additional authors not shown)

Abstract: We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3… ▽ More We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time. △ Less

Submitted 27 May, 2020; originally announced May 2020.

Comments: Accepted for publication by ApJL

Journal ref: ApJL 896 L29 (2020)

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

Authors: Che-Yu Chen, Erica A. Behrens, Jasmin E. Washington, Laura M. Fissel, Rachel K. Friesen, Zhi-Yun Li, Jaime E. Pineda, Adam Ginsburg, Helen Kirk, Samantha Scibelli, Felipe Alves, Elena Redaelli, Paola Caselli, Anna Punanova, James Di Francesco, Erik Rosolowsky, Stella S. R. Offner, Peter G. Martin, Ana Chacón-Tanarro, Hope H. -H. Chen, Michael C. -Y. Chen, Jared Keown, Youngmin Seo, Yancy Shirley, Hector G. Arce , et al. (4 additional authors not shown)

Abstract: The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular cl… ▽ More The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comments: 18 pages, 11 figures, accepted for publication in MNRAS

A Galactic-scale gas wave in the Solar Neighborhood

Authors: João Alves, Catherine Zucker, Alyssa A. Goodman, Joshua S. Speagle, Stefan Meingast, Thomas Robitaille, Douglas P. Finkbeiner, Edward F. Schlafly, Gregory M. Green

Abstract: For the past 150 years, the prevailing view of the local Interstellar Medium (ISM) was based on a peculiarity known as the Gould's Belt, an expanding ring of young stars, gas, and dust, tilted about 20$^\circ$ to the Galactic plane. Still, the physical relation between local gas clouds has remained practically unknown because the distance accuracy to clouds is of the same order or larger than thei… ▽ More For the past 150 years, the prevailing view of the local Interstellar Medium (ISM) was based on a peculiarity known as the Gould's Belt, an expanding ring of young stars, gas, and dust, tilted about 20$^\circ$ to the Galactic plane. Still, the physical relation between local gas clouds has remained practically unknown because the distance accuracy to clouds is of the same order or larger than their sizes. With the advent of large photometric surveys and the Gaia satellite astrometric survey this situation has changed. Here we report the 3-D structure of all local cloud complexes. We find a narrow and coherent 2.7 kpc arrangement of dense gas in the Solar neighborhood that contains many of the clouds thought to be associated with the Gould Belt. This finding is inconsistent with the notion that these clouds are part of a ring, disputing the Gould Belt model. The new structure comprises the majority of nearby star-forming regions, has an aspect ratio of about 1:20, and contains about 3 million solar masses of gas. Remarkably, the new structure appears to be undulating and its 3-D distribution is well described by a damped sinusoidal wave on the plane of the Milky Way, with an average period of about 2 kpc and a maximum amplitude of about 160 pc. Our results represent a first step in the revision of the local gas distribution and Galactic structure and offer a new, broader context to studies on the transformation of molecular gas into stars. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: Published in Nature on 7 January 2020. For data, interactive visualizations, and more information see: http://tinyurl.com/radwave (https://sites.google.com/cfa.harvard.edu/radcliffewave/)

Journal ref: Nature, 7 January 2020

Constraints on the Emission of Gamma Rays from M31 with HAWC

Authors: HAWC Collaboration, A. Albert, R. Alfaro, C. Alvarez, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti, J. Cotzomi, S. Coutiño de León, E. De la Fuente, C. de León, S. Dichiara, B. L. Dingus, M. A. DuVernois, K. Engel, C. Espinoza , et al. (60 additional authors not shown)

Abstract: Cosmic rays, along with stellar radiation and magnetic fields, are known to make up a significant fraction of the energy density of galaxies such as the Milky Way. When cosmic rays interact in the interstellar medium, they produce gamma-ray emission which provides an important indication of how the cosmic rays propagate. Gamma rays from the Andromeda Galaxy (M31), located 785 kpc away, provide a u… ▽ More Cosmic rays, along with stellar radiation and magnetic fields, are known to make up a significant fraction of the energy density of galaxies such as the Milky Way. When cosmic rays interact in the interstellar medium, they produce gamma-ray emission which provides an important indication of how the cosmic rays propagate. Gamma rays from the Andromeda Galaxy (M31), located 785 kpc away, provide a unique opportunity to study cosmic-ray acceleration and diffusion in a galaxy with a structure and evolution very similar to the Milky Way. Using 33 months of data from the High Altitude Water Cherenkov Observatory, we search for TeV gamma rays from the galactic plane of M31. We also investigate past and present evidence of galactic activity in M31 by searching for Fermi Bubble-like structures above and below the galactic nucleus. No significant gamma-ray emission is observed, so we use the null result to compute upper limits on the energy density of cosmic rays $>10$ TeV in M31. The computed upper limits are approximately ten times higher than expected from the extrapolation of the Fermi LAT results. △ Less

Submitted 25 February, 2020; v1 submitted 13 January, 2020; originally announced January 2020.

Comments: Accepted for publication in ApJ

A Compendium of Distances to Molecular Clouds in the Star Formation Handbook

Authors: Catherine Zucker, Joshua S. Speagle, Edward F. Schlafly, Gregory M. Green, Douglas P. Finkbeiner, Alyssa Goodman, João Alves

Abstract: Accurate distances to local molecular clouds are critical for understanding the star and planet formation process, yet distance measurements are often obtained inhomogeneously on a cloud-by-cloud basis. We have recently developed a method which combines stellar photometric data with Gaia DR2 parallax measurements in a Bayesian framework to infer the distances of nearby dust clouds to a typical acc… ▽ More Accurate distances to local molecular clouds are critical for understanding the star and planet formation process, yet distance measurements are often obtained inhomogeneously on a cloud-by-cloud basis. We have recently developed a method which combines stellar photometric data with Gaia DR2 parallax measurements in a Bayesian framework to infer the distances of nearby dust clouds to a typical accuracy of $\sim5\%$. After refining the technique to target lower latitudes and incorporating deep optical data from DECam in the southern Galactic plane, we have derived a catalog of distances to molecular clouds in Reipurth (2008, Star Formation Handbook, vols I and II) which contains a large fraction of the molecular material in the solar neighborhood. Comparison with distances derived from maser parallax measurements towards the same clouds shows our method produces consistent distances with $\lesssim10\%$ scatter for clouds across our entire distance spectrum (150 pc $-$ 2.5 kpc). We hope this catalog of homogeneous distances will serve as a baseline for future work. △ Less

Submitted 2 January, 2020; originally announced January 2020.

Comments: Accepted for publication in Astronomy & Astrophysics. Full catalog of distances available at https://doi.org/10.7910/DVN/07L7YZ or the CDS

Journal ref: A&A 633, A51 (2020)

Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) -- Full Data Release

Authors: Ian W. Stephens, Tyler L. Bourke, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, John J. Tobin, Héctor G. Arce, Sarah I. Sadavoy, Eduard I. Vorobyov, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Mark A. Gurwell, Alyssa A. Goodman

Abstract: We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1… ▽ More We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $μ$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: Accepted to ApJS

Constraints on Lorentz invariance violation from HAWC observations of gamma rays above 100 TeV

Authors: HAWC Collaboration, A. Albert, R. Alfaro, C. Alvarez, J. R. Angeles Camacho, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, V. Baghmanyan, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti, J. Cotzomi, S. Coutiño de León, E. De la Fuente, C. de León, B. L. Dingus, M. A. DuVernois, J. C. Díaz-Vélez , et al. (68 additional authors not shown)

Abstract: Due to the high energies and long distances to the sources, astrophysical observations provide a unique opportunity to test possible signatures of Lorentz invariance violation (LIV). Superluminal LIV enables the decay of photons at high energy. The High Altitude Water Cherenkov (HAWC) Observatory is among the most sensitive gamma-ray instruments currently operating above 10 TeV. HAWC finds evidenc… ▽ More Due to the high energies and long distances to the sources, astrophysical observations provide a unique opportunity to test possible signatures of Lorentz invariance violation (LIV). Superluminal LIV enables the decay of photons at high energy. The High Altitude Water Cherenkov (HAWC) Observatory is among the most sensitive gamma-ray instruments currently operating above 10 TeV. HAWC finds evidence of 100 TeV photon emission from at least four astrophysical sources. These observations exclude, for the strongest of the limits set, the LIV energy scale to $2.2\times10^{31}$ eV, over 1800 times the Planck energy and an improvement of 1 to 2 orders of magnitude over previous limits. △ Less

Submitted 25 March, 2020; v1 submitted 18 November, 2019; originally announced November 2019.

Comments: Accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 124, 131101 (2020)

Constraining the Local Burst Rate Density of Primordial Black Holes with HAWC

Authors: A. Albert, R. Alfaro, C. Alvarez, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, V. Baghmanyan, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti, J. Cotzomi, E. De la Fuente, C. de León, B. L. Dingus, M. A. DuVernois, J. C. Díaz-Vélez, R. W. Ellsworth, K. L. Engel, C. Espinoza , et al. (62 additional authors not shown)

Abstract: Primordial Black Holes (PBHs) may have been created by density fluctuations in the early Universe and could be as massive as $> 10^9$ solar masses or as small as the Planck mass. It has been postulated that a black hole has a temperature inversely-proportional to its mass and will thermally emit all species of fundamental particles via Hawking Radiation. PBHs with initial masses of… ▽ More Primordial Black Holes (PBHs) may have been created by density fluctuations in the early Universe and could be as massive as $> 10^9$ solar masses or as small as the Planck mass. It has been postulated that a black hole has a temperature inversely-proportional to its mass and will thermally emit all species of fundamental particles via Hawking Radiation. PBHs with initial masses of $\sim 5 \times 10^{14}$ g (approximately one gigaton) should be expiring today with bursts of high-energy gamma radiation in the GeV--TeV energy range. The High Altitude Water Cherenkov (HAWC) Observatory is sensitive to gamma rays with energies of $\sim$300 GeV to past 100 TeV, which corresponds to the high end of the PBH burst spectrum. With its large instantaneous field-of-view of $\sim 2$ sr and a duty cycle over 95%, the HAWC Observatory is well suited to perform an all-sky search for PBH bursts. We conducted a search using 959 days of HAWC data and exclude the local PBH burst rate density above $3400~\mathrm{pc^{-3}~yr^{-1}}$ at 99% confidence, the strongest limit on the local PBH burst rate density from any existing electromagnetic measurement. △ Less

Submitted 17 February, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: Corresponding authors: K.L. Engel & A. Peisker. 13 pages, 5 figures

MSC Class: 85-05

Synthetic Large-Scale Galactic Filaments -- on their Formation, Physical Properties, and Resemblance to Observations

Authors: Catherine Zucker, Rowan Smith, Alyssa Goodman

Abstract: Using a population of large-scale filaments extracted from an AREPO simulation of a Milky Way-like galaxy, we seek to understand the extent to which observed large-scale filament properties (with lengths $\gtrsim 100$ pc) can be explained by galactic dynamics alone. From an observer's perspective in the disk of the galaxy, we identify filaments forming purely due to galactic dynamics, without the… ▽ More Using a population of large-scale filaments extracted from an AREPO simulation of a Milky Way-like galaxy, we seek to understand the extent to which observed large-scale filament properties (with lengths $\gtrsim 100$ pc) can be explained by galactic dynamics alone. From an observer's perspective in the disk of the galaxy, we identify filaments forming purely due to galactic dynamics, without the effects of feedback or local self-gravity. We find that large-scale Galactic filaments are intrinsically rare, and we estimate that at maximum approximately one filament per $\rm kpc^{2}$ should be identified in projection, when viewed from the direction of our Sun in the Milky Way. In this idealized scenario, we find filaments in both the arm and interarm regions, and hypothesize that the former may be due to gas compression in the spiral-potential wells, with the latter due to differential rotation. Using the same analysis pipeline applied previously to observations, we analyze the physical properties of large-scale Galactic filaments, and quantify their sensitivity to projection effects and galactic environment (i.e. whether they lie in the arm or interarm regions). We find that observed "Giant Molecular Filaments" are consistent with being non-self-gravitating structures dominated by galactic dynamics. Straighter, narrower, and denser "Bone-like" filaments, like the paradigmatic Nessie filament, have similar column densities, velocity gradients, and Galactic plane heights ($z\approx$ 0 pc) to those in our simple model, but additional physical effects (such as feedback and self-gravity) must be invoked to explain their lengths and widths. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 887, Issue 2, article id. 186, 17 pp. (2019)

Terrestrial and Martian Heat Flow Limits on Dark Matter

Authors: Joseph Bramante, Andrew Buchanan, Alan Goodman, Eesha Lodhi

Abstract: If dark matter is efficiently captured by a planet, energy released in its annihilation can exceed that planet's total heat output. Building on prior work, we treat Earth's composition and dark matter capture in detail and present improved limits on dark matter-nucleon scattering cross sections for dark matter masses ranging from 0.1 to $10^{10}$ GeV. We also extend Earth limits by applying the sa… ▽ More If dark matter is efficiently captured by a planet, energy released in its annihilation can exceed that planet's total heat output. Building on prior work, we treat Earth's composition and dark matter capture in detail and present improved limits on dark matter-nucleon scattering cross sections for dark matter masses ranging from 0.1 to $10^{10}$ GeV. We also extend Earth limits by applying the same treatment to Mars. The scope of dark matter models considered is expanded to include spin-dependent nuclear interactions including isospin-independent, proton only, and neutron only interactions. We find that Earth and Mars heating bounds are alleviated for dark matter s-wave self-annihilation cross sections $\lesssim 10^{-44}~{\rm cm^2}$. △ Less

Submitted 20 January, 2020; v1 submitted 25 September, 2019; originally announced September 2019.

Comments: 20 pages, 18 figures, PRD version

Journal ref: Phys. Rev. D 101, 043001 (2020)

Multiple Galactic Sources with Emission Above 56 TeV Detected by HAWC

Authors: HAWC Collaboration, A. U. Abeysekara, A. Albert, R. Alfaro, J. R. Angeles Camacho, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, V. Baghmanyan, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti, J. Cotzomi, S. Coutiño de León, E. De la Fuente, C. de León, S. Dichiara, B. L. Dingus, M. A. DuVernois , et al. (77 additional authors not shown)

Abstract: We present the first catalog of gamma-ray sources emitting above 56 and 100 TeV with data from the High Altitude Water Cherenkov (HAWC) Observatory, a wide field-of-view observatory capable of detecting gamma rays up to a few hundred TeV. Nine sources are observed above 56 TeV, all of which are likely Galactic in origin. Three sources continue emitting past 100 TeV, making this the highest-energy… ▽ More We present the first catalog of gamma-ray sources emitting above 56 and 100 TeV with data from the High Altitude Water Cherenkov (HAWC) Observatory, a wide field-of-view observatory capable of detecting gamma rays up to a few hundred TeV. Nine sources are observed above 56 TeV, all of which are likely Galactic in origin. Three sources continue emitting past 100 TeV, making this the highest-energy gamma-ray source catalog to date. We report the integral flux of each of these objects. We also report spectra for three highest-energy sources and discuss the possibility that they are PeVatrons. △ Less

Submitted 9 January, 2020; v1 submitted 18 September, 2019; originally announced September 2019.

Comments: Accepted by Physical Review Letters

Journal ref: Phys. Rev. Lett. 124, 021102 (2020)

HAWC Contributions to the 36th International Cosmic Ray Conference (ICRC2019)

Authors: A. U. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, J. D. Álvarez, J. R. Angeles Camacho, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, V. Baghmanyan, A. S. Barber, J. Becerra Gonzalez, E. Belmont-Moreno, S. Y. BenZvi, D. Berley, J. Braun, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti12, J. Cotzomi, S. Coutiño de León , et al. (105 additional authors not shown)

Abstract: List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA. List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: List of proceedings from the HAWC Collaboration presented at ICRC2019. Corrected typos in the index of the previous version. Follow the "HTML" link to access the list

Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures

Authors: Hope How-Huan Chen, Jaime E. Pineda, Stella S. R. Offner, Alyssa A. Goodman, Andreas Burkert, Rachel K. Friesen, Erik Rosolowsky, Samantha Scibelli, Yancy Shirley

Abstract: We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that… ▽ More We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the "net rotational motions" of the droplets. We find a ratio between rotational and gravitational energies, $β$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $β\sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field. △ Less

Submitted 18 October, 2019; v1 submitted 12 August, 2019; originally announced August 2019.

Comments: accepted to ApJ in Oct, 2019; in press

The Southern Wide-Field Gamma-Ray Observatory (SWGO): A Next-Generation Ground-Based Survey Instrument for VHE Gamma-Ray Astronomy

Authors: P. Abreu, A. Albert, R. Alfaro, C. Alvarez, R. Arceo, P. Assis, F. Barao, J. Bazo, J. F. Beacom, J. Bellido, S. BenZvi, T. Bretz, C. Brisbois, A. M. Brown, F. Brun, M. Buscemi, K. S. Caballero Mora, P. Camarri, A. Carramiñana, S. Casanova, A. Chiavassa, R. Conceição, G. Cotter, P. Cristofari, S. Dasso , et al. (73 additional authors not shown)

Abstract: We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to und… ▽ More We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to understand extreme astrophysical phenomena throughout the universe. A detailed description of science topics addressed by SWGO is available in the science case white paper [1]. The development of SWGO will draw on extensive experience within the community in designing, constructing, and successfully operating wide-field instruments using observations of extensive air showers. The detector will consist of a compact inner array of particle detection units surrounded by a sparser outer array. A key advantage of the design of SWGO is that it can be constructed using current, already proven technology. We estimate a construction cost of 54M USD and a cost of 7.5M USD for 5 years of operation, with an anticipated US contribution of 20M USD ensuring that the US will be a driving force for the SWGO effort. The recently formed SWGO collaboration will conduct site selection and detector optimization studies prior to construction, with full operations foreseen to begin in 2026. Throughout this document, references to science white papers submitted to the Astro2020 Decadal Survey with particular relevance to the key science goals of SWGO, which include unveiling Galactic particle accelerators [2-10], exploring the dynamic universe [11-21], and probing physics beyond the Standard Model [22-25], are highlighted in red boldface. △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: Astro2020 APC White Paper

Quantifying allowable motion to achieve safe dose escalation in pancreatic SBRT

Authors: Yijun Ding, Warren G. Campbell, Moyed Miften, Yevgeniy Vinogradskiy, Karyn A. Goodman, Tracey E. Schefter, Bernard L. Jones

Abstract: Tumor motion plays a key role in the safe delivery of Stereotactic Body Radiotherapy (SBRT) for pancreatic cancer. The purpose of this study was to use tumor motion data measured in patients to establish limits on motion magnitude for safe delivery of pancreatic SBRT. Using 91 sets of pancreatic tumor motion data measured in patients, we calculated motion-convolved dose for 25 pancreatic cancer pa… ▽ More Tumor motion plays a key role in the safe delivery of Stereotactic Body Radiotherapy (SBRT) for pancreatic cancer. The purpose of this study was to use tumor motion data measured in patients to establish limits on motion magnitude for safe delivery of pancreatic SBRT. Using 91 sets of pancreatic tumor motion data measured in patients, we calculated motion-convolved dose for 25 pancreatic cancer patients, and established the maximum amount of motion allowable while satisfying error thresholds on key dose metrics. In our patient cohort, the mean [min-max] allowable motion for 33/40/50 Gy to the PTV was 11.9 [6.3-22.4], 10.4 [5.2-19.1] and 9.0 [4.2-16.0] mm, respectively. Maximum allowable motion decreased as dose was escalated, and was smaller in patients with larger tumors. The effects of motion on pancreatic SBRT are highly variable between patients and there is potential to allow more motion in certain patients, even in dose-escalated scenarios. In our dataset, a conservative limit of 6.3 mm would ensure safe treatment of all patients treated to 33 Gy in 5 fractions. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: 9 pages, 5 figures, supplement included. Submitted to Practical Radiation Oncology

Measurement of the Crab Nebula Spectrum Past 100 TeV with HAWC

Authors: HAWC Collaboration, A. U. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, J. D. Álvarez, J. R. Angeles Camacho, R. Acero, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, V. Baghmanyan, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Cabellero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti, J. Cotzomi, S. Coutiño de León, E. De la Fuente, C. de León , et al. (80 additional authors not shown)

Abstract: We present TeV gamma-ray observations of the Crab Nebula, the standard reference source in ground-based gamma-ray astronomy, using data from the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory. In this analysis we use two independent energy-estimation methods that utilize extensive air shower variables such as the core position, shower angle, and shower lateral energy distribution. In c… ▽ More We present TeV gamma-ray observations of the Crab Nebula, the standard reference source in ground-based gamma-ray astronomy, using data from the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory. In this analysis we use two independent energy-estimation methods that utilize extensive air shower variables such as the core position, shower angle, and shower lateral energy distribution. In contrast, the previously published HAWC energy spectrum roughly estimated the shower energy with only the number of photomultipliers triggered. This new methodology yields a much improved energy resolution over the previous analysis and extends HAWC's ability to accurately measure gamma-ray energies well beyond 100 TeV. The energy spectrum of the Crab Nebula is well fit to a log parabola shape $\left(\frac{dN}{dE} = φ_0 \left(E/\textrm{7 TeV}\right)^{-α-β\ln\left(E/\textrm{7 TeV}\right)}\right)$ with emission up to at least 100 TeV. For the first estimator, a ground parameter that utilizes fits to the lateral distribution function to measure the charge density 40 meters from the shower axis, the best-fit values are $φ_o$=(2.35$\pm$0.04$^{+0.20}_{-0.21}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $α$=2.79$\pm$0.02$^{+0.01}_{-0.03}$, and $β$=0.10$\pm$0.01$^{+0.01}_{-0.03}$. For the second estimator, a neural network which uses the charge distribution in annuli around the core and other variables, these values are $φ_o$=(2.31$\pm$0.02$^{+0.32}_{-0.17}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $α$=2.73$\pm$0.02$^{+0.03}_{-0.02}$, and $β$=0.06$\pm$0.01$\pm$0.02. The first set of uncertainties are statistical; the second set are systematic. Both methods yield compatible results. These measurements are the highest-energy observation of a gamma-ray source to date. △ Less

Submitted 17 September, 2019; v1 submitted 29 May, 2019; originally announced May 2019.

Comments: published in ApJ

Journal ref: The Astrophysical Journal, Volume 881, Number 2 (2019)

Efficient methanol production on the dark side of a prestellar core

Authors: Jorma Harju, Jaime E. Pineda, Anton I. Vasyunin, Paola Caselli, Stella S. R. Offner, Alyssa A. Goodman, Mika Juvela, Olli Sipilae, Alexandre Faure, Romane Le Gal, Pierre Hily-Blant, Joao Alves, Luca Bizzocchi, Andreas Burkert, Hope Chen, Rachel K. Friesen, Rolf Guesten, Philip C. Myers, Anna Punanova, Claire Rist, Erik Rosolowsky, Stephan Schlemmer, Yancy Shirley, Silvia Spezzano, Charlotte Vastel , et al. (1 additional authors not shown)

Abstract: We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particul… ▽ More We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity. △ Less

Submitted 4 May, 2020; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: Accepted for publication in the Astrophysical Journal

Pulsars in a Bubble? Following Electron Diffusion in the Galaxy with TeV Gamma Rays

Authors: H. Fleischhack, A. Albert, C. Alvarez, R. Arceo, J. F. Beacom, R. Bird, C. A. Brisbois, K. S. Caballero-Mora, A. Carraminana, S. Casanova, P. Cristofari, P. Coppi, B. L. Dingus, M. A. DuVernois, K. L. Engel, J. A. Goodman, T. Greenshaw, J. P. Harding, B. Hona, P. H. Huentemeyer, H. Li, T. Linden, K. Malone, J. Martinez-Castro, M. A. Mostafa , et al. (10 additional authors not shown)

Abstract: TeV Halos, extended regions of TeV gamma-ray emission around middle-aged pulsars, have recently been established as a new source class in gamma-ray astronomy. These halos have been attributed to relativistic electrons and positrons that have left the acceleration region close to the pulsar and are diffusing in the surrounding medium. Measuring the morphology of TeV Halos enables for the first ti… ▽ More TeV Halos, extended regions of TeV gamma-ray emission around middle-aged pulsars, have recently been established as a new source class in gamma-ray astronomy. These halos have been attributed to relativistic electrons and positrons that have left the acceleration region close to the pulsar and are diffusing in the surrounding medium. Measuring the morphology of TeV Halos enables for the first time a direct measurement of the electron diffusion on scales of tens of parsecs. There are hints that the presence of relativistic particles affects the diffusion rate in the pulsars' surroundings. Understanding electron diffusion is necessary to constrain the origins of the apparent `excess' of cosmic-ray positrons at tens of GeV. TeV Halos can also be used to find mis-aligned pulsars, as well as study certain properties of the Galaxy's pulsar population. Future VHE gamma-ray instruments will detect more of those TeV Halos and determine how much pulsars contribute to the observed cosmic-ray electron and positron fluxes, and how they affect diffusion in their environments. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Whitepaper for the US Astro2020 decadal survey

Astro2020 Science White Paper: The Next Decade of Astroinformatics and Astrostatistics

Authors: A. Siemiginowska, G. Eadie, I. Czekala, E. Feigelson, E. B. Ford, V. Kashyap, M. Kuhn, T. Loredo, M. Ntampaka, A. Stevens, A. Avelino, K. Borne, T. Budavari, B. Burkhart, J. Cisewski-Kehe, F. Civano, I. Chilingarian, D. A. van Dyk, G. Fabbiano, D. P. Finkbeiner, D. Foreman-Mackey, P. Freeman, A. Fruscione, A. A. Goodman, M. Graham , et al. (27 additional authors not shown)

Abstract: Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New met… ▽ More Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers

The Role of Machine Learning in the Next Decade of Cosmology

Authors: Michelle Ntampaka, Camille Avestruz, Steven Boada, Joao Caldeira, Jessi Cisewski-Kehe, Rosanne Di Stefano, Cora Dvorkin, August E. Evrard, Arya Farahi, Doug Finkbeiner, Shy Genel, Alyssa Goodman, Andy Goulding, Shirley Ho, Arthur Kosowsky, Paul La Plante, Francois Lanusse, Michelle Lochner, Rachel Mandelbaum, Daisuke Nagai, Jeffrey A. Newman, Brian Nord, J. E. G. Peek, Austin Peel, Barnabas Poczos , et al. (5 additional authors not shown)

Abstract: In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster an… ▽ More In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors. △ Less

Submitted 14 January, 2021; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: Submitted to the Astro2020 call for science white papers

Science Case for a Wide Field-of-View Very-High-Energy Gamma-Ray Observatory in the Southern Hemisphere

Authors: A. Albert, R. Alfaro, H. Ashkar, C. Alvarez, J. Álvarez, J. C. Arteaga-Velázquez, H. A. Ayala Solares, R. Arceo, J. A. Bellido, S. BenZvi, T. Bretz, C. A. Brisbois, A. M. Brown, F. Brun, K. S. Caballero-Mora, A. Carosi, A. Carramiñana, S. Casanova, P. M. Chadwick, G. Cotter, S. Coutiño De León, P. Cristofari, S. Dasso, E. de la Fuente, B. L. Dingus , et al. (78 additional authors not shown)

Abstract: We outline the science motivation for SGSO, the Southern Gamma-Ray Survey Observatory. SGSO will be a next-generation wide field-of-view gamma-ray survey instrument, sensitive to gamma-rays in the energy range from 100 GeV to hundreds of TeV. Its science topics include unveiling galactic and extragalactic particle accelerators, monitoring the transient sky at very high energies, probing particle p… ▽ More We outline the science motivation for SGSO, the Southern Gamma-Ray Survey Observatory. SGSO will be a next-generation wide field-of-view gamma-ray survey instrument, sensitive to gamma-rays in the energy range from 100 GeV to hundreds of TeV. Its science topics include unveiling galactic and extragalactic particle accelerators, monitoring the transient sky at very high energies, probing particle physics beyond the Standard Model, and the characterization of the cosmic ray flux. SGSO will consist of an air shower detector array, located in South America. Due to its location and large field of view, SGSO will be complementary to other current and planned gamma-ray observatories such as HAWC, LHAASO, and CTA. △ Less

Submitted 22 February, 2019; originally announced February 2019.

A Large Catalog of Accurate Distances to Local Molecular Clouds: The Gaia DR2 Edition

Authors: Catherine Zucker, Joshua S. Speagle, Edward F. Schlafly, Gregory M. Green, Douglas P. Finkbeiner, Alyssa A. Goodman, João Alves

Abstract: We present a uniform catalog of accurate distances to local molecular clouds informed by the Gaia DR2 data release. Our methodology builds on that of Schlafly et al. (2014). First, we infer the distance and extinction to stars along sightlines towards the clouds using optical and near-infrared photometry. When available, we incorporate knowledge of the stellar distances obtained from Gaia DR2 para… ▽ More We present a uniform catalog of accurate distances to local molecular clouds informed by the Gaia DR2 data release. Our methodology builds on that of Schlafly et al. (2014). First, we infer the distance and extinction to stars along sightlines towards the clouds using optical and near-infrared photometry. When available, we incorporate knowledge of the stellar distances obtained from Gaia DR2 parallax measurements. We model these per-star distance-extinction estimates as being caused by a dust screen with a 2-D morphology derived from Planck at an unknown distance, which we then fit for using a nested sampling algorithm. We provide updated distances to the Schlafly et al. (2014) sightlines towards the Dame et al. (2001) and Magnani et al. (1985) clouds, finding good agreement with the earlier work. For a subset of 27 clouds, we construct interactive pixelated distance maps to further study detailed cloud structure, and find several clouds which display clear distance gradients and/or are comprised of multiple components. We use these maps to determine robust average distances to these clouds. The characteristic combined uncertainty on our distances is approximately 5-6%, though this can be higher for clouds at farther distances, due to the limitations of our single-cloud model. △ Less

Submitted 19 April, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: Accepted for publication in The Astrophysical Journal. Many figures are interactive; click the link in the figure captions of the PDF to access, or see the online version of the published article

Journal ref: The Astrophysical Journal, Volume 879, Issue 2, article id. 125, 20 pp. (2019)

MAGIC and Fermi-LAT gamma-ray results on unassociated HAWC sources

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak, P. Colin , et al. (318 additional authors not shown)

Abstract: The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the fir… ▽ More The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail. △ Less

Submitted 13 January, 2019; originally announced January 2019.

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

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 485, Issue 1, May 2019, Pages 356-366

All-Sky Measurement of the Anisotropy of Cosmic Rays at 10 TeV and Mapping of the Local Interstellar Magnetic Field

Authors: HAWC Collaboration, A. U. Abeysekara, R. Alfaro, C. Alvarez, J. D. Álvarez, R. Arceo, J. C. Arteaga-Velázquez, D. Avila Rojas, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti, J. Cotzomi, J. C. Díaz-Vélez, C. De León, E. De la Fuente, S. Dichiara, M. A. DuVernois, C. Espinoza, D. W. Fiorino, H. Fleischhack, N. Fraija , et al. (382 additional authors not shown)

Abstract: We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the HAWC and IceCube observatories in the Northern and Southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and holds a key to probe into the propagation pr… ▽ More We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the HAWC and IceCube observatories in the Northern and Southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and holds a key to probe into the propagation properties of TeV cosmic rays through our local interstellar medium and the interaction between the interstellar and heliospheric magnetic fields. From the map we determine the horizontal dipole components of the anisotropy $δ_{0h} = 9.16 \times 10^{-4}$ and $δ_{6h} = 7.25 \times 10^{-4}~(\pm0.04 \times 10^{-4})$. In addition, we infer the direction ($229.2\pm 3.5^\circ$ RA , $11.4\pm 3.0^\circ$ Dec.) of the interstellar magnetic field from the boundary between large scale excess and deficit regions from which we estimate the missing corresponding vertical dipole component of the large scale anisotropy to be $δ_N \sim -3.97 ^{+1.0}_{-2.0} \times 10^{-4}$. △ Less

Submitted 24 January, 2019; v1 submitted 13 December, 2018; originally announced December 2018.

Comments: Accepted by The Astrophysical Journal

Journal ref: ApJ, 871, (2019) 96

Searching for Dark Matter Sub-structure with HAWC

Authors: A. U. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, R. Arceo, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, A. Carramiñana, S. Casanova, J. Cotzomi, S. Coutiño de León, C. De León, E. De la Fuente, S. Dichiara, B. L. Dingus, M. A. DuVernois, J. C. Díaz-Vélez, K. Engel, C. Espinoza, H. Fleischhack , et al. (53 additional authors not shown)

Abstract: Simulations of dark matter show a discrepancy between the expected number of Galactic dark matter sub-halos and how many have been optically observed. Some of these unseen satellites may exist as dark dwarf galaxies: sub-halos like dwarf galaxies with no luminous counterpart. Assuming WIMP dark matter, it may be possible to detect these unseen sub-halos from gamma-ray signals originating from dark… ▽ More Simulations of dark matter show a discrepancy between the expected number of Galactic dark matter sub-halos and how many have been optically observed. Some of these unseen satellites may exist as dark dwarf galaxies: sub-halos like dwarf galaxies with no luminous counterpart. Assuming WIMP dark matter, it may be possible to detect these unseen sub-halos from gamma-ray signals originating from dark matter annihilation. The High Altitude Water Cherenkov Observatory (HAWC) is a very high energy (500 GeV to 100 TeV) gamma ray detector with a wide field-of-view and near continuous duty cycle, making HAWC ideal for unbiased sky surveys. We perform such a search for gamma ray signals from dark dwarfs in the Milky Way halo. We perform a targeted search of HAWC gamma-ray sources which have no known association with lower-energy counterparts, based on an unbiased search of the entire sky. With no sources found to strongly prefer dark matter models, we calculate the ability of HAWC to observe dark dwarfs. We also compute the HAWC sensitivity to potential future detections for a given model of dark matter substructure. △ Less

Submitted 20 June, 2019; v1 submitted 28 November, 2018; originally announced November 2018.

Comments: 20 pages, 10 figures. submitted to JCAP

Suppressing Diffusion-Mediated Exciton Annihilation in 2D Semiconductors Using the Dielectric Environment

Authors: Aaron J. Goodman, Der-Hsien Lien, Geun Ho Ahn, Leo L. Spiegel, Matin Amani, Adam P. Willard, Ali Javey, William A. Tisdale

Abstract: Atomically thin semiconductors such as monolayer MoS2 and WS2 exhibit nonlinear exciton-exciton annihilation at notably low excitation densities (below ~10 excitons/um2 in MoS2). Here, we show that the density threshold at which annihilation occurs can be tuned by changing the underlying substrate. When the supporting substrate is changed from SiO2 to Al2O3 or SrTiO3, the rate constant for second-… ▽ More Atomically thin semiconductors such as monolayer MoS2 and WS2 exhibit nonlinear exciton-exciton annihilation at notably low excitation densities (below ~10 excitons/um2 in MoS2). Here, we show that the density threshold at which annihilation occurs can be tuned by changing the underlying substrate. When the supporting substrate is changed from SiO2 to Al2O3 or SrTiO3, the rate constant for second-order exciton-exciton annihilation, k_XX [cm2/s], is reduced by one or two orders of magnitude, respectively. Using transient photoluminescence microscopy, we measure the effective room-temperature exciton diffusion coefficient in chemical-treated MoS2 to be D = 0.06 +/- 0.01 cm2/s, corresponding to a diffusion length of LD = 350 nm for an exciton lifetime of τ = 20 ns, which is independent of the substrate. These results, together with numerical simulations, suggest that the effective exciton-exciton annihilation radius monotonically decreases with increasing refractive index of the underlying substrate. Exciton-exciton annihilation limits the overall efficiency of 2D semiconductor devices operating at high exciton densities; the ability to tune these interactions via the dielectric environment is an important step toward more efficient optoelectronic technologies featuring atomically thin materials. △ Less

Submitted 2 November, 2018; originally announced November 2018.

Very high energy particle acceleration powered by the jets of the microquasar SS 433

Authors: HAWC Collaboration, A. U. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, J. D. Álvarez, R. Arceo, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, M. Castillo, U. Cotti, J. Cotzomi, S. Coutiño de León, C. De León, E. De la Fuente, J. C. Díaz-Vélez, S. Dichiara , et al. (82 additional authors not shown)

Abstract: SS 433 is a binary system containing a supergiant star that is overflowing its Roche lobe with matter accreting onto a compact object (either a black hole or neutron star). Two jets of ionized matter with a bulk velocity of $\sim0.26c$ extend from the binary, perpendicular to the line of sight, and terminate inside W50, a supernova remnant that is being distorted by the jets. SS 433 differs from o… ▽ More SS 433 is a binary system containing a supergiant star that is overflowing its Roche lobe with matter accreting onto a compact object (either a black hole or neutron star). Two jets of ionized matter with a bulk velocity of $\sim0.26c$ extend from the binary, perpendicular to the line of sight, and terminate inside W50, a supernova remnant that is being distorted by the jets. SS 433 differs from other microquasars in that the accretion is believed to be super-Eddington, and the luminosity of the system is $\sim10^{40}$ erg s$^{-1}$. The lobes of W50 in which the jets terminate, about 40 pc from the central source, are expected to accelerate charged particles, and indeed radio and X-ray emission consistent with electron synchrotron emission in a magnetic field have been observed. At higher energies (>100 GeV), the particle fluxes of $γ$ rays from X-ray hotspots around SS 433 have been reported as flux upper limits. In this energy regime, it has been unclear whether the emission is dominated by electrons that are interacting with photons from the cosmic microwave background through inverse-Compton scattering or by protons interacting with the ambient gas. Here we report TeV $γ$-ray observations of the SS 433/W50 system where the lobes are spatially resolved. The TeV emission is localized to structures in the lobes, far from the center of the system where the jets are formed. We have measured photon energies of at least 25 TeV, and these are certainly not Doppler boosted, because of the viewing geometry. We conclude that the emission from radio to TeV energies is consistent with a single population of electrons with energies extending to at least hundreds of TeV in a magnetic field of $\sim16$~micro-Gauss. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: Preprint version of Nature paper. Contacts: S. BenZvi, B. Dingus, K. Fang, C.D. Rho , H. Zhang, H. Zhou

Journal ref: Nature 562 (2018), 82-85

Droplets I: Pressure-Dominated Sub-0.1 pc Coherent Structures in L1688 and B18

Authors: Hope How-Huan Chen, Jaime E. Pineda, Alyssa A. Goodman, Andreas Burkert, Stella S. R. Offner, Rachel K. Friesen, Philip C. Myers, Felipe Alves, Hector G. Arce, Paola Caselli, Ana Chacon-Tanarro, Michael Chun-Yuan Chen, James Di Francesco, Adam Ginsburg, Jared Keown, Helen Kirk, Peter G. Martin, Christopher Matzner, Anna Punanova, Elena Redaelli, Erik Rosolowsky, Samantha Scibelli, Young Min Seo, Yancy Shirley, Ayushi Singh

Abstract: We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transiti… ▽ More We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process. △ Less

Submitted 15 May, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: Accepted by ApJ in April, 2019

Journal ref: 2019ApJ...877...93C

Histogram of oriented gradients: a technique for the study of molecular cloud formation

Authors: J. D. Soler, H. Beuther, M. Rugel, Y. Wang, P. C. Clark, S. C. O. Glover, P. F. Goldsmith, M. Heyer, L. D. Anderson, A. Goodman, Th. Henning, J. Kainulainen, R. S. Klessen, S. N. Longmore, N. M. McClure-Griffiths, K. M. Menten, J. C. Mottram, J. Ott, S. E. Ragan, R. J. Smith, J. S. Urquhart, F. Bigiel, P. Hennebelle, N. Roy, P. Schilke

Abstract: We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial corre… ▽ More We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels. We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight. We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall. △ Less

Submitted 2 January, 2019; v1 submitted 21 September, 2018; originally announced September 2018.

Comments: 32 pages, 36 figures. Accepted for publication at A&A (28DEC2018)

Report number: AA/2018/34300

Journal ref: A&A 622, A166 (2019)

VERITAS and Fermi-LAT observations of new HAWC sources

Authors: VERITAS Collaboration, A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, V. Bugaev, A. J. Chromey, M. P. Connolly, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, M. Hutten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson , et al. (259 additional authors not shown)

Abstract: The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detect… ▽ More The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Comments: Accepted for publication in the ApJ, Corresponding author: Nahee Park (VERITAS Collaboration), John W. Hewitt (Fermi-LAT Collaboration), Ignacio Taboada (HAWC Collaboration), 30 pages, 12 figures

Journal ref: ApJ 866 (2018) no.1, 24

3D shape of Orion A from Gaia DR2

Authors: Josefa E. Grossschedl, Joao Alves, Stefan Meingast, Christine Ackerl, Joana Ascenso, Herve Bouy, Andreas Burkert, Jan Forbrich, Verena Fuernkranz, Alyssa Goodman, Alvaro Hacar, Gabor Herbst-Kiss, Charles J. Lada, Irati Larreina, Kieran Leschinski, Marco Lombardi, Andre Moitinho, Daniel Mortimer, Eleonora Zari

Abstract: We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhan… ▽ More We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhanced star-forming (bent) Head, and a lower density and star-formation quieter $\sim$75 pc long Tail. The true extent of Orion A is not the projected $\sim$40 pc but $\sim$90 pc, making it by far the largest molecular cloud in the local neighborhood. Its aspect ratio ($\sim$30:1) and high column-density fraction ($\sim45\%$) make it similar to large-scale Milky Way filaments ("bones"), despite its distance to the galactic mid-plane being an order of magnitude larger than typically found for these structures. △ Less

Submitted 17 August, 2018; originally announced August 2018.

Comments: 9 pages, 5 figures, 4 tables, accepted for publication in A&A

Journal ref: A&A 619, A106 (2018)

Constraints on Spin-Dependent Dark Matter Scattering with Long-Lived Mediators from TeV Observations of the Sun with HAWC

Authors: A. Albert, R. Alfaro, C. Alvarez, R. Arceo, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, T. Capistràn, A. Carramiñana, S. Casanova, M. Castillo, J. Cotzomi, S. Coutiño de León, C. De León, E. De la Fuente, S. Dichiara, B. L. Dingus, M. A. DuVernois, J. C. Díaz-Vélez, K. Engel, O. Enríquez-Rivera , et al. (69 additional authors not shown)

Abstract: We analyze the Sun as a source for the indirect detection of dark matter through a search for gamma rays from the solar disk. Capture of dark matter by elastic interactions with the solar nuclei followed by annihilation to long-lived mediators can produce a detectable gamma-ray flux. We search three years of data from the High Altitude Water Cherenkov Observatory and find no statistically signific… ▽ More We analyze the Sun as a source for the indirect detection of dark matter through a search for gamma rays from the solar disk. Capture of dark matter by elastic interactions with the solar nuclei followed by annihilation to long-lived mediators can produce a detectable gamma-ray flux. We search three years of data from the High Altitude Water Cherenkov Observatory and find no statistically significant detection of TeV gamma-ray emission from the Sun. Using this, we constrain the spin-dependent elastic scattering cross section of dark matter with protons for dark matter masses above 1 TeV, assuming an unstable mediator with a favorable lifetime. The results complement constraints obtained from Fermi-LAT observations of the Sun and together cover WIMP masses between 4 GeV and $10^6$ GeV. The cross section constraints for mediator decays to gamma rays can be as strong as $\sim10^{-45}$ cm$^{-2}$, which is more than four orders of magnitude stronger than current direct-detection experiments for 1 TeV dark matter mass. The cross-section constraints at higher masses are even better, nearly 7 orders of magnitude better than the current direct-detection constraints for 100 TeV dark matter mass. This demonstration of sensitivity encourages detailed development of theoretical models in light of these powerful new constraints. △ Less

Submitted 15 November, 2018; v1 submitted 16 August, 2018; originally announced August 2018.

Comments: 11 pages, 4 figures. See also companion paper 1808.05620. Accepted for publication in Physical Review D

Report number: MIT-CTP/5038

Journal ref: Phys. Rev. D 98, 123012 (2018)

First HAWC Observations of the Sun Constrain Steady TeV Gamma-Ray Emission

Authors: A. Albert, R. Alfaro, C. Alvarez, R. Arceo, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, E. Belmont-Moreno, S. Y. BenZvi, C. Brisbois, K. S. Caballero-Mora, T. Capistràn, A. Carramiñana, S. Casanova, M. Castillo, J. Cotzomi, S. Coutiño de León, C. De León, E. De la Fuente, S. Dichiara, B. L. Dingus, M. A. DuVernois, J. C. Díaz-Vélez, K. Engel, O. Enríquez-Rivera , et al. (70 additional authors not shown)

Abstract: Steady gamma-ray emission up to at least 200 GeV has been detected from the solar disk in the Fermi-LAT data, with the brightest, hardest emission occurring during solar minimum. The likely cause is hadronic cosmic rays undergoing collisions in the Sun's atmosphere after being redirected from ingoing to outgoing in magnetic fields, though the exact mechanism is not understood. An important new tes… ▽ More Steady gamma-ray emission up to at least 200 GeV has been detected from the solar disk in the Fermi-LAT data, with the brightest, hardest emission occurring during solar minimum. The likely cause is hadronic cosmic rays undergoing collisions in the Sun's atmosphere after being redirected from ingoing to outgoing in magnetic fields, though the exact mechanism is not understood. An important new test of the gamma-ray production mechanism will follow from observations at higher energies. Only the High Altitude Water Cherenkov (HAWC) Observatory has the required sensitivity to effectively probe the Sun in the TeV range. Using three years of HAWC data from November 2014 to December 2017, just prior to the solar minimum, we search for 1--100 TeV gamma rays from the solar disk. No evidence of a signal is observed, and we set strong upper limits on the flux at a few $10^{-12}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. Our limit, which is the most constraining result on TeV gamma rays from the Sun, is $\sim10\%$ of the theoretical maximum flux (based on a model where all incoming cosmic rays produce outgoing photons), which in turn is comparable to the Fermi-LAT data near 100 GeV. The prospects for a first TeV detection of the Sun by HAWC are especially high during solar minimum, which began in early 2018. △ Less

Submitted 1 November, 2018; v1 submitted 16 August, 2018; originally announced August 2018.

Comments: 14 pages, 6 figures. See also companion paper 1808.05624. Accepted for publication in Physical Review D

Report number: MIT-CTP/5037

Journal ref: Phys. Rev. D 98, 123011 (2018)

Multi-Messenger Astrophysics: Harnessing the Data Revolution

Authors: Gabrielle Allen, Warren Anderson, Erik Blaufuss, Joshua S. Bloom, Patrick Brady, Sarah Burke-Spolaor, S. Bradley Cenko, Andrew Connolly, Peter Couvares, Derek Fox, Avishay Gal-Yam, Suvi Gezari, Alyssa Goodman, Darren Grant, Paul Groot, James Guillochon, Chad Hanna, David W. Hogg, Kelly Holley-Bockelmann, D. Andrew Howell, David Kaplan, Erik Katsavounidis, Marek Kowalski, Luis Lehner, Daniel Muthukrishna , et al. (5 additional authors not shown)

Abstract: The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugura… ▽ More The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe. △ Less

Submitted 12 July, 2018; originally announced July 2018.

Comments: workshop summary, 20 pages, by the participants of the Cyberinfrastructure for Multi-Messenger Astrophysics NSF-funded Workshop Held on 23-24 May 2018 at the University of Maryland

Mass Assembly of Stellar Systems and their Evolution with the SMA - 1.3 mm Subcompact Data Release

Authors: Ian W. Stephens, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, Tyler L. Bourke, Eduard I. Vorobyov, John J. Tobin, Sarah I. Sadavoy, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Alyssa A. Goodman, Héctor G. Arce, Mark Gurwell

Abstract: We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we ca… ▽ More We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES . △ Less

Submitted 19 June, 2018; originally announced June 2018.

Comments: Accepted ApJS

Ultrafast Charge Transfer at a Quantum Dot/2D Materials Interface Probed by Second Harmonic Generation

Authors: Aaron J. Goodman, Nabeel S. Dahod, William A. Tisdale

Abstract: Hybrid quantum dot (QD) / transition metal dichalcogenide (TMD) heterostructures are attractive components of next generation optoelectronic devices, which take advantage of the spectral tunability of QDs and the charge and exciton transport properties of TMDs. Here, we demonstrate tunable electronic coupling between CdSe QDs and monolayer WS$_2$ using variable length alkanethiol ligands on the QD… ▽ More Hybrid quantum dot (QD) / transition metal dichalcogenide (TMD) heterostructures are attractive components of next generation optoelectronic devices, which take advantage of the spectral tunability of QDs and the charge and exciton transport properties of TMDs. Here, we demonstrate tunable electronic coupling between CdSe QDs and monolayer WS$_2$ using variable length alkanethiol ligands on the QD surface. Using femtosecond time-resolved second harmonic generation (SHG) microscopy, we show that electron transfer from photoexcited CdSe QDs to single-layer WS$_2$ occurs on ultrafast (50 fs - 1 ps) timescales. Moreover, in the samples exhibiting the fastest charge transfer rates ($\leq$ 50 fs) we observed oscillations in the time-domain signal corresponding to an acoustic phonon mode of the donor QD, which coherently modulates the SHG response of the underlying WS$_2$ layer. These results reveal surprisingly strong electronic coupling at the QD/TMD interface and demonstrate the usefulness of time-resolved SHG for exploring ultrafast electronic-vibrational dynamics in TMD heterostructures. △ Less

Submitted 10 June, 2018; originally announced June 2018.

Journal ref: J. Phys. Chem. Lett. 9 (2018), 4227-4232

Phase-Modulated Degenerate Parametric Amplification Microscopy

Authors: Yunan Gao, Aaron J. Goodman, Pin-Chun Shen, Jing Kong, William A. Tisdale

Abstract: Second-order nonlinear optical interactions, including second harmonic generation (SHG) and sum-frequency generation (SFG), can reveal a wealth of information about chemical, electronic, and vibrational dynamics at the nanoscale. Here, we demonstrate a powerful and flexible new approach, called phase-modulated degenerate parametric amplification (DPA). The technique, which allows for facile retrie… ▽ More Second-order nonlinear optical interactions, including second harmonic generation (SHG) and sum-frequency generation (SFG), can reveal a wealth of information about chemical, electronic, and vibrational dynamics at the nanoscale. Here, we demonstrate a powerful and flexible new approach, called phase-modulated degenerate parametric amplification (DPA). The technique, which allows for facile retrieval of both the amplitude and phase of the second-order nonlinear optical response, has many advantages over conventional or heterodyne-detected SHG, including the flexibility to detect the signal at either the second harmonic or fundamental field wavelength. We demonstrate the capabilities of this approach by imaging multi-grain flakes of single-layer MoS2. We identify the absolute crystal orientation of each MoS2 domain and resolve grain boundaries with high signal contrast and sub-diffraction-limited spatial resolution. This robust all-optical method can be used to characterize structure and dynamics in organic and inorganic systems, including biological tissue, soft materials, and metal and semiconductor nanostructures, and is particularly well-suited for imaging in media that are absorptive or highly scattering to visible and ultraviolet light. △ Less

Submitted 9 June, 2018; originally announced June 2018.

Journal ref: Nano Letters (2018)

Dense gas kinematics and a narrow filament in the Orion A OMC1 region using NH3

Authors: Kristina Monsch, Jaime E. Pineda, Hauyu Baobab Liu, Catherine Zucker, Hope How-Huan Chen, Kate Pattle, Stella S. R. Offner, James Di Francesco, Adam Ginsburg, Barbara Ercolano, Héctor G. Arce, Rachel Friesen, Helen Kirk, Paola Caselli, Alyssa A. Goodman

Abstract: We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL reg… ▽ More We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments. △ Less

Submitted 18 December, 2018; v1 submitted 5 June, 2018; originally announced June 2018.

Comments: Accepted for publication in ApJ. The combined data cubes of the NH3 (1,1) and (2,2) transitions as well as the resulting parameter maps are provided as FITS-files on Harvard Dataverse (https://doi.org/10.7910/DVN/QLD7TC); updated references

Journal ref: ApJ 861, 77 (2018)

New constraints on oscillation parameters from $ν_e$ appearance and $ν_μ$ disappearance in the NOvA experiment

Authors: M. A. Acero, P. Adamson, L. Aliaga T. Alion, V. Allakhverdian, N. Anfimov A. Antoshkin, E. Arrieta-Diaz, A. Aurisano A. Back, C. Backhouse, M. Baird N. Balashov, B. A. Bambah, K. Bays B. Behera, S. Bending, R. Bernstein V. Bhatnagar, B. Bhuyan, J. Bian T. Blackburn, J. Blair, A. Bolshakova P. Bour, C. Bromberg, J. Brown N. Buchanan, A. Butkevich, V. Bychkov M. Campbell, T. J. Carroll, E. Catano-Mur A. Cedeno, S. Childress, B. C. Choudhary B. Chowdhury , et al. (104 additional authors not shown)

Abstract: We present updated results from the NOvA experiment for $ν_μ\rightarrowν_μ$ and $ν_μ\rightarrowν_e$ oscillations from an exposure of $8.85\times10^{20}$ protons on target, which represents an increase of 46% compared to our previous publication. The results utilize significant improvements in both the simulations and analysis of the data. A joint fit to the data for $ν_μ$ disappearance and $ν_e$ a… ▽ More We present updated results from the NOvA experiment for $ν_μ\rightarrowν_μ$ and $ν_μ\rightarrowν_e$ oscillations from an exposure of $8.85\times10^{20}$ protons on target, which represents an increase of 46% compared to our previous publication. The results utilize significant improvements in both the simulations and analysis of the data. A joint fit to the data for $ν_μ$ disappearance and $ν_e$ appearance gives the best fit point as normal mass hierarchy, $Δm^2_{32} = 2.44\times 10^{-3}{{\rm eV}^2}/c^4$, $\sin^2θ_{23} = 0.56$, and $δ_{CP} = 1.21π$. The 68.3% confidence intervals in the normal mass hierarchy are $Δm^2_{32} \in [2.37,2.52]\times 10^{-3}{{\rm eV}^2}/c^4$, $\sin^2θ_{23} \in [0.43,0.51] \cup [0.52,0.60]$, and $δ_{CP} \in [0,0.12π] \cup [0.91π,2π]$. The inverted mass hierarchy is disfavored at the 95% confidence level for all choices of the other oscillation parameters. △ Less

Submitted 6 August, 2018; v1 submitted 31 May, 2018; originally announced June 2018.

Journal ref: Phys. Rev. D 98, 032012 (2018)

New Thinking on, and with, Data Visualization

Authors: Alyssa A. Goodman, Michelle A. Borkin, Thomas P. Robitaille

Abstract: As the complexity and volume of datasets have increased along with the capabilities of modular, open-source, easy-to-implement, visualization tools, scientists' need for, and appreciation of, data visualization has risen too. Until recently, scientists thought of the "explanatory" graphics created at a research project's conclusion as "pretty pictures" needed only for journal publication or public… ▽ More As the complexity and volume of datasets have increased along with the capabilities of modular, open-source, easy-to-implement, visualization tools, scientists' need for, and appreciation of, data visualization has risen too. Until recently, scientists thought of the "explanatory" graphics created at a research project's conclusion as "pretty pictures" needed only for journal publication or public outreach. The plots and displays produced during a research project - often intended only for experts - were thought of as a separate category, what we here call "exploratory" visualization. In this view, discovery comes from exploratory visualization, and explanatory visualization is just for communication. Our aim in this paper is to spark conversation amongst scientists, computer scientists, outreach professionals, educators, and graphics and perception experts about how to foster flexible data visualization practices that can facilitate discovery and communication at the same time. We present an example of a new finding made using the glue visualization environment to demonstrate how the border between explanatory and exploratory visualization is easily traversed. The linked-view principles as well as the actual code in glue are easily adapted to astronomy, medicine, and geographical information science - all fields where combining, visualizing, and analyzing several high-dimensional datasets yields insight. Whether or not scientists can use such a flexible "undisciplined" environment to its fullest potential without special training remains to be seen. We conclude with suggestions for improving the training of scientists in visualization practices, and of computer scientists in the iterative, non-workflow-like, ways in which modern science is carried out. △ Less

Submitted 29 May, 2018; originally announced May 2018.

Comments: Submitted as an invited "Perspectives" Paper for PNAS, in conjunction with the 2018 Sackler Colloquium

Observation of Anisotropy of TeV Cosmic Rays with Two Years of HAWC

Authors: A. U. Abeysekara, R. Alfaro, C. Alvarez, J. D. Alvarez, R. Arceo, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, A. Becerril, E. Belmont-Moreno, S. Y. BenZvi, A. Bernal, J. Braun, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, M. Castillo, U. Cotti, J. Cotzomi, C. De León, E. De la Fuente, R. Diaz Hernandez, S. Dichiara, B. L. Dingus , et al. (78 additional authors not shown)

Abstract: After two years of operation, the High-Altitude Water Cherenkov (HAWC) Observatory has analyzed the TeV cosmic-ray sky over an energy range between $2.0$ and $72.8$ TeV. The HAWC detector is a ground-based air-shower array located at high altitude in the state of Puebla, Mexico. Using 300 light-tight water tanks, it collects the Cherenkov light from the particles of extensive air showers from prim… ▽ More After two years of operation, the High-Altitude Water Cherenkov (HAWC) Observatory has analyzed the TeV cosmic-ray sky over an energy range between $2.0$ and $72.8$ TeV. The HAWC detector is a ground-based air-shower array located at high altitude in the state of Puebla, Mexico. Using 300 light-tight water tanks, it collects the Cherenkov light from the particles of extensive air showers from primary gamma rays and cosmic rays. This detection method allows for uninterrupted observation of the entire overhead sky (2~sr instantaneous, 8.5~sr integrated) in the energy range from a few TeV to hundreds of TeV. Like other detectors in the northern and southern hemisphere, HAWC observes an energy-dependent anisotropy in the arrival direction distribution of cosmic rays. The observed cosmic-ray anisotropy is dominated by a dipole moment with phase $α\approx40^{\circ}$ and amplitude that slowly rises in relative intensity from $8\times10^{-4}$ at 2 TeV to $14\times10^{-4}$ around 30.3 TeV, above which the dipole decreases in strength. A significant large-scale ($>60^{\circ}$ in angular extent) signal is also observed in the quadrupole and octupole moments, and significant small-scale features are also present, with locations and shapes consistent with previous observations. Compared to previous measurements in this energy range, the HAWC cosmic-ray sky maps improve on the energy resolution and fit precision of the anisotropy. These data can be used in an effort to better constrain local cosmic-ray accelerators and the intervening magnetic fields. △ Less

Submitted 20 July, 2018; v1 submitted 4 May, 2018; originally announced May 2018.

Comments: 22 pages, 14 figures, 2 tables, submission to ApJ

Relative Alignment Between the Magnetic Field and Molecular Gas Structure in the Vela C Giant Molecular Cloud using Low and High Density Tracers

Authors: Laura M. Fissel, Peter A. R. Ade, Francesco E. Angilè, Peter Ashton, Steven J. Benton, Che-Yu Chen, Maria Cunningham, Mark J. Devlin, Bradley Dober, Rachel Friesen, Yasuo Fukui, Nicholas Galitzki, Natalie N. Gandilo, Alyssa Goodman, Claire-Elise Green, Paul Jones, Jeffrey Klein, Patrick King, Andrei L. Korotkov, Zhi-Yun Li, Vicki Lowe, Peter G. Martin, Tristan G. Matthews, Lorenzo Moncelsi, Fumitaka Nakamura , et al. (15 additional authors not shown)

Abstract: We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$μ$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low den… ▽ More We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$μ$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud. △ Less

Submitted 2 April, 2019; v1 submitted 24 April, 2018; originally announced April 2018.

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

Search for Dark Matter Gamma-ray Emission from the Andromeda Galaxy with the High-Altitude Water Cherenkov Observatory

Authors: A. Albert, R. Alfaro, C. Alvarez, J. D. Alvarez, R. Arceo, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, A. Becerril, E. Belmont-Moreno, S. Y. BenZvi, A. Bernal, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, M. Castillo, U. Cotti, J. Cotzomi, S. Coutiño de León, C. De León, S. Dichiara, B. L. Dingus, M. A. DuVernois , et al. (76 additional authors not shown)

Abstract: The Andromeda Galaxy (M31) is a nearby ($\sim$780 kpc) galaxy similar to our own Milky Way. Observational evidence suggests that it resides in a large halo of dark matter (DM), making it a good target for DM searches. We present a search for gamma rays from M31 using 1017 days of data from the High Altitude Water Cherenkov (HAWC) Observatory. With its wide field of view and constant monitoring, HA… ▽ More The Andromeda Galaxy (M31) is a nearby ($\sim$780 kpc) galaxy similar to our own Milky Way. Observational evidence suggests that it resides in a large halo of dark matter (DM), making it a good target for DM searches. We present a search for gamma rays from M31 using 1017 days of data from the High Altitude Water Cherenkov (HAWC) Observatory. With its wide field of view and constant monitoring, HAWC is well-suited to search for DM in extended targets like M31. No DM annihilation or decay signal was detected for DM masses from 1 to 100 TeV in the $b\bar{b}$, $t\bar{t}$, $τ^{+}τ^{-}$, $μ^{+}μ^{-}$, and $W^{+}W^{-}$ channels. Therefore we present limits on those processes. Our limits nicely complement the existing body of DM limits from other targets and instruments. Specifically the DM decay limits from our benchmark model are the most constraining for DM masses from 25 TeV to 100 TeV in the $b\bar{b}, t\bar{t}$ and $μ^{+}μ{-}$ channels. In addition to DM-specific limits, we also calculate general gamma-ray flux limits for M31 in 5 energy bins from 1 TeV to 100 TeV. △ Less

Submitted 13 March, 2019; v1 submitted 2 April, 2018; originally announced April 2018.

Comments: published in JCAP. Figure 6 (bottom) has been corrected

The CARMA-NRO Orion Survey

Authors: Shuo Kong, Héctor G. Arce, Jesse R. Feddersen, John M. Carpenter, Fumitaka Nakamura, Yoshito Shimajiri, Andrea Isella, Volker Ossenkopf-Okada, Anneila I. Sargent, Álvaro Sánchez-Monge, Sümeyye T. Suri, Jens Kauffmann, Thushara Pillai, Jaime E. Pineda, Jin Koda, John Bally, Dariusz C. Lis, Paolo Padoan, Ralf Klessen, Steve Mairs, Alyssa Goodman, Paul Goldsmith, Peregrine McGehee, Peter Schilke, Peter J. Teuben , et al. (13 additional authors not shown)

Abstract: We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in… ▽ More We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $Δ$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion △ Less

Submitted 30 March, 2018; originally announced March 2018.

Comments: 46 pages, 28 figures, 2 tables, accepted by ApJS

Mapping Distances Across the Perseus Molecular Cloud Using CO Observations, Stellar Photometry, and Gaia DR2 Parallax Measurements

Authors: Catherine Zucker, Edward F. Schlafly, Joshua S. Speagle, Gregory M. Green, Stephen K. N. Portillo, Douglas P. Finkbeiner, Alyssa A. Goodman

Abstract: We present a new technique to determine distances to major star-forming regions across the Perseus Molecular Cloud, using a combination of stellar photometry, astrometric data, and $\rm ^{12} CO$ spectral-line maps. Incorporating the Gaia DR2 parallax measurements when available, we start by inferring the distance and reddening to stars from their Pan-STARRS1 and 2MASS photometry, based on a techn… ▽ More We present a new technique to determine distances to major star-forming regions across the Perseus Molecular Cloud, using a combination of stellar photometry, astrometric data, and $\rm ^{12} CO$ spectral-line maps. Incorporating the Gaia DR2 parallax measurements when available, we start by inferring the distance and reddening to stars from their Pan-STARRS1 and 2MASS photometry, based on a technique presented in Green et al. 2014; Green et al. 2015 and implemented in their 3D "Bayestar" dust map of three-quarters of the sky. We then refine the Green et al. technique by using the velocity slices of a CO spectral cube as dust templates and modeling the cumulative distribution of dust along the line of sight towards these stars as a linear combination of the emission in the slices. Using a nested sampling algorithm, we fit these per-star distance-reddening measurements to find the distances to the CO velocity slices towards each star-forming region. This results in distance estimates explicitly tied to the velocity structure of the molecular gas. We determine distances to the B5, IC348, B1, NGC1333, L1448, and L1451 star-forming regions and find that individual clouds are located between $\approx 275-300$ pc, with typical combined uncertainties of $\approx 5\%$. We find that the velocity gradient across Perseus corresponds to a distance gradient of about 25 pc, with the eastern portion of the cloud farther away than the western portion. We determine an average distance to the complex of $294\pm 17$ pc, about 60 pc higher than the distance derived to the western portion of the cloud using parallax measurements of water masers associated with young stellar objects. The method we present is not limited to the Perseus Complex, but may be applied anywhere on the sky with adequate CO data in the pursuit of more accurate 3D maps of molecular clouds in the solar neighborhood and beyond. △ Less

Submitted 17 October, 2018; v1 submitted 23 March, 2018; originally announced March 2018.

Comments: Accepted for publication in The Astrophysical Journal

Constraining the $\bar{p}/p$ Ratio in TeV Cosmic Rays with Observations of the Moon Shadow by HAWC

Authors: A. U. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, R. Arceo, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, E. Belmont-Moreno, S. Y. BenZvi, J. Braun, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, M. Castillo, U. Cotti, J. Cotzomi, S. Coutiño de León, C. De León, E. D la Fuentem, R. Diaz Hernandez, S. Dichiara, B. L. Dingus , et al. (76 additional authors not shown)

Abstract: An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed towards the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a func… ▽ More An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed towards the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a function of the energy and charge of the cosmic rays. The displacement provides a natural tool for momentum/charge discrimination that can be used to study the composition of cosmic rays. Using 33 months of data comprising more than 80 billion cosmic rays measured by the High Altitude Water Cherenkov (HAWC) observatory, we have analyzed the Moon shadow to search for TeV antiprotons in cosmic rays. We present our first upper limits on the $\bar{p}/p$ fraction, which in the absence of any direct measurements, provide the tightest available constraints of $\sim1\%$ on the antiproton fraction for energies between 1 and 10 TeV. △ Less

Submitted 22 April, 2018; v1 submitted 24 February, 2018; originally announced February 2018.

Comments: 10 pages, 5 figures. Accepted by Physical Review D

Journal ref: Phys. Rev. D 97, 102005 (2018)

Learning multiagent coordination in the absence of communication channels

Authors: Aaron Goodman

Abstract: In this work, we develop a reinforcement learning protocol for a multiagent coordination task in a discrete state and action space: an iterated prisoner's dilemma game extended into a team based, winner-take all tournament, which forces the agents to collude in order to maximize their reward. By disallowing extra communication channels, the agents are forced to embed their coordination strategy in… ▽ More In this work, we develop a reinforcement learning protocol for a multiagent coordination task in a discrete state and action space: an iterated prisoner's dilemma game extended into a team based, winner-take all tournament, which forces the agents to collude in order to maximize their reward. By disallowing extra communication channels, the agents are forced to embed their coordination strategy into their actions in the prisoner's dilemma game. We develop a representation of the iterated prisoners dilemma that makes it amenable to Q-learning. We find that the reinforcement learning strategy is able to consistently train agents that can win the winner take all iterated prisoners dilemma tournament. By using a game with discrete state and action space, we are able to better analyze and understand both the dynamics and the communication protocols that are established between the agents. We find that the agents adapt a number of interesting behaviors, such as the formation of benevolent dictators, that minimize inequality of scores. We also find that the agents settle on a remarkably consistent symbology in their actions, such that agents from independent trials are able to collude with each other without further training. △ Less

Submitted 14 June, 2018; v1 submitted 16 February, 2018; originally announced February 2018.

Estimating city-level travel patterns using street imagery: a case study of using Google Street View in Britain

Authors: Rahul Goel, Leandro M. T. Garcia, Anna Goodman, Rob Johnson, Rachel Aldred, Manoradhan Murugesan, Soren Brage, Kavi Bhalla, James Woodcock

Abstract: Street imagery is a promising big data source providing current and historical images in more than 100 countries. Previous studies used this data to audit built environment features. Here we explore a novel application, using Google Street View (GSV) to predict travel patterns at the city level. We sampled 34 cities in Great Britain. In each city, we accessed GSV images from 1000 random locations… ▽ More Street imagery is a promising big data source providing current and historical images in more than 100 countries. Previous studies used this data to audit built environment features. Here we explore a novel application, using Google Street View (GSV) to predict travel patterns at the city level. We sampled 34 cities in Great Britain. In each city, we accessed GSV images from 1000 random locations from years overlapping with the 2011 Census and the 2011-2013 Active People Survey (APS). We manually annotated images into seven categories of road users. We developed regression models with the counts of images of road users as predictors. Outcomes included Census-reported commute shares of four modes (walking plus public transport, cycling, motorcycle, and car), and APS-reported past-month participation in walking and cycling. In bivariate analyses, we found high correlations between GSV counts of cyclists (GSV-cyclists) and cycle commute mode share (r=0.92) and past-month cycling (r=0.90). Likewise, GSV-pedestrians was moderately correlated with past-month walking for transport (r=0.46), GSV-motorcycles was moderately correlated with commute share of motorcycles (r=0.44), and GSV-buses was highly correlated with commute share of walking plus public transport (r=0.81). GSV-car was not correlated with car commute mode share (r=-0.12). However, in multivariable regression models, all mode shares were predicted well. Cross-validation analyses showed good prediction performance for all the outcomes except past-month walking. Street imagery is a promising new big data source to predict urban mobility patterns. Further testing across multiple settings is warranted both for cross-sectional and longitudinal assessments. △ Less

Submitted 8 February, 2018; originally announced February 2018.

Comments: Paper submitted for peer review. 7 figures. 3 Tables