-
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.
-
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.
-
The cosmic ray flux measured by the Telescope Array Low Energy Extension
(TALE) exhibits three spectral features: the knee, the dip in the $10^{16}$ eV
decade, and the second knee. Here the spectrum has been measured for the first
time using fluorescence telescopes, which provide a calorimetric,
model-independent result. The spectrum appears to be a rigidity-dependent
cutoff sequence, where the knee is made by the hydrogen and helium portions of
the composition, the dip comes from the reduction in composition from helium to
metals, the rise to the second knee occurs due to intermediate range nuclei,
and the second knee is the iron knee.
-
We report on a measurement of the cosmic ray energy spectrum by the Telescope
Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air
fluorescence detector is also sensitive to the Cherenkov light produced by
shower particles. Low energy cosmic rays, in the PeV energy range, are
detectable by TALE as "Cherenkov Events". Using these events, we measure the
energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100
PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This
allows for the extension of the measurement to energies greater than a few EeV.
In this paper, we will describe the detector, explain the technique, and
present results from a measurement of the spectrum using $\sim 1000$ hours of
observation. The observed spectrum shows a clear steepening near $10^{17.1}$
eV, along with an ankle-like structure at $10^{16.2}$ eV. These features
present important constraints on galactic cosmic rays origin and propagation
models. The feature at $10^{17.1}$ eV may also mark the end of the galactic
cosmic rays flux and the start of the transition to extra-galactic sources.
-
An intermediate-scale energy spectrum anisotropy has been found in the
arrival directions of ultra-high energy cosmic rays of energies above
$10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the
Telescope Array surface detector. A relative energy distribution test is done
comparing events inside oversampled spherical caps of equal exposure, to those
outside, using the Poisson likelihood ratio. The center of maximum significance
is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies
$10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$
eV. The post-trial probability of this energy anisotropy, appearing by chance
anywhere on an isotropic sky, is found by Monte Carlo simulation to be
$9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).
-
The energy of the ultrahigh energy spectral cutoff was measured, integrating
over the northern hemisphere sky, by the Telescope Array (TA) collaboration, to
be $10^{19.78 \pm 0.06}$ eV, in agreement with the High Resolution Fly's Eye
(HiRes) experiment, whereas the Pierre Auger experiment, integrating over the
southern hemisphere sky, measured the cutoff to be at 10$^{19.62 \pm 0.02}$ eV.
An 11% energy scale difference between the TA and Auger does not account for
this difference. However, in comparing the spectra of the Telescope Array and
Pierre Auger experiments in the band of declination common to both experiments
($-15.7^{\circ} < \delta < 24.8^{\circ}$) we have found agreement in the energy
of the spectral cutoff. While the Auger result is essentially unchanged, the TA
cutoff energy has changed to $10^{19.59 \pm 0.06}$ eV. In this paper we argue
that this is an astrophysical effect.
-
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.
-
Dynamic and thermal processes regulate the structure of the multi-phase
interstellar medium (ISM), and ultimately establish how galaxies evolve through
star formation. Thus, to constrain ISM models and better understand the
interplay of these processes, it is of great interest to measure the thermal
pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158
$\mu$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH
sample, we have measured thermal pressures in 534 predominantly atomic regions
with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal
pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a
sub-sample of regions with conditions similar to those of the diffuse, neutral
gas in the Galactic plane, we find thermal pressures that follow a log-normal
distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$.
These results are consistent with thermal pressure measurements using other
observational methods. We find that $P_{\rm th}$ increases with radiation field
strength and star formation activity, as expected from the close link between
the heating of the gas and the star formation rate. Our thermal pressure
measurements fall in the regime where a two-phase ISM with cold and warm
neutral medium could exist in pressure equilibrium. Finally, we find that the
midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical
weight of the overlying ISM, consistent with results from hydrodynamical
simulations of self-regulated star formation in galactic disks.
-
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.
-
The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized
gas in the interstellar medium. By combining data from 21 galaxies selected
from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141
spatially resolved regions with a typical size of ~1 kiloparsec, with
observations of both [NII] far-infrared lines. We measure [NII] 122/205 line
ratios in the ~0.6-6 range, which corresponds to electron gas densities
$n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in
the electron density within individual galaxies can be as a high as a factor of
~50, frequently with strong radial gradients. We find that $n_e$ increases as a
function of infrared color, dust-weighted mean starlight intensity, and star
formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII]
transitions is related to the ionizing photon flux, we investigate their
reliability as tracers of the star formation rate (SFR). We derive relations
between the [NII] emission and SFR in the low-density limit and in the case of
a log-normal distribution of densities. The scatter in the correlation between
[NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of
the $n_e$ distribution. For regions with $n_e$ close to or higher than the
[NII] line critical densities, the low-density limit [NII]-based SFR
calibration systematically underestimates the SFR since [NII] emission is
collisionally quenched. Finally, we investigate the relation between [NII]
emission, SFR, and $n_e$ by comparing our observations to predictions from the
MAPPINGS-III code.
-
Molecular dynamics simulations of aqueous NaCl, KCl, NaI, and KI solutions
are used to study the effects of salts on the properties of the liquid/vapor
interface. The simulations use the models which include both charge transfer
and polarization effects. Pairing and the formation of larger ion clusters
occurs both in the bulk and surface region, with a decreased tendency to form
larger clusters near the interface. An analysis of the roughness of the surface
reveals that the chloride salts, which have less tendency to be near the
surface, have a roughness that is less than pure water, while the iodide salts,
which have a greater surface affinity, have a larger roughness. This suggests
that ions away from the surface and ions near the surface affect the interface
in opposite ways.
-
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.
-
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.
-
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.
-
Many of the transformative processes in the Universe have taken place in
regions obscured by dust, and are best studied with far-IR spectroscopy. We
present the Cryogenic-Aperture Large Infrared-Submillimeter Telescope
Observatory (CALISTO), a 5-meter class, space-borne telescope actively cooled
to 4 K, emphasizing moderate-resolution spectroscopy in the crucial 35 to 600
micron band. CALISTO will enable NASA and the world to study the rise of heavy
elements in the Universe's first billion years, chart star formation and black
hole growth in dust-obscured galaxies through cosmic time, and conduct a census
of forming planetary systems in our region of the Galaxy. CALISTO will
capitalize on rapid progress in both format and sensitivity of far-IR
detectors. Arrays with a total count of a few 100,000 detector pixels will form
the heart of a suite of imaging spectrometers in which each detector reaches
the photon background limit.
This document contains a large overview paper on CALISTO, as well as six 2-3
page scientific white papers, all prepared in response to NASA's Cosmic Origins
Program Analysis Group (COPAG's) request for input on future mission concepts.
The Far-IR Science Interest Group will meet from 3-5 June 2015 with the
intention of reaching consensus on the architecture for the Far-IR Surveyor
mission. This white paper describes one of the architectures to be considered
by the community. One or more companion papers will describe alternative
architectures.
-
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 manner in which gas accretes and orbits within circumnuclear rings has
direct implications for the star formation process. In particular, gas may be
compressed and shocked at the inflow points, resulting in bursts of star
formation at these locations. Afterwards the gas and young stars move together
through the ring. In addition, star formation may occur throughout the ring, if
and when the gas reaches sufficient density to collapse under gravity. These
two scenarios for star formation in rings are often referred to as the `pearls
on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS
observations, obtained as part of the KINGFISH Open Time Key Program, along
with archival Spitzer and ground-based observations from the SINGS Legacy
project, to investigate the heating and cooling of the interstellar medium in
the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved
estimates of the stellar FUV flux available for heating, with the gas and dust
cooling derived from the FIR continuum and line emission, we show that while
star formation is indeed dominant at the inflow points in NGC 4736, additional
star formation is needed to balance the gas heating and cooling throughout the
ring. This additional component most likely arises from the general increase in
gas density in the ring over its lifetime. Our data provide strong evidence,
therefore, for a combination of the two paradigms for star formation in the
ring in NGC4736.
-
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.
-
We have searched for intermediate-scale anisotropy in the arrival directions
of ultrahigh-energy cosmic rays with energies above 57~EeV in the northern sky
using data collected over a 5 year period by the surface detector of the
Telescope Array experiment. We report on a cluster of events that we call the
hotspot, found by oversampling using 20$^\circ$-radius circles. The hotspot has
a Li-Ma statistical significance of 5.1$\sigma$, and is centered at
R.A.=146.7$^{\circ}$, Dec.=43.2$^{\circ}$. The position of the hotspot is about
19$^{\circ}$ off of the supergalactic plane. The probability of a cluster of
events of 5.1$\sigma$ significance, appearing by chance in an isotropic
cosmic-ray sky, is estimated to be 3.7$\times$10$^{-4}$ (3.4$\sigma$).
-
We 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.
-
Motived by recent ISM studies, we present high quality MMT and Gemini
spectroscopic observations of H II regions in the nearby spiral galaxies NGC
628 and NGC 2403 in order to measure their chemical abundance gradients. Using
long-slit and multi-object mask optical spectroscopy, we obtained measurements
of the temperature sensitive auroral lines [O III] {\lambda}4363 and/or [N II]
{\lambda}5755 at a strength of 4{\sigma} or greater in 11 H II regions in NGC
628 and 7 regions in NGC 2403. These observations allow us, for the first time,
to derive an oxygen abundance gradient in NGC 628 based solely on "direct"
oxygen abundances of H II regions: 12 + log(O/H) = (8.43+/-0.03) +
(-0.017+/-0.002) x Rg (dex/kpc), with a dispersion in log(O/H) of {\sigma} =
0.10 dex, from 14 regions with a radial coverage of ~2-19 kpc. This is a
significantly shallower slope than found by previous "strong-line" abundance
studies. In NGC 2403, we derive an oxygen abundance gradient of 12 + log(O/H) =
(8.48+/-0.04) + (-0.032+/-0.007) x Rg (dex/kpc), with a dispersion in log(O/H)
of {\sigma} = 0.07 dex, from 7 H II with a radial coverage of ~1-10 kpc.
Additionally, we measure the N, S, Ne, and Ar abundances. We find the N/O
ratio decreases with increasing radius for the inner disk, but reaches a
plateau past R25 in NGC 628. NGC 2403 also has a negative N/O gradient with
radius, but we do not sample the outer disk of the galaxy past R25 and so do
not see evidence for a plateau. This bi-modal pattern measured for NGC 628
indicates dominant contributions from secondary nitrogen inside of the R25
transition and dominantly primary nitrogen farther out. As expected for
{\alpha}-process elements, S/O, Ne/O, and Ar/O are consistent with constant
values over a range in oxygen abundance.
-
The metal content of a galaxy is a key property for distinguishing between
viable galaxy evolutionary scenarios, and it strongly influences many of the
physical processes in the interstellar medium. An absolute and robust
determination of extragalactic metallicities is essential in constraining
models of chemical enrichment and chemical evolution, however, current gas
phase abundance determinations from optical fine-structure lines are uncertain
to 0.8 dex as conversion of these optical line fluxes to abundances is strongly
dependent on the electron temperature of the ionized gas. In contrast, the
far-IR emission lines can be used to derive an O++ abundance that is relatively
insensitive to temperature, while the ratio of the optical to far-IR lines
provides a consistent temperature to be used in the derivation of an O$^+$
abundance. We present observations of the [O III] 88 \micron\ fine-structure
line in NGC 628 that were obtained as part of the KINGFISH program. These data
are combined with optical IFU data to derive oxygen abundances for seven HII
regions. We find the abundance of these regions to all lie between the high and
low values of strong line calibrations and in agreement with estimates that
assume temperature fluctuations are present in the HII regions.
-
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.