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.

Observations of large-scale radio emissions prove the existence of shock accelerated cosmic-ray electrons in galaxy clusters, while the lack of detected $\gamma$-rays limits the acceleration of cosmic-ray protons in galaxy clusters. This challenges our understanding of how diffusive shock acceleration works. In this work, we couple the most updated recipes for shock acceleration in the intracluster medium to state-of-the-art magneto-hydrodynamical simulations of massive galaxy clusters. Furthermore, we use passive tracer particles to follow the evolution of accelerated cosmic-rays. We show that when the interplay between magnetic field topology and the feedback from accelerated cosmic rays is taken into account, the latest developments of particle acceleration theory give results which are compatible with observational constraints.

Although a merging galaxy cluster is a useful laboratory to study many important astrophysical processes and fundamental physics, only limited interpretations are possible without careful analysis of the merger scenario. However, the study is demanding because a thorough comparison of multi-wavelength observations with detailed numerical simulations is required. In this paper, we present such a study for the off-axis binary merger A115. The system possesses a number of remarkable observational features, but no convincing merger scenario, explaining the shape and location of the radio relic in harmony with the orientation of the cometary X-ray tails, has been presented. Our hydrodynamical simulation, with adaptive mesh refinement, suggests that the cometary X-ray tail of A115 might be a slingshot tail and can arise $\sim0.3$ Gyrs after the impact and before the two subclusters reach their apocenters. This scenario can predict the location and orientation of the giant radio relic, which is parallel to the northern X-ray tail. In addition, our study indicates that diffusive shock acceleration alone cannot generate the observed radio power unless aided by stronger magnetic fields and/or more significant presence of fossil electrons.

Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

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 examine the cosmic-ray protons (CRp) accelerated at collisionless shocks in galaxy clusters using cosmological structure formation simulations. We find that in the intracluster medium (ICM) within the virial radius of simulated clusters, only $\sim7$\% of shock kinetic energy flux is dissipated by the shocks that are expected to accelerate CRp, that is, supercritical, quasi-parallel ($Q_\parallel$) shocks with sonic Mach number $M_s\ge2.25$. The rest is dissipated at subcritical shocks and quasi-perpendicular shocks, both of which may not accelerate CRp. Adopting the diffusive shock acceleration (DSA) model recently presented in Ryu et al. (2019), we quantify the DSA of CRp in simulated clusters. The average fraction of the shock kinetic energy transferred to CRp via DSA is assessed at $\sim(1-2)\times10^{-4}$. We also examine the energization of CRp through reacceleration using a model based on the test-particle solution. Assuming that the ICM plasma passes through shocks three times on average through the history of the universe and that CRp are reaccelerated only at supercritical $Q_\parallel$-shocks, the CRp spectrum flattens by $\sim0.05-0.1$ in slope and the total amount of CRp energy increases by $\sim40-80$\% from reacceleration. We then estimate diffuse $\gamma$-ray and neutrino emissions, resulting from inelastic collisions between CRp and thermal protons. The predicted $\gamma$-ray emissions from simulated clusters lie mostly below the upper limits set by Fermi-LAT for observed clusters. The neutrino fluxes towards nearby clusters would be $\lesssim10^{-4}$ of the IceCube flux at $E_{\nu}=1$ PeV and $\lesssim10^{-6}$ of the atmospheric neutrino flux in the energy range of $E_{\nu}\leq1$ TeV.

Radio relics in galaxy clusters are giant diffuse synchrotron sources powered in cluster outskirts by merger shocks. Although the relic-shock connection has been consolidated in recent years by a number of observations, the details of the mechanisms leading to the formation of relativistic particles in this environment are still not well understood. The diffusive shock acceleration (DSA) theory is a commonly adopted scenario to explain the origin of cosmic rays at astrophysical shocks, including those in radio relics in galaxy clusters. However, in a few specific cases it has been shown that the energy dissipated by cluster shocks is not enough to reproduce the luminosity of the relics via DSA of thermal particles. Studies based on samples of radio relics are required to further address this limitation of the mechanism. In this paper, we focus on ten well-studied radio relics with underlying shocks observed in the X-rays and calculate the electron acceleration efficiency of these shocks that is necessary to reproduce the observed radio luminosity of the relics. We find that in general the standard DSA cannot explain the origin of the relics if electrons are accelerated from the thermal pool with an efficiency significantly smaller than 10%. Our results show that other mechanisms, such as shock re-acceleration of supra-thermal seed electrons or a modification of standard DSA, are required to explain the formation of radio relics.

The preacceleration of electrons through reflection and shock drift acceleration (SDA) is essential for the diffusive shock acceleration (DSA) of nonthermal electrons in collisionless shocks. Previous studies suggested that, in weak quasi-perpendicular ($Q_\perp$) shocks in the high-$\beta$ ($\beta=P_{\rm gas}/P_{\rm B}$) intracluster medium (ICM), the temperature anisotropy due to SDA-reflected electrons can drive the electron firehose instability, which excites oblique nonpropagating waves in the shock foot. In this paper, we investigate, through a linear analysis and particle-in-cell (PIC) simulations, the firehose instabilities driven by an electron temperature anisotropy (ETAFI) and also by a drifting electron beam (EBFI) in $\beta\sim100$ ICM plasmas. The EBFI should be more relevant in describing the self-excitation of upstream waves in $Q_\perp$-shocks, since backstreaming electrons in the shock foot behave more like an electron beam rather than an anisotropic bi-Maxwellian population. We find that the basic properties of the two instabilities, such as the growth rate, $\gamma$, and the wavenumber of fast-growing oblique modes are similar in the ICM environment, with one exception; while the waves excited by the ETAFI are nonpropagating ($\omega_r=0$), those excited by the EBFI have a non-zero frequency ($\omega_r\neq0$). However, the frequency is small with $\omega_r<\gamma$. Thus, we conclude that the interpretation of previous studies for the nature of upstream waves based on the ETAFI remains valid in $Q_\perp$-shocks in the ICM.

We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.14 \pm 0.01$. We compute the shock strength using density jumps to conclude that the Mach number of the merging subcluster is small ($M \leq 1.15_{-0.09}^{+0.14}$). We also present pseudo-density, projected temperature, pseudo-pressure, and pseudo-entropy maps. Based on the pseudo-entropy map we conclude that the cluster is undergoing a mild merger, consistent with the small Mach number. On the other hand, radio relics extend over Mpc scale in the A3411-12 system, which strongly suggests that a population of energetic electrons already existed over extended regions of the cluster.

A high concentration of ultra-high-energy cosmic ray (UHECR) events, called a hotspot, was reported by the Telescope Array (TA) experiment, but its origin still remains unsolved. One of the obstacles is that there is no astronomical object, which could be the source, behind the TA hotpot. In an effort to understand the origin of the TA hotspot, we suggested a model based on the magnetized cosmic web structure. The UHECRs were produced from sources in the Virgo cluster and were initially confined by cluster magnetic fields for a certain period. Next, some of them preferentially escaped to and propagated along filaments. Eventually, they were scattered by filament magnetic fields, and come to us. To examine the model, we followed the propagation trajectories of UHE protons in a simulated universe with clusters, filaments, and voids, by employing a number of models for cosmic magnetic fields. In this study, we present some of the initial results, such as the ratio between the particles directly escaping from the clusters to the voids and particles escaping from the clusters to the filaments. We also discuss the feasibility of our model for the origin of the hotspot by examining the trajectories of the UHE protons.

The existence of microgauss magnetic fields in galaxy clusters have been established through observations of synchrotron radiation and Faraday rotation. They are conjectured to be generated via small-scale dynamo by turbulent flow motions in the intracluster medium (ICM). Some of giant radio relics, on the other hand, show the structures of synchrotron polarization vectors, organized over the scales of $\sim$ Mpc, challenging the turbulence origin of cluster magnetic fields. Unlike turbulence in the interstellar medium, turbulence in the ICM is subsonic. And it is driven sporadically in highly stratified backgrounds, when major mergers occur during the hierarchical formation of clusters. To investigate quantitatively the characteristics of turbulence dynamo in such ICM environment, we performed a set of turbulence simulations using a high-order-accurate, magnetohydrodynamic (MHD) code. We find that turbulence dynamo could generate the cluster magnetic fields up to the observed level from the primordial seed fields of $10^{-15}$ G or so within the age of the universe, if the MHD description of the ICM could be extended down to $\sim$ kpc scales. However, highly organized structures of polarization vectors, such as those observed in the Sausage relic, are difficult to be reproduced by the shock compression of turbulence-generated magnetic fields. This implies that the modeling of giant radio relics may require the pre-existing magnetic fields organized over $\sim$ Mpc scales.

Low sonic Mach number shocks form in the intracluster medium (ICM) during the formation of the large-scale structure of the universe. Nonthermal cosmic-ray (CR) protons are expected to be accelerated via diffusive shock acceleration (DSA) in those ICM shocks, although observational evidence for the $\gamma$-ray emission of hadronic origin from galaxy clusters has yet to be established. Considering the results obtained from recent plasma simulations, we improve the analytic test-particle DSA model for weak quasi-parallel ($Q_\parallel$) shocks, previously suggested by \citet{kang2010}. In the model CR spectrum, the transition from the postshock thermal to CR populations occurs at the injection momentum, $p_{\rm inj}$, above which protons can undergo the full DSA process. As the shock energy is transferred to CR protons, the postshock gas temperature should decrease accordingly and the subshock strength weakens due to the dynamical feed of the CR pressure to the shock structure. This results in the reduction of the injection fraction, although the postshock CR pressure approaches an asymptotic value when the CR spectrum extends to the relativistic regime. Our new DSA model self-consistently accounts for such behaviors and adopts better estimations for $p_{\rm inj}$. With our model DSA spectrum, the CR acceleration efficiency ranges $\eta\sim10^{-3}-0.01$ for supercritical, $Q_\parallel$-shocks with sonic Mach number $2.25\lesssim M_{\rm s}\lesssim5$ in the ICM. Based on \citet{ha2018b}, on the other hand, we argue that proton acceleration would be negligible in subcritical shocks with $M_{\rm s}<2.25$.

Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$\beta$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encounters and accretion, along with violent feedback from galaxies and relativistic jets from active galactic nuclei, drive winds, gravity waves, turbulence and shocks within the ICM. Those dynamics, in turn, generate cluster-scale magnetic fields and accelerate and mediate the transport of high-energy charged particles. Kinetic-scale, collective plasma processes define the basic character and fundamental signatures of these ICM phenomena, which are observed primarily by X-ray and radio astronomers.

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.

Shock waves play an important role in turbulent astrophysical media by compressing the gas and dissipating the turbulent energy into the thermal energy. We here study shocks in magnetohydrodynamic turbulence using high-resolution simulations. Turbulent Mach numbers of $\mathcal{M}_\text{turb}=0.5-7$ and initial magnetic fields of plasma beta $\beta_0=0.1 - 10$ are considered, targeting turbulences in interstellar and intracluster media. Specifically, we present the statistics of fast and slow shocks, such as the distribution of shock Mach numbers ($M_\text{s}$) and the energy dissipation at shocks, based on refined methodologies for their quantifications. While most shocks form with low $M_\text{s}$, strong shocks follow exponentially decreasing distributions of $M_\text{s}$. More shocks appear for larger $\mathcal{M}_\text{turb}$ and larger $\beta_0$. Fast shock populations dominate over slow shocks if $\beta_0\gg1$, but substantial populations of slow shocks develop in the cases of $\beta\lesssim1$, i.e., strong background fields. The shock dissipation of turbulent energy occurs preferentially at fast shocks with $M_\text{s}\lesssim$ a few to several, and the dissipation at strong shocks shows exponentially decreasing functions of $M_\text{s}$. The energy dissipation at shocks, normalized to the energy injection, $\epsilon_\text{shock}/\epsilon_\text{inj}$, is estimated to be in the range of $\sim0.1-0.5$, except for the case of $\mathcal{M}_\text{turb}=0.5$ and $\beta_0=0.1$ where the shock dissipation is negligible. The fraction decreases with $\mathcal{M}_\text{turb}$; it is close to $\sim0.4-0.6$ for $\mathcal{M}_\text{turb}=0.5$, while it is $\sim0.1-0.25$ for $\mathcal{M}_\text{turb}=7$. The rest of the turbulent energy is expected to dissipate through the turbulent cascade.

Giant radio relics in the outskirts of galaxy clusters are known to be lit up by the relativistic electrons produced via diffusive shock acceleration (DSA) in shocks with low sonic Mach numbers, $M_{\rm s}\lesssim3$. The particle acceleration at these collisionless shocks critically depends on the kinetic plasma processes that govern the injection to DSA. Here, we study the preacceleration of suprathermal electrons in weak, quasi-perpendicular ($Q_\perp$) shocks in the hot, high-$\beta$ ($\beta = P_{\rm gas}/P_{\rm B}$) intracluster medium (ICM) through two-dimensional particle-in-cell simulations. \citet{guo2014a,guo2014b} showed that in high-$\beta$ $Q_\perp$-shocks, some of incoming electrons could be reflected upstream and gain energy via shock drift acceleration (SDA). The temperature anisotropy due to the SDA-energized electrons then induces the electron firehose instability (EFI), and oblique waves are generated, leading to a Fermi-like process and multiple cycles of SDA in the preshock region. We find that such electron preacceleration is effective only in shocks above a critical Mach number $M_{\rm ef}^*\approx2.3$. This means that in ICM plasmas, $Q_\perp$-shocks with $M_{\rm s}\lesssim2.3$ may not efficiently accelerate electrons. We also find that even in $Q_\perp$-shocks with $M_{\rm s}\gtrsim2.3$, electrons may not reach high enough energies to be injected to the full Fermi-I process of DSA, because long-wavelength waves are not developed via the EFI alone. Our results indicate that additional electron preaccelerations are required for DSA in ICM shocks, and the presence of fossil relativistic electrons in the shock upstream region may be necessary to explain observed radio relics.

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.

Ultra-high-energy cosmic rays (UHECRs) are known to come from outside of our Galaxy, but their origin still remains unknown. The Telescope Array (TA) experiment recently identified a high concentration in the arrival directions of UHECRs with energies above $5.7 \times 10^{19} eV$, called hotspot. We here report the presence of filaments of galaxies, connected to the Virgo Cluster, on the sky around the hotspot, and a statistically significant correlation between hotspot events and the filaments. With 5-year TA data, the maximum significance of binomial statistics for the correlation is estimated to be 6.1 $\sigma$ at correlation angle 3.4 degree. The probability that the above significance appears by chance is $\sim 2.0 \times 10^{-8}$ (5.6 $\sigma$). Based on this finding, we suggest a model for the origin of TA hotspot UHECRs; they are produced at sources in the Virgo Cluster, and escape to and propagate along filaments, before they are scattered toward us. This picture requires the filament magnetic fields of strength $\gtrsim 20$ nG, which need to be confirmed in future observations.

Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a fifth order weighted essentially non-oscillatory scheme for constrained-transport magnetohydrodynamics into the code WOMBAT. We establish the correctness of our implementation with an extensive number tests. We find that the fifth order scheme performs as accurately as a common second order scheme at half the resolution. We argue that for a given solution quality the new scheme is more computationally efficient than lower order schemes in three dimensions. We also establish the performance characteristics of the solver in the WOMBAT framework. Our implementation fully vectorizes using flattened arrays in thread-local memory. It performs at about 0.6 Million zones per second per node on Intel Broadwell. We present scaling tests of the code up to 98 thousand cores on the Cray XC40 machine "Hazel Hen", with a sustained performance of about 5 percent of peak at scale.

In galaxy clusters, modern radio interferometers observe non-thermal radio sources with unprecedented spatial and spectral resolution. For the first time, the new data allows to infer the structure of the intra-cluster magnetic fields on small scales via Faraday tomography. This leap forward demands new numerical models for the amplification of magnetic fields in cosmic structure formation - the cosmological magnetic dynamo. Here we present a novel numerical approach to astrophyiscal MHD simulations aimed to resolve this small-scale dynamo in future cosmological simulations. As a first step, we implement a fifth order WENO scheme in the new code WOMBAT. We show that this scheme doubles the effective resolution of the simulation and is thus less expensive than common second order schemes. WOMBAT uses a novel approach to parallelization and load balancing developed in collaboration with performance engineers at Cray Inc. This will allow us scale simulation to the exaflop regime and achieve kpc resolution in future cosmological simulations of galaxy clusters. Here we demonstrate the excellent scaling properties of the code and argue that resolved simulations of the cosmological small scale dynamo within the whole virial radius are possible in the next years.

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.

Despite progress in understanding radio relics, there are still open questions regarding the underlying particle acceleration mechanisms. In this paper we present deep 1--4 GHz VLA observations of CIZA\,J2242.8+5301 ($z=0.1921$), a double radio relic cluster characterized by small projection on the plane of the sky. Our VLA observations reveal, for the first time, the complex morphology of the diffuse sources and the filamentary structure of the northern relic. We discover new faint diffuse radio emission extending north of the main northern relic. Our Mach number estimates for the northern and southern relics, based on the radio spectral index map obtained using the VLA observations and existing LOFAR and GMRT data, are consistent with previous radio and X-ray studies ($\mathcal{M}_{\rm RN}=2.58\pm0.17$ and $\mathcal{M}_{\rm RS}=2.10\pm0.08$). However, color-color diagrams and modelings suggest a flatter injection spectral index than the one obtained from the spectral index map, indicating that projection effects might be not entirely negligible. The southern relic consists of five "arms". Embedded in it, we find a tailed radio galaxy which seems to be connected to the relic. A spectral index flattening, where the radio tail connects to the relic, is also measured. We propose that the southern relic may trace AGN fossil electrons that are re-accelerated at a shock, with an estimated strength of $\mathcal{M}=2.4$. High-resolution mapping of other tailed radio galaxies also supports a scenario where AGN fossil electrons are revived by the merger event and could be related to the formation of some diffuse cluster radio emission.

Collisionless shocks with low sonic Mach numbers, $M_{\rm s} \lesssim 4$, are expected to accelerate cosmic ray (CR) protons via diffusive shock acceleration (DSA) in the intracluster medium (ICM). However, observational evidence for CR protons in the ICM has yet to be established. Performing particle-in-cell simulations, we study the injection of protons into DSA and the early development of a nonthermal particle population in weak shocks in high $\beta$ ($\approx 100$) plasmas. Reflection of incident protons, self-excitation of plasma waves via CR-driven instabilities, and multiple cycles of shock drift acceleration are essential to the early acceleration of CR protons in supercritical quasi-parallel shocks. We find that only in ICM shocks with $M_{\rm s} \gtrsim M_{\rm s}^*\approx 2.25$, a sufficient fraction of incoming protons are reflected by the overshoot in the shock electric potential and magnetic mirror at locally perpendicular magnetic fields, leading to efficient excitation of magnetic waves via CR streaming instabilities and the injection into the DSA process. Since a significant fraction of ICM shocks have $M_{\rm s} < M_{\rm s}^*$, CR proton acceleration in the ICM might be less efficient than previously expected. This may explain why the diffuse gamma-ray emission from galaxy clusters due to proton-proton collisions has not been detected so far.

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.

Massive stars blow powerful stellar winds throughout their evolutionary stages from the main sequence to Wolf-Rayet phases. The amount of mechanical energy deposited in the interstellar medium by the wind from a massive star can be comparable to the explosion energy of a core-collapse supernova that detonates at the end of its life. In this study, we estimate the kinetic energy deposition by massive stars in our Galaxy by considering the integrated Galactic initial mass function and modeling the stellar wind luminosity. The mass loss rate and terminal velocity of stellar winds during the main sequence, red supergiant, and Wolf-Rayet stages are estimated by adopting theoretical calculations and observational data published in the literature. We find that the total stellar wind luminosity due to all massive stars in the Galaxy is about $L_w\approx 1.1\times 10^{41}$ erg/s, which is about 1/4 of the power of supernova explosions, $L_{SN} \approx 4.8\times 10^{41}$ erg/s. If we assume that $\sim 1-10$ % of the wind luminosity could be converted to Galactic cosmic rays (GCRs) through collisonless shocks such as termination shocks in stellar bubbles and superbubbles, colliding-wind shocks in binaries, and bow-shocks of massive runaway stars, stellar winds might be expected to make a significant contribution to GCR production, though lower than that of supernova remnants.

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

Non-detection of $\gamma$-ray emission from galaxy clusters has challenged diffusive shock acceleration (DSA) of cosmic-ray (CR) protons at weak collisionless shocks that are expected to form in the intracluster medium. As an effort to address this problem, we here explore possible roles of Alfv\'en waves self-excited via resonant streaming instability during the CR acceleration at parallel shocks. The mean drift of Alfv\'en waves may either increase or decrease the scattering center compression ratio, depending on the postshock cross-helicity, leading to either flatter or steeper CR spectra. We first examine such effects at planar shocks, based on the transport of Alfv\'en waves in the small amplitude limit. For the shock parameters relevant to cluster shocks, Alfv\'enic drift flattens the CR spectrum slightly, resulting in a small increase of the CR acceleration efficiency, $\eta$. We then consider two additional, physically motivated cases: (1) postshock waves are isotropized via MHD and plasma processes across the shock transition and (2) postshock waves contain only forward waves propagating along with the flow due to a possible gradient of CR pressure behind the shock. In these cases, Alfv\'enic drift could reduce $\eta$ by as much as a factor of 5 for weak cluster shocks. For the canonical parameters adopted here, we suggest $\eta\sim10^{-4}-10^{-2}$ for shocks with sonic Mach number $M_{\rm s}\approx2-3$. The possible reduction of $\eta$ may help ease the tension between non-detection of $\gamma$-rays from galaxy clusters and DSA predictions.

X-ray shocks and radio relics detected in the cluster outskirts are commonly interpreted as shocks induced by mergers of sub-clumps. We study the properties of merger shocks in merging galaxy clusters, using a set of cosmological simulations for the large-scale structure formation of the universe. As a representative case, we here focus on the simulated clusters that undergo almost head-on collisions with mass ratio $\sim2$. Due to the turbulent nature of the intracluster medium, shock surfaces are not smooth, but composed of shocks with different Mach numbers. As the merger shocks expand outward from the core to the outskirts, the average Mach number, $\left<M_s\right>$, increases in time. We suggest that the shocks propagating along the merger axis could be manifested as X-ray shocks and/or radio relics. The kinetic energy through the shocks, $F_\phi$, peaks at $\sim1$ Gyr after their initial launching, or at $\sim1-2$ Mpc from the core. Because of the Mach number dependent model adopted here for the cosmic ray (CR) acceleration efficiency, their CR-energy-weighted Mach number is higher with $\left< M_s \right>_{\rm CR}\sim3-4$, compared to the kinetic-energy-weighted Mach number, $\left<M_s\right>_{\phi}\sim2-3$. Most energetic shocks are to be found ahead of the lighter dark matter (DM) clump, while the heavier DM clump is located in the opposite side of clusters. Although our study is limited to the merger case considered, the results such as the means and variations of shock properties and their time evolution could be compared with the observed characteristics of merger shocks, constraining interpretations of relevant observations.

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.

A recent study by Posselt et al. (2017) reported the deepest X-ray investigation of the Geminga pulsar wind nebula (PWN) by using \emph{Chandra X-ray Observatory}. In comparison with previous studies of this system, a number of new findings have been reported and we found these suggest the possible variabilities in various components of this PWN. This motivates us to carry out a dedicated search for the morphological and spectral variations of this complex nebula. We have discovered variabilities on timescales from a few days to a few months from different components of the nebula. The fastest change occurred in the circumstellar environment at a rate of 80 per cent of the speed of light. One of the most spectacular results is the wiggling of a half light-year long tail as an extension of the jet, which is significantly bent by the ram pressure. The jet wiggling occurred at a rate of about 20 per cent of the speed of light. This twisted structure can possibly be a result of a propagating torsional Alf\`{v}en wave. We have also found evidence of spectral hardening along this tail for a period of about nine months.

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.

According to structure formation simulations, weak shocks with typical Mach number, $M_{\rm s}\lesssim 3$, are expected to form in merging galaxy clusters. The presence of such shocks has been indicated by X-ray and radio observations of many merging clusters. In particular, diffuse radio sources known as radio relics could be explained by synchrotron-emitting electrons accelerated via diffusive shock acceleration (Fermi I) at quasi-perpendicular shocks. Here we also consider possible roles of stochastic acceleration (Fermi II) by compressive MHD turbulence downstream of the shock. Then we explore a puzzling discrepancy that for some radio relics, the shock Mach number inferred from the radio spectral index is substantially larger than that estimated from X-ray observations. This problem could be understood, if shock surfaces associated with radio relics consist of multiple shocks with different strengths.In that case, X-ray observations tend to pick up the part of shocks with lower Mach numbers and higher kinetic energy flux, while radio emissions come preferentially from the part of shocks with higher Mach numbers and higher cosmic ray (CR) production. We also show that the Fermi I reacceleration model with preexisting fossil electrons supplemented by Fermi II acceleration due to postshock turbulence could reproduce observed profiles of radio flux densities and integrated radio spectra of two giant radio relics. This study demonstrates the CR electrons can be accelerated at collisionless shocks in galaxy clusters just like supernova remnant shock in the interstellar medium and interplanetary shocks in the solar wind.

Faraday tomography allows astronomers to probe the distribution of magnetic field along the line of sight (LOS), but that can be achieved only after Faraday spectrum is interpreted. However, the interpretation is not straightforward, mainly because Faraday spectrum is complicated due to turbulent magnetic field; it ruins the one-to-one relation between the Faraday depth and the physical depth, and appears as many small-scale features in Faraday spectrum. In this paper, employing "simple toy models" for the magnetic field, we describe numerically as well as analytically the characteristic properties of Faraday spectrum. We show that Faraday spectrum along "multiple loss" can be used to extract the global properties of magnetic field. Specifically, considering face-on spiral galaxies and modeling turbulent magnetic field as a random field with single coherence length, we numerically calculate Faraday spectrum along a number of LOSs and its shape-characterizing parameters, that is, the moments. When multiple LOSs cover a region of $\gtrsim (10\ {\rm coherence\ length)^2}$, the shape of Faraday spectrum becomes smooth and the shape-characterizing parameters are well specified. With the Faraday spectrum constructed as a sum of Gaussian functions with different means and variances, we analytically show that the parameters are expressed in terms of the regular and turbulent components of LOS magnetic field and the coherence length. We also consider the turbulent magnetic field modeled with power-law spectrum, and study how the magnetic field is revealed in Faraday spectrum. Our work suggests a way toward obtaining the information of magnetic field from Faraday tomography study.

Although many of the observed properties of giant radio relics detected in the outskirts of galaxy clusters can be explained by relativistic electrons accelerated at merger-driven shocks, significant puzzles remain. In the case of the so-called Toothbrush relic, the shock Mach number estimated from X-ray observations ($M_{\rm X}\approx1.2-1.5$) is substantially weaker than that inferred from the radio spectral index ($M_{\rm rad}\approx2.8$).Toward understanding such a discrepancy, we here consider the following diffusive shock acceleration (DSA) models:(1) weak-shock models with $M_{\rm s}\lesssim 2$ and a preexisting population of cosmic-ray electrons (CRe) with a flat energy spectrum,and (2) strong-shock models with $M_{\rm s}\approx3$ and either shock-generated suprathermal electrons or preexisting fossil CRe. We calculate the synchrotron emission from the accelerated CRe, following the time evolution of the electron DSA, and subsequent radiative cooling and postshock turbulent acceleration (TA). We find that both models could reproduce reasonably well the observed integrated radio spectrum of the Toothbrush relic, but the observed broad transverse profile requires the stochastic acceleration by downstream turbulence, which we label "turbulent acceleration" or TA to distinguish it from DSA. Moreover, to account for the almost uniform radio spectral index profile along the length of the relic, the weak-shock models require a preshock region over 400~kpc with a uniform population of preexisting CRe with a high cutoff energy ($\gtrsim 40$ GeV). Due to the short cooling time, it is challenging to explain the origin of such energetic electrons. Therefore, we suggest the strong-shock models with low-energy seed CRe ($\lesssim 150$~MeV) are preferred for the radio observations of this relic.

The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics.

Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $\mu$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization of the Universe.

On the largest scales, the Universe consists of voids and filaments making up the cosmic web. Galaxy clusters are located at the knots in this web, at the intersection of filaments. Clusters grow through accretion from these large-scale filaments and by mergers with other clusters and groups. In a growing number of galaxy clusters, elongated Mpc-size radio sources have been found, so-called radio relics. These relics are thought to trace relativistic electrons in the intracluster plasma accelerated by low-Mach number collisionless shocks generated by cluster-cluster merger events. A long-standing problem is how low-Mach number shocks can accelerate electrons so efficiently to explain the observed radio relics. Here we report on the discovery of a direct connection between a radio relic and a radio galaxy in the merging galaxy cluster Abell 3411-3412. This discovery indicates that fossil relativistic electrons from active galactic nuclei are re-accelerated at cluster shocks. It also implies that radio galaxies play an important role in governing the non-thermal component of the intracluster medium in merging clusters.

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.

Turbulence is a key ingredient for the evolution of the intracluster medium, whose properties can be predicted with high resolution numerical simulations. We present initial results on the generation of solenoidal and compressive turbulence in the intracluster medium during the formation of a small-size cluster using highly resolved, non-radiative cosmological simulations, with a refined monitoring in time. In this first of a series of papers, we closely look at one simulated cluster whose formation was distinguished by a merger around $z \sim 0.3$. We separate laminar gas motions, turbulence and shocks with dedicated filtering strategies and distinguish the solenoidal and compressive components of the gas flows using Hodge-Helmholtz decomposition. Solenoidal turbulence dominates the dissipation of turbulent motions ($\sim 95\%$) in the central cluster volume at all epochs. The dissipation via compressive modes is found to be more important ($\sim 30 \%$ of the total) only at large radii ($\geq 0.5 ~r_{\rm vir}$) and close to merger events. We show that enstrophy (vorticity squared) is good proxy of solenoidal turbulence. All terms ruling the evolution of enstrophy (i.e. baroclinic, compressive, stretching and advective terms) are found to be significant, but in amounts that vary with time and location. Two important trends for the growth of enstrophy in our simulation are identified: first, enstrophy is continuously accreted into the cluster from the outside, and most of that accreted enstrophy is generated near the outer accretion shocks by baroclinic and compressive processes. Second, in the cluster interior vortex stretching is dominant, although the other terms also contribute substantially.

We simulate shock-free and shocked viscous accretion flow onto a black hole in a two dimensional cylindrical geometry, where initial conditions were chosen from analytical solutions. The simulation code used the Lagrangian Total Variation Diminishing (LTVD) and remap routine, which enabled us to attain high accuracy in capturing shocks and to handle the angular momentum distribution correctly. Inviscid shock-free accretion disk solution produced a thick disk structure, while the viscous shock-free solution attained a Bondi-like structure, but in either case, no jet activity nor any QPO-like activity developed. The steady state shocked solution in the inviscid, as well as, in the viscous regime, matched theoretical predictions well. However, increasing viscosity renders the accretion shock unstable. Large amplitude shock oscillation is accompanied by intermittent, transient inner multiple shocks. Such oscillation of the inner part of disk is interpreted as the source of QPO in hard X-rays observed in micro-quasars. Strong shock oscillation induces strong episodic jet emission. The jets also showed existence of shocks, which are produced as one shell hits the preceding one. The periodicity of jets and shock oscillation were similar. The jets for higher viscosity parameter are evidently stronger and faster.

We examine the proposal that the dispersion measures (DMs) and Faraday rotation measures (RMs) of extragalactic linearly-polarized fast radio bursts (FRBs) can be used to probe the intergalactic magnetic field (IGMF) in filaments of galaxies. The DM through the cosmic web is dominated by contributions from the warm-hot intergalactic medium (WHIM) in filaments and from the gas in voids. On the other hand, RM is induced mostly by the hot medium in galaxy clusters, and only a fraction of it is produced in the WHIM. We show that if one excludes FRBs whose sightlines pass through galaxy clusters, the line-of-sight strength of the IGMF in filaments, $B_{||}$, is approximately $C(\langle 1+z \rangle/f_{DM})(RM/DM)$, where $C$ is a known constant. Here, the redshift of the FRB is not required to be known; $f_{DM}$ is the fraction of total DM due to the WHIM, while $\langle 1+z \rangle$ is the redshift of interevening gas weighted by the WHIM gas density, both of which can be evaluated for a given cosmology model solely from the DM of an FRB. Using data on structure formation simulations and a model IGMF, we show that $C(\langle 1+z \rangle/f_{DM})(RM/DM)$ closely reproduces the density-weighted line-of-sight strength of the IGMF in filaments of the large-scale structure.

Most of the observed features of radio {\it gischt} relics such as spectral steepening across the relic width and power-law-like integrated spectrum can be adequately explained by diffusive shock acceleration (DSA) model, in which relativistic electrons are (re-)accelerated at shock waves induced in the intracluster medium. However, the steep spectral curvature in the integrated spectrum above $\sim 2$ GHz detected in some radio relics such as the Sausage relic in cluster CIZA J2242.8+5301 may not be interpreted by simple radiative cooling of postshock electrons. In order to understand such steepening, we here consider a model in which a spherical shock sweeps through and then exits out of a finite-size cloud with fossil relativistic electrons. The ensuing integrated radio spectrum is expected to steepen much more than predicted for aging postshock electrons, since the re-acceleration stops after the cloud-crossing time. Using DSA simulations that are intended to reproduce radio observations of the Sausage relic, we show that both the integrated radio spectrum and the surface brightness profile can be fitted reasonably well, if a shock of speed, $u_s \sim 2.5-2.8\times 10^3 \kms$, and sonic Mach number, $M_s \sim 2.7-3.0$, traverses a fossil cloud for $\sim 45$ Myr and the postshock electrons cool further for another $\sim 10$ Myr. This attempt illustrates that steep curved spectra of some radio gischt relics could be modeled by adjusting the shape of the fossil electron spectrum and adopting the specific configuration of the fossil cloud.

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.

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.

Radio relics detected in the outskirts of galaxy clusters are thought to trace radio-emitting relativistic electrons accelerated at cosmological shocks. In this study, using the cosmological hydrodynamic simulation data for the large-scale structure formation and adopting a diffusive shock acceleration (DSA) model for the production of cosmic-ray (CR) electrons, we construct mock radio and X-ray maps of simulated galaxy clusters that are projected in the sky plane. Various properties of shocks and radio relics, including the shock Mach number, radio spectral index and luminosity are extracted from the synthetic maps and compared with observations. A substantial fraction of radio and X-ray shocks identified in these maps involve multiple shock surfaces along line of sights (LoSs), and the morphology of shock distributions in the maps depends on the projection direction. Among multiple shocks in a given LoS, radio observations tend to pick up stronger shocks with flatter radio spectra, while X-ray observations preferentially select weaker shocks with larger kinetic energy flux. As a result, the shock Mach numbers and locations derived from radio and X-ray observations could differ from each other in some cases. We also find that the distributions of the spectral index and radio power of the synthetic radio relics are somewhat inconsistent with those of observed real relics; a bit more radio relics have been observed closer to the cluster core and with steeper spectral indices. We suggest the inconsistency could be explained, if very weak shocks with $M_s \la 2$ accelerate CR electrons more efficiently, compared to the DSA model adopted here.

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 analyze high resolution simulations of compressible, MHD turbulence with properties resembling conditions in galaxy clusters. The flow is driven to turbulence Mach number $\mathcal{M}_t \sim 1/2$ in an isothermal medium with an initially very weak, uniform seed magnetic field ($\beta = P_g/P_B = 10^6$). Since cluster turbulence is likely to result from a mix of sheared (solenoidal) and compressive forcing processes, we examine the distinct turbulence properties for both cases. In one set of simulations velocity forcing is entirely solenoidal ($\nabla\cdot \delta {\vec u} = 0$), while in the other it is entirely compressive ($\nabla\times \delta {\vec u} = 0$). Both cases develop a mixture of solenoidal and compressive turbulent motions, since each generates the other. The development of compressive turbulent motions leads to shocks, even when the turbulence is solenoidally forced and subsonic. Shocks, in turn, produce and amplify vorticity, which is especially important in compressively forced turbulence. To clarify those processes we include a pair of appendices that look in detail at vorticity evolution in association with shocks. From our simulation analyses we find that magnetic fields amplified to near saturation levels in predominantly solenoidal turbulence can actually enhance vorticity on small scales by concentrating and stabilizing shear. The properties, evolution rates and relative contributions of the kinetic and magnetic turbulent elements depend strongly on the character of the forcing. Specifically, shocks are stronger, but vorticity evolution and magnetic field amplification are slower and weaker when the turbulence is compressively forced. We identify a simple relation to estimate characteristic shock strengths in terms of the turbulence Mach number and the character of the forcing. Our results will be helpful in understanding flow motions in galaxy clusters.

In order to understand certain observed features of arc-like giant radio relics such as the rareness, uniform surface brightness, and curved integrated spectra, we explore a diffusive shock acceleration (DSA) model for radio relics in which a spherical shock impinges on a magnetized cloud containing fossil relativistic electrons. Toward this end, we perform DSA simulations of spherical shocks with the parameters relevant for the Sausage radio relic in cluster CIZA J2242.8+5301, and calculate the ensuing radio synchrotron emission from re-accelerated electrons. Three types of fossil electron populations are considered: a delta-function like population with the shock injection momentum, a power-law distribution, and a power-law with an exponential cutoff. The surface brightness profile of radio-emitting postshock region and the volume-integrated radio spectrum are calculated and compared with observations. We find that the observed width of the Sausage relic can be explained reasonably well by shocks with speed $u_s \sim 3\times 10^3 \kms$ and sonic Mach number $M_s \sim 3$. These shocks produce curved radio spectra that steepen gradually over $(0.1-10) \nu_{\rm br}$ with break frequency $ \nu_{\rm br}\sim 1$ GHz, if the duration of electron acceleration is $\sim 60 - 80$ Myr. However, the abrupt increase of spectral index above $\sim 1.5$ GHz observed in the Sausage relic seems to indicate that additional physical processes, other than radiative losses, operate for electrons with $\gamma_e \gtrsim 10^4$.

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.

Deep surveys with the SKA1-MID array offer for the first time the opportunity to systematically explore the polarization properties of the microJy source population. Our knowledge of the polarized sky approaching these levels is still very limited. In total intensity the population will be dominated by star-forming and normal galaxies to intermediate redshifts ($z \sim1-2$), and low-luminosity AGN to high redshift. The polarized emission from these objects is a powerful probe of their intrinsic magnetic fields and of their magnetic environments. For redshift of order 1 and above the broad bandwidth of the mid-bands span the Faraday thick and thin regimes allowing study of the intrinsic polarization properties of these objects as well as depolarization from embedded and foreground plasmas. The deep field polarization images will provide Rotation Measures data with very high solid angle density allowing a sensitive statistical analysis of the angular variation of RM on critical arc-minute scales from a magnetic component of Large Scale Structure of the Universe.

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 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 Faraday dispersion function (FDF), which can be derived from an observed polarization spec- trum by Faraday rotation measure synthesis, is a profile of polarized emissions as a function of Faraday depth. We study intrinsic FDFs along sight lines through face-on, Milky-Way-like galaxies by means of a sophisticated galactic model incorporating 3D MHD turbulence, and investigate how much the FDF contains information intrinsically. Since the FDF reflects distributions of thermal and cosmic- ray electrons as well as magnetic fields, it has been expected that the FDF could be a new probe to examine internal structures of galaxies. We, however, find that an intrinsic FDF along a sight line through a galaxy is very complicated, depending significantly on actual configurations of turbulence. We perform 800 realizations of turbulence, and find no universal shape of the FDF even if we fix the global parameters of the model. We calculate the probability distribution functions of the standard deviation, skewness, and kurtosis of FDFs and compare them for models with different global pa- rameters. Our models predict that the presence of vertical magnetic fields and large scale-height of cosmic-ray electrons tend to make the standard deviation relatively large. Contrastingly, differences in skewness and kurtosis are relatively less significant.

Rotation measure (RM) grids of extragalactic radio sources have been widely used for studying cosmic magnetism. But their potential for exploring the intergalactic magnetic field (IGMF) in filaments of galaxies is unclear, since other Faraday-rotation media such as the radio source itself, intervening galaxies, and the interstellar medium of our Galaxy are all significant contributors. We study statistical techniques for discriminating the Faraday rotation of filaments from other sources of Faraday rotation in future large-scale surveys of radio polarization. We consider a 30 degree times 30 degree field-of-view toward the south Galactic pole, while varying the number of sources detected in both present and future observations. We select sources located at high redshifts and toward which depolarization and optical absorption systems are not observed, so as to reduce the RM contributions from the sources and intervening galaxies. It is found that a high-pass filter can satisfactorily reduce the RM contribution from the Galaxy, since the angular scale of this component toward high Galactic latitudes would be much larger than that expected for the IGMF. Present observations do not yet provide a sufficient source density to be able to estimate the RM of filaments. However, from the proposed approach with forthcoming surveys, we predict significant residuals of RM that should be ascribable to filaments. The predicted structure of the IGMF down to scales of 0.1 degree should be observable with data from the SKA, if we achieve selections of sources toward which sight lines do not contain intervening galaxies and RM errors are less than a few rad m-2.

We consider a phenomenological model for the thermal leakage injection in the diffusive shock acceleration (DSA) process, in which suprathermal protons and electrons near the shock transition zone are assumed to have the so-called $\kappa$-distributions produced by interactions of background thermal particles with pre-existing and/or self-excited plasma/MHD waves or turbulence. The $\kappa$-distribution has a power-law tail, instead of an exponential cutoff, well above the thermal peak momentum. So there are a larger number of potential seed particles with momentum, above that required for participation in the DSA process. As a result, the injection fraction for the $\kappa$-distribution depends on the shock Mach number much less severely compared to that for the Maxwellian distribution. Thus, the existence of $\kappa$-like suprathermal tails at shocks would ease the problem of extremely low injection fractions, especially for electrons and especially at weak shocks such as those found in the intracluster medium. We suggest that the injection fraction for protons ranges $10^{-4}-10^{-3}$ for a $\kappa$-distribution with $10 < \kappa_p < 30$ at quasi-parallel shocks, while the injection fraction for electrons becomes $10^{-6}-10^{-5}$ for a $\kappa$-distribution with $\kappa_e < 2$ at quasi-perpendicular shocks. For such $\kappa$ values the ratio of cosmic ray electrons to protons naturally becomes $K_{e/p}\sim 10^{-3}-10^{-2}$, which is required to explain the observed ratio for Galactic cosmic rays.

The outskirts of galaxy clusters are continuously disturbed by mergers and gas infall along filaments, which in turn induce turbulent flow motions and shock waves. We examine the properties of shocks that form within $r_{200}$ in sample galaxy clusters from structure formation simulations. While most of these shocks are weak and inefficient accelerators of cosmic rays (CRs), there are a number of strong, energetic shocks which can produce large amounts of CR protons via diffusive shock acceleration. We show that the energetic shocks reside mostly in the outskirts and a substantial fraction of them are induced by infall of the warm-hot intergalactic medium from filaments. As a result, the radial profile of the CR pressure in the intracluster medium is expected to be broad, dropping off more slowly than that of the gas pressure, and might be even temporarily inverted, peaking in the outskirts. The volume-integrated momentum spectrum of CR protons inside $r_{200}$ has the power-law slope of $4.25 - 4.5$, indicating that the average Mach number of the shocks of main CR production is in the range of $\left< M_s \right>_{\rm CR} \approx 3 - 4$. We suggest that some radio relics with relatively flat radio spectrum could be explained by primary electrons accelerated by energetic infall shocks with $M_s >~ 3$ induced in the cluster outskirts.

Unveiling the intergalactic magnetic field (IGMF) in filaments of galaxies is a very important and challenging subject in modern astronomy. In order to probe the IGMF from rotation measures (RMs) of extragalactic radio sources, we need to separate RMs due to other origins such as the source, intervening galaxies, and our Galaxy. In this paper, we discuss observational strategies for the separation by means of Faraday tomography (Faraday RM Synthesis). We consider an observation of a single radio source such as a radio galaxy or a quasar viewed through the Galaxy and the cosmic web. We then compare the observation with another observation of a neighbor source with a small angular separation. Our simulations with simple models of the sources suggest that it would be not easy to detect the RM due to the IGMF of order ~ 1 rad/m/m, an expected value for the IGMF through a single filament. Contrary to it, we find that the RM of at least ~10 rad/m/m could be detected with the SKA or its pathfinders/precursors, if we achieve selections of ideal sources. These results would be improved if we incorporate decomposition techniques such as RMCLEAN and QU-fitting. We discuss feasibility of the strategies for cases with complex Galactic emissions as well as with effects of observational noise and radio frequency interferences.

Wu et al. (2009) suggested an empirical relation between the magnetic field strength along the line of sight (LOS) and the dispersion of Faraday rotation measure (RM) distribution in turbulent media with root-mean-square sonic Mach number $M_s \simeq 1$. In this paper, we extend the work by incorporating the Mach number dependence. Media with $\sim 0.5 < M_s < \sim 2$ are considered to cover the Mach number range of the warm ionized medium (WIM) of our Galaxy. Three-dimensional, magnetohydrodynamic isothermal turbulence simulations with solenoidal forcing are used. We suggest a new relation among the LOS magnetic field strength, the dispersion of RM distribution, and the Mach number, which approximately represents the relation for Alfv\'enic perturbations. In addition, we suggest a relation between the Mach number and the dispersion of log-normal distribution of emission measure (EM), which is basically the relation for the Mach number and the density dispersion. The relations could be used for a quick and rough estimation of the LOS magnetic field strength in the turbulent WIM.

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.

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.

The thermal state of the fluid is governed by the ratio of the thermal and the rest energy. This brings the composition of the fluid into the picture. Although, fluid composed of lighter particles (e.g: electron-positron pair plasma) at same temperature, is more relativistic compared to fluids with finite baryon loading, but this is not necessarily true when baryon poor transonic fluid are compared with each other. It can be shown that the transonic pair-fluid is the least relativistic. This has far reaching consequences on accreting flows around compact objects and are expected to have similar effect on relativistic outflows and explosive events as well.

We investigate the capability of ongoing radio telescopes for probing Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) in the large-scale structure of the universe which is expected to be of order $O(1) {\rm rad/m^2}$. We consider polarization observations of a compact radio source such as quasars behind a diffuse source such as the Galaxy, and calculate Stokes parameters $Q$ and $U$ assuming a simple model of the Faraday dispersion functions with Gaussian shape. Then, we perform the Fisher analysis to estimate the expected errors in the model parameters from QU-fitting of polarization intensity, accounting for sensitivities and frequency bands of Australian Square Kilometer Array Pathfinder, Low Frequency Array, and the Giant Meterwave Radio Telescope. Finally, we examine the condition on the source intensities which are required to detect the IGMF. Our analysis indicates that the QU-fitting is promising for forthcoming wideband polarimetry to explore RM due to the IGMF in filaments of galaxies.

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.