Radiative transfer describes the propagation of electromagnetic radiation through an interacting medium. This process is often simulated by the use of the Monte Carlo method, which involves the probabilistic determination and tracking of simulated photon packages. In the regime of high optical depths, this approach encounters difficulties since a proper representation of the various physical processes can only be achieved by considering high numbers of simulated photon packages. As a consequence, the demand for computation time rises accordingly and thus practically puts a limit on the optical depth of models that can be simulated. Here we present a method that aims to solve the problem of high optical depths in dusty media, which relies solely on the use of unbiased Monte Carlo radiative transfer. For that end, we identified and precalculated repeatedly occuring and simulated processes, stored their outcome in a multidimensional cumulative distribution function, and immediately replaced the basic Monte Carlo transfer during a simulation by that outcome. During the precalculation, we generated emission spectra as well as deposited energy distributions of photon packages traveling from the center of a sphere to its rim. We carried out a performance test of the method to confirm its validity and gain a boost in computation speed by up to three orders of magnitude. We then applied the method to a simple model of a viscously heated circumstellar disk, and we discuss the necessity of finding a solution for the optical depth problem with regard to a proper temperature calculation. We find that the impact of an incorrect treatment of photon packages in highly optically thick regions extents even to optically thin regions, thus, changing the overall observational appearance of the disk.

Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.

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.

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.

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

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.

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.

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

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.