We present a multi-wavelength study of the active galactic nucleus in the nearby ($D=14.1$ Mpc) low mass galaxy IC 750, which has circumnuclear 22 GHz water maser emission. The masers trace a nearly edge-on, warped disk $\sim$0.2 pc in diameter, coincident with the compact nuclear X-ray source which lies at the base of the $\sim$kpc-scale extended X-ray emission. The position-velocity structure of the maser emission indicates the central black hole (BH) has a mass less than $1.4 \times 10^5~M_\odot$. Keplerian rotation curves fitted to these data yield enclosed masses between $4.1 \times 10^4~M_\odot$ and $1.4 \times 10^5~M_\odot$, with a mode of $7.2 \times 10^4~M_\odot$. Fitting the optical spectrum, we measure a nuclear stellar velocity dispersion $\sigma_* = 110.7^{+12.1}_{-13.4}$~{\rm km~s}$^{-1}.$ From near-infrared photometry, we fit a bulge mass of $(7.3 \pm 2.7) \times 10^8~M_\odot$ and a stellar mass of $1.4 \times 10^{10}~M_\odot$. The mass upper limit of the intermediate mass black hole in IC 750 falls roughly two orders of magnitude below the $M_{\rm BH}-\sigma_*$ relation and roughly one order of magnitude below the $M_{\rm BH}-M_{\rm Bulge}$ and $M_{\rm BH}-M_*$ relations -- larger than the relations' intrinsic scatters of (0.58 $\pm$ 0.09) dex, 0.69 dex, and (0.65 $\pm$ 0.09) dex, respectively. These offsets could be due to larger scatter at the low mass end of these relations. Alternatively, black hole growth is intrinsically inefficient in galaxies with low bulge and/or stellar masses, which causes the black holes to be under-massive relative to their hosts, as predicted by some galaxy evolution simulations.

We present a new investigation of the dependence of H2O maser detection rates and properties on the mid-IR AGN luminosity, L_AGN, and the obscuring column density, N_H, based on mid-IR and hard X-ray photometry. Based on spectral energy distribution fitting that allows for decomposition of the black hole accretion and star-formation components in the mid-infrared, we show that the megamaser (disk maser) detection rate increases sharply for galaxies with 12 micron AGN luminosity L^{AGN}_{12 micron} greater than 10^42 erg/s, from ~<3%(~<2%) to ~12%(~5%). By using the ratio of the observed X-ray to mid-IR AGN luminosity as an indicator of N_H, we also find that the megamaser (disk maser) detection rates are boosted to 15%(7%) and 20%(9%) for galaxies with N_H >= 10^23 cm^{-2} and N_H >= 10^{24} cm^{-2}, respectively. Combining these column density cuts with a constraint for high L^{AGN}_{12 micron} (>=10^42 erg/s) predicts further increases in the megamaser (disk maser) detection rates to 19%(8%) and 27%(14%), revealing unprecedented potential boosts of the megamaser and disk maser detection rates by a factor of 7-15 relative to the current rates, depending on the chosen sample selection criteria. A noteworthy aspect of these new predictions is that the completeness rates are only compromised mildly, with the rates remaining at the level of ~95%(~50%) for sources with N_H >= 10^{23} cm^{-2} (N_H >= 10^24 cm^-2). Applying these selection methods to the current X-ray AGN surveys predicts the detection of >~15 new megamaser disks.

We have constructed an all-sky catalog of optical AGNs with $z < 0.09$, based on optical spectroscopy, from the parent sample of galaxies in the 2MASS Redshift Survey (2MRS), a near-complete census of the nearby universe. Our catalog consists of 1929 broad line AGNs, and 6562 narrow line AGNs which satisfy the \citet{Kauffmann03} criteria, of which 3607 also satisfy the \citet{Kewley01} criteria. We also report emission line widths, fluxes, flux errors, and signal-to-noise ratios of all the galaxies in our spectroscopic sample, allowing users to customize the selection criteria. Although we uniformly processed the spectra of galaxies from a homogeneous parent sample, inhomogeneities persist due to the differences in the quality of the obtained spectra, taken with different instruments, and the unavailability of spectra for $\sim$20\% of the galaxies. We quantify how the differences in spectral quality affect not only the AGN detection rates but also broad line to narrow line AGN ratios. We find that the inhomogeneities primarily stem from the continuum signal-to-noise (S/N) in the spectra near the emission lines of interest. We fit for the AGN fraction as a function of continuum S/N and assign AGN likelihoods to galaxies which were not identified as AGNs using the available spectra. This correction results in a catalog suitable for statistical studies. This work also paves the way for a truly homogeneous and complete nearby AGN catalog by identifying galaxies whose AGN status needs to be verified with higher quality spectra, quantifying the spectral quality necessary to do so.

Water megamasers from circumnuclear disks in galaxy centers provide the most accurate measurements of supermassive black hole masses and uniquely probe the sub-parsec accretion processes. At the same time, these systems offer independent crucial constraints of the Hubble Constant in the nearby universe, and thus, the arguably best single constraint on the nature of dark energy. The chances of finding these golden standards are however abysmally low, at an overall =< 3% for any level of water maser emission detected at 22 GHz, and =< 1% for those exhibiting the disk-like configuration. We provide here a thorough summary of the current state of the detection of water megamaser disks, along with a novel investigation of the likelihood of increasing their detection rates based on a multivariate parameter analysis of optical and mid-infrared photometric properties of the largest database of galaxies surveyed for the 22 GHz emission. We find that galaxies with water megamaser emission tend to associate with strong emission in all WISE mid-infrared wavelengths, with the strongest enhancement in the W4 band, at 22 micron, as well as with previously proposed and newly found indicators of AGN strength in the mid-infrared, such as red W1-W2 and W1-W4 colors, and the integrated mid-infrared luminosity of the host galaxy. These trends offer a potential boost of the megamaser detection rates to 6 - 15%, or a factor of 2-8 relative to the current rates, depending on the chosen sample selection criteria, while fostering real chances for discovering >= 20 new megamaser disks.

We have conducted a study to quantify the systematic differences resulting from using different stellar population models in optical spectroscopic identification of type II AGN. We examined the different AGN detection fractions of 7069 nearby galaxies (z <= 0.09) with SDSS DR8 spectra when using the Bruzual & Charlot (2003, BC03), Vazdekis et al. (2010, MILES), and solar metallicity Maraston and Stromback (2011) (MS11solar) stellar population models. The line fluxes obtained using BC03 and MS11solar are publicly available from SDSS data releases. We find that the BC03 templates result in systematically higher BPT line ratios and consequently higher AGN fractions and the MS11solar templates result in systematically lower line ratios and AGN fractions compared with the MILES templates. Using MILES as the standard, BC03 results in 25% "false positives" and MS11solar results in 22% "false negatives" when using the Kewley et al. (2001a) boundary for AGN identification. The fraction of galaxies whose AGN identification changes for different templates is luminosity dependent, ranging from a few percent for L[OIII]5007 >= 10^40 erg s-1 and increasing to ~ 50% for L[OIII]5007 <= 10^38 erg s-1. These results suggest that template choice should be accounted for when using and comparing the AGN and emission line fluxes from different catalogs.

We present the second and final data release of the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). Data are presented for 63 new galaxies not included in the first data release, and we provide 2D emission-line fitting products for the full S7 sample of 131 galaxies. The S7 uses the WiFeS instrument on the ANU 2.3m telescope to obtain spectra with a spectral resolution of R = 7000 in the red (540 - 700 nm) and R = 3000 in the blue (350 - 570 nm), over an integral field of 25x38 arcsec^2 with 1x1 arcsec^2 spatial pixels. The S7 contains both the largest sample of active galaxies and the highest spectral resolution of any comparable integral field survey to date. The emission-line fitting products include line fluxes, velocities and velocity dispersions across the WiFeS field of view, and an artificial neural network has been used to determine the optimal number of Gaussian kinematic components for emission-lines in each spaxel. Broad Balmer lines are subtracted from the spectra of nuclear spatial pixels in Seyfert 1 galaxies before fitting the narrow lines. We bin nuclear spectra and measure reddening-corrected nuclear fluxes of strong narrow lines for each galaxy. The nuclear spectra are classified on optical diagnostic diagrams, where the strength of the coronal line [FeVII]6087 is shown to be correlated with [OIII]/Hbeta. Maps revealing gas excitation and kinematics are included for the entire sample, and we provide notes on the newly observed objects.

The most rapidly evolving regions of galaxies often display complex optical spectra with emission lines excited by massive stars, shocks and accretion onto supermassive black holes. Standard calibrations (such as for the star formation rate) cannot be applied to such mixed spectra. In this paper we isolate the contributions of star formation, shock excitation and active galactic nucleus (AGN) activity to the emission line luminosities of individual spatially resolved regions across the central 3 $\times$ 3 kpc$^2$ region of the active barred spiral galaxy NGC$\sim$613. The star formation rate and AGN luminosity calculated from the decomposed emission line maps are in close agreement with independent estimates from data at other wavelengths. The star formation component traces the B-band stellar continuum emission, and the AGN component forms an ionization cone which is aligned with the nuclear radio jet. The optical line emission associated with shock excitation is cospatial with strong $H_2$ and [Fe II] emission and with regions of high ionized gas velocity dispersion ($\sigma > 100$ km s$^{-1}$). The shock component also traces the outer boundary of the AGN ionization cone and may therefore be produced by outflowing material interacting with the surrounding interstellar medium. Our decomposition method makes it possible to determine the properties of star formation, shock excitation and AGN activity from optical spectra, without contamination from other ionization mechanisms.

The optical spectra of Seyfert galaxies are often dominated by emission lines excited by both star formation and AGN activity. Standard calibrations (such as for the star formation rate) are not applicable to such composite (mixed) spectra. In this paper, we describe how integral field data can be used to spectrally and spatially separate emission associated with star formation from emission associated with accretion onto an active galactic nucleus (AGN). We demonstrate our method using integral field data for two AGN host galaxies (NGC 5728 and NGC 7679) from the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). The spectra of NGC 5728 and NGC 7679 form clear sequences of AGN fraction on standard emission line ratio diagnostic diagrams. We show that the emission line luminosities of the majority (> 85 per cent) of spectra along each AGN fraction sequence can be reproduced by linear superpositions of the emission line luminosities of one AGN dominated spectrum and one star formation dominated spectrum. We separate the Halpha, Hbeta, [N II]$\lambda$6583, [S II]$\lambda \lambda$6716, 6731, [O III]$\lambda$5007 and [O II]$\lambda \lambda$3726, 3729 luminosities of every spaxel into contributions from star formation and AGN activity. The decomposed emission line images are used to derive the star formation rates and AGN bolometric luminosities for NGC 5728 and NGC 7679. Our calculated values are mostly consistent with independent estimates from data at other wavelengths. The recovered star forming and AGN components also have distinct spatial distributions which trace structures seen in high resolution imaging of the galaxies, providing independent confirmation that our decomposition has been successful.

We investigate the relative significance of radiation pressure and gas pressure in the extended narrow line regions (ENLRs) of four Seyfert galaxies from the integral field Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). We demonstrate that there exist two distinct types of starburst-AGN mixing curves on standard emission line diagnostic diagrams which reflect the balance between gas pressure and radiation pressure in the ENLR. In two of the galaxies the ENLR is radiation pressure dominated throughout and the ionization parameter remains constant (log U ~ 0). In the other two galaxies radiation pressure is initially important, but gas pressure becomes dominant as the ionization parameter in the ENLR decreases from log U ~ 0 to -3.4 <= log U <= -3.2. Where radiation pressure is dominant, the AGN regulates the density of the interstellar medium on kpc scales and may therefore have a direct impact on star formation activity and/or the incidence of outflows in the host galaxy to scales far beyond the zone of influence of the black hole. We find that both radiation pressure dominated and gas pressure dominated ENLRs are dynamically active with evidence for outflows, indicating that radiation pressure may be an important source of AGN feedback even when it is not dominant over the entire ENLR.

Aims. Study the connection between the masing disk and obscuring torus in Seyfert 2 galaxies. Methods. We present a uniform X-ray spectral analysis of the high energy properties of 14 nearby megamaser Active Galactic Nuclei observed by NuSTAR. We use a simple analytical model to localize the maser disk and understand its connection with the torus by combining NuSTAR spectral parameters with available physical quantities from VLBI mapping. Results. Most of the sources analyzed are heavily obscured, showing a column density in excess of $\sim 10^{23}$ cm$^{-2}$. In particular, $79\%$ are Compton-thick ($N_{\rm H} > 1.5 \times 10^{24}$ cm$^{-2}$). Using column densities measured by NuSTAR, with the assumption that the torus is the extension of the maser disk, and further assuming a reasonable density profile, the torus dimensions can be predicted. They are found to be consistent with mid-IR interferometry parsec-scale observations of Circinus and NGC 1068. In this picture, the maser disk is intimately connected to the inner part of the torus. It is probably made of a large number of molecular clouds connecting the torus and the outer part of the accretion disk, giving rise to a thin disk rotating in most cases in Keplerian or sub-Keplerian motion. This toy model explains the established close connection between water megamaser emission and nuclear obscuration as a geometric effect.

Studies discerning whether there is a significant correlation between UHECR arrival directions and optical AGN are hampered by the lack of a uniformly selected and complete all-sky optical AGN catalog. To remedy this, we are preparing such a catalog based on the 2MASS Redshift Survey (2MRS), a spectroscopic sample of $\sim 44,500$ galaxies complete to a K magnitude of 11.75 over 91% of the sky. We have analyzed the available optical spectra of these 2MRS galaxies ($\sim 80$% of the galaxies), in order to identify the AGN amongst them with uniform criteria. We present a first-stage release of the AGN catalog for the southern sky, based on spectra from the 6dF Galaxy survey and CTIO telescope. Providing a comparably uniform and complete catalog for the northern sky is more challenging because the spectra for the northern galaxies were taken with different instruments.

The Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7) is a targeted survey probing the narrow-line regions (NLRs) of a representative sample of ~140 nearby (z<0.02) Seyfert galaxies by means of optical integral field spectroscopy. The survey is based on a homogeneous data set observed using the Wide Field Spectrograph WiFeS. The data provide a 25x38 arcsec$^2$ field-of-view around the galaxy centre at typically ~1.5 arcsec spatial resolution and cover a wavelength range between ~3400 - 7100 $\AA$ at spectral resolutions of ~100 km s$^{-1}$ and ~50 km s$^{-1}$ in the blue and red parts, respectively. The survey is primarily designed to study gas excitation and star formation around AGN, with a special focus on the shape of the AGN ionising continuum, the interaction between radio jets and the NLR gas, and the nature of nuclear LINER emission. We provide an overview of the current status of S7-based results and present new results for NGC 6300.

We present Wide Field Spectrograph (WiFeS) integral field spectroscopy and HST FOS spectroscopy for the LINER galaxy NGC 1052. We infer the presence of a turbulent accretion flow forming a small-scale accretion disk. We find a large-scale outflow and ionisation cone along the minor axis of the galaxy. Part of this outflow region is photoionised by the AGN, and shares properties with the ENLR of Seyfert galaxies, but the inner ($R \lesssim 1.0$~arcsec) accretion disk and the region around the radio jet appear shock excited. The emission line properties can be modelled by a "double shock" model in which the accretion flow first passes through an accretion shock in the presence of a hard X-ray radiation, and the accretion disk is then processed through a cocoon shock driven by the overpressure of the radio jets. This model explains the observation of two distinct densities ($\sim10^4$ and $\sim10^6$ cm$^{-3}$), and provides a good fit to the observed emission line spectrum. We derive estimates for the velocities of the two shock components and their mixing fractions, the black hole mass, the accretion rate needed to sustain the LINER emission and derive an estimate for the jet power. Our emission line model is remarkably robust against variation of input parameters, and so offers a generic explanation for the excitation of LINER galaxies, including those of spiral type such as NGC 3031 (M81).

Here we describe the \emph{Siding Spring Southern Seyfert Spectroscopic Snapshot Survey} (S7) and present results on 64 galaxies drawn from the first data release. The S7 uses the Wide Field Spectrograph (WiFeS) mounted on the ANU 2.3m telescope located at the Siding Spring Observatory to deliver an integral field of $38\times25$~ arcsec at a spectral resolution of $R=7000$ in the red ($530-710$nm), and $R=3000$ in the blue ($340-560$nm). {From these data cubes we have extracted the Narrow Line Region (NLR) spectra from a 4 arc sec aperture centred on the nucleus. We also determine the H$\beta$ and [OIII]~$\lambda$5007 fluxes in the narrow lines, the nuclear reddening, the reddening-corrected relative intensities of the observed emission lines, and the H$\beta$ and \lOIII\ luminosities {determined from spectra for which the stellar continuum has been removed.} We present a set of images of the galaxies in [OIII]~$\lambda$5007, [NII]~$\lambda$6584 and H$\alpha$ which serve to delineate the spatial extent of the extended narrow line region (ENLR) and {\bf also to} reveal the structure and morphology of the surrounding \HII\ regions. Finally, we provide a preliminary discussion of those Seyfert~1 and Seyfert~2 galaxies which display coronal emission lines in order to explore the origin of these lines.

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.

Here we present the first results from the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7) which aims to investigate the physics of ~140 radio-detected southern active Galaxies with z<0.02 through Integral Field Spectroscopy using the Wide Field Spectrograph (WiFeS). This instrument provides data cubes of the central 38 x 25 arc sec. of the target galaxies in the waveband 340-710nm with the unusually high resolution of R=7000 in the red (530-710nm), and R=3000 in the blue (340-560nm). These data provide the morphology, kinematics and the excitation structure of the extended narrow-line region, probe relationships with the black hole characteristics and the host galaxy, measures host galaxy abundance gradients and the determination of nuclear abundances from the HII regions. From photoionisation modelling, we may determine the shape of the ionising spectrum of the AGN, discover whether AGN metallicities differ from nuclear abundances determined from HII regions, and probe grain destruction in the vicinity of the AGN. Here we present some preliminary results and modelling of both Seyfert galaxies observed as part of the survey.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

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.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

We report on Chandra X-ray observations of possible-AGNs which have been correlated with Ultra-high Energy Cosmic Rays (UHECRs) observed by the Pierre Auger Collaboration. Combining our X-ray observations with optical observations, we conclude that one-third of the 21 Veron-Cetty Veron (VCV) galaxies correlating with UHECRs in the first Auger data-release are actually not AGNs. We review existing optical observations of the 20 VCV galaxies correlating with UHECRs in the second Auger data-release and determine that three of them are not AGNs and two are uncertain. Overall, of the 57 published UHECRs with |b|>10 degrees, 22 or 23 correlate with true AGNs using the Auger correlation parameters. We also measured the X-ray luminosity of ESO139-G12 to complete the determination of the bolometric luminosities of AGNs correlating with UHECRs in the first data-set. Apart from two candidate sources which require further observation, we determined bolometric luminosities for the candidate galaxies of the second dataset. We find that only two of the total of 69 published UHECRs correlate with AGNs (IC5135 and IC4329a) which are powerful enough in their steady-state to accelerate protons to the observed energies of their correlated UHECRs. The GZK expectation is that about 45% of the sources of UHECRs above 60 EeV should be contained within the z<0.018 volume defined by the Auger scan analysis, so an observed level of 30-50% correlation with weak AGNs is compatible with the suggestion that AGNs experience transient high-luminosity states during which they accelerate UHECRs.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.