Showing 1–50 of 161 results for author: Smith, J D

Isotropy of cosmic rays beyond $10^{20}$ eV favors their heavy mass composition

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the resul… ▽ More We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the results are consistent with a relatively heavy injected composition at E ~ 10 EeV that becomes lighter up to E ~ 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extra-galactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 8 pages, 3 figures, accepted for publication in PRL

Mass composition of ultra-high energy cosmic rays from distribution of their arrival directions with the Telescope Array

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale struc… ▽ More We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extra-galactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields (GMF). The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E > 100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as $2 \times 10^{-5}$ Mpc$^{-3}$, that is the conservative lower limit for the source number density. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

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

Observation of Declination Dependence in the Cosmic Ray Energy Spectrum

Authors: The Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, J. W. Belz, D. R. Bergman, I. Buckland, W. Campbell, B. G. Cheon, K. Endo, A. Fedynitch, T. Fujii, K. Fujisue, K. Fujita, M. Fukushima, G. Furlich, Z. Gerber, N. Globus, W. Hanlon, N. Hayashida, H. He, K. Hibino, R. Higuchi, D. Ikeda, T. Ishii , et al. (101 additional authors not shown)

Abstract: We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements fr… ▽ More We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements from different observatories introduces the issue of possible systematic differences between detectors and analyses, we validate the methodology of the comparison by examining the region of the sky where the apertures of the two observatories overlap. Although the spectra differ in this region, we find that there is only a $1.8σ$ difference between the spectrum measurements when anisotropic regions are removed and a fiducial cut in the aperture is applied. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 8 pages, 6 figures

JWST Observations of Starbursts: Cold Clouds and Plumes Launching in the M82 Outflow

Authors: Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Serena A. Cronin, Laura A. Lopez, J. D. Smith, Danielle A. Berg, Sebastian Lopez, Sylvain Veilleux, Paul P. van der Werf, Torsten Böker, Leindert A. Boogaard, Laura Lenkić, Simon C. O. Glover, Vicente Villanueva, Divakara Mayya, Thomas S. -Y. Lai, Daniel A. Dale, Kimberly L. Emig, Fabian Walter, Monica Relaño , et al. (6 additional authors not shown)

Abstract: In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority (… ▽ More In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$α$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$α$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

Disentangling the co-evolution of galaxies and supermassive black holes with PRIMA

Authors: L. Bisigello, C. Gruppioni, A. Bolatto, L. Ciesla, A. Pope, L. Armus, L., J. D. Smith, R. Somerville, L. Y. A. Yung, R. J. Wright, C. M. Bradford, J. Glenn, A. Feltre

Abstract: The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-I… ▽ More The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-IR provide unmatched diagnostic power to determine the properties of gas and dust, measure gas-phase metallicities and map cold galactic outflows in even the most obscured galaxies. By combining multi-band photometric surveys with low and high-resolution far-IR spectroscopy, the PRobe far-Infrared Mission for Astrophysics (PRIMA), a concept for a far-IR, 1.8m-diameter, cryogenically cooled observatory, will revolutionise the field of galaxy evolution by taking advantage of this IR toolkit to find and study dusty galaxies across galactic time. In this work, we make use of the phenomenological simulation SPRITZ and the Santa Cruz semi-analytical model to describe how a moderately deep multi-band PRIMA photometric survey can easily reach beyond previous IR missions to detect and study galaxies down to $10^{11}\,L_{\odot}$ beyond cosmic noon and at least up to z=4, even in the absence of gravitational lensing. By decomposing the spectral energy distribution (SED) of these photometrically selected galaxies, we show that PRIMA can be used to accurately measure the relative AGN power, the mass fraction contributed by polycyclic aromatic hydrocarbon (PAH) and the total IR luminosity. At the same time, spectroscopic follow up with PRIMA will allow to trace both the star formation and black hole accretion rates (SFR, BHAR), the gas phase metallicities and the mass outflow rates of cold gas in hundreds to thousands of individual galaxies to z=2. △ Less

Submitted 3 June, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: 16 pages, 17 figures, accepted for publication in A&A

Spectroscopic Constraints on the Mid-Infrared Attenuation Curve: I -- Attenuation Model using PAH Emission

Authors: Thomas S. -Y. Lai, J. D. T. Smith, Els Peeters, Henrik W. W. Spoon, Shunsuke Baba, Masatoshi Imanishi, Takao Nakagawa

Abstract: We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9… ▽ More We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9.7 um and 18 um silicate absorption features make MIR spectral modeling and interpretation challenging, particularly for heavily obscured galaxies. By adopting the basic premise that PAH bands have relatively consistent intrinsic ratios within dusty starbursting galaxies, we can, for the first time, empirically determine the overall shape of the MIR attenuation curve by measuring the differential attenuation at specific PAH wavelengths. Our attenuation model shows PAH emission in most (U)LIRGs is unambiguously subjected to attenuation, and we find strong evidence that PAH bands undergo differential attenuation as obscuration increases. Compared to pre-existing results, the MIR attenuation curve derived from the model favors relatively gray continuum absorption from 3-8 $μ$m and silicate features with intermediate strength at 9.7 um but with stronger than typical 18 um opacity. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 15 pages, 5 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

The Impact of an AGN on PAH Emission in Galaxies: the Case of Ring Galaxy NGC 4138

Authors: G. P. Donnelly, J. D. T. Smith, B. T. Draine, A. Togi, T. S. -Y. Lai, L. Armus, D. A. Dale, V. Charmandaris

Abstract: We present a focused study of radially-resolved varying PAH emission in the low-luminosity AGN-host NGC 4138 using deep Spitzer/IRS spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138,… ▽ More We present a focused study of radially-resolved varying PAH emission in the low-luminosity AGN-host NGC 4138 using deep Spitzer/IRS spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138, and we use this model to predict the corresponding changes in the PAH emission spectrum. We find that PAH band ratios are strongly variable with radius in this galaxy with short-to-long wavelength band ratios peaking in the starburst ring. The changing mix of starlight appears to have a considerable effect on the trends in these band ratios, and our radiation model predicts the shapes of these trends. However, the amplitude of observed variation is ~2.5 times larger than predicted for some ratios. A cutoff of small grains in the PAH size distribution, as has been suggested for AGN, together with changes in PAH ionization fraction could explain the behavior of the shorter bands, but this model fails to reproduce longer band behaviors. Additionally, we find that short-to-long wavelength PAH band ratios increase slightly within ~270pc of the center, suggesting that the AGN may directly influence PAH emission there. △ Less

Submitted 23 April, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 19 pages, 13 figures, published in ApJ. Updated values in radiation model to the same as the published version

JWST Observations of Starbursts: Polycyclic Aromatic Hydrocarbon Emission at the Base of the M 82 Galactic Wind

Authors: Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, Simon C. O. Glover, Karin M. Sandstrom, Evan D. Skillman, John Chisholm, Vicente Villanueva , et al. (15 additional authors not shown)

Abstract: We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find promine… ▽ More We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas. △ Less

Submitted 21 April, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

TEMPLATES: Direct Abundance Constraints for Two Lensed Lyman-Break Galaxies

Authors: Brian Welch, Grace M. Olivier, Taylor A. Hutchison, Jane R. Rigby, Danielle A. Berg, Manuel Aravena, Matthew B. Bayliss, Jack E. Birkin, Scott C. Chapman, Håkon Dahle, Gourav Khullar, Keunho J. Kim, Guillaume Mahler, Matthew A. Malkan, Desika Narayanan, Kedar A. Phadke, Keren Sharon, J. D. T. Smith, Manuel Solimano, Justin S. Spilker, Joaquin D. Viera, David Vizgan

Abstract: Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$λλ$7320,7330 auroral emission lines, and set an upper li… ▽ More Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$λλ$7320,7330 auroral emission lines, and set an upper limit for the [NII]$\lambda5755$ line. For the second galaxy, SGAS1226$+$21 ($z = 2.925$), we do not detect any auroral lines, and report upper limits. With these measurements and upper limits, we constrain the electron temperatures in different ionization zones within both of these galaxies. For SGAS1723$+$34, where auroral lines are detected, we calculate direct oxygen and nitrogen abundances, finding an N/O ratio consistent with observations of nearby ($z\sim 0$) galaxies. These observations highlight the potent combination of JWST and gravitational lensing to measure faint emission lines in individual distant galaxies and to directly study the chemical abundance patterns in those galaxies. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Submitted to the Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program. 15 pages, 6 figures. Comments welcome

Radio Jet Feedback on the Inner Disk of Virgo Spiral Galaxy Messier 58

Authors: Patrick M. Ogle, Ivan E. Lopez, Victoria Reynaldi, Aditya Togi, R. Michael Rich, Javier Roman, Osmin Caceres, Zhuofu, Li, Grant Donnelly, J. D. T. Smith, Philip N. Appleton, Lauranne Lanz

Abstract: Spitzer spectral maps reveal a disk of highly luminous, warm (>150 K) H2 in the center of the massive spiral galaxy Messier 58, which hosts a radio-loud AGN. The inner 2.6 kpc of the galaxy appears to be overrun by shocks from the radio jet cocoon. Gemini NIRI imaging of the H2 1-0 S(1) emission line, ALMA CO 2-1, and HST multiband imagery indicate that much of the molecular gas is shocked in-situ… ▽ More Spitzer spectral maps reveal a disk of highly luminous, warm (>150 K) H2 in the center of the massive spiral galaxy Messier 58, which hosts a radio-loud AGN. The inner 2.6 kpc of the galaxy appears to be overrun by shocks from the radio jet cocoon. Gemini NIRI imaging of the H2 1-0 S(1) emission line, ALMA CO 2-1, and HST multiband imagery indicate that much of the molecular gas is shocked in-situ, corresponding to lanes of dusty molecular gas that spiral towards the galaxy nucleus. The CO 2-1 and ionized gas kinematics are highly disturbed, with velocity dispersion up to 300 km/s. Dissipation of the associated kinetic energy and turbulence, likely injected into the ISM by radio-jet driven outflows, may power the observed molecular and ionized gas emission from the inner disk. The PAH fraction and composition in the inner disk appear to be normal, in spite of the jet and AGN activity. The PAH ratios are consistent with excitation by the interstellar radiation field from old stars in the bulge, with no contribution from star formation. The phenomenon of jet-shocked H2 may substantially reduce star formation and help to regulate the stellar mass of the inner disk and supermassive black hole in this otherwise normal spiral galaxy. Similarly strong H2 emission is found at the centers of several nearby spiral and lenticular galaxies with massive bulges and radio-loud AGN. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 18 pages, 12 figures, 1 Table, accepted for publication in the Astrophysical Journal

Can massive stars form in low mass clouds?

Authors: Jamie D. Smith, Sarah E. Jaffa, Martin G. H. Krause

Abstract: The conditions required for massive star formation are debated, particularly whether massive stars must form in conjunction with massive clusters. Some authors have advanced the view that stars of any mass (below the total cluster mass) can form in clusters of any mass with some probability (random sampling). Others pointed out that the scatter in the determinations of the most massive star mass f… ▽ More The conditions required for massive star formation are debated, particularly whether massive stars must form in conjunction with massive clusters. Some authors have advanced the view that stars of any mass (below the total cluster mass) can form in clusters of any mass with some probability (random sampling). Others pointed out that the scatter in the determinations of the most massive star mass for a given cluster mass was consistent with the measurement error, such that the mass of the most massive star was determined by the total cluster mass (optimal sampling). Here we investigate the relation between cluster mass (M\textsubscript{ecl}) and the maximum stellar mass (M\textsubscript{max}) using a suite of SPH simulations. Varying cloud mass and turbulence random seed results in a range of cluster masses which we compare with their respective maximum star masses. We find that more massive clusters will have, on average, higher mass stars with this trend being steeper at lower cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.31}$ for $M\textsubscript{ecl}<500M\,_{\odot}$) and flattening at higher cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.11}$ for $M\textsubscript{ecl}>500M\,_{\odot}$). This rules out purely stochastic star formation in our simulations. Significant scatter in the maximum masses with identical initial conditions also rules out the possibility that the relation is purely deterministic (that is that a given cluster mass will result in a specific maximum stellar mass). In conclusion our simulations disagree with both random and optimal sampling of the initial mass function. △ Less

Submitted 28 September, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted in MNRAS

JWST's TEMPLATES for Star Formation: The First Resolved Gas-Phase Metallicity Maps of Dust-Obscured Star-Forming Galaxies at $z$ $\sim$ 4

Authors: Jack E. Birkin, Taylor A. Hutchison, Brian Welch, Justin S. Spilker, Manuel Aravena, Matthew B. Bayliss, Jared Cathey, Scott C. Chapman, Anthony H. Gonzalez, Gayathri Gururajan, Christopher C. Hayward, Gourav Khullar, Keunho J. Kim, Guillaume Mahler, Matthew A. Malkan, Desika Narayanan, Grace M. Olivier, Kedar A. Phadke, Cassie Reuter, Jane R. Rigby, J. D. T. Smith, Manuel Solimano, Nikolaus Sulzenauer, Joaquin D. Vieira, David Vizgan , et al. (1 additional authors not shown)

Abstract: We present the first spatially resolved maps of gas-phase metallicity for two dust-obscured star-forming galaxies (DSFGs) at $z\sim$ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H$α$ and [NII] emission lines. Empirical optical line calibrations are used to determine that the sources are globally enriched to near-solar levels… ▽ More We present the first spatially resolved maps of gas-phase metallicity for two dust-obscured star-forming galaxies (DSFGs) at $z\sim$ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H$α$ and [NII] emission lines. Empirical optical line calibrations are used to determine that the sources are globally enriched to near-solar levels. While one source shows elevated [NII]/H$α$ ratios and broad H$α$ emission consistent with the presence of an AGN in a $\gtrsim$1kpc region, we argue that both systems have already undergone significant metal enrichment as a result of their extremely high star formation rates. Utilizing ALMA rest-frame 380$μ$m continuum and [CI]($^3$P$_2$-$^3$P$_1$) line maps we compare the spatial variation of the metallicity and gas-to-dust ratio in the two galaxies, finding the two properties to be anticorrelated on highly resolved spatial scales, consistent with various literature studies of $z\sim$ 0 galaxies. The data are indicative of the enormous potential of JWST to probe the enrichment of the interstellar medium on $\sim$kpc scales in extremely dust-obscured systems at $z\sim$ 4 and beyond. △ Less

Submitted 28 August, 2023; v1 submitted 19 July, 2023; originally announced July 2023.

Comments: 14 pages, 6 figures, accepted for publication in ApJ

Spatial variations in aromatic hydrocarbon emission in a dust-rich galaxy

Authors: Justin S. Spilker, Kedar A. Phadke, Manuel Aravena, Melanie Archipley, Matthew B. Bayliss, Jack E. Birkin, Matthieu Bethermin, James Burgoyne, Jared Cathey, Scott C. Chapman, Hakon Dahle, Anthony H. Gonzalez, Gayathri Gururajan, Christopher C. Hayward, Yashar D. Hezaveh, Ryley Hill, Taylor A. Hutchison, Keunho J. Kim, Seonwoo Kim, David Law, Ronan Legin, Matthew A. Malkan, Daniel P. Marrone, Eric J. Murphy, Desika Narayanan , et al. (13 additional authors not shown)

Abstract: Dust grains absorb half of the radiation emitted by stars throughout the history of the universe, re-emitting this energy at infrared wavelengths. Polycyclic aromatic hydrocarbons (PAHs) are large organic molecules that trace millimeter-size dust grains and regulate the cooling of the interstellar gas within galaxies. Observations of PAH features in very distant galaxies have been difficult due to… ▽ More Dust grains absorb half of the radiation emitted by stars throughout the history of the universe, re-emitting this energy at infrared wavelengths. Polycyclic aromatic hydrocarbons (PAHs) are large organic molecules that trace millimeter-size dust grains and regulate the cooling of the interstellar gas within galaxies. Observations of PAH features in very distant galaxies have been difficult due to the limited sensitivity and wavelength coverage of previous infrared telescopes. Here we present JWST observations that detect the 3.3um PAH feature in a galaxy observed less than 1.5 billion years after the Big Bang. The high equivalent width of the PAH feature indicates that star formation, rather than black hole accretion, dominates the infrared emission throughout the galaxy. The light from PAH molecules, large dust grains, and stars and hot dust are spatially distinct from one another, leading to order-of-magnitude variations in the PAH equivalent width and the ratio of PAH to total infrared luminosity across the galaxy. The spatial variations we observe suggest either a physical offset between the PAHs and large dust grains or wide variations in the local ultraviolet radiation field. Our observations demonstrate that differences in the emission from PAH molecules and large dust grains are a complex result of localized processes within early galaxies. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: Published in Nature 5 June 2023 at https://www.nature.com/articles/s41586-023-05998-6. MIRI MRS reduction notebook is available at https://github.com/jwst-templates

The Physical Drivers and Observational Tracers of CO-to-H2 Conversion Factor Variations in Nearby Barred Galaxy Centers

Authors: Yu-Hsuan Teng, Karin M. Sandstrom, Jiayi Sun, Munan Gong, Alberto D. Bolatto, I-Da Chiang, Adam K. Leroy, Antonio Usero, Simon C. O. Glover, Ralf S. Klessen, Daizhong Liu, Miguel Querejeta, Eva Schinnerer, Frank Bigiel, Yixian Cao, Melanie Chevance, Cosima Eibensteiner, Kathryn Grasha, Frank P. Israel, Eric J. Murphy, Lukas Neumann, Hsi-An Pan, Francesca Pinna, Mattia C. Sormani, J. D. T. Smith , et al. (2 additional authors not shown)

Abstract: The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $α_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NG… ▽ More The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $α_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NGC 3627 and NGC 4321 tracing $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales. Our multi-line modeling and Bayesian likelihood analysis of these datasets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of $α_\rm{CO}$. The central 300 pc nuclei in both galaxies show strong enhancement of temperature $T_\rm{k}>100$ K and density $n_\rm{H_2}>10^3$ cm$^{-3}$. Assuming a CO-to-H$_2$ abundance of $3\times10^{-4}$, we derive 4-15 times lower $α_\rm{CO}$ than the Galactic value across our maps, which agrees well with previous kpc-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of $α_\rm{CO}$ with low-J $^{12}$CO optical depths ($τ_\rm{CO}$), as well as an anti-correlation with $T_\rm{k}$. The $τ_\rm{CO}$ correlation explains most of the $α_\rm{CO}$ variation in the three galaxy centers, whereas changes in $T_\rm{k}$ influence $α_\rm{CO}$ to second order. Overall, the observed line width and $^{12}$CO/$^{13}$CO 2-1 line ratio correlate with $τ_\rm{CO}$ variation in these centers, and thus they are useful observational indicators for $α_\rm{CO}$ variation. We also test current simulation-based $α_\rm{CO}$ prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted for publication in ApJ; 30 pages of main text + 3 appendices

A Framework for Modeling Polycyclic Aromatic Hydrocarbon Emission in Galaxy Evolution Simulations

Authors: Desika Narayanan, J. D. Smith, Brandon Hensley, Qi Li, Chia-Yu Hu, Karin Sandstrom, Paul Torrey, Mark Vogelsberger, Federico Marinacci, Laura Sales

Abstract: We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain size distributions, and ionization states with a new on-the-fly model for dust evolution in hydrodynamic galaxy simulations.… ▽ More We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain size distributions, and ionization states with a new on-the-fly model for dust evolution in hydrodynamic galaxy simulations. We apply these models to 3 idealized arepo galaxy evolution simulations within the smuggle physics framework. We use these simulations to develop numerical experiments investigating the buildup of PAH masses and luminosities in galaxies in idealized analogs of the Milky Way, a dwarf galaxy, and starburst disk. Our main results follow. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to efficiently shatter dust grains, driving up the fraction of ultra small grains. At the same time, in our model large radiation fields per unit gas density convert aliphatic grains into aromatics. The fraction of dust grains in the form of PAHs (q_PAH) can be understood as a consequence of these processes, and in our model PAHs form primarily from interstellar processing (shattering) of larger grains rather than from the growth of smaller grains. We find that the hardness of the radiation field plays a larger role than variations in the grain size distribution in setting the total integrated PAH luminosities, though cosmological simulations are necessary to fully investigate the complex interplay of processes that drive PAH band luminosities in galaxies. Finally, we highlight feature PAH strength variations, cautioning against the usage of emission templates with constant feature strength ratios. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: Submitted to ApJ; comments welcome

GOALS-JWST: Tracing AGN Feedback on the Star-Forming ISM in NGC 7469

Authors: Thomas S. -Y. Lai, Lee Armus, Vivian U, Tanio Diaz-Santos, Kirsten L. Larson, Aaron Evans, Matthew A. Malkan, Philip Appleton, Jeff Rich, Francisco Muller-Sanchez, Hanae Inami, Thomas Bohn, Jed McKinney, Luke Finnerty, David R. Law, Sean Linden, Anne M. Medling, George C. Privon, Yiqing Song, Sabrina Stierwalt, Paul P. van der Werf, Loreto Barcos-Muñoz, J. D. T. Smith, Aditya Togi, Susanne Aalto , et al. (12 additional authors not shown)

Abstract: We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study AGN feedback in the local Universe. We take advantage of the high spatial/spectral… ▽ More We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study AGN feedback in the local Universe. We take advantage of the high spatial/spectral resolution of JWST/MIRI to isolate the star-forming regions surrounding the central active nucleus and study the properties of the dust and warm molecular gas on ~100 pc scales. The starburst ring exhibits prominent Polycyclic Aromatic Hydrocarbon (PAH) emission, with grain sizes and ionization states varying by only ~30%, and a total star formation rate of $\rm 10 - 30 \ M_\odot$/yr derived from fine structure and recombination emission lines. Using pure rotational lines of H2, we detect 1.2$\times$10$^{7} \rm \ M_\odot$ of warm molecular gas at a temperature higher than 200 K in the ring. All PAH bands get significantly weaker towards the central source, where larger and possibly more ionized grains dominate the emission. However, the bulk of the dust and molecular gas in the ring appears unaffected by the ionizing radiation or the outflowing wind from the AGN. These observations highlight the power of JWST to probe the inner regions of dusty, rapidly evolving galaxies for signatures of feedback and inform models that seek to explain the co-evolution of supermassive black holes and their hosts. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 13 pages, 4 figures, 2 tables, Submitted to ApJL

Star Cluster Formation in Clouds with Externally Driven Turbulence

Authors: Jamie D. Smith, James E. Dale, Sarah E. Jaffa, Martin G. H. Krause

Abstract: Star clusters are known to be formed in turbulent molecular clouds. How turbulence is driven in molecular clouds and what effect this has on star formation is still unclear. We compare a simulation setup with turbulent driving everywhere in a periodic box with a setup where turbulence is only driven around the outside of the box. We analyse the resulting gas distribution, kinematics, and the popul… ▽ More Star clusters are known to be formed in turbulent molecular clouds. How turbulence is driven in molecular clouds and what effect this has on star formation is still unclear. We compare a simulation setup with turbulent driving everywhere in a periodic box with a setup where turbulence is only driven around the outside of the box. We analyse the resulting gas distribution, kinematics, and the population of stars that are formed from the cloud. Both setups successfully produce a turbulent velocity field with a power law structure function, the externally driven cloud has a more central, monolithic, clump, while the fully driven cloud has many smaller, more dispersed, clumps. The star formation follows the cloud morphology producing large clusters, with high star forming efficiency in the externally driven simulations and sparse individual star formation with much lower star formation efficiency in the fully driven case. We conclude that the externally driven method, which resembles a Global Hierarchical Collapse (GHC) scenario, produces star clusters that more closely match with observations. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: 8 pages, 6 figures, accepted by MNRAS

First High-speed Video Camera Observations of a Lightning Flash Associated with a Downward Terrestrial Gamma-ray Flash

Authors: R. U. Abbasi, M. M. F. Saba, J. W. Belz, P. R. Krehbiel, W. Rison, N. Kieu, D. R. da Silva, Dan Rodeheffer, M. A. Stanley, J. Remington, J. Mazich, R. LeVon, K. Smout, A. Petrizze, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii , et al. (127 additional authors not shown)

Abstract: In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field… ▽ More In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field fast antenna, and the National Lightning Detection Network. The cloud-to-ground flash associated with the observed TGF was formed by a fast downward leader followed by a very intense return stroke peak current of -154 kA. The TGF occurred while the downward leader was below cloud base, and even when it was halfway in its propagation to ground. The suite of gamma-ray and lightning instruments, timing resolution, and source proximity offer us detailed information and therefore a unique look at the TGF phenomena. △ Less

Submitted 9 August, 2023; v1 submitted 10 May, 2022; originally announced May 2022.

Journal ref: Geophysical Research Letters, 50, e2023GL102958 (2023)

The State of the Molecular Gas in Post-Starburst Galaxies

Authors: K. Decker French, Adam Smercina, Kate Rowlands, Akshat Tripathi, Ann I. Zabludoff, J. D. Smith, Desika Narayanan, Yujin Yang, Yancy Shirley, Katey Alatalo

Abstract: The molecular gas in galaxies traces both the fuel for star formation and the processes that can enhance or suppress star formation. Observations of the molecular gas state can thus point to when and why galaxies stop forming stars. In this study, we present ALMA observations of the molecular gas in galaxies evolving through the post-starburst phase. These galaxies have low current star formation… ▽ More The molecular gas in galaxies traces both the fuel for star formation and the processes that can enhance or suppress star formation. Observations of the molecular gas state can thus point to when and why galaxies stop forming stars. In this study, we present ALMA observations of the molecular gas in galaxies evolving through the post-starburst phase. These galaxies have low current star formation rates, regardless of the SFR tracer used, with recent starbursts ending within the last 600 Myr. We present CO (3-2) observations for three post-starburst galaxies, and dense gas HCN/HCO+/HNC (1-0) observations for six (four new) post-starburst galaxies. The post-starbursts have low excitation traced by the CO spectral line energy distribution (SLED) up to CO (3-2), more similar to early-type than starburst galaxies. The low excitation indicates that lower density rather than high temperatures may suppress star formation during the post-starburst phase. One galaxy displays a blueshifted outflow traced by CO (3-2). MaNGA observations show that the ionized gas velocity is disturbed relative to the stellar velocity field, with a blueshifted component aligned with the molecular gas outflow, suggestive of a multiphase outflow. Low ratios of HCO+/CO, indicating low fractions of dense molecular gas relative to the total molecular gas, are seen throughout post-starburst phase, except for the youngest post-starburst galaxy considered here. These observations indicate that the impact of any feedback or quenching processes may be limited to low excitation and weak outflows in the cold molecular gas during the post-starburst phase. △ Less

Submitted 12 December, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Comments: 33 pages, 24 figures, 4 tables, ApJ in press

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

Authors: The ANTARES collaboration, A. Albert, S. Alves, M. André, M. Anghinolfi, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, B. Belhorma, M. Bendahman, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzaş, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (1025 additional authors not shown)

Abstract: For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for corre… ▽ More For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs. △ Less

Submitted 23 August, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 39 pages, 7 figures, 4 tables; updated source files including xml authorlist

Report number: FERMILAB-PUB-22-033-AD-PPD-SCD-TD

Journal ref: ApJ 934 164 (2022)

Direct Far-Infrared Metal Abundances (FIRA) I: M101

Authors: C. Lamarche, J. D. Smith, K. Kreckel, S. T. Linden, N. S. J. Rogers, E. Skillman, D. Berg, E. Murphy, R. Pogge, G. P. Donnelly, R. Kennicutt Jr., A. Bolatto, K. Croxall, B. Groves, C. Ferkinhoff

Abstract: Accurately determining gas-phase metal-abundances within galaxies is critical as metals strongly affect the physics of the interstellar medium (ISM). To date, the vast majority of widely-used gas-phase abundance-indicators rely on emission from bright optical-lines, whose emissivities are highly sensitive to the electron temperature. Alternatively, direct-abundance methods exist that measure the t… ▽ More Accurately determining gas-phase metal-abundances within galaxies is critical as metals strongly affect the physics of the interstellar medium (ISM). To date, the vast majority of widely-used gas-phase abundance-indicators rely on emission from bright optical-lines, whose emissivities are highly sensitive to the electron temperature. Alternatively, direct-abundance methods exist that measure the temperature of the emitting gas directly, though these methods usually require challenging observations of highly-excited auroral lines. Low-lying far-infrared (FIR) fine-structure lines are largely insensitive to electron temperature and thus provide an attractive alternative to optically-derived abundances. Here, we introduce the far-infrared abundances (FIRA) project, which employs these FIR transitions, together with both radio free-free emission and hydrogen recombination-lines, to derive direct, absolute gas-phase oxygen-abundances. Our first target is M101, a nearby spiral-galaxy with a relatively steep abundance gradient. Our results are consistent with the O$^{++}$ electron-temperatures and absolute oxygen-abundances derived using optical direct-abundance methods by the CHemical Abundance Of Spirals (CHAOS) program, with a small difference ($\sim$ 1.5$σ$) in the radial abundance-gradients derived by the FIR/free-free-normalized vs. CHAOS/direct-abundance techniques. This initial result demonstrates the validity of the FIRA methodology $-$ with the promise of determining absolute metal-abundances within dusty star-forming galaxies, both locally and at high redshift. △ Less

Submitted 6 December, 2021; v1 submitted 19 November, 2021; originally announced November 2021.

Comments: 28 pages, 11 figures, 5 tables, accepted for publication in the Astrophysical Journal

Observation of Variations in Cosmic Ray Single Count Rates During Thunderstorms and Implications for Large-Scale Electric Field Changes

Authors: R. U. Abbasi, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, M. Hayashi , et al. (140 additional authors not shown)

Abstract: We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of t… ▽ More We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of the 700 km$^{2}$ detector, without dealing with the limitation of narrow exposure in time and space using balloons and aircraft detectors. In this work, variations in the cosmic ray intensity (single count rate) using the TASD, were studied and found to be on average at the $\sim(0.5-1)\%$ and up to 2\% level. These observations were found to be both in excess and in deficit. They were also found to be correlated with lightning in addition to thunderstorms. These variations lasted for tens of minutes; their footprint on the ground ranged from 6 to 24 km in diameter and moved in the same direction as the thunderstorm. With the use of simple electric field models inside the cloud and between cloud to ground, the observed variations in the cosmic ray single count rate were recreated using CORSIKA simulations. Depending on the electric field model used and the direction of the electric field in that model, the electric field magnitude that reproduces the observed low-energy cosmic ray single count rate variations was found to be approximately between 0.2-0.4 GV. This in turn allows us to get a reasonable insight on the electric field and its effect on cosmic ray air showers inside thunderstorms. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

Authors: Yu-Hsuan Teng, Karin M. Sandstrom, Jiayi Sun, Adam K. Leroy, L. Clifton Johnson, Alberto D. Bolatto, J. M. Diederik Kruijssen, Andreas Schruba, Antonio Usero, Ashley T. Barnes, Frank Bigiel, Guillermo A. Blanc, Brent Groves, Frank P. Israel, Daizhong Liu, Erik Rosolowsky, Eva Schinnerer, J. D. Smith, Fabian Walter

Abstract: The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred sp… ▽ More The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line radiative transfer modeling and a Bayesian likelihood analysis, we infer the H$_2$ density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of $2{-}3\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$1 kpc and a high temperature of 30$-$60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H$_2$ abundance of $3\times10^{-4}$, our analysis yields $α_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $α_\rm{CO} < 0.1$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted $α_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $α_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows. △ Less

Submitted 28 December, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 31 pages, 17 figures; accepted for publication in ApJ; font size of text and figures adjusted

Indications of a Cosmic Ray Source in the Perseus-Pisces Supercluster

Authors: Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon , et al. (135 additional authors not shown)

Abstract: The Telescope Array Collaboration has observed an excess of events with $E \ge 10^{19.4} ~{\rm eV}$ in the data which is centered at (RA, dec) = ($19^\circ$, $35^\circ$). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about $70 ~{\rm Mpc}$ distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A… ▽ More The Telescope Array Collaboration has observed an excess of events with $E \ge 10^{19.4} ~{\rm eV}$ in the data which is centered at (RA, dec) = ($19^\circ$, $35^\circ$). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about $70 ~{\rm Mpc}$ distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A Li-Ma oversampling analysis with $20^\circ$-radius circles indicates an excess in the arrival direction of events with a local significance of about 4 standard deviations. The probability of having such excess close to the PPSC by chance is estimated to be 3.5 standard deviations. This result indicates that a cosmic ray source likely exists in that supercluster. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 8 pages, 4 figures, 1 table

Milky Way Mid-Infrared Spitzer Spectroscopic Extinction Curves: Continuum and Silicate Features

Authors: Karl D. Gordon, Karl A. Misselt, Jeroen Bouwman, Geoffrey C. Clayton, Marjorie Decleir, Dean C. Hines, Yvonne Pendleton, George Rieke, J. D. T. Smith, D. C. B. Whittet

Abstract: We measured the mid-infrared (MIR) extinction using Spitzer photometry and spectroscopy (3.6--37 micron) for a sample of Milky Way sightlines (mostly) having measured ultraviolet extinction curves. We used the pair method to determine the MIR extinction that we then fit with a power law for the continuum and modified Drude profiles for the silicate features. We derived 16 extinction curves having… ▽ More We measured the mid-infrared (MIR) extinction using Spitzer photometry and spectroscopy (3.6--37 micron) for a sample of Milky Way sightlines (mostly) having measured ultraviolet extinction curves. We used the pair method to determine the MIR extinction that we then fit with a power law for the continuum and modified Drude profiles for the silicate features. We derived 16 extinction curves having a range of A(V) (1.8-5.5) and R(V) values (2.4-4.3). Our sample includes two dense sightlines that have 3 micron ice feature detections and weak 2175 A bumps. The average A(lambda)/A(V) diffuse sightline extinction curve we calculate is lower than most previous literature measurements. This agrees better with literature diffuse dust grain models, though it is somewhat higher. The 10 micron silicate feature does not correlate with the 2175 A bump, for the first time providing direct observational confirmation that these two features arise from different grain populations. The strength of the 10 micron silicate feature varies by $\sim$2.5 and is not correlated with A(V) or R(V). It is well fit by a modified Drude profile with strong correlations seen between the central wavelength, width, and asymmetry. We do not detect other features with limits in A(lambda)/A(V) units of 0.0026 (5--10 micron), 0.004 (10--20 micron), and 0.008 (20-40 micron). We find that the standard prescription of estimating R(V) from C times E(K_s-V)/E(B-V) has C = -1.14 and a scatter of $\sim$7%. Using the IRAC 5.6 micron band instead of K_s gives C = -1.03 and the least scatter of $\sim$3\%. △ Less

Submitted 12 May, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: 27 pages, 18 figures, ApJ accepted

Characterizing the Multi-Phase Origin of [CII] Emission in M101 and NGC 6946 with Velocity Resolved Spectroscopy

Authors: Elizabeth Tarantino, Alberto D. Bolatto, Rodrigo Herrera-Camus, Andrew I. Harris, Mark Wolfire, Christof Buchbender, Kevin V. Croxall, Daniel A. Dale, Brent Groves, Rebecca C. Levy, Denise Riquelme, J. D. T. Smith, Jürgen Stutzki

Abstract: The [CII] fine-structure transition at 158 micron is frequently the brightest far-infrared line in galaxies. Due to its low ionization potential, C+ can trace the ionized, atomic, and molecular phases of the ISM. We present velocity resolved [CII] and [NII] pointed observations from SOFIA/GREAT on ~500 pc scales in the nearby galaxies M101 and NGC 6946 and investigate the multi-phase origin of [CI… ▽ More The [CII] fine-structure transition at 158 micron is frequently the brightest far-infrared line in galaxies. Due to its low ionization potential, C+ can trace the ionized, atomic, and molecular phases of the ISM. We present velocity resolved [CII] and [NII] pointed observations from SOFIA/GREAT on ~500 pc scales in the nearby galaxies M101 and NGC 6946 and investigate the multi-phase origin of [CII] emission over a range of environments. We show that ionized gas makes a negligible contribution to the [CII] emission in these positions using [NII] observations. We spectrally decompose the [CII] emission into components associated with the molecular and atomic phases using existing CO(2-1) and HI data and show that a peak signal-to-noise ratio of 10-15 is necessary for a reliable decomposition. In general, we find that in our pointings greater than or equal to 50% of the [CII] emission arises from the atomic phase, with no strong dependence on star formation rate, metallicity, or galactocentric radius. We do find a difference between pointings in these two galaxies, where locations in NGC 6946 tend to have larger fractions of [CII] emission associated with the molecular phase than in M101. We also find a weak but consistent trend for fainter [CII] emission to exhibit a larger contribution from the atomic medium. We compute the thermal pressure of the cold neutral medium through the [CII] cooling function and find log(P_th/k)=3.8-4.6 [K cm^-3], a value slightly higher than similar determinations, likely because our observations are biased towards star-forming regions. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 23 pages, 11 figures, 4 tables, accepted to The Astrophysical Journal

Great Observatories: The Past and Future of Panchromatic Astrophysics

Authors: L. Armus, S. T. Megeath, L. Corrales, M. Marengo, A. Kirkpatrick, J. D. Smith, M. Meyer, S. Gezari, R. P. Kraft, S. McCandliss, S. Tuttle, M. Elvis, M. Bentz, B. Binder, F. Civano, D. Dragomir, C. Espaillat, S. Finkelstein, D. B. Fox, M. Greenhouse, E. Hamden, J. Kauffmann, G. Khullar, J. Lazio, J. Lee , et al. (8 additional authors not shown)

Abstract: NASA's Great Observatories have opened up the electromagnetic spectrum from space, providing sustained access to wavelengths not accessible from the ground. Together, Hubble, Compton, Chandra, and Spitzer have provided the scientific community with an agile and powerful suite of telescopes with which to attack broad scientific questions, and react to a rapidly changing scientific landscape. As the… ▽ More NASA's Great Observatories have opened up the electromagnetic spectrum from space, providing sustained access to wavelengths not accessible from the ground. Together, Hubble, Compton, Chandra, and Spitzer have provided the scientific community with an agile and powerful suite of telescopes with which to attack broad scientific questions, and react to a rapidly changing scientific landscape. As the existing Great Observatories age, or are decommissioned, community access to these wavelengths will diminish, with an accompanying loss of scientific capability. This report, commissioned by the NASA Cosmic Origins, Physics of the Cosmos and Exoplanet Exploration Program Analysis Groups (PAGs), analyzes the importance of multi-wavelength observations from space during the epoch of the Great Observatories, providing examples that span a broad range of astrophysical investigations. △ Less

Submitted 31 March, 2021; originally announced April 2021.

Comments: A report by the Great Observatories Science Analysis Group (SAG-10), commissioned by NASA's Cosmic Origins, Physics of the Cosmos, and Exoplanet Exploration Program Analysis Groups. 87 pages, 23 figures

Surface detectors of the TAx4 experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino , et al. (124 additional authors not shown)

Abstract: Telescope Array (TA) is the largest ultrahigh energy cosmic-ray (UHECR) observatory in the Northern Hemisphere. It explores the origin of UHECRs by measuring their energy spectrum, arrival-direction distribution, and mass composition using a surface detector (SD) array covering approximately 700 km$^2$ and fluorescence detector (FD) stations. TA has found evidence for a cluster of cosmic rays with… ▽ More Telescope Array (TA) is the largest ultrahigh energy cosmic-ray (UHECR) observatory in the Northern Hemisphere. It explores the origin of UHECRs by measuring their energy spectrum, arrival-direction distribution, and mass composition using a surface detector (SD) array covering approximately 700 km$^2$ and fluorescence detector (FD) stations. TA has found evidence for a cluster of cosmic rays with energies greater than 57 EeV. In order to confirm this evidence with more data, it is necessary to increase the data collection rate.We have begun building an expansion of TA that we call TAx4. In this paper, we explain the motivation, design, technical features, and expected performance of the TAx4 SD. We also present TAx4's current status and examples of the data that have already been collected. △ Less

Submitted 1 March, 2021; originally announced March 2021.

Comments: 26 pages, 17 figures, submitted to Nuclear Inst. and Methods in Physics Research, A

The Case for Thermalization as a Contributor to the [CII] Deficit

Authors: Jessica Sutter, Daniel A. Dale, Karin Sandstrom, J. D. T. Smith, Alberto Bolatto, Mederic Boquien, Daniela Calzetti, Kevin V. Croxall, Ilse De Looze, Maude Galametz, Brent A. Groves, George Helou, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Eric W. Pelligrini, Christine Wilson, Mark G. Wolfire

Abstract: The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission fr… ▽ More The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission from 18 Local Volume galaxies has been split based on conditions within the interstellar medium where it originated. This is completed using the Key Insights in Nearby Galaxies: a Far-Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) surveys and the wide-range of wavelength information, from UV to far-infrared emission lines, available for a selection of star-forming regions within these samples. By comparing these subdivided [CII] emissions to isolated infrared emission and other properties, we find that the thermalization (collisional de-excitation) of the [CII] line in HII regions plays a significant role in the deficit observed in our sample. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: To be published in MNRAS, 8 pages, 6 figures

All the PAHs: an AKARI-Spitzer Cross Archival Spectroscopic Survey of Aromatic Emission in Galaxies

Authors: Thomas S. -Y. Lai, J. D. T. Smith, Shunsuke Baba, Henrik W. W. Spoon, Masatoshi Imanishi

Abstract: We present a large sample of 2.5-38 $μm$ galaxy spectra drawn from a cross-archival comparison in the AKARI-Spitzer Extragalactic Spectral Survey (ASESS), and investigate a subset of 113 star-forming galaxies with prominent polycyclic aromatic hydrocarbon (PAH) emission spanning a wide range of star formation properties. With AKARI's extended 2.5-5 $μm$ wavelength coverage, we self-consistently mo… ▽ More We present a large sample of 2.5-38 $μm$ galaxy spectra drawn from a cross-archival comparison in the AKARI-Spitzer Extragalactic Spectral Survey (ASESS), and investigate a subset of 113 star-forming galaxies with prominent polycyclic aromatic hydrocarbon (PAH) emission spanning a wide range of star formation properties. With AKARI's extended 2.5-5 $μm$ wavelength coverage, we self-consistently model for the first time all PAH emission bands using a modified version of PAHFIT. We find $L_\mathrm{PAH\,3.3}/L_\mathrm{IR}$ $\sim$ 0.1% and the 3.3 $μm$ PAH feature contributes $\sim$1.5-3% to the total PAH power -- somewhat less than earlier dust models have assumed. We establish a calibration between 3.3 $μm$ PAH emission and star formation rate, but also find regimes where it loses reliability, including at high luminosity and low metallicity. The 3.4 $μm$ aliphatic emission and a broad plateau feature centered at 3.47 $μm$ are also modeled. As the shortest wavelength PAH feature, 3.3 $μm$ is susceptible to attenuation, leading to a factor of $\sim$ 3 differences in the inferred star formation rate at high obscuration with different assumed attenuation geometries. Surprisingly, $L_\mathrm{PAH\, 3.3} / L_\mathrm{Σ\, PAH}$ shows no sign of decline at high luminosities, and the low metallicity dwarf galaxy II Zw 40 exhibits an unusually strong 3.3 $μm$ band; both results suggest either that the smallest PAHs are better able to survive under intense radiation fields than presumed, or that PAH emission is shifted to shorter wavelengths in intense and high energy radiation environments. A photometric surrogate for 3.3 $μm$ PAH luminosity using JWST/NIRCam is provided and found to be highly reliable at low redshift. △ Less

Submitted 10 October, 2020; originally announced October 2020.

Comments: 29 pages, 18 figures, accepted in ApJ

Observations of the Origin of Downward Terrestrial Gamma-Ray Flashes

Authors: J. W. Belz, P. R. Krehbiel, J. Remington, M. A. Stanley, R. U. Abbasi, R. LeVon, W. Rison, D. Rodeheffer, the Telescope Array Scientific Collaboration, :, T. Abu-Zayyad, M. Allen, E. Barcikowski, D. R. Bergman, S. A. Blake, M. Byrne, R. Cady, B. G. Cheon, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich , et al. (116 additional authors not shown)

Abstract: In this paper we report the first close, high-resolution observations of downward-directed terrestrial gamma-ray flashes (TGFs) detected by the large-area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (I… ▽ More In this paper we report the first close, high-resolution observations of downward-directed terrestrial gamma-ray flashes (TGFs) detected by the large-area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (IBPs) in the first few milliseconds of negative cloud-to-ground and low-altitude intracloud flashes, and that the IBPs are produced by a newly-identified streamer-based discharge process called fast negative breakdown. The observations indicate the relativistic runaway electron avalanches (RREAs) responsible for producing the TGFs are initiated by embedded spark-like transient conducting events (TCEs) within the fast streamer system, and potentially also by individual fast streamers themselves. The TCEs are inferred to be the cause of impulsive sub-pulses that are characteristic features of classic IBP sferics. Additional development of the avalanches would be facilitated by the enhanced electric field ahead of the advancing front of the fast negative breakdown. In addition to showing the nature of IBPs and their enigmatic sub-pulses, the observations also provide a possible explanation for the unsolved question of how the streamer to leader transition occurs during the initial negative breakdown, namely as a result of strong currents flowing in the final stage of successive IBPs, extending backward through both the IBP itself and the negative streamer breakdown preceding the IBP. △ Less

Submitted 12 October, 2020; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: Typo fixed and reference added. Manuscript is 36 pages. Supplemental Information is 42 pages. This paper is to be published in the Journal of Geophysical Research: Atmospheres. Online data repository: Open Science Framework DOI: 10.17605/OSF.IO/Z3XDA

Search for Large-scale Anisotropy on Arrival Directions of Ultra-high-energy Cosmic Rays Observed with the Telescope Array Experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino , et al. (121 additional authors not shown)

Abstract: Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is… ▽ More Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger. △ Less

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

Comments: 6 pages, 3 figures, 1 table, Proofed title. Added journal reference and DOI

Journal ref: The Astrophysical Journal Letters 898, L28 (2020)

Measurement of the Proton-Air Cross Section with Telescope Array's Black Rock Mesa and Long Ridge Fluorescence Detectors, and Surface Array in Hybrid Mode

Authors: R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino, R. Higuchi , et al. (120 additional authors not shown)

Abstract: Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab.… ▽ More Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $σ^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $σ^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $σ^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. D 102, 062004 (2020)

Evidence for a Supergalactic Structure of Magnetic Deflection Multiplets of Ultra-High Energy Cosmic Rays

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino , et al. (119 additional authors not shown)

Abstract: Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and st… ▽ More Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$σ$. The ten years of data post-trial significance is 4.1$σ$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented. △ Less

Submitted 2 July, 2020; v1 submitted 14 May, 2020; originally announced May 2020.

Using [CII] 158 micron Emission From Isolated ISM Phases as a Star-Formation Rate Indicator

Authors: Jessica Sutter, Daniel A. Dale, Kevin V. Croxall, Eric W. Pelligrini, J. D. T. Smith, Phillip N. Appleton, Pedro Beirao, Alberto D. Bolatto, Daniela Calzetti, Alison Crocker, Ilse De Looze, Bruce Draine, Maud Galametz, Brent A. Groves, George Helou, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Helene Roussel, Mark G. Wolfire

Abstract: The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that origi… ▽ More The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 19 pages, 8 figures, 8 tables

Search for Ultra-High-Energy Neutrinos with the Telescope Array Surface Detector

Authors: R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino, K. Honda , et al. (112 additional authors not shown)

Abstract: We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the p… ▽ More We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the properties of the shower front and the lateral distribution function. △ Less

Submitted 12 May, 2020; v1 submitted 9 May, 2019; originally announced May 2019.

Comments: 10 pages, 4 figures, accepted to JETP

Journal ref: JETP Vol. 158 (8(2)) (2020)

Search for point sources of ultra-high energy photons with the Telescope Array surface detector

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. diMatteo, T. Fujii, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (114 additional authors not shown)

Abstract: The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (… ▽ More The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper. △ Less

Submitted 9 March, 2020; v1 submitted 30 March, 2019; originally announced April 2019.

Comments: accepted to MNRAS, 11 pages, 4 figures, 2 tables; results in text-file format are supplemented to paper source

Report number: INR-TH-2019-005

Journal ref: MNRAS 492 (2020), 3984

Simultaneous Measurements of Star Formation and Supermassive Black Hole Growth in Galaxies

Authors: Alexandra Pope, Lee Armus, Eric Murphy, Susanne Aalto, David Alexander, Philip Appleton, Amy Barger, Matt Bradford, Peter Capak, Caitlin Casey, Vassilis Charmandaris, Ranga Chary, Asantha Cooray, Jim Condon, Tanio Diaz Santos, Mark Dickinson, Duncan Farrah, Carl Ferkinhoff, Norman Grogin, Ryan Hickox, Allison Kirkpatrick, Kohno Kotaro, Allison Matthews, Desika Narayanan, Dominik Riechers , et al. (7 additional authors not shown)

Abstract: Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star for… ▽ More Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White paper submitted to Astro2020 Decadal Survey

Constraints on the diffuse photon flux with energies above $10^{18}$ eV using the surface detector of the Telescope Array experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (118 additional authors not shown)

Abstract: We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the wav… ▽ More We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the waveforms sensitive to the shower muon content. A total number of two photon candidates found in the search is fully compatible with the expected background. The $95\%\,$CL limits on the diffuse flux of the photons with energies greater than $10^{18.0}$, $10^{18.5}$, $10^{19.0}$, $10^{19.5}$ and $10^{20.0}$ eV are set at the level of $0.067$, $0.012$, $0.0036$, $0.0013$, $0.0013~\mbox{km}^{-2}\mbox{yr}^{-1}\mbox{sr}^{-1}$ correspondingly. △ Less

Submitted 19 March, 2019; v1 submitted 9 November, 2018; originally announced November 2018.

Comments: 7 pages, 2 figures

Report number: INR-TH-2018-026

Journal ref: Astroparticle Physics, 110 (2019) 8-14

Probing the High-Redshift Universe with SPICA: Toward the Epoch of Reionization and Beyond

Authors: E. Egami, S. Gallerani, R. Schneider, A. Pallottini, L. Vallini, E. Sobacchi, A. Ferrara, S. Bianchi, M. Bocchio, S. Marassi, L. Armus, L. Spinoglio, A. W. Blain, M. Bradford, D. L. Clements, H. Dannerbauer, J. A. Fernández-Ontiveros, E. González-Alfonso, M. J. Griffin, C. Gruppioni, H. Kaneda, K. Kohno, S. C. Madden, H. Matsuhara, P. Najarro , et al. (14 additional authors not shown)

Abstract: With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper… ▽ More With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $μ$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI. △ Less

Submitted 25 September, 2018; v1 submitted 20 September, 2018; originally announced September 2018.

Comments: 24 pages, 8 figures, accepted for publication in PASA (SPICA special issue)

Comprehensive comparison of models for spectral energy distributions from 0.1 micron to 1 mm of nearby star-forming galaxies

Authors: L. K. Hunt, I. De Looze, M. Boquien, R. Nikutta, A. Rossi, S. Bianchi, D. A. Dale, G. L. Granato, R. C. Kennicutt, L. Silva, L. Ciesla, M. Relano, S. Viaene, B. Brandl, D. Calzetti, K. V. Croxall, B. T. Draine, M. Galametz, K. D. Gordon, B. A. Groves, G. Helou, R. Herrera-Camus, J. L. Hinz, J. Koda, S. Salim , et al. (4 additional authors not shown)

Abstract: We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of… ▽ More We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR. △ Less

Submitted 13 November, 2018; v1 submitted 11 September, 2018; originally announced September 2018.

Comments: 42 pages, 23 figures, 8 tables, 2 appendices: updated to v2, A&A, accepted

Journal ref: A&A 621, A51 (2019)

Testing a reported correlation between arrival directions of ultrahigh-energy cosmic rays and a flux pattern from nearby starburst galaxies using Telescope Array data

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (117 additional authors not shown)

Abstract: The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did n… ▽ More The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did not optimize the parameter values but kept them fixed to the best-fit values found by Auger, namely 9.7% for the anisotropic fraction of cosmic rays assumed to originate from the SBGs in the list and 12.9° for the angular scale of the correlations. The energy threshold we adopted is 43 EeV, corresponding to 39 EeV in Auger when taking into account the energy-scale difference between two experiments. We find that the TA data is compatible with isotropy to within 1.1σ and with the Auger result to within 1.4σ, meaning that it is not capable to discriminate between these two hypotheses. △ Less

Submitted 22 October, 2018; v1 submitted 5 September, 2018; originally announced September 2018.

Comments: More accurate title, fixed an inconsistency in the declination range used, with a few figures slightly changing as a result, other minor revisions

Journal ref: R. U. Abbasi et al 2018 ApJL 867 L27

Infrared molecular hydrogen lines in GRB host galaxies

Authors: K. Wiersema, A. Togi, D. Watson, L. Christensen, J. P. U. Fynbo, B. P. Gompertz, A. B. Higgins, A. J. Levan, S. R. Oates, S. Schulze, J. D. T. Smith, E. R. Stanway, R. L. C. Starling, D. Steeghs, N. R. Tanvir

Abstract: Molecular species, most frequently H_2, are present in a small, but growing, number of gamma-ray burst (GRB) afterglow spectra at redshifts z~2-3, detected through their rest-frame UV absorption lines. In rare cases, lines of vibrationally excited states of H_2 can be detected in the same spectra. The connection between afterglow line-of-sight absorption properties of molecular (and atomic) gas, a… ▽ More Molecular species, most frequently H_2, are present in a small, but growing, number of gamma-ray burst (GRB) afterglow spectra at redshifts z~2-3, detected through their rest-frame UV absorption lines. In rare cases, lines of vibrationally excited states of H_2 can be detected in the same spectra. The connection between afterglow line-of-sight absorption properties of molecular (and atomic) gas, and the observed behaviour in emission of similar sources at low redshift, is an important test of the suitability of GRB afterglows as general probes of conditions in star formation regions at high redshift. Recently, emission lines of carbon monoxide have been detected in a small sample of GRB host galaxies, at sub-mm wavelengths, but no searches for H_2 in emission have been reported yet. In this paper we perform an exploratory search for rest-frame K band rotation-vibrational transitions of H_2 in emission, observable only in the lowest redshift GRB hosts (z<0.22). Searching the data of four host galaxies, we detect a single significant rotation-vibrational H_2 line candidate, in the host of GRB 031203. Re-analysis of Spitzer mid-infrared spectra of the same GRB host gives a single low significance rotational line candidate. The (limits on) line flux ratios are consistent with those of blue compact dwarf galaxies in the literature. New instrumentation, in particular on the JWST and the ELT, can facilitate a major increase in our understanding of the H_2 properties of nearby GRB hosts, and the relation to H_2 absorption in GRBs at higher redshift. △ Less

Submitted 28 August, 2018; originally announced August 2018.

Comments: 7 pages, 4 figures, MNRAS accepted

Mass composition of ultra-high-energy cosmic rays with the Telescope Array Surface Detector Data

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (118 additional authors not shown)

Abstract: The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events in… ▽ More The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events induced by the primary protons and iron. An average atomic mass of UHECR is presented for energies $10^{18.0}-10^{20.0}\ \mbox{eV}$. The average atomic mass of primary particles shows no significant energy dependence and corresponds to $\langle \ln A \rangle = 2.0 \pm 0.1 (stat.) \pm 0.44 (syst.)$. The result is compared to the mass composition obtained by the Telescope Array with $\mbox{X}_{\mbox{max}}$ technique along with the results of other experiments. Possible systematic errors of the method are discussed. △ Less

Submitted 24 January, 2019; v1 submitted 10 August, 2018; originally announced August 2018.

Comments: 11 pages, 9 figures

Journal ref: Phys. Rev. D 99, 022002 (2019)

Why Post-Starburst Galaxies are Now Quiescent

Authors: K. Decker French, Ann I. Zabludoff, Ilsang Yoon, Yancy Shirley, Yujin Yang, Adam Smercina, J. D. Smith, Desika Narayanan

Abstract: Post-starburst or "E+A" galaxies are rapidly transitioning from star-forming to quiescence. While the current star formation rate of post-starbursts is already at the level of early type galaxies, we recently discovered that many have large CO-traced molecular gas reservoirs consistent with normal star forming galaxies. These observations raise the question of why these galaxies have such low star… ▽ More Post-starburst or "E+A" galaxies are rapidly transitioning from star-forming to quiescence. While the current star formation rate of post-starbursts is already at the level of early type galaxies, we recently discovered that many have large CO-traced molecular gas reservoirs consistent with normal star forming galaxies. These observations raise the question of why these galaxies have such low star formation rates. Here we present an ALMA search for the denser gas traced by HCN (1--0) and HCO+ (1--0) in two CO-luminous, quiescent post-starburst galaxies. Intriguingly, we fail to detect either molecule. The upper limits are consistent with the low star formation rates and with early-type galaxies. The HCN/CO luminosity ratio upper limits are low compared to star-forming and even many early type galaxies. This implied low dense gas mass fraction explains the low star formation rates relative to the CO-traced molecular gas and suggests the state of the gas in post-starburst galaxies is unusual, with some mechanism inhibiting its collapse to denser states. We conclude that post-starbursts galaxies are now quiescent because little dense gas is available, in contrast to the significant CO-traced lower density gas reservoirs that still remain. △ Less

Submitted 30 May, 2018; originally announced May 2018.

Comments: accepted for publication in ApJ

Study of muons from ultra-high energy cosmic ray air showers measured with the Telescope Array experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. Di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (117 additional authors not shown)

Abstract: One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured numbe… ▽ More One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high $\it{muon \, purity}$. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of $\sim$65% at a lateral distance of 2000 m from the shower core is $1.72 \pm 0.10{\rm (stat.)} \pm 0.37 {\rm (syst.)}$ times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a $1.67 \pm 0.10 \pm 0.36$ excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances. △ Less

Submitted 11 April, 2018; originally announced April 2018.

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

The Knee and the Second Knee of the Cosmic-Ray Energy Spectrum

Authors: T. Abu-Zayyad, D. Ivanov, C. C. H. Jui, J. H. Kim, J. N. Matthews, J. D. Smith, S. B. Thomas, G. B. Thomson, Z. Zundel

Abstract: The cosmic ray flux measured by the Telescope Array Low Energy Extension (TALE) exhibits three spectral features: the knee, the dip in the $10^{16}$ eV decade, and the second knee. Here the spectrum has been measured for the first time using fluorescence telescopes, which provide a calorimetric, model-independent result. The spectrum appears to be a rigidity-dependent cutoff sequence, where the kn… ▽ More The cosmic ray flux measured by the Telescope Array Low Energy Extension (TALE) exhibits three spectral features: the knee, the dip in the $10^{16}$ eV decade, and the second knee. Here the spectrum has been measured for the first time using fluorescence telescopes, which provide a calorimetric, model-independent result. The spectrum appears to be a rigidity-dependent cutoff sequence, where the knee is made by the hydrogen and helium portions of the composition, the dip comes from the reduction in composition from helium to metals, the rise to the second knee occurs due to intermediate range nuclei, and the second knee is the iron knee. △ Less

Submitted 19 March, 2018; originally announced March 2018.

Comments: 10 pages, 5 figures, submitted to Astroparticle physics

The Cosmic-Ray Energy Spectrum between 2 PeV and 2 EeV Observed with the TALE detector in monocular mode

Authors: R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. Di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino, K. Honda , et al. (116 additional authors not shown)

Abstract: We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a lo… ▽ More We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources. △ Less

Submitted 3 March, 2018; originally announced March 2018.

Comments: 26 pages, 24 figures, submitted to ApJ

Journal ref: ApJ 865 74 (2018)

Evidence of Intermediate-Scale Energy Spectrum Anisotropy of Cosmic Rays E$\geq$10$^{19.2}$ eV with the Telescope Array Surface Detector

Authors: R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. Di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino, K. Honda , et al. (115 additional authors not shown)

Abstract: An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson l… ▽ More An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$σ_{global}$). △ Less

Submitted 17 February, 2018; v1 submitted 14 February, 2018; originally announced February 2018.

After The Fall: The Dust and Gas in E+A Post-Starburst Galaxies

Authors: A. Smercina, J. D. T. Smith, D. A. Dale, K. D. French, K. V. Croxall, S. Zhukovska, A. Togi, E. F. Bell, A. F. Crocker, B. T. Draine, T. H. Jarrett, C. Tremonti, Yujin Yang, A. I. Zabludoff

Abstract: The traditional picture of post-starburst galaxies as dust- and gas-poor merger remnants, rapidly transitioning to quiescence, has been recently challenged. Unexpected detections of a significant ISM in many post-starbursts raise important questions. Are they truly quiescent and, if so, what mechanisms inhibit further star formation? What processes dominate their ISM energetics? We present an infr… ▽ More The traditional picture of post-starburst galaxies as dust- and gas-poor merger remnants, rapidly transitioning to quiescence, has been recently challenged. Unexpected detections of a significant ISM in many post-starbursts raise important questions. Are they truly quiescent and, if so, what mechanisms inhibit further star formation? What processes dominate their ISM energetics? We present an infrared spectroscopic and photometric survey of 33 SDSS-selected E+A post-starbursts, aimed at resolving these questions. We find compact, warm dust reservoirs with high PAH abundances, and total gas and dust masses significantly higher than expected from stellar recycling alone. Both PAH/TIR and dust-to-burst stellar mass ratios are seen to decrease with post-burst age, indicative of the accumulating effects of dust destruction and an incipient transition to hot, early-type ISM properties. Their infrared spectral properties are unique, with dominant PAH emission, very weak nebular lines, unusually strong H$_{2}$ rotational emission, and deep ${\rm [C\, II]}$ deficits. There is substantial scatter among SFR indicators, and both PAH and TIR luminosities provide overestimates. Even as potential upper limits, all tracers show that the SFR has typically experienced a more than two order-of-magnitude decline since the starburst, and that the SFR is considerably lower than expected given both their stellar masses and molecular gas densities. These results paint a coherent picture of systems in which star formation was, indeed, rapidly truncated, but in which the ISM was $\textit{not}$ completely expelled, and is instead supported against collapse by latent or continued injection of turbulent or mechanical heating. The resulting aging burst populations provide a "high-soft" radiation field which seemingly dominates the E+As' unusual ISM energetics. △ Less

Submitted 6 June, 2018; v1 submitted 13 February, 2018; originally announced February 2018.

Comments: 44 pages, 21 figures, Appendix Table 2 updated