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$\sigma$. The ten years of data post-trial significance is 4.1$\sigma$. 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.

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.

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, $\sigma^{\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. $\sigma^{\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 $\sigma^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.

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.

We measure the gas disc thicknesses of the edge-on galaxy NGC 4013 and the less edge-on galaxies (NGC 4157 and 5907) using CO (CARMA/OVRO) and/or HI (EVLA) observations. We also estimate the scale heights of stars and/or the star formation rate (SFR) for our sample of five galaxies using Spitzer IR data (3.6 $\mu$m and 24 $\mu$m). We derive the average volume densities of the gas and the SFR using the measured scale heights along with radial surface density profiles. Using the volume density that is more physically relevant to the SFR than the surface density, we investigate the existence of a volumetric star formation law (SFL), how the volumetric SFL is different from the surface-density SFL, and how the gas pressure regulates the SFR based on our galaxy sample. We find that the volumetric and surface SFLs in terms of the total gas have significantly different slopes, while the volumetric and surface SFLs in terms of the molecular gas do not show any noticeable difference. The volumetric SFL for the total gas has a flatter power-law slope of 1.26 with a smaller scatter of 0.19 dex compared to the slope (2.05) and the scatter (0.25 dex) of the surface SFL. The molecular gas SFLs have similar slopes of 0.78 (volume density) and 0.77 (surface density) with the same rms scatter. We show that the interstellar gas pressure is strongly correlated with the SFR but find no significant difference between the correlations based on the volume and surface densities.

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.

We combine Hubble Space Telescope (HST) Paschen $\beta$ (Pa$\beta$) imaging with ground-based, previously published H$\alpha$ maps to estimate the attenuation affecting H$\alpha$, A(H$\alpha$), across the nearby, face-on galaxies NGC 5194 and NGC 6946. We estimate A(H$\alpha$) in ~ 2,000 independent 2" ~75 pc diameter apertures in each galaxy, spanning out to a galactocentric radius of almost 10 kpc. In both galaxies, A(H$\alpha$) drops with radius, with a bright, high attenuation inner region, though in detail the profiles differ between the two galaxies. Regions with the highest attenuation-corrected H$\alpha$ luminosity show the highest attenuation, but the observed H$\alpha$ luminosity of a region is not a good predictor of attenuation in our data. Consistent with much previous work, the IR-to-H$\alpha$ color does a good job of predicting A(H$\alpha$). We calculate the best-fit empirical coefficients for use combining H$\alpha$ with 8, 12, 24, 70, or 100 $\mu$m to correct for attenuation. These agree well with previous work but we also measure significant scatter around each of these linear relations. The local atomic plus molecular gas column density, N(H), also predicts A(H$\alpha$) well. We show that a screen with magnitude ~ 0.2 times the expected for a Milky Way gas-to-dust value does a reasonable job of explaining A(H$\alpha$) as a function of N(H). This could be expected if only ~ 40% of gas and dust directly overlap regions of H$\alpha$ emission.

High-resolution ($<$0.1pc) ALMA observations of the 30Dor-10 molecular cloud 15pc north of R136 are presented. The $^{12}$CO 2-1 emission morphology contains clumps near the locations of known mid-infrared massive protostars, as well as a series of parsec-long filaments oriented almost directly towards R136. There is elevated kinetic energy (linewidths at a given size scale) in 30Dor-10 compared to other LMC and Galactic star formation regions, consistent with large scale energy injection to the region. Analysis of the cloud substructures is performed by segmenting emission into disjoint approximately round "cores" using clumpfind, by considering the hierarchical structures defined by isointensity contours using dendrograms, and by segmenting into disjoint long thin "filaments" using Filfinder. Identified filaments have widths $\sim$0.1pc. The inferred balance between gravity and kinematic motions depends on the segmentation method: Entire objects identified with clumpfind are consistent with free-fall collapse or virial equilibrium with moderate external pressure, whereas many dendrogram-identified parts of hierarchical structures have higher mass surface densities $\Sigma_{LTE}$ than if gravitational and kinetic energies were in balance. Filaments have line masses that vary widely compared to the critical line mass calculated assuming thermal and nonthermal support. Velocity gradients in the region do not show any strong evidence for accretion of mass along filaments. The upper end of the "core" mass distribution is consistent with a power-law with the same slope as the stellar initial mass function.

We present an empirical relation between the cold gas surface density ($\Sigma_{\rm gas}$) and the optical extinction (${\rm A_V}$) in a sample of 103 galaxies from the Extragalactic Database for Galaxy Evolution (EDGE) survey. This survey provides CARMA interferometric CO observations for 126 galaxies included in the Calar Alto Legacy Integral Field Area (CALIFA) survey. The matched, spatially resolved nature of these data sets allows us to derive the $\Sigma_{\rm gas}$-${\rm A_V}$ relation on global, radial, and kpc (spaxel) scales. We determine ${\rm A_V}$ from the Balmer decrement (H$\alpha$/H$\beta$). We find that the best fit for this relation is $\Sigma_{\rm gas} ({\rm M_\odot pc^{-2}})\sim~26~\times~ {\rm A_V}({\rm mag})$, and that it does not depend on the spatial scale used for the fit. However, the scatter in the fits increases as we probe smaller spatial scales, reflecting the complex relative spatial distributions of stars, gas, and dust. We investigate the $\Sigma_{\rm gas}$/ ${\rm A_V}$ ratio on radial and spaxel scales as a function of ${\rm EW(H\alpha)}$. We find that at larger values of ${\rm EW(H\alpha)}$ (i.e., actively star-forming regions) this ratio tend to converge to the value expected for dust-star mixed geometries ($\sim$ 30 $\mathrm{M_{\odot} \,pc^{-2}\,mag^{-1}}$). On radial scales, we do not find a significant relation between the $\Sigma_{\rm gas}$/${\rm A_V}$ ratio and the ionized gas metallicity. We contrast our estimates of $\Sigma_{\rm gas}$ using ${\rm A_V}$ with compilations in the literature of the gas fraction on global and radial scales as well as with well known scaling relations such as the radial star-formation law and the $\Sigma_{\rm gas}$-$\Sigma_*$ relation. These tests show that optical extinction is a reliable proxy for estimating $\Sigma_{\rm gas}$ in the absence of direct sub/millimeter observations of the cold gas.

Cold atomic hydrogen clouds are the precursors of molecular clouds. Due to self-absorption, the opacity of cold atomic hydrogen may be high, and this gas may constitute an important mass component of the interstellar medium (ISM). Atomic hydrogen gas can be cooled to temperatures much lower than found in the cold neutral medium (CNM) through collisions with molecular hydrogen. In this paper, we search for HI Narrow Self-Absorption (HINSA) features in the Large Magellanic Cloud (LMC) as an indicator of such cold HI clouds, and use the results to quantify atomic masses and atomic-to-molecular gas ratio. Our search for HINSA features was conducted towards molecular clouds in the LMC using the ATCA+Parkes HI survey and the MAGMA CO survey. HINSA features are prevalent in the surveyed sightlines. This is the first detection of HINSA in an external galaxy. The HINSA-HI/$\rm{H}_{2}$ ratio in the LMC varies from 0.5\e{-3} to 3.4\e{-3} (68\% interval), with a mean value of $(1.31 \pm 0.03)$\e{-3}, after correcting for the effect of foreground HI gas. This is similar to the Milky Way value and indicates that similar fractions of cold gas exist in the LMC and the Milky Way, despite their differing metallicities, dust content and radiation fields. The low ratio also confirms that, as with the Milky Way, the formation timescale of molecular clouds is short. The ratio shows no radial gradient, unlike the case for stellar metallicity. No correlation is found between our results and those from previous HI absorption studies of the LMC.

We present a comparative study of the size-line width relation for substructures within six molecular clouds in the Large Magellanic Cloud (LMC) mapped with the Atacama Large Millimeter/submillimeter Array (ALMA). Our sample extends our previous study, which compared a Planck detected cold cloud in the outskirts of the LMC with the 30 Doradus molecular cloud and found the typical line width for 1 pc radius structures to be 5 times larger in 30 Doradus. By observing clouds with intermediate levels of star formation activity, we find evidence that line width at a given size increases with increasing local and cloud-scale 8${\mu}$m intensity. At the same time, line width at a given size appears to independently correlate with measures of mass surface density. Our results suggest that both virial-like motions due to gravity and local energy injection by star formation feedback play important roles in determining intracloud dynamics.

We present the first results from the Small Magellanic Cloud portion of a new Australia Telescope Compact Array (ATCA) HI absorption survey of both of the Magellanic Clouds, comprising over 800 hours of observations. Our new HI absorption line data allow us to measure the temperature and fraction of cold neutral gas in a low metallicity environment. We observed 22 separate fields, targeting a total of 55 continuum sources against 37 of which we detected HI absorption; from this we measure a column density weighted mean average spin temperature of $<T_{s}>=150$ K. Splitting the spectra into individual absorption line features, we estimate the temperatures of different gas components and find an average cold gas temperature of $\sim{30}$ K for this sample, lower than the average of $\sim{40}$ K in the Milky Way. The HI appears to be evenly distributed throughout the SMC and we detect absorption in $67\%$ of the lines of sight in our sample, including some outside the main body of the galaxy ($N_{\text{HI}}>2\times{10^{21}}$ cm$^{-2}$). The optical depth and temperature of the cold neutral atomic gas shows no strong trend with location spatially or in velocity. Despite the low metallicity environment, we find an average cold gas fraction of $\sim{20\%}$, not dissimilar from that of the Milky Way.

We present a multilinear analysis to determine the significant predictors of star formation in galaxies using the combined EDGE-CALIFA sample of galaxies. We analyze 1845 kpc-scale lines of sight across 39 galaxies with molecular line emission measurements from EDGE combined with optical IFU data drawn from CALIFA. We use the Least Absolute Shrinkage and Selection Operator (LASSO) to identify significant factors in predicting star formation rates. We find that the local star formation rate surface density is increased by higher molecular gas surface densities and stellar surface densities. In contrast, we see lower star formation rates in systems with older stellar populations, higher gas- and stellar-phase metallicities and larger galaxy masses. We also find a significant increase in star formation rate with galactocentric radius normalized by the disk scale length, which suggests additional parameters regulating star formation rate not explored in this study.

We investigate the prevalence, properties, and kinematics of extraplanar diffuse ionized gas (eDIG) in a sample of 25 edge-on galaxies selected from the CALIFA survey. We measure ionized gas scale heights from ${\rm H\alpha}$ and find that 90% have measurable scale heights with a median of $0.8^{+0.7}_{-0.4}$ kpc. From the ${\rm H\alpha}$ kinematics, we find that 60% of galaxies show a decrease in the rotation velocity as a function of height above the midplane. This lag is characteristic of eDIG, and we measure a median lag of 21 km s$^{-1}$ kpc$^{-1}$ which is comparable to lags measured in the literature. We also investigate variations in the lag with radius. $\rm H{\small I}$ lags have been reported to systematically decrease with galactocentric radius. We find both increasing and decreasing ionized gas lags with radius, as well as a large number of galaxies consistent with no radial lag variation, and investigate these results in the context of internal and external origins for the lagging ionized gas. We confirm that the ${\rm [S{\small II}]}$/${\rm H\alpha}$ and ${\rm [N{\small II}]}$/${\rm H\alpha}$ line ratios increase with height above the midplane as is characteristic of eDIG. The ionization of the eDIG is dominated by star-forming complexes (leaky ${\rm H{\small II}}$ regions). We conclude that the lagging ionized gas is turbulent ejected gas likely resulting from star formation activity in the disk as opposed to gas in the stellar thick disk or bulge. This is further evidence for the eDIG being a product of stellar feedback and for the pervasiveness of this WIM-like phase in many local star-forming galaxies.

We perform binary evolution calculations on helium star - carbon-oxygen white dwarf (CO WD) binaries using the stellar evolution code MESA. This single degenerate channel may contribute significantly to thermonuclear supernovae at short delay times. We examine the thermal-timescale mass transfer from a 1.1 - 2.0 $M_{\odot}$ helium star to a 0.90 - 1.05 $M_{\odot}$ CO WD for initial orbital periods in the range 0.05 - 1 day. Systems in this range may produce a thermonuclear supernova, helium novae, a helium star - oxygen-neon WD binary, or a detached double CO WD binary. Our time-dependent calculations that resolve the stellar structures of both binary components allow accurate distinction between the eventual formation of a thermonuclear supernova (via central ignition of carbon burning) and that of an ONe WD (in the case of off-center ignition). Furthermore, we investigate the effect of a slow WD wind which implies a specific angular momentum loss from the binary that is larger than typically assumed. We find that this does not significantly alter the region of parameter space over which systems evolve toward thermonuclear supernovae. Our determination of the correspondence between initial binary parameters and the final outcome informs population synthesis studies of the contribution of the helium donor channel to thermonuclear supernovae. In addition, we constrain the orbital properties and observable stellar properties of the progenitor binaries of thermonuclear supernovae and helium novae.

There is wide consensus that galaxy outflows are one of the most important processes determining the evolution of galaxies through cosmic time, for example playing a key role in shaping the galaxy mass function. Our understanding of outflows and their drivers, however, is in its infancy --- this is particularly true for the cold (neutral atomic and molecular) phases of outflows, which present observational and modeling challenges. Here we outline several key open questions, briefly discussing the requirements of the observations necessary to make progress, and the relevance of several existing and planned facilities. It is clear that galaxy outflows, and particularly cold outflows, will remain a topic of active research for the next decade and beyond.

We describe a next major frontier in observational studies of galaxy evolution and star formation: linking the physical conditions in the cold, star-forming interstellar medium to host galaxy and local environment. The integrated gas content of galaxies has been surveyed extensively over the last decades. The link between environment and cold gas density, turbulence, excitation, dynamical state, and chemical makeup remain far less well understood. We know that these properties do vary dramatically and theoretical work posits a strong connection between the state of the gas, its ability to form stars, and the impact of stellar feedback. A next major step in the field will be to use sensitive cm-, mm-, and submm-wave spectroscopy and high resolution spectroscopic imaging to survey the state of cold gas across the whole local galaxy population. Such observations have pushed the capabilities of the current generation of telescopes. We highlight three critical elements for progress in the next decade: (1) robust support and aggressive development of ALMA, (2) the deployment of very large heterodyne receiver arrays on single dish telescopes, and (3) development of a new interferometric array that dramatically improves on current capabilities at cm- and mm-wavelengths (~ 1-115 GHz).

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.

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on available rheological measurements[4], that solid layers of N2 ice approximately greater than 1 km thick should convect for estimated present-day heat flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-km-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of N2-ice viscosity implies that the SP ice layer convects in the so-called sluggish lid regime[5], a unique convective mode heretofore not definitively observed in the Solar System. Average surface horizontal velocities of a few cm/yr imply surface transport or renewal times of ~500,000 years, well under the 10 Myr upper limit crater retention age for Sputnik Planum[2]. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help explain the high albedos of some of them.

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.

The aim of this study is to characterize the distribution and basic properties of the natal gas associated with high-mass young stellar objects (YSOs) in isolated environments in the Large Magellanic Cloud (LMC). High-mass stars usually form in Giant Molecular Clouds (GMCs) as part of a young stellar cluster, but some OB stars are observed far from GMCs. By examining the spatial coincidence between the high-mass YSOs and 12CO (J = 1-0) emission detected by NANTEN and Mopra observations, we selected ten high-mass YSOs that are located away from any of the NANTEN clouds but are detected by the Mopra pointed observations. The ALMA observations revealed that a compact molecular cloud whose mass is a few thousand solar masses or smaller is associated with the high-mass YSOs, which indicates that these compact clouds are the sites of high-mass star formation. The high-density and high-temperature throughout the clouds are explained by the severe photodissociation of CO due to the lower metallicity than in the Galaxy. The star formation efficiency ranges from several to as high as ~ 40%, indicating efficient star formation in these environments. The enhanced turbulence may be a cause of the efficient star formation therein, as judged from the gas velocity information and the association with the lower density gas.

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{\deg} 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{\sigma} and with the Auger result to within 1.4{\sigma}, meaning that it is not capable to discriminate between these two hypotheses.

Gas density is widely believed to play a governing role in star formation. However, the exact role of density in setting the star formation rate remains debated. We also lack a general theory that explains how the gas density distribution in galaxies is set. The primary factor preventing the resolution of these issues is the limited number of observations of the gas density distribution across diverse environments. Centimeter- and millimeter-wave spectroscopy offer the most promising way forward in this field, but the key density-sensitive transitions are faint compared to the capabilities of current telescopes. In this chapter, we describe how a next-generation Very Large Array (ngVLA) represents the natural next step forward in this sensitivity-limited field. Such a facility would provide a crucial link between the `Milky Way' and `Extragalactic' views of star formation and dramatically advance our understanding of the drive and role of gas density in galaxies, building on current results from ALMA, NOEMA, the Green Bank Telescope, and other current facilities working in this area.

Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved $^{12}$CO(J=1-0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and H$\alpha$ rotation curves using the same geometric parameters out to $\gtrsim$1 $R_e$. We find that $\sim$75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence supported disk. Using observations of the H$\gamma$ transition also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and H$\alpha$ rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.

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.

Deriving circular velocities of galaxies from stellar kinematics can provide an estimate of their total dynamical mass, provided a contribution from the velocity dispersion of the stars is taken into account. Molecular gas (e.g., CO) on the other hand, is a dynamically cold tracer and hence acts as an independent circular velocity estimate without needing such a correction. In this paper we test the underlying assumptions of three commonly used dynamical models, deriving circular velocities from stellar kinematics of 54 galaxies (S0-Sd) that have observations of both stellar kinematics from the CALIFA survey, and CO kinematics from the EDGE survey. We test the Asymmetric Drift Correction (ADC) method, as well as Jeans, and Schwarzschild models. The three methods each reproduce the CO circular velocity at 1Re to within 10%. All three methods show larger scatter (up to 20%) in the inner regions (R < 0.4Re) which may be due to an increasingly spherical mass distribution (which is not captured by the thin disk assumption in ADC), or non-constant stellar M/L ratios (for both the JAM and Schwarzschild models). This homogeneous analysis of stellar and gaseous kinematics validates that all three models can recover Mdyn at 1Re to better than 20%, but users should be mindful of scatter in the inner regions where some assumptions may break down.

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.

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$$\sigma_{global}$).

The energy of the ultrahigh energy spectral cutoff was measured, integrating over the northern hemisphere sky, by the Telescope Array (TA) collaboration, to be $10^{19.78 \pm 0.06}$ eV, in agreement with the High Resolution Fly's Eye (HiRes) experiment, whereas the Pierre Auger experiment, integrating over the southern hemisphere sky, measured the cutoff to be at 10$^{19.62 \pm 0.02}$ eV. An 11% energy scale difference between the TA and Auger does not account for this difference. However, in comparing the spectra of the Telescope Array and Pierre Auger experiments in the band of declination common to both experiments ($-15.7^{\circ} < \delta < 24.8^{\circ}$) we have found agreement in the energy of the spectral cutoff. While the Auger result is essentially unchanged, the TA cutoff energy has changed to $10^{19.59 \pm 0.06}$ eV. In this paper we argue that this is an astrophysical effect.

We present the molecular cloud properties of N55 in the Large Magellanic Cloud using $^{12}$CO(1-0) and $^{13}$CO(1-0) observations obtained with Atacama Large Millimeter Array. We have done a detailed study of molecular gas properties, to understand how the cloud properties of N55 differ from Galactic clouds. Most CO emission appears clumpy in N55, and molecular cores that have YSOs show larger linewidths and masses. The massive clumps are associated with high and intermediate mass YSOs. The clump masses are determined by local thermodynamic equilibrium and virial analysis of the $^{12}$CO and $^{13}$CO emissions. These mass estimates lead to the conclusion that, (a) the clumps are in self-gravitational virial equilibrium, and (b) the $^{12}$CO(1-0)-to-H$_2$ conversion factor, X$_{\rm CO}$, is 6.5$\times$10$^{20}$cm$^{-2}$(K km s$^{-1}$)$^{-1}$. This CO-to-H$_2$ conversion factor for N55 clumps is measured at a spatial scale of $\sim$0.67 pc, which is about two times higher than the X$_{\rm CO}$ value of Orion cloud at a similar spatial scale. The core mass function of N55 clearly show a turnover below 200M$_{\odot}$, separating the low-mass end from the high-mass end. The low-mass end of the $^{12}$CO mass spectrum is fitted with a power law of index 0.5$\pm$0.1, while for $^{13}$CO it is fitted with a power law index 0.6$\pm$0.2. In the high-mass end, the core mass spectrum is fitted with a power index of 2.0$\pm$0.3 for $^{12}$CO, and with 2.5$\pm$0.4 for $^{13}$CO. This power-law behavior of the core mass function in N55 is consistent with many Galactic clouds.

We present a kpc-scale analysis of the relationship between the molecular depletion time ($\tau_\mathrm{dep}^\mathrm{mol}$) and the orbital time ($\tau_\mathrm{orb}$) across the field of 39 face-on local galaxies, selected from the EDGE-CALIFA sample. We find that, on average, 5% of the available molecular gas is converted into stars per orbital time, or $\tau_\mathrm{dep}^\mathrm{mol}\sim20\tau_\mathrm{orb}$. The resolved relation shows a scatter of $\sim0.5$ dex. The scatter is ascribable to galaxies of different morphologies that follow different $\tau_\mathrm{dep}^\mathrm{mol}-\tau_\mathrm{orb}$ relations which decrease in steepness from early- to late-types. The morphologies appear to be linked with the star formation rate surface density, the molecular depletion time, and the orbital time, but they do not correlate with the molecular gas content of the galaxies in our sample. We speculate that in our molecular gas rich, early-type galaxies, the morphological quenching (in particular the disc stabilization via shear), rather than the absence of molecular gas, is the main factor responsible for their current inefficient star formation.

We present high-resolution (sub-parsec) observations of a giant molecular cloud in the nearest star-forming galaxy, the Large Magellanic Cloud. ALMA Band 6 observations trace the bulk of the molecular gas in $^{12}$CO(2-1) and high column density regions in $^{13}$CO(2-1). Our target is a quiescent cloud (PGCC G282.98-32.40, which we refer to as the "Planck cold cloud" or PCC) in the southern outskirts of the galaxy where star-formation activity is very low and largely confined to one location. We decompose the cloud into structures using a dendrogram and apply an identical analysis to matched-resolution cubes of the 30 Doradus molecular cloud (located near intense star formation) for comparison. Structures in the PCC exhibit roughly 10 times lower surface density and 5 times lower velocity dispersion than comparably sized structures in 30 Dor, underscoring the non-universality of molecular cloud properties. In both clouds, structures with relatively higher surface density lie closer to simple virial equilibrium, whereas lower surface density structures tend to exhibit super-virial line widths. In the PCC, relatively high line widths are found in the vicinity of an infrared source whose properties are consistent with a luminous young stellar object. More generally, we find that the smallest resolved structures ("leaves") of the dendrogram span close to the full range of line widths observed across all scales. As a result, while the bulk of the kinetic energy is found on the largest scales, the small-scale energetics tend to be dominated by only a few structures, leading to substantial scatter in observed size-linewidth relationships.

We present results from the EDGE survey, a spatially resolved CO(1-0) follow-up to CALIFA, an optical Integral Field Unit (IFU) survey of local galaxies. By combining the data products of EDGE and CALIFA, we study the variation in molecular gas depletion time ($\tau_{\rm dep}$) on kiloparsec scales in 52 galaxies. We divide each galaxy into two parts: the center, defined as the region within $0.1 \ R_{25}$, and the disk, defined as the region between $0.1$ and $0.7 \ R_{25}$. We find that 14 galaxies show a shorter $\tau_{\rm dep}$ ($\sim 1$ Gyr) in the center relative to that in the disk ($\tau_{\rm dep} \sim 2.4$ Gyrs), which means the central region in those galaxies is more efficient at forming stars per unit molecular gas mass. This finding implies that the centers with shorter $\tau_{\rm dep}$ resemble the intermediate regime between galactic disks and starburst galaxies. Furthermore, the central drop in $\tau_{\rm dep}$ is correlated with a central increase in the stellar surface density, suggesting that a shorter $\tau_{\rm dep}$ is associated with molecular gas compression by the stellar gravitational potential. We argue that varying the CO-to-H$_2$ conversion factor only exaggerates the central drop of $\tau_{\rm dep}$.

We present a $^{13}\mathrm{CO} (J = 1 \rightarrow 0)$ mapping survey of 12 nearby galaxies from the CARMA STING sample. The line intensity ratio $\mathcal{R} \equiv I[^{12}\mathrm{CO} (J = 1 \rightarrow 0)]/I[^{13}\mathrm{CO} (J = 1 \rightarrow 0)]$ is derived to study the variations in molecular gas properties. For 11 galaxies where it can be measured with high significance, the spatially resolved $\mathcal{R}$ on (sub-)kiloparsec scales varies by up to a factor of 3--5 within a galaxy. Lower $\mathcal{R}$ values are usually found in regions with weaker $^{12}\rm CO$. We attribute this apparent trend to a bias against measuring large $\mathcal{R}$ values when $^{12}\rm CO$ is weak. Limiting our analysis to the $^{12}\rm CO$ bright regions that are less biased, we do not find $\mathcal{R}$ on (sub)kpc scales correlate with galactocentric distance, velocity dispersion or the star formation rate. The lack of correlation between SFR and $\mathcal{R}$ indicates that the CO optical depth is not sensitive to stellar energy input, or that any such sensitivity is easily masked by other factors. Extending the analysis to all regions with $\rm ^{12}CO$ emission by spectral stacking, we find that 5 out of 11 galaxies show higher stacked $\mathcal{R}$ for galactocentric radii of $\gtrsim 1$ kpc and $\Sigma_{\mathrm{SFR}} \lesssim 0.1 \rm \ M_{sun} yr^{-1} kpc^{-2}$, which could result from a greater contribution from diffuse gas. Moreover, significant galaxy-to-galaxy variations are found in $\mathcal{R}$, but the global $\mathcal{R}$ does not strongly depend on dust temperature, inclination, or metallicity of the galaxy.

We present interferometric CO observations made with the Combined Array for Millimeter-wave Astronomy (CARMA) of galaxies from the Extragalactic Database for Galaxy Evolution survey (EDGE). These galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) sample, mapped with optical integral field spectroscopy. EDGE provides good quality CO data (3$\sigma$ sensitivity $\Sigma_{\rm mol}$ $\sim$ 11 M$_\odot$ pc$^{-2}$ before inclination correction, resolution $\sim1.4$ kpc) for 126 galaxies, constituting the largest interferometric CO survey of galaxies in the nearby universe. We describe the survey, the data characteristics, the data products, and present initial science results. We find that the exponential scale-lengths of the molecular, stellar, and star-forming disks are approximately equal, and galaxies that are more compact in molecular gas than in stars tend to show signs of interaction. We characterize the molecular to stellar ratio as a function of Hubble type and stellar mass, present preliminary results on the resolved relations between the molecular gas, stars, and star formation rate, and discuss the dependence of the resolved molecular depletion time on stellar surface density, nebular extinction, and gas metallicity. EDGE provides a key dataset to address outstanding topics regarding gas and its role in star formation and galaxy evolution, which will be publicly available on completion of the quality assessment.

The Telescope Array (TA) experiment is located in the western desert of Utah, USA, and observes ultra high energy cosmic rays (UHECRs) in the Northern hemisphere. At the highest energies, $E>10$~EeV, the shape of cosmic ray energy spectrum may carry an imprint of the source density distribution along the line of sight different in different directions of the sky. In this study, we search for such directional variations in the shape of the energy spectrum using events observed with the Telescope Array's surface detector. We divide the TA field of view into two nearly equal-exposure regions: the "on-source" region which we define as $\pm 30^\circ$ of the supergalactic plane containing mostly nearby structures, and the complementary "off-source" region where the sources are further away on average. We compare the UHECR spectra in these regions by fitting them to the broken power law and comparing the resulting parameters. We find that the off-source spectrum has an earlier break at highest energies. The chance probability to obtain such or larger difference in statistically equivalent distributions is estimated as $6.2\pm1.1\times10^{-4}$ ($3.2\sigma$) by a Monte-Carlo simulation. The observed difference in spectra is in a reasonable quantitative agreement with a simplified model that assumes that the UHECR sources trace the galaxy distribution from the 2MRS catalogue, primary particles are protons and the magnetic deflections can be neglected.

We present and analyze deep Herschel/HIFI observations of the [CII] 158um, [CI] 609um, and [CI] 370um lines towards 54 lines-of-sight (LOS) in the Large and Small Magellanic clouds. These observations are used to determine the physical conditions of the line--emitting gas, which we use to study the transition from atomic to molecular gas and from C^+ to C^0 to CO in their low metallicity environments. We trace gas with molecular fractions in the range 0.1<f(H2)<1, between those in the diffuse H2 gas detected by UV absorption (f(H2)<0.2) and well shielded regions in which hydrogen is essentially completely molecular. The C^0 and CO column densities are only measurable in regions with molecular fractions f(H2)>0.45 in both the LMC and SMC. Ionized carbon is the dominant gas-phase form of this element that is associated with molecular gas, with C^0 and CO representing a small fraction, implying that most (89% in the LMC and 77% in the SMC) of the molecular gas in our sample is CO-dark H2. The mean X_CO conversion factors in our LMC and SMC sample are larger than the value typically found in the Milky Way. When applying a correction based on the filling factor of the CO emission, we find that the values of X_CO in the LMC and SMC are closer to that in the Milky Way. The observed [CII] intensity in our sample represents about 1% of the total far-infrared intensity from the LOSs observed in both Magellanic Clouds.

CO observations allow estimations of the gas content of molecular clouds, which trace the reservoir of cold gas fuelling star formation, as well as to determine extinction via H$_2$ column density, N(H$_2$). Here, we studied millimetric and optical properties at 26 supernovae (SNe) locations of different types in a sample of 23 nearby galaxies by combining molecular $^{12}$C$^{16}$O (J = 1 $\rightarrow$ 0) resolved maps from the EDGE survey and optical Integral Field Spectroscopy from the CALIFA survey. We found an even clearer separation between type II and type Ibc SNe in terms of molecular gas than what we found in the optical using H$\alpha$ emission as a proxy for current SF rate, which reinforces the fact that SNe Ibc are more associated with SF-environments. While A$_V$ at SN locations is similar for SNe II and SNe Ibc, and higher compared to SNe Ia, N(H$_2$) is significantly higher for SNe Ibc than for SNe II and SNe Ia. When compared to alternative extinction estimations directly from SN photometry and spectroscopy, we find that our SNe Ibc have also redder color excess but showed standard Na I D absorption pseudo-equivalent widths ($\sim$1 \AA). In some cases we find no extinction when estimated from the environment, but high amounts of extinction when measured from SN observations, which suggests that circumstellar material or dust sublimation may be playing a role. This work serves as a benchmark for future studies combining last generation millimeter and optical IFS instruments to reveal the local environmental properties of extragalactic SNe.

Cosmic rays in the energy range $10^{18.0}$ - $10^{18.5}$ eV are thought to have a light, probably protonic, composition. To study their origin one can search for anisotropy in their arrival directions. Extragalactic cosmic rays should be isotropic, but galactic cosmic rays of this type should be seen mostly along the galactic plane, and there should be a shortage of events coming from directions near the galactic anticenter. This is due to the fact that, under the influence of the galactic magnetic field, the transition from ballistic to diffusive behavior is well advanced, and this qualitative picture persists over the whole energy range. Guided by models of the galactic magnetic field that indicate that the enhancement along the galactic plane should have a standard deviation of about 20$^\circ$ in galactic latitude, and the deficit in the galactic anticenter direction should have a standard deviation of about 50$^\circ$ in galactic longitude, we use the data of the Telescope Array surface detector in $10^{18.0}$ to $10^{18.5}$ eV energy range to search for these effects. The data are isotropic. Neither an enhancement along the galactic plane nor a deficit in the galactic anticenter direction is found. Using these data we place an upper limit on the fraction of EeV cosmic rays of galactic origin at 1.3% at 95% confidence level.

We present encube $-$ a qualitative, quantitative and comparative visualisation and analysis system, with application to high-resolution, immersive three-dimensional environments and desktop displays. encube extends previous comparative visualisation systems by considering: 1) the integration of comparative visualisation and analysis into a unified system; 2) the documentation of the discovery process; and 3) an approach that enables scientists to continue the research process once back at their desktop. Our solution enables tablets, smartphones or laptops to be used as interaction units for manipulating, organising, and querying data. We highlight the modularity of encube, allowing additional functionalities to be included as required. Additionally, our approach supports a high level of collaboration within the physical environment. We show how our implementation of encube operates in a large-scale, hybrid visualisation and supercomputing environment using the CAVE2 at Monash University, and on a local desktop, making it a versatile solution. We discuss how our approach can help accelerate the discovery rate in a variety of research scenarios.

Observations suggest that there is a significant fraction of O-stars in the field of the Milky Way that appear to have formed in isolation or in low mass clusters ($<$100 $M_\odot$). The existence of these high-mass stars that apparently formed in the field challenges the generally accepted paradigm, which requires star formation to occur in clustered environments. In order to understand the physical conditions for the formation of these stars, it is necessary to observe isolated high-mass stars while they are still forming. With the $Hubble$ $Space$ $Telescope$, we observe the seven most isolated massive ($>$8 $M_\odot$) young stellar objects (MYSOs) in the Large Magellanic Cloud (LMC). The observations show that while these MYSOs are remote from other MYSOs, OB associations, and even from known giant molecular clouds, they are actually not isolated at all. Imaging reveals $\sim$100 to several hundred pre--main-sequence (PMS) stars in the vicinity of each MYSO. These previously undetected PMS stars form prominent compact clusters around the MYSOs, and in most cases they are also distributed sparsely across the observed regions. Contrary to what previous high-mass field star studies show, these observations suggest that high-mass stars may not be able to form in clusters with masses less than 100 $M_\odot$. If these MYSOs are indeed the best candidates for isolated high-mass star formation, then the lack of isolation is at odds with random sampling of the IMF. Moreover, while isolated MYSOs may not exist, we find evidence that isolated clusters containing O-stars can exist, which in itself is rare.

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

We discuss the thermal pressures ($n_H T$) in predominantly cold, neutral interstellar gas in the Magellanic Clouds, derived from analyses of the fine-structure excitation of neutral carbon, as seen in high-resolution HST/STIS spectra of seven diverse sight lines in the LMC and SMC. Detailed fits to the line profiles of the absorption from C I, C I*, and C I** yield consistent column densities for the 3--6 C I multiplets detected in each sight line. In the LMC and SMC, $N$(C I$_{\rm tot}$) is consistent with Galactic trends versus $N$(Na I) and $N$(CH), but is slightly lower versus $N$(K I) and $N$(H$_2$). As for $N$(Na I) and $N$(K I), $N$(C I$_{\rm tot}$) is generally significantly lower, for a given $N$(H$_{\rm tot}$), in the LMC and (especially) in the SMC, compared to the local Galactic relationship. For the LMC and SMC components with well determined column densities for C I, C I*, and C I**, the derived thermal pressures are typically factors of a few higher than the values found for most cold, neutral clouds in the Galactic ISM. Such differences are consistent with the predictions of models for clouds in systems (like the LMC and SMC) that are characterized by lower metallicities, lower dust-to-gas ratios, and enhanced radiation fields -- where higher pressures are required for stable cold, neutral clouds. The pressures may be further enhanced by energetic activity (e.g., due to stellar winds, star formation, and/or supernova remnants) in several of the regions probed by these sight lines. Comparisons are made with the C I observed in some quasar absorption-line systems.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

We present an updated status of the EDGE project, which is a survey of 125 local galaxies in the $^{12}$CO($1-0$) and $^{13}$CO($1-0$) lines. We combine the molecular data of the EDGE survey with the stellar and ionized gas maps of the CALIFA survey to give a comprehensive view of the dependence of the star formation efficiency, or equivalently, the molecular gas depletion time, on various local environments, such as the stellar surface density, metallicity, and radius from the galaxy center. This study will provide insight into the parameters that drive the star formation efficiency in galaxies at $z \sim 0$.

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

This white paper discusses how a "next-generation" Very Large Array (ngVLA) operating in the frequency range 1-116 GHz could be a groundbreaking tool to study the detailed astrophysics of the "matter cycle" in the Milky Way and other galaxies. If optimized for high brightness sensitivity, the ngVLA would bring detailed microwave spectroscopy and modeling of the full radio spectral energy distribution into regular use as survey tools at resolutions of 0.1- 1 arcseconds. This wavelength range includes powerful diagnostics of density, excitation, and chemistry in the cold ISM, as well as multiple tracers of the rate of recent star formation, the magnetic field, shocks, and properties of the ionized ISM. We highlight design considerations that would make this facility revolutionary in this area, the foremost of which is a large amount of collecting area on ~km-length baselines. We also emphasize the strong case for harnessing the large proposed collecting area of the ngVLA for very long baseline applications as part of the core design. This would allow measurements of the three dimensional space motions of galaxies to beyond the Local Group and mapping of the Milky Way out to the far side of the disk. Finally, we discuss the gains from the proposed combination of very high resolution and sensitivity to thermal emission, which include observing the feeding of black holes and resolving forming protoclusters.

In this work we are reporting on the measurement of the proton-air inelastic cross section $\sigma^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $\sigma^{\rm inel}_{\rm p-air}$ the proton-proton cross section $\sigma_{\rm p-p}$ value is also determined at $\sqrt{s} = 95_{-8}^{+5}$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is the hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector collected over five years. The value of the $\sigma^{\rm inel}_{\rm p-air}$ is found to be equal to $567.0 \pm 70.5 [{\rm Stat.}] ^{+29}_{-25} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] _{-17}^{+19} [{\rm Sys.}] $mb.

We report high spatial resolution observations of Giant Molecular Clouds (GMCs) in the nearby spiral galaxies NGC 6946, M101 and NGC 628 obtained with the CARMA telescope. We observed CO(1-0) over regions with active star formation, and higher resolution observations of CO(2-1) have allowed us to resolve some of the largest GMCs. Using a Bayesian fitting approach, we generate scaling relations for the sizes, line widths, and virial masses of the structures identified in this work. We do not find evidence for a tight power law relation between size and line width, although the limited dynamic range in cloud size remains a clear issue in our analysis. Additionally, we use a Bayesian approach to analyze the Kennicutt-Schmidt (K-S) relation for the identified structures. We find that the distribution of slopes are broadly distributed, mainly due to the limited dynamic range of our measured H2 mass surface density, and being most consistent with super-linear relations. On the other hand, when we use the Bayesian approach to analyze the K-S relation for a uniform grid, the distributions of slopes is consistent with sub-linear relations. On-arm regions tend to have higher star formation rates than inter-arm regions. As in NGC 6946, in M101 we find regions where the star formation efficiency shows marked peaks at specific galoctocentric radii. On the other hand, the distribution of SFE in NGC 628 is more contiguous. We hypothesize that differences in the distribution of SFE may be indicative of different processes driving the spiral structure.

We present the first detection of HCO$^+$ absorption in the Magellanic System. Using the Australia Telescope Compact Array (ATCA), we observed 9 extragalactic radio continuum sources behind the Magellanic System and detected HCO$^+$ absorption towards one source located behind the leading edge of the Magellanic Bridge. The detection is located at LSR velocity of $v=214.0 \pm 0.4\rm\,km\,s^{-1}$, with a full width at half maximum of $\Delta v=4.5\pm 1.0\rm\,km\,s^{-1}$ and optical depth of $\tau(\rm HCO^+)=0.10\pm 0.02$. Although there is abundant neutral hydrogen (HI) surrounding the sightline in position-velocity space, at the exact location of the absorber the HI column density is low, $<10^{20}\rm\,cm^{-2}$, and there is little evidence for dust or CO emission from Planck observations. While the origin and survival of molecules in such a diffuse environment remains unclear, dynamical events such as HI flows and cloud collisions in this interacting system likely play an important role.

To test the theoretical understanding that finding bright CO emission depends primarily on dust shielding, we investigate the relationship between CO emission ($I_{\rm CO}$) and the amount of dust (estimated from IR emission and expressed as "$A_V$") across the Large Magellanic Cloud, the Small Magellanic Cloud, and the Milky Way. We show that at our common resolution of 10 pc scales, $I_{\rm CO}$ given a fixed line-of-sight $A_V$ is similar across all three systems despite the difference in metallicity. We find some evidence for a secondary dependence of $I_{\rm CO}$ on radiation field; in the LMC, $I_{\rm CO}$ at a given $A_V$ is smaller in regions of high $T_{\rm dust}$, perhaps because of an increased photodissociating radiation field. We suggest a simple but useful picture in which the CO-to-H$_2$ conversion factor (\xco) depends on two separable factors: (1) the distribution of gas column densities, which maps to an extinction distribution via a dust-to-gas ratio; and (2) the dependence of $I_{\rm CO}$ on $A_V$. Assuming that the probability distribution function (PDF) of local Milky Way clouds is universal, this approach predicts a dependence of \xco\ on $Z$ between $Z^{-1}$ and $Z^{-2}$ above about a third solar metallicity. Below this metallicity, CO emerges from only the high column density parts of the cloud and so depends very sensitively on the adopted PDF and the H$_2$/{\sc Hi} prescription. The PDF of low metallicity clouds is thus of considerable interest and the uncertainty associated with even an ideal prescription for \xco\ at very low metallicity will be large.

We report on the search for steady point-like sources of neutral particles around 10$^{18}$ eV between 2008 May and 2013 May with the scintillator surface detector of the Telescope Array experiment. We found overall no significant point-like excess above 0.5 EeV in the northern sky. Subsequently, we also searched for coincidence with the Fermi bright Galactic sources. No significant coincidence was found within the statistical uncertainty. Hence, we set an upper limit on the neutron flux that corresponds to an averaged flux of 0.07 km$^{-2}$ yr$^{-1}$ for $E>1$ EeV in the northern sky at the 95% confidence level. This is the most stringent flux upper limit in a northern sky survey assuming point-like sources. The upper limit at the 95% confidence level on the neutron flux from Cygnus X-3 is also set to 0.2 km$^{-2}$ yr$^{-1}$ for $E>0.5$ EeV. This is an order of magnitude lower than previous flux measurements.

The spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.

Previous measurements of the composition of Ultra-High Energy Cosmic Rays(UHECRs) made by the High Resolution Fly's Eye(HiRes) and Pierre Auger Observatory(PAO) are seemingly contradictory, but utilize different detection methods, as HiRes was a stereo detector and PAO is a hybrid detector. The five year Telescope Array(TA) Middle Drum hybrid composition measurement is similar in some, but not all, respects in methodology to PAO, and good agreement is evident between data and a light, largely protonic, composition when comparing the measurements to predictions obtained with the QGSJetII-03 and QGSJet-01c models. These models are also in agreement with previous HiRes stereo measurements, confirming the equivalence of the stereo and hybrid methods. The data is incompatible with a pure iron composition, for all models examined, over the available range of energies. The elongation rate and mean values of Xmax are in good agreement with Pierre Auger Observatory data. This analysis is presented using two methods: data cuts using simple geometrical variables and a new pattern recognition technique.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

We present a study of the vertical structure of the gaseous and stellar disks in a sample of edge-on galaxies (NGC 4157, 4565, and 5907) using BIMA/CARMA 12CO (J = 1 --> 0), VLA H I, and Spitzer 3.6 micron data. In order to take into account projection effects when we measure the disk thickness as a function of radius, we first obtain the inclination by modeling the radio data. Using the measurement of the disk thicknesses and the derived radial profiles of gas and stars, we estimate the corresponding volume densities and vertical velocity dispersions. Both stellar and gas disks have smoothly varying scale heights and velocity dispersions, contrary to assumptions of previous studies. Using the velocity dispersions, we find that the gravitational instability parameter Q follows a fairly uniform profile with radius and is greater than or equal to 1 across the star forming disk. The star formation law has a slope that is significantly different from those found in more face-on galaxy studies, both in deprojected and pixel-by-pixel plots. Midplane gas pressure based on the varying scale heights and velocity dispersions appears to roughly hold a power-law correlation with the midplane volume density ratio.

We have searched for intermediate-scale anisotropy in the arrival directions of ultrahigh-energy cosmic rays with energies above 57~EeV in the northern sky using data collected over a 5 year period by the surface detector of the Telescope Array experiment. We report on a cluster of events that we call the hotspot, found by oversampling using 20$^\circ$-radius circles. The hotspot has a Li-Ma statistical significance of 5.1$\sigma$, and is centered at R.A.=146.7$^{\circ}$, Dec.=43.2$^{\circ}$. The position of the hotspot is about 19$^{\circ}$ off of the supergalactic plane. The probability of a cluster of events of 5.1$\sigma$ significance, appearing by chance in an isotropic cosmic-ray sky, is estimated to be 3.7$\times$10$^{-4}$ (3.4$\sigma$).

We present parsec scale interferometric maps of HCN(1-0) and HCO$^{+}$(1-0) emission from dense gas in the star-forming region 30 Doradus, obtained using the Australia Telescope Compact Array. This extreme star-forming region, located in the Large Magellanic Cloud (LMC), is characterized by a very intense ultraviolet ionizing radiation field and sub-solar metallicity, both of which are expected to impact molecular cloud structure. We detect 13 bright, dense clumps within the 30 Doradus-10 giant molecular cloud. Some of the clumps are aligned along a filamentary structure with a characteristic spacing that is consistent with formation via the varicose fluid instability. Our analysis shows that the filament is gravitationally unstable and collapsing to form stars. There is a good correlation between HCO$^{+}$ emission in the filament and signatures of recent star formation activity including H$_{2}$O masers and young stellar objects (YSOs). YSOs seem to continue along the same direction of the filament toward the massive compact star cluster R136 in the southwest. We present detailed comparisons of clump properities (masses, linewidths, sizes) in 30Dor-10 to those in other star forming regions of the LMC (N159, N113, N105, N44). Our analysis shows that the 30Dor-10 clumps have similar mass but wider linewidths and similar HCN/HCO$^{+}$ (1-0) line ratios as clumps detected in other LMC star-forming regions. Our results suggest that the dense molecular gas clumps in the interior of 30Dor-10 are well-shielded against the intense ionizing field that is present in the 30Doradus region.

The extragalactic X-ray binary IC 10 X-1 has attracted attention as it is possibly the host of the most massive stellar-mass black-hole (BH) known to date. Here we consider all available observational constraints and construct its evolutionary history up to the instant just before the formation of the BH. Our analysis accounts for the simplest possible history that includes three evolutionary phases: binary orbital dynamics at core collapse, common envelope (CE) evolution, and evolution of the BH--helium star binary progenitor of the observed system. We derive the complete set of constraints on the progenitor system at various evolutionary stages. Specifically: right before the core collapse event, we find the mass of the BH immediate progenitor to be > 31 Msun (at 95% of confidence, same hereafter). The magnitude of the natal kick imparted to the BH is constrained to be < 130 km/s. Furthermore, we find that the "enthalpy" formalism recently suggested by Ivanova & Chaichenets is able to explain the existence of IC 10 X-1 without the need of invoking unreasonably high CE efficiencies. With this physically motivated formalism, we find that the CE efficiency required to explain the system is in the range of 0.6--1.

Supernovae are Nature's high-energy, high density laboratory experiments, reaching densities in excess of nuclear densities and temperatures above 10MeV. Astronomers have built up a suite of diagnostics to study these supernovae. If we can utilize these diagnostics, and tie them together with a theoretical understanding of supernova physics, we can use these cosmic explosions to study the nature of matter at these extreme densities and temperatures. Capitalizing on these diagnostics will require understanding a wide range of additional physics. Here we review the diagnostics and the physics needed to use them to learn about the supernova engine, and ultimate nuclear physics.

Although the details of the core-collapse supernova mechanism are not fully understood, it is generally accepted that the energy released in the collapse produces a shock that disrupts the star and produces the explosion. Some of the stellar material does not receive enough energy to escape the potential well of the newly formed neutron star and it falls back on to the core. This fallback plays an important role in setting compact remnant masses and has been invoked to explain a number of observed phenomena: late-time neutrino emission, r-process element production, peculiar supernovae, and long-duration gamma-ray bursts. To examine the link between fallback and these observed phenomena, we must understand the nature of fallback and currently, multiple mechanisms have been proposed that lead to different characteristics of the fallback. In this paper, we model a series of three-dimensional explosions by following the fallback for two different massive progenitors to study both the nature of the fallback and the role of asymmetric explosions on this fallback. One of the most-cited mechanisms behind fallback argues that it is caused by the reverse shock produced when the outward moving supernova shock decelerates in the relatively flat density profile of the stellar envelope. We show that this mechanism plays a minor, if any, role in driving fallback and most supernova fallback is produced by prompt mechanisms.

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

Star formation does not occur until the onset of gravitational collapse inside giant molecular clouds. However, the conditions that initiate cloud collapse and regulate the star formation process remain poorly understood. Local processes such as turbulence and magnetic fields can act to promote or prevent collapse. On larger scales, the galactic potential can also influence cloud stability and is traditionally assessed by the tidal and shear effects. In this paper, we examine the stability of giant molecular clouds (GMCs) in the Large Magellanic Cloud (LMC) against shear and the galactic tide using CO data from the Magellanic Mopra Assessment (MAGMA) and rotation curve data from the literature. We calculate the tidal acceleration experienced by individual GMCs and determine the minimum cloud mass required for tidal stability. We also calculate the shear parameter, which is a measure of a clouds susceptibility to disruption via shearing forces in the galactic disk. We examine whether there are correlations between the properties and star forming activity of GMCs and their stability against shear and tidal disruption. We find that the GMCs are in approximate tidal balance in the LMC, and that shear is unlikely to affect their further evolution. GMCs with masses close to the minimal stable mass against tidal disruption are not unusual in terms of their mass, location, or CO brightness, but we note that GMCs with large velocity dispersion tend to be more sensitive to tidal instability. We also note that GMCs with smaller radii, which represent the majority of our sample, tend to more strongly resist tidal and shear disruption. Our results demonstrate that star formation in the LMC is not inhibited by to tidal or shear instability.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

We investigate the correlation between CO and HI emission in 18 nearby galaxies from the CARMA Survey Toward IR-Bright Nearby Galaxies (STING) at sub-kpc and kpc scales. Our sample, spanning a wide range in stellar mass and metallicity, reveals evidence for a metallicity dependence of the HI column density measured in regions exhibiting CO emission. Such a dependence is predicted by the equilibrium model of McKee & Krumholz, which balances H_2 formation and dissociation. The observed HI column density is often smaller than predicted by the model, an effect we attribute to unresolved clumping, although values close to the model prediction are also seen. We do not observe HI column densities much larger than predicted, as might be expected were there a diffuse HI component that did not contribute to H_2 shielding. We also find that the H_2 column density inferred from CO correlates strongly with the stellar surface density, suggesting that the local supply of molecular gas is tightly regulated by the stellar disk.