We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least in most part from AGN activity. The AGN fractional luminosities are higher than those in either sub-millimeter galaxies, or AGN selected via other mid-infrared criteria. Adopting the $G$, M$_{20}$ and $A$ morphological parameters, together with traditional classification boundaries, infers that three quarters of the sample as mergers. Our sample do not, however, show any correlation between the considered morphological parameters and either infrared luminosity or AGN fractional luminosity. Moreover, their asymmetries and effective radii are distributed identically to those of massive galaxies at $z\sim2$. We conclude that our sample is not preferentially associated with mergers, though a significant merger fraction is still plausible. Instead, we propose that our sample are examples of the massive galaxy population at $z\sim2$ that harbor a briefly luminous, "flickering" AGN, and in which the $G$ and M$_{20}$ values have been perturbed, due to either the AGN, and/or the earliest formation stages of a bulge in an inside-out manner. Furthermore, we find that the mass assembly of the central black holes in our sample leads the mass assembly of any bulge component. Finally, we speculate that our sample represent a small fraction of the immediate antecedents of compact star-forming galaxies at $z\sim2$.

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.

We analyze the three-year SDSS-II Superernova (SN) Survey data and identify a sample of 1070 photometric SN Ia candidates based on their multi-band light curve data. This sample consists of SN candidates with no spectroscopic confirmation, with a subset of 210 candidates having spectroscopic redshifts of their host galaxies measured, while the remaining 860 candidates are purely photometric in their identification. We describe a method for estimating the efficiency and purity of photometric SN Ia classification when spectroscopic confirmation of only a limited sample is available, and demonstrate that SN Ia candidates from SDSS-II can be identified photometrically with ~91% efficiency and with a contamination of ~6%. Although this is the largest uniform sample of SN candidates to date for studying photometric identification, we find that a larger spectroscopic sample of contaminating sources is required to obtain a better characterization of the background events. A Hubble diagram using SN candidates with no spectroscopic confirmation, but with host galaxy spectroscopic redshifts, yields a distance modulus dispersion that is only ~20 - 40% larger than that of the spectroscopically-confirmed SN Ia sample alone with no significant bias. A Hubble diagram with purely photometric classification and redshift-distance measurements, however, exhibit biases that require further investigation for precision cosmology.

A common problem in ultra-high energy cosmic ray physics is the comparison of energy spectra. The question is whether the spectra from two experiments or two regions of the sky agree within their statistical and systematic uncertainties. We develop a method to directly compare energy spectra for ultra-high energy cosmic rays from two different regions of the sky in the same experiment without reliance on agreement with a theoretical model of the energy spectra. The consistency between the two spectra is expressed in terms of a Bayes factor, defined here as the ratio of the likelihood of the two-parent source hypothesis to the likelihood of the one-parent source hypothesis. Unlike other methods, for example chi^2 tests, the Bayes factor allows for the calculation of the posterior odds ratio and correctly accounts for non-Gaussian uncertainties. The latter is particularly important at the highest energies, where the number of events is very small.

We present the first principal component analysis (PCA) applied to a sample of 119 Spitzer Infrared Spectrograph (IRS) spectra of local ultraluminous infrared galaxies (ULIRGs) at z<0.35. The purpose of this study is to objectively and uniquely characterise the local ULIRG population using all information contained in the observed spectra. We have derived the first three principal components (PCs) from the covariance matrix of our dataset which account for over 90% of the variance. The first PC is characterised by dust temperatures and the geometry of the mix of source and dust. The second PC is a pure star formation component. The third PC represents an anti-correlation between star formation activity and a rising AGN. Using the first three PCs, we are able to accurately reconstruct most of the spectra in our sample. Our work shows that there are several factors that are important in characterising the ULIRG population, dust temperature, geometry, star formation intensity, AGN contribution, etc. We also make comparison between PCA and other diagnostics such as ratio of the 6.2 microns PAH emission feature to the 9.7 micron silicate absorption depth and other observables such as optical spectral type.

We present mid-infrared spectra of six FeLoBAL QSOs at 1<z<1.8, taken with the Spitzer space telescope. The spectra span a range of shapes, from hot dust dominated AGN with silicate emission at 9.7 microns, to moderately obscured starbursts with strong Polycyclic Aromatic Hydrocarbon (PAH) emission. The spectrum of one object, SDSS 1214-0001, shows the most prominent PAHs yet seen in any QSO at any redshift, implying that the starburst dominates the mid-IR emission with an associated star formation rate of order 2700 solar masses per year. With the caveats that our sample is small and not robustly selected, we combine our mid-IR spectral diagnostics with previous observations to propose that FeLoBAL QSOs are at least largely comprised of systems in which (a) a merger driven starburst is ending, (b) a luminous AGN is in the last stages of burning through its surrounding dust, and (c) which we may be viewing over a restricted line of sight range.

Upcoming large-scale ground- and space- based supernova surveys will face a challenge identifying supernova candidates largely without the use of spectroscopy. Over the past several years, a number of supernova identification schemes have been proposed that rely on photometric information only. Some of these schemes use color-color or color-magnitude diagrams; others simply fit supernova data to models. Both of these approaches suffer a number of drawbacks partially addressed in the so-called Bayesian-based supernova classification techniques. However, Bayesian techniques are also problematic in that they typically require that the supernova candidate be one of a known set of supernova types. This presents a number of problems, the most obvious of which is that there are bound to be objects that do not conform to any presently known model in large supernova candidate samples. We propose a new photometric classification scheme that uses a Bayes factor based on color in order to identify supernovae by type. This method does not require knowledge of the complete set of possible astronomical objects that could mimic a supernova signal. Further, as a Bayesian approach, it accounts for all systematic and statistical uncertainties of the measurements in a single step. To illustrate the use of the technique, we apply it to a simulated dataset for a possible future large-scale space-based Joint Dark Energy Mission and demonstrate how it could be used to identify Type Ia supernovae. The method's utility in pre-selecting and ranking supernova candidates for possible spectroscopic follow-up -- i.e., its usage as a supernova trigger -- will be briefly discussed.

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

Technical reports on operations and features of the Pierre Auger Observatory, including ongoing and planned enhancements and the status of the future northern hemisphere portion of the Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the correlations of ultra-high energy cosmic ray directions with extra-Galactic objects, of general anisotropy, of photons and neutrinos, and of other astrophysical effects, with the Pierre Auger Observatory. Contributions to the 31st ICRC, Lodz, Poland, July 2009.

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory.

Reports on the atmospheric monitoring, calibration, and other operating systems of the Pierre Auger Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers. Submissions to the 31st ICRC, Lodz, Poland (July 2009).

We apply methods from Bayesian inferencing and graph theory to a dataset of 102 mid-infrared spectra, and archival data from the optical to the millimeter, to construct an evolutionary paradigm for z<0.4 infrared-luminous galaxies (ULIRGs). We propose that the ULIRG lifecycle consists of three phases. The first phase lasts from the initial encounter until approximately coalescence. It is characterized by homogeneous mid-IR spectral shapes, and IR emission mainly from star formation, with a contribution from an AGN in some cases. At the end of this phase, a ULIRG enters one of two evolutionary paths depending on the dynamics of the merger, the available quantities of gas, and the masses of the black holes in the progenitors. On one branch, the contributions from the starburst and the AGN to the total IR luminosity decline and increase respectively. The IR spectral shapes are heterogeneous, likely due to feedback from AGN-driven winds. Some objects go through a brief QSO phase at the end. On the other branch, the decline of the starburst relative to the AGN is less pronounced, and few or no objects go through a QSO phase. We show that the 11.2 micron PAH feature is a remarkably good diagnostic of evolutionary phase, and identify six ULIRGs that may be archetypes of key stages in this lifecycle.

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

Searches for statistically significant correlations between arrival directions of ultra-high energy cosmic rays and classes of astrophysical objects are common in astroparticle physics. We present a method to test potential correlation signals of a priori unknown strength and evaluate their statistical significance sequentially, i.e., after each incoming new event in a running experiment. The method can be applied to data taken after the test has concluded, allowing for further monitoring of the signal significance. It adheres to the likelihood principle and rigorously accounts for our ignorance of the signal strength.

We present a new measurement of the volumetric rate of Type Ia supernova up to a redshift of 1.7, using the Hubble Space Telescope (HST) GOODS data combined with an additional HST dataset covering the North GOODS field collected in 2004. We employ a novel technique that does not require spectroscopic data for identifying Type Ia supernovae (although spectroscopic measurements of redshifts are used for over half the sample); instead we employ a Bayesian approach using only photometric data to calculate the probability that an object is a Type Ia supernova. This Bayesian technique can easily be modified to incorporate improved priors on supernova properties, and it is well-suited for future high-statistics supernovae searches in which spectroscopic follow up of all candidates will be impractical. Here, the method is validated on both ground- and space-based supernova data having some spectroscopic follow up. We combine our volumetric rate measurements with low redshift supernova data, and fit to a number of possible models for the evolution of the Type Ia supernova rate as a function of redshift. The data do not distinguish between a flat rate at redshift > 0.5 and a previously proposed model, in which the Type Ia rate peaks at redshift >1 due to a significant delay from star-formation to the supernova explosion. Except for the highest redshifts, where the signal to noise ratio is generally too low to apply this technique, this approach yields smaller or comparable uncertainties than previous work.

Air fluorescence detectors measure the energy of ultra-high energy cosmic rays by collecting fluorescence light emitted from nitrogen molecules along the extensive air shower cascade. To ensure a reliable energy determination, the light signal needs to be corrected for atmospheric effects, which not only attenuate the signal, but also produce a non-negligible background component due to scattered Cherenkov light and multiple-scattered light. The correction requires regular measurements of the aerosol attenuation length and the aerosol phase function, defined as the probability of light scattered in a given direction. At the Pierre Auger Observatory in Malargue, Argentina, the phase function is measured on an hourly basis using two Aerosol Phase Function (APF) light sources. These sources direct a UV light beam across the field of view of the fluorescence detectors; the phase function can be extracted from the image of the shots in the fluorescence detector cameras. This paper describes the design, current status, standard operation procedure, and performance of the APF system at the Pierre Auger Observatory.

The air fluorescence detectors (FDs) of the Pierre Auger Observatory are vital for the determination of the air shower energy scale. To compensate for variations in atmospheric conditions that affect the energy measurement, the Observatory operates an array of monitoring instruments to record hourly atmospheric conditions across the detector site, an area exceeding 3,000 square km. This paper presents results from four instruments used to characterize the aerosol component of the atmosphere: the Central Laser Facility (CLF), which provides the FDs with calibrated laser shots; the scanning backscatter lidars, which operate at three FD sites; the Aerosol Phase Function monitors (APFs), which measure the aerosol scattering cross section at two FD locations; and the Horizontal Attenuation Monitor (HAM), which measures the wavelength dependence of aerosol attenuation.

FRAM - F/(Ph)otometric Robotic Atmospheric Monitor is the latest addition to the atmospheric monitoring instruments of the Pierre Auger Observatory. An optical telescope equipped with CCD camera and photometer, it automatically observes a set of selected standard stars and a calibrated terrestrial source. Primarily, the wavelength dependence of the attenuation is derived and the comparison between its vertical values (for stars) and horizontal values (for the terrestrial source) is made. Further, the integral vertical aerosol optical depth can be obtained. A secondary program of the instrument, the detection of optical counterparts of gamma-ray bursts, has already proven successful. The hardware setup, software system, data taking procedures, and first analysis results are described in this paper.

The Pierre Auger Observatory in Malargue, Argentina, is designed to study the origin of ultrahigh energy cosmic rays with energies above 10^18 eV. The energy calibration of the detector is based on a system of four air fluorescence detectors. To obtain reliable calorimetric information from the fluorescence stations, the atmospheric conditions at the experiment's site need to be monitored continuously during operation. One of the components of the observatory's atmospheric monitoring system is a set of four elastic backscatter lidar stations, one station at each of the fluorescence detector sites. This paper describes the design, current status, standard operation procedure, and performance of the lidar system of the Pierre Auger Observatory.

This paper presents a novel method for determining the probability that a supernova candidate belongs to a known supernova type (such as Ia, Ibc, IIL, \emph{etc.}), using its photometric information alone. It is validated with Monte Carlo, and both space- and ground- based data. We examine the application of the method to well-sampled as well as poorly sampled supernova light curves and investigate to what extent the best currently available supernova models can be used for typing supernova candidates. Central to the method is the assumption that a supernova candidate belongs to a group of objects that can be modeled; we therefore discuss possible ways of removing anomalous or less well understood events from the sample. This method is particularly advantageous for analyses where the purity of the supernova sample is of the essence, or for those where it is important to know the number of the supernova candidates of a certain type (\emph{e.g.}, in supernova rate studies).

The current measurements of the cosmic ray energy spectrum at ultra-high energies ($\text{E}>10^{19}$ eV) are characterized by large systematic errors and poor statistics. In addition, the experimental results of the two experiments with the largest published data sets, AGASA and HiRes, appear to be inconsistent with each other, with AGASA seeing an unabated continuation of the energy spectrum even at energies beyond the GZK cutoff energy at $10^{19.6}$ eV. Given the importance of the related astrophysical questions regarding the unknown origin of these highly energetic particles, it is crucial that the extent to which these measurements disagree be well understood. Here we evaluate the consistency of the two measurements for the first time with a model-independent method that accounts for the large statistical and systematic errors of current measurements. We further compare the AGASA and HiRes spectra with the recently presented Auger spectrum. The method directly compares two measurements, bypassing the introduction of theoretical models for the shape of the energy spectrum. The inconsistency between the observations is expressed in terms of a Bayes Factor, a standard statistic defined as the ratio of a separate parent source hypothesis to a single parent source hypothesis. Application to the data shows that the two-parent hypothesis is disfavored. We expand the method to allow comparisons between an experimental flux and that predicted by any model.