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.

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.

Context. EE Cep is one of few eclipsing binary systems with a dark, dusty disk around an invisible object similar to {\epsilon} Aur. The system is characterized by grey and asymmetric eclipses every 5.6 yr, with a significant variation in their photometric depth, ranging from ~ 0 m .5 to ~ 2 m .0. Aims. The main aim of the observational campaign of the EE Cep eclipse in 2014 was to test the model of disk precession (Galan et al. 2012). We expected that this eclipse would be one of the deepest with a depth of ~ 2 m .0. Methods. We collected multicolor observations from almost 30 instruments located in Europe and North America. This photometric data covers 243 nights during and around the eclipse. We also analyse the low- and high-resolution spectra from several instruments. Results. The eclipse was shallow with a depth of 0 m .71 in V-band. The multicolor photometry illustrates small color changes during the eclipse with a total amplitude of order ~ +0 m . 15 in B-I color index. The linear ephemeris for this system is updated by including new times of minima, measured from the three most recent eclipses at epochs E = 9, 10 and 11. New spectroscopic observations were acquired, covering orbital phases around the eclipse, which were not observed in the past and increased the data sample, filling some gaps and giving a better insight into the evolution of the H {\alpha} and NaI spectral line profiles during the primary eclipse. Conclusions. The eclipse of EE Cep in 2014 was shallower than expected 0 m .71 instead of ~ 2 m . 0. This means that our model of disk precession needs revision.

The past two decades have seen dramatic progress in our knowledge of the population of young stars of age $< 200\,$Myr that lie within $150\,$pc of the Sun. These nearby, young stars, most of which are found in loose, comoving groups, provide the opportunity to explore (among many other things) the dissolution of stellar clusters and their diffusion into the field star population. Here, we exploit the combination of astrometric and photometric data from {\it Gaia} and photometric data from GALEX (UV) and 2MASS (near-IR) in an attempt to identify additional nearby, young, late-type stars. Specifically, we present a sample of 146 GALEX UV-selected late-type (predominantly K-type) field stars with {\it Gaia}-based distances $< 125\,$pc (based on {\it Gaia} Data Release 1) that have isochronal ages $< 80\,$Myr even if equal-components binaries. We investigate the spectroscopic and kinematic properties of this sample. Despite their young isochronal ages, only $\sim 10$ per cent) of stars among this sample can be confidently associated with established nearby, young moving groups (MGs). These candidate MG members include 5 stars newly identified in this study. The vast majority of our sample of 146 nearby young star candidates have anomalous kinematics relative to the known MGs. These stars may hence represent a previously unrecognised population of young stars that has recently mixed into the older field star population. We discuss the implications and caveats of such a hypothesis---including the intriguing fact that, in addition to their non-young-star-like kinematics, the majority of the UV-selected, isochronally young field stars within $50\,$pc appear surprisingly X-ray faint.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

Gaia is currently revolutionizing modern astronomy. However, much of the Galactic plane, center and the spiral arm regions are obscured by interstellar extinction, rendering them inaccessible because Gaia is an optical instrument. An all-sky near infrared (NIR) space observatory operating in the optical NIR, separated in time from the original Gaia would provide microarcsecond NIR astrometry and millimag photometry to penetrate obscured regions unraveling the internal dynamics of the Galaxy.

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 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 past two decades have seen dramatic progress in our knowledge of the population of young stars of age <200 Myr that lie within 150 pc of the Sun. These nearby, young stars, most of which are found in loose, comoving groups, provide the opportunity to explore (among many other things) the dissolution of stellar clusters and their diffusion into the field star population. In the age of Gaia, this potential can now be fully exploited. We have identified, and are now investigating, a sample of nearly 400 Galex UV-selected late-type (K and early-M) field stars with Gaia-based distances <120 pc and isochronal ages <=80 Myr (even if binaries). Only a small percentage (<10%) of stars among this (kinematically unbiased) sample can be confidently associated with established nearby, young moving groups (NYMGs). The majority display anomalous kinematics, relative to the known NYMGs. These stars may hence represent a previously unrecognized population of young stars that has recently mixed into the older field star population. We discuss the implications and caveats of such a hypothesis---including the intriguing fact that, in addition to their non-young-star-like kinematics, the majority of the UV-selected, isochronally young field stars within 50 pc appear surprisingly X-ray faint.

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.

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.

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.

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}$).

BANYAN $\Sigma$ is a new Bayesian algorithm to identify members of young stellar associations within 150 pc of the Sun. It includes 27 young associations with ages in the range ~1-800 Myr, modelled with multivariate Gaussians in 6-dimensional XYZUVW space. It is the first such multi-association classification tool to include the nearest sub-groups of the Sco-Cen OB star-forming region, the IC 2602, IC 2391, Pleiades and Platais 8 clusters, and the $\rho$ Ophiuchi, Corona Australis, and Taurus star-formation regions. A model of field stars is built from a mixture of multivariate Gaussians based on the Besan\c{c}on Galactic model. The algorithm can derive membership probabilities for objects with only sky coordinates and proper motion, but can also include parallax and radial velocity measurements, as well as spectrophotometric distance constraints from sequences in color-magnitude or spectral type-magnitude diagrams. BANYAN $\Sigma$ benefits from an analytical solution to the Bayesian marginalization integrals that makes it more accurate and significantly faster than its predecessor BANYAN II. A contamination versus hit rate analysis is presented and demonstrates that BANYAN $\Sigma$ achieves a better classification performance than other moving group classification tools, especially in terms of cross-contamination between young associations. An updated list of bona fide members in the 27 young associations, augmented by the Gaia-DR1 release, are presented. This new tool will make it possible to analyze large data sets such as the upcoming Gaia-DR2 to identify new young stars. IDL and Python versions of BANYAN $\Sigma$ are made available with this publication. (shortened)

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 spectroscopy of emission lines for 81 Seyfert 1 and 104 Seyfert 2 galaxies in the IRAS 12$\mu$m galaxy sample. We analyzed the emission-line luminosity functions, reddening, and other gas diagnostics. The narrow-line regions (NLR) of Sy1 and 2 galaxies do not significantly differ from each other in most of these diagnostics. Combining the H$\alpha$/H$\beta$ ratio with a new reddening indicator-the [SII]6720/[OII]3727 ratio, we find the average $E(B-V)=0.49\pm0.35$ for Sy1s and $0.52\pm0.26$ for Sy2s. The NLR of Sy1 galaxies has only marginally higher ionization than the Sy2s. Our sample includes 22 Sy1.9s and 1.8s. In their narrow lines, these low-luminosity Seyferts are more similar to the Sy2s than the Sy1s. We construct a BPT diagram, and include the Sy1.8s and 1.9s. They overlap the region occupied by the Sy2s. The C IV equivalent width correlates more strongly with [O III]/H$\beta$ than with UV luminosity. The Sy1 and Sy2 luminosity functions of [OII]3727 and [OIII]5007 are indistinguishable. Unlike the LF's of Seyfert galaxies measured by SDSS, ours are nearly flat at low L. The larger number of faint Sloan "AGN" is attributable to their inclusion of weakly emitting LINERs and H II+AGN "composite" nuclei, which do not meet our classification criteria for Seyferts. An Appendix investigates which emission line luminosities provide the most reliable measures of the total non-stellar luminosity. The hard X-ray or near-ultraviolet continuum luminosity can be crudely predicted from either the [O III]5007 luminosity, or the combination of [O III]+H$\beta$, or [N II]+H$\alpha$ lines, with a scatter of $\pm\,4$ times for the Sy1s and $\pm\,10$ times for the Sy2s. The latter two hybrid (NLR+BLR) indicators have the advantage of predicting the same HX luminosity independent of Seyfert type.

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.

The Galex Nearby Young Star Survey (GALNYSS) has yielded a sample of $\sim$2000 UV-selected objects that are candidate nearby ($D \stackrel{<}{\sim}$150 pc), young (age $\sim$10--100 Myr), late-type stars. Here, we evaluate the distances and ages of the subsample of (19) GALNYSS stars with Gaia Data Release 1 (DR1) parallax distances $D \le 120$ pc. The overall youth of these 19 mid-K to early-M stars is readily apparent from their positions relative to the loci of main sequence stars and giants in Gaia-based color-magnitude and color-color diagrams constructed for all Galex- and WISE-detected stars with parallax measurements included in DR1. The isochronal ages of all 19 stars lie in the range $\sim$10--100 Myr. Comparison with Li-based age estimates indicates a handful of these stars may be young main-sequence binaries rather than pre-main sequence stars. Nine of the 19 objects have not previously been considered as nearby, young stars, and all but one of these are found at declinations north of $+$30$^\circ$. The Gaia DR1 results presented here indicate that the GALNYSS sample includes several hundred nearby, young stars, a substantial fraction of which have not been previously recognized as having ages $\stackrel{<}{\sim}$100 Myr.

TYC 8241 2652 1 is a young star that showed a strong mid-infrared (mid-IR, 8-25 mu) excess in all observations before 2008 consistent with a dusty disk. Between 2008 and 2010 the mid-IR luminosity of this system dropped dramatically by at least a factor of 30 suggesting a loss of dust mass of an order of magnitude or more. We aim to constrain possible models including removal of disk material by stellar activity processes, the presence of a binary companion, or other explanations suggested in the literature. We present new X-ray observations, optical spectroscopy, near-IR interferometry, and mid-IR photometry of this system to constrain its parameters and further explore the cause of the dust mass loss. In X-rays TYC 8241 2652 1 has all properties expected from a young star: Its luminosity is in the saturation regime and the abundance pattern shows enhancement of O/Fe. The photospheric Ha line is filled with a weak emission feature, indicating chromospheric activity consistent with the observed level of coronal emission. Interferometry does not detect a companion and sets upper limits on the companion mass of 0.2, 0.35, 0.1 and 0.05 M_sun at projected physical separations of 0.1-4 AU,4-5 AU, 5-10 AU, and 10-30 AU, respectively (assuming a distance of 120.9 pc). Our mid-IR measurements, the first of the system since 2012, are consistent with the depleted dust level seen after 2009. The new data confirms that stellar activity is unlikely to destroy the dust in the disk and shows that scenarios where either TYC 8241 2652 1 heats the disk of a binary companion or a potential companion heats the disk of TYC 8241 2652 1 are unlikely.

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.

An increasing number of observations have shown that gaseous debris discs are not an exception. However, until now we only knew of cases around A stars. Here we present the first detection of 12CO (2-1) disc emission around an F star, HD 181327, obtained with ALMA observations at 1.3 mm. The continuum and CO emission are resolved into an axisymmetric disc with ring-like morphology. Using a Markov chain Monte Carlo method coupled with radiative transfer calculations we study the dust and CO mass distribution. We find the dust is distributed in a ring with a radius of 86.0 +- 0.4 AU and a radial width of 23.2 +- 1.0 AU. At this frequency the ring radius is smaller than in the optical, revealing grain size segregation expected due to radiation pressure. We also report on the detection of low level continuum emission beyond the main ring out to ~200 AU. We model the CO emission in the non-LTE regime and we find that the CO is co-located with the dust, with a total CO gas mass ranging between 1.2x10^-6 Mearth and 2.9x10^-6 Mearth, depending on the gas kinetic temperature and collisional partners densities. The CO densities and location suggest a secondary origin, i.e. released from icy planetesimals in the ring. We derive a CO+CO2 cometary composition that is consistent with Solar system comets. Due to the low gas densities it is unlikely that the gas is shaping the dust distribution.

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.

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.

X-ray observations of pre-main sequence (pre-MS) stars of M-type probe coronal emission and offer a means to investigate magnetic activity at the stellar-substellar boundary. Recent observations of main sequence (MS) stars at this boundary display a decrease in fractional X-ray luminosity ($L_{X}$/$L_{bol}$) by almost two orders of magnitude for spectral types M7 and later. We investigate magnetic activity and search for a decrease in X-ray emission in the pre-MS progenitors of these MS stars. We present XMM-Newton X-ray observations and preliminary results for ~10 nearby (30-70 pc), very low mass pre-MS stars in the relatively unexplored age range of 10-30 Myr. We compare the fractional X-ray luminosities of these 10-30 Myr old stars to younger (1-3 Myr) pre-MS brown dwarfs and find no dependence on spectral type or age suggesting that X-ray activity declines at an age later than ~30 Myr in these very low-mass stars.

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.

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.

Low-mass, pre-main sequence stars possess intense high-energy radiation fields as a result of their strong stellar magnetic activity. This stellar UV and X-ray radiation may have a profound impact on the lifetimes of protoplanetary disks. We aim to constrain the X-ray-induced photoevaporation rates of protoplanetary disks orbiting low-mass stars by analyzing serendipitous XMM-Newton and Chandra X-ray observations of candidate nearby (D $<$ 100 pc), young (age $<$ 100 Myr) M stars identified in the GALEX Nearby Young-Star Survey (GALNYSS).

We have carried out an ALMA survey of 15 confirmed or candidate low-mass (<0.2M$_\odot$) members of the TW Hya Association (TWA) with the goal of detecting molecular gas in the form of CO emission, as well as providing constraints on continuum emission due to cold dust. Our targets have spectral types of M4-L0 and hence represent the extreme low end of the TWA's mass function. Our ALMA survey has yielded detections of 1.3mm continuum emission around 4 systems (TWA 30B, 32, 33, & 34), suggesting the presence of cold dust grains. All continuum sources are unresolved. TWA 34 further shows 12CO(2-1) emission whose velocity structure is indicative of Keplerian rotation. Among the sample of known ~7-10 Myr-old star/disk systems, TWA 34, which lies just ~50 pc from Earth, is the lowest mass star thus far identified as harboring cold molecular gas in an orbiting disk.

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.

Circumstellar dust disks have been observed around many nearby stars. However, many stars are part of binary or multiple stellar systems. A natural question arises regarding the presence and properties of such disks in systems with more than one star. To address this, we consider a sample of 449 systems (spectral types A-M) observed with the Herschel Space Observatory as part of the DEBRIS program. We have examined the stellar multiplicity of this sample by gathering information from the literature and performing an adaptive optics imaging survey at Lick Observatory. Five new companions were revealed with our program. In total, we identify 188 (42%) binary or multiple star systems. The multiplicity of the sample is examined with regards to the detection of circumstellar disks for stars of spectral types AFGK. In general, disks are less commonly detected around binaries than single stars, though the disk frequency is comparable among A stars regardless of multiplicity. However, this sample reveals the period distribution of disk-bearing binaries is consistent with that of non-disk binaries and with comparison field samples. We find that the properties of disks in binary systems are not statistically different from those around single stars. Although the frequency of disk-bearing FGK binaries may be lower than in single star systems, the processes behind disk formation and the characteristics of these disks are comparable among both populations.

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 study of the planet-debris disk connection can shed light on the formation and evolution of planetary systems, and may help predict the presence of planets around stars with certain disk characteristics. In preliminary analyses of the Herschel DEBRIS and DUNES surveys, Wyatt et al. (2012) and Marshall et al. (2014) identified a tentative correlation between debris and low-mass planets. Here we use the cleanest possible sample out these surveys to assess the presence of such a correlation, discarding stars without known ages, with ages < 1 Gyr and with binary companions <100 AU, to rule out possible correlations due to effects other than planet presence. In our sample of 204 FGK stars, we do not find evidence that debris disks are more common or more dusty around stars harboring high-mass or low-mass planets compared to a control sample without identified planets, nor that debris disks are more or less common (or more or less dusty) around stars harboring multiple planets compared to single-planet systems. Diverse dynamical histories may account for the lack of correlations. The data show the correlation between the presence of high-mass planets and stellar metallicity, but no correlation between the presence of low-mass planets or debris and stellar metallicity. Comparing the observed cumulative distribution of fractional luminosity to those expected from a Gaussian distribution, we find that a distribution centered on the Solar system's value fits well the data, while one centered at 10 times this value can be rejected. This is of interest in the context of future terrestrial planet characterization because it indicates that there are good prospects for finding a large number of debris disk systems (i.e. with evidence of harboring the building blocks of planets) with exozodiacal emission low enough to be appropriate targets for an ATLAST-type mission to search for biosignatures.

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 have conducted the first comprehensive mm-wave molecular emission line surveys of the evolved circumstellar disks orbiting the nearby T Tauri stars TW Hya and V4046 Sgr AB. Both disks are known to retain significant residual gaseous components, despite the advanced ages of their host stars. Our unbiased broad-band radio spectral surveys of the TW Hya and V4046 Sgr disks were performed with the Atacama Pathfinder Experiment (APEX) 12 meter telescope and are intended to yield a complete census of bright molecular emission lines in the range 275-357 GHz (1.1-0.85 mm). We find that lines of 12CO, 13CO, HCN, CN, and C2H, all of which lie in the higher-frequency range, constitute the strongest molecular emission from both disks in the spectral region surveyed. The molecule C2H is detected here for the first time in both disks, as is CS in the TW Hya disk. The survey results also include the first measurements of the full suite of hyperfine transitions of CN N=3-2 and C2H N=4-3 in both disks. Modeling of these CN and C2H hyperfine complexes in the spectrum of TW Hya indicates that the emission from both species is optically thick and may originate from very cold disk regions. It furthermore appears that the fractional abundances of CN and C2H are significantly enhanced in these evolved protoplanetary disks relative to the fractional abundances of the same molecules in the environments of deeply embedded protostars.

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$).

The UKIDSS Galactic Plane Survey (GPS) began in 2005 as a 7 year effort to survey ~1800 square degrees of the northern Galactic plane in the J, H, and K passbands. The survey included a second epoch of K band data, with a baseline of 2 to 8 years, for the purpose of investigating variability and measuring proper motions. We have calculated proper motions for 167 Million sources in a 900 square degree area located at l > 60 degrees in order to search for new high proper motion objects. Visual inspection has verified 617 high proper motion sources (> 200 mas/yr) down to K=17, of which 153 are new discoveries. Among these we have a new spectroscopically confirmed T5 dwarf, an additional T dwarf with estimated type T6, 13 new L dwarf candidates, and two new common proper motion systems containing ultracool dwarf candidates. We provide improved proper motions for an additional 12 high proper motion stars that were independently discovered in the WISE dataset during the course of this investigation.

We present spectra of an M-type, binary star system (LDS 5606) that belongs to the nearby ~20 Myr old beta Pictoris moving group. Both stars are very dusty; the dustier member displays optical emission lines from eight elements indicative of ongoing mass accretion. The spectra of both stars contain oxygen forbidden line emission at 6302 and 5579 A, consistent with a recent model of far ultraviolet photodissociation of OH molecules in a circumstellar disk. These are the oldest dwarf stars presently known to display such a phenomenon. The spectral energy distribution of the dustier star indicates substantial quantities of dust as hot as 900 K, and its fractional infrared luminosity (L_IR)/L_bol) is almost as large as that of the main sequence record holder, V488 Per. The LDS 5606 binary joins a remarkable group of very dusty, old, T Tauri stars that belong to widely separated multiple systems.

We report the identification of a new wide separation binary (LDS 5606) in the ~20 Myr-old beta Pic moving group. This M5+M5 pair has a projected separation of 26'', or ~1700 AU at a distance of 65 pc. Both stars host warm circumstellar disks and many strong hydrogen and helium emission lines. Spectroscopic observations reveal signatures of youth for both stars and on-going mass accretion in the primary. The properties of LDS 5606 make it an older analog to the ~8 Myr TWA 30 system, which is also composed of a pair of widely separated mid-M dwarfs, each hosting their own warm circumstellar disks. LDS 5606 joins a rather exclusive club of only 3 other known stellar systems where both members of a binary, far from any molecular cloud, are orbited by detected circumstellar disks.

We present basic observational data and association membership analysis for 45 young and active low-mass stellar systems from the ongoing RECONS photometry and astrometry program at the Cerro Tololo Inter-American Observatory. Most of these systems have saturated X-ray emission (log(Lx/Lbol) > -3.5) based on X-ray fluxes from the ROSAT All-Sky Survey, and many are significantly more luminous than main-sequence stars of comparable color. We present parallaxes and proper motions, Johnson-Kron-Cousins VRI photometry, and multiplicity observations from the CTIOPI program on the CTIO 0.9m telescope. To this we add low-resolution optical spectroscopy and line measurements from the CTIO 1.5m telescope, and interferometric binary measurements from the Hubble Space Telescope Fine Guidance Sensors. We also incorporate data from published sources: JHKs photometry from the 2MASS point source catalog; X-ray data from the ROSAT All-Sky Survey; and radial velocities from literature sources. Within the sample of 45 systems, we identify 21 candidate low-mass pre-main-sequence members of nearby associations, including members of beta Pictoris, TW Hydrae, Argus, AB Doradus, two ambiguous 30 Myr old systems, and one object that may be a member of the Ursa Major moving group. Of the 21 candidate young systems, 14 are newly identified as a result of this work, and six of those are within 25 parsecs of the Sun.

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.

Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model.

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 describe a method that exploits data from the GALEX ultraviolet and WISE and 2MASS infrared source catalogs, combined with proper motions and empirical pre-main sequence isochrones, to identify candidate nearby, young, low-mass stars. Applying our method across the full GALEX- covered sky, we identify 2031 mostly M-type stars that, for an assumed age of 10 (100) Myr, all lie within ~150 (~90) pc of Earth. The distribution of M spectral subclasses among these ~2000 candidate young stars peaks sharply in the range M3-M4; these subtypes constitute 50% of the sample, consistent with studies of the M star population in the immediate solar neighborhood. We focus on a subset of 58 of these candidate young M stars in the vicinity of the Tucana-Horologium Association. Only 20 of these 58 candidates were detected in the ROSAT All-Sky X-ray Survey -- reflecting the greater sensitivity of GALEX for purposes of identifying active nearby, young stars, particularly for stars of type M4 and later. Based on statistical analysis of the kinematics and/or spectroscopic followup of these 58 M stars, we find that 50% (29 stars) indeed have properties consistent with Tuc-Hor membership, while 12 are potential new members of the Columba Association, and two may be AB Dor moving group members. Hence, ~75% of our initial subsample of 58 candidates are likely members of young (age ~10-40 Myr) stellar moving groups within 100 pc, verifying that the stellar color- and kinematics-based selection algorithms described here can be used to efficiently isolate nearby, young, low-mass objects from among the field star population. Future studies will focus on characterizing additional subsamples selected from among this list of candidate nearby, young M stars.

We search for correlations between positions of extragalactic objects and arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) with primary energy $E \ge 40$ EeV as observed by the surface detector array of the Telescope Array (TA) experiment during the first 40 months of operation. We examined several public astronomical object catalogs, including the Veron-Cetty and Veron catalog of active galactic nuclei. We counted the number of TA events correlated with objects in each catalog as a function of three parameters: the maximum angular separation between a TA event and an object, the minimum energy of the events, and the maximum redshift of the objects. We determine combinations of these parameters which maximize the correlations, and calculate the chance probabilities of having the same levels of correlations from an isotropic distribution of UHECR arrival directions. No statistically significant correlations are found when penalties for scanning over the above parameters and for searching in several catalogs are taken into account.

We measure the spectrum of cosmic rays with energies greater than $10^{18.2}$ eV with the Fluorescence Detectors (FDs) and the Surface Detectors (SDs) of the Telescope Array Experiment using the data taken in our first 2.3-year observation from May 27 2008 to September 7 2010. A hybrid air shower reconstruction technique is employed to improve accuracies in determination of arrival directions and primary energies of cosmic rays using both FD and SD data. The energy spectrum presented here is in agreement with our previously published spectra and the HiRes results.