A Monte Carlo program for the simulation of electromagnetic cascades initiated by high-energy photons and electrons interacting with extragalactic background light (EBL) is presented. Pair production and inverse Compton scattering on EBL photons as well as synchrotron losses and deflections of the charged component in extragalactic magnetic fields (EGMF) are included in the simulation. Weighted sampling of the cascade development is applied to reduce the number of secondary particles and to speed up computations. As final result, the simulation procedure provides the energy, the observation angle, and the time delay of secondary cascade particles at the present epoch. Possible applications are the study of TeV blazars and the influence of the EGMF on their spectra or the calculation of the contribution from ultrahigh energy cosmic rays or dark matter to the diffuse extragalactic gamma-ray background. As an illustration, we present results for deflections and time-delays relevant for the derivation of limits on the EGMF.

The large neutrino fluxes emitted with a distinct flavor hierarchy from core-collapse supernovae (SNe) during the post-bounce accretion phase, offer the best opportunity to detect effects from neutrino flavor oscillations. We perform a dedicated study of the SN neutrino flavor evolution during the accretion phase, using results from recent neutrino radiation hydrodynamics simulations. In contrast to what expected in the presence of only neutrino-neutrino interactions, we find that the multi-angle effects associated with the dense ordinary matter suppress collective oscillations. This is related to the high matter densities during the accretion phase in core-collapse SNe of massive iron-core progenitors. The matter suppression implies that neutrino oscillations will start outside the neutrino transport region and therefore will have a negligible impact on the neutrino heating and the explosion dynamics. Furthermore, the possible detection of the next galactic SN neutrino signal from the accretion phase, based on the usual Mikheyev- Smirnov-Wolfenstein effect in the SN mantle and Earth matter effects, can reveal the neutrino mass hierarchy in the case that the mixing angle $\theta_{13}$ is not very small.

The usual description of self-induced neutrino flavor conversions in core collapse supernovae (SNe) is based on the dominance of the neutrino density n_nu over the net electron density n_e. However, this condition is not met during the post-bounce accretion phase, when the dense matter in a SN is piled up above the neutrinosphere. As recently pointed-out, a dominant matter term in the anisotropic SN environment would dephase the flavor evolution for neutrinos traveling on different trajectories, challenging the occurrence of the collective behavior in the dense neutrino gas. Using the results from recent long term simulations of core-collapse SN explosions, based on three flavor Boltzmann neutrino transport in spherical symmetry, we find that both the situations of complete matter suppression (when n_e >> n_nu) and matter-induced decoherence (when n_e \gtrsim n_nu) of flavor conversions are realized during the accretion phase. The matter suppression at high densities prevents any possible impact of the neutrino oscillations on the neutrino heating and hence on the dynamics of the explosion. Furthermore, it changes the interpretation of the Earth matter effect on the SN neutrino signal during the accretion phase, allowing the possibility of the neutrino mass hierarchy discrimination at not too small values of the leptonic mixing angle \theta_{13} (i.e. \sin^2{\theta}_{13} \gtrsim 10^{-3}).

We calculate the energy spectra of cosmic rays (CR) and their secondaries produced in a supernova remnant (SNR), taking into account the time-dependence of the SNR shock. We model the trajectories of charged particles as a random walk with a prescribed diffusion coefficient, accelerating the particles at each shock crossing. Secondary production by CRs colliding with gas is included as a Monte Carlo process. We find that SNRs produce less antimatter than suggested previously: The positron/electron ratio and the antiproton/proton ratio are a few percent and few $\times 10^{-5}$, respectively. Moreover, the obtained positron/electron ratio decreases with energy, while the antiproton/proton ratio rises at most by a factor of two above 10 GeV.

High energy photons from blazars can initiate electromagnetic pair cascades interacting with the extragalactic photon background. The charged component of such cascades is deflected and delayed by extragalactic magnetic fields (EGMF), reducing thereby the observed point-like flux and leading potentially to multi degree images in the GeV energy range. We calculate the fluence of 1ES 0229+200 as seen by Fermi-LAT for different EGMF profiles using a Monte Carlo simulation for the cascade development. The non-observation of 1ES 0229+200 by Fermi-LAT suggests that the EGMF fills at least 60% of space with fields stronger than {\cal O}(10^{-16}-10^{-15})G for life times of TeV activity of {\cal O}(10^2-10^4)yr. Thus the (non-) observation of GeV extensions around TeV blazars probes the EGMF in voids and puts strong constraints on the origin of EGMFs: Either EGMFs were generated in a space filling manner (e.g. primordially) or EGMFs produced locally (e.g. by galaxies) have to be efficiently transported to fill a significant volume fraction, as e.g. by galactic outflows.

It was recently realized that three-flavor effects could peculiarly modify the development of spectral splits induced by collective oscillations, for supernova neutrinos emitted during the cooling phase of a protoneutron star. We systematically explore this case, explaining how the impact of these three-flavor effects depends on the ordering of the neutrino masses. In inverted mass hierarchy, the solar mass splitting gives rise to instabilities in regions of the (anti)neutrino energy spectra that were otherwise stable under the leading two-flavor evolution governed by the atmospheric mass splitting and by the 1-3 mixing angle. As a consequence, the high-energy spectral splits found in the electron (anti)neutrino spectra disappear, and are transferred to other flavors. Imperfect adiabaticity leads to smearing of spectral swap features. In normal mass hierarchy, the three-flavor and the two-flavor instabilities act in the same region of the neutrino energy spectrum, leading to only minor departures from the two-flavor treatment.

We calculate the energy spectra of cosmic rays (CR) and their secondaries produced in a supernova remnant (SNR), taking into account the time-dependence of the SNR shock. We model the trajectories of charged particles as a random walk with a prescribed diffusioncoefficient, accelerating the particles at each shock crossing. Secondary production by CRs colliding with gas is included as a Monte Carlo process. We find that SNRs produce less antimatter than suggested previously: The positron/electron ratio and the antiproton/proton ratio are a few percent and few $\times 10^{-5}$, respectively. Both ratios do not rise with energy.

TeV gamma-rays have been observed from blazars as well as from radio galaxies like M87 and Cen A. In leptonic models, gamma-rays above the pair production threshold can escape from the ultra-relativistic jet, since large Lorentz factors reduce the background photon densities compared to those required for isotropic emission. Here we discuss an alternative scenario, where VHE photons are generated as secondaries from UHECR interaction in the AGN core. We show that TeV gamma-rays can escape from the core despite large IR and UV backgrounds. For the special case of Cen A, we study if the various existing observations from the far infra-red to the UHE range can be reconciled within this picture.

We consider the effect of non-standard neutrino interactions (NSI, for short) on the propagation of neutrinos through the supernova (SN) envelope within a three-neutrino framework and taking into account the presence of a neutrino background. We find that for given NSI parameters, with strength generically denoted by $\varepsilon_{ij}$, neutrino evolution exhibits a significant time dependence. For $|\varepsilon_{\tau\tau}|\gtrsim$ $10^{-3}$ the neutrino survival probability may become sensitive to the $\theta_{23}$ octant and the sign of $\varepsilon_{\tau\tau}$. In particular, if $\varepsilon_{\tau\tau}\gtrsim 10^{-2}$ an internal $I$-resonance may arise independently of the matter density. For typical values found in SN simulations this takes place in the same dense-neutrino region above the neutrinosphere where collective effects occur, in particular during the synchronization regime. This resonance may lead to an exchange of the neutrino fluxes entering the bipolar regime. The main consequences are (i) bipolar conversion taking place for normal neutrino mass hierarchy and (ii) a transformation of the flux of low-energy $\nu_e$, instead of the usual spectral swap.

High energy photons from blazars interact within tens of kpc with the extragalactic photon background, initiating electromagnetic pair cascades. The charged component of such cascades is deflected by extragalactic magnetic fields (EGMF), leading to halos even around initially point-like sources. We calculate the intensity profile of the resulting secondary high-energy photons for different assumptions on the initial source spectrum and the strength of the EGMF, employing also fields found earlier in a constrained simulation of structure formation including MHD processes. We find that the observation of halos around blazars like Mrk~180 probes an interesting range of EGMF strengths and acceleration models: In particular, blazar halos test if the photon energy spectrum at the source extends beyond \sim 100 TeV and how anisotropic this high energy component is emitted.

We calculated for the nearest active galactic nucleus (AGN), Centaurus A (Cen A), the flux of high energy cosmic rays and of accompanying secondary photons and neutrinos expected from hadronic interactions in the source. We used as two basic models for the generation of ultrahigh energy cosmic rays (UHECR) shock acceleration in the radio jet and acceleration in the regular electromagnetic field close to the core of the AGN, normalizing the UHECR flux to the observations of the Auger experiment. Here we compare the previously obtained photon fluxes with the recent data reported by the Fermi LAT and H.E.S.S. collaborations. In the case of the core model, we find good agreement both for the predicted spectral shape and the overall normalization between our earlier results and the H.E.S.S. observations for a primary proton energy $dN/dE\propto E^{-\alpha}$ with $\alpha\sim 2$ or smaller. A broken-power law with high-energy part $\alpha=-2.7$ leads to photon fluxes in excess of the Fermi measurements. The energy spectrum of the photon fluxes obtained by us for the jet scenario is in all cases at variance with the H.E.S.S. and Fermi observations.

We give a very brief overview of collective effects in neutrino oscillations in core collapse supernovae where refractive effects of neutrinos on themselves can considerably modify flavor oscillations, with possible repercussions for future supernova neutrino detection. We discuss synchronized and bipolar oscillations, the role of energy and angular neutrino modes, as well as three-flavor effects. We close with a short summary and some open questions.

We calculate for the nearest active galactic nucleus (AGN), Centaurus A, the flux of high energy cosmic rays and of accompanying secondary photons and neutrinos expected from hadronic interactions in the source. We use as two basic models for the generation of ultrahigh energy cosmic rays (UHECR) shock acceleration in the radio jet and acceleration in the regular electromagnetic field close to the core of the AGN. While scattering on photons dominates in scenarios with acceleration close to the core, scattering on gas becomes more important if acceleration takes place along the jet. Normalizing the UHECR flux from Centaurus A to the observations of the Auger experiment, the neutrino flux may be marginally observable in a 1 km$^3$ neutrino telescope, if a steep UHECR flux $\d N/\d E\propto E^{-\alpha}$ with $\alpha=2.7$ extends down to $10^{17}$ eV. The associated photon flux is close to or exceeds the observational data of atmospheric Cherenkov and $\gamma$-ray telescopes for $\alpha\gsim 2$. In particular, we find that already present data favour either a softer UHECR injection spectrum than $\alpha=2.7$ for Centaurus A or a lower UHECR flux than expected from the normalization to the Auger observations.

For neutrinos streaming from a supernova (SN) core, dense matter suppresses self-induced flavor transformations if the electron density n_e significantly exceeds the neutrino density n_nu in the conversion region. If n_e is comparable to n_nu one finds multi-angle decoherence, whereas the standard self-induced transformation behavior requires that in the transformation region n_nu is safely above n_e. This condition need not be satisfied in the early phase after supernova core bounce. Our new multi-angle effect is a subtle consequence of neutrinos traveling on different trajectories when streaming from a source that is not point-like.

We study three-flavor collective neutrino transformations in the dense-neutrino region above the neutrino sphere of a supernova core. We find that two-flavor conversions driven by the atmospheric mass difference and the 13-mixing angle capture the full effect if one neglects the second-order difference between the muon and tau neutrino refractive index. Including this "mu-tau matter term" provides a resonance at a density of approximately 3 x 10^7 g cm^-3 that typically causes significant modifications of the overall electron neutrino and antineutrino survival probabilities. This effect is surprisingly sensitive to deviations from maximal 23-mixing, being different for each octant.

In the dense-neutrino region at 50-400 km above the neutrino sphere in a supernova, neutrino-neutrino interactions cause large flavor transformations. We study when the multi-angle nature of the neutrino trajectories leads to flavor decoherence between different angular modes. We consider a two-flavor mixing scenario between nu_e and another flavor nu_x and assume the usual hierarchy F(nu_e)>F{antinu_e)>F(nu_x)=F(antinu_x) for the number fluxes. We define epsilon=(F(nu_e)-F(antinu_e))/(F(antinu_e)-F(antinu_x)) as a measure for the deleptonization flux which is the one crucial parameter. The transition between the quasi single-angle behavior and multi-angle decoherence is abrupt as a function of epsilon. For typical choices of other parameters, multi-angle decoherence is suppressed for epsilon>0.3, but a much smaller asymmetry suffices if the neutrino mass hierarchy is normal and the mixing angle small. The critical epsilon depends logarithmically on the neutrino luminosity. In a realistic supernova scenario, the deleptonization flux is probably enough to suppress multi-angle decoherence.

We calculate the yield of high energy neutrinos produced in astrophysical sources for arbitrary interaction depths $\tau_0$ and magnetic field strengths $B$. We take into account energy loss processes like synchrotron radiation and diffusion of charged particles in turbulent magnetic fields as well as the scattering of secondaries on background photons and the direct production of charm neutrinos. Meson-photon interactions are simulated with an extended version of the SOPHIA model. Diffusion leads to an increased path-length before protons leave the source of size R_s and therefore magnetized sources lose their transparency below the energy $E\sim 10^{18}{\rm eV} (R_s/{\rm pc}) (B/{\rm mG}) \tau_0^{1/\alpha}$, with $\alpha=1/3$ and 1 for Kolmogorov and Bohm diffusion, respectively. Moreover, the neutrino flux is suppressed above the energy where synchrotron energy losses become important for charged particles. As a consequence, the energy spectrum and the flavor composition of neutrinos are strongly modified both at low and high energies even for sources with $\tau_0\lsim 1$.

We review how a high-statistics observation of the neutrino signal from a future galactic core-collapse supernova (SN) may be used to discriminate between different neutrino mixing scenarios. Most SN neutrinos are emitted in the accretion and cooling phase, during which the flavor-dependent differences of the emitted neutrino spectra are small and rather uncertain. Therefore the discrimination between neutrino mixing scenarios using these neutrinos should rely on observables independent of the SN neutrino spectra. We discuss two complementary methods that allow for the positive identification of the mass hierarchy without knowledge of the emitted neutrino fluxes, provided that the 13-mixing angle is large, $\sin^2\theta_{13}\gg 10^{-5}$. These two approaches are the observation of modulations in the neutrino spectra by Earth matter effects or by the passage of shock waves through the SN envelope. If the value of the 13-mixing angle is unknown, using additionally the information encoded in the prompt neutronization $\nu_e$ burst--a robust feature found in all modern SN simulations--can be sufficient to fix both the neutrino hierarchy and to decide whether $\theta_{13}$ is ``small'' or ``large.''

We calculate the yield of high energy neutrinos produced in astrophysical sources with negligible magnetic fields varying their interaction depth from nearly transparent to opaque. We take into account the scattering of secondaries on background photons as well as the direct production of neutrinos in decays of charm mesons. If multiple scattering of nucleons becomes important, the neutrino spectra from meson and muon decays are strongly modified with respect to transparent sources. Characteristic for neutrino sources containing photons as scattering targets is a strong energy-dependence of the ratio $R^0$ of $\nu_\mu$ and $\nu_e$ fluxes at the sources, ranging from $R^0=\phi_\mu/\phi_e\sim 0$ below threshold to $R^0\sim 4$ close to the energy where the decay length of charged pions and kaons equals their interaction length on target photons. Above this energy, the neutrino flux is strongly suppressed and depends mainly on charm production.

One of the robust features found in simulations of core-collapse supernovae (SNe) is the prompt neutronization burst, i.e. the first $\sim 25$ milliseconds after bounce when the SN emits with very high luminosity mainly $\nu_e$ neutrinos. We examine the dependence of this burst on variations in the input of current SN models and find that recent improvements of the electron capture rates as well as uncertainties in the nuclear equation of state or a variation of the progenitor mass have only little effect on the signature of the neutronization peak in a megaton water Cherenkov detector for different neutrino mixing schemes. We show that exploiting the time-structure of the neutronization peak allows one to identify the case of a normal mass hierarchy and large 13-mixing angle $\theta_{13}$, where the peak is absent. The robustness of the predicted total event number in the neutronization burst makes a measurement of the distance to the SN feasible with a precision of about 5%, even in the likely case that the SN is optically obscured.

A few seconds after bounce in a core-collapse supernova, the shock wave passes the density region corresponding to resonant neutrino oscillations with the ``atmospheric'' neutrino mass difference. The transient violation of the adiabaticity condition manifests itself in an observable modulation of the neutrino signal from a future galactic supernova. In addition to the shock wave propagation effects that were previously studied, a reverse shock forms when the supersonically expanding neutrino-driven wind collides with the slower earlier supernova ejecta. This implies that for some period the neutrinos pass two subsequent density discontinuities, giving rise to a ``double dip'' feature in the average neutrino energy as a function of time. We study this effect both analytically and numerically and find that it allows one to trace the positions of the forward and reverse shocks. We show that the energy dependent neutrino conversion probabilities allow one to detect oscillations even if the energy spectra of different neutrino flavors are the same as long as the fluxes differ. These features are observable in the \bar\nu_e signal for an inverted and in the \nu_e signal for a normal neutrino mass hierarchy, provided the 13-mixing angle is ``large'' (sin^2\theta_{13}\gg 10^{-5}).

The Earth matter effects on supernova (SN) neutrinos can be identified at a single detector through peaks in the Fourier transform of their ``inverse energy'' spectrum. The positions of these peaks are independent of the SN models and therefore the peaks can be used as a robust signature of the Earth matter effects, which in turn can distinguish between different neutrino mixing scenarios. Whereas only one genuine peak is observable when the neutrinos traverse only the Earth mantle, traversing also the core gives rise to multiple peaks. We calculate the strengths and positions of these peaks analytically and explore their features at a large scintillation detector as well as at a megaton water Cherenkov detector through Monte Carlo simulations. We propose a simple algorithm to identify the peaks in the actual data and quantify the chances of a peak identification as a function of the location of the SN in the sky.

A future galactic SN can be located several hours before the optical explosion through the MeV-neutrino burst, exploiting the directionality of $\nu$-$e$-scattering in a water Cherenkov detector such as Super-Kamiokande. We study the statistical efficiency of different methods for extracting the SN direction and identify a simple approach that is nearly optimal, yet independent of the exact SN neutrino spectra. We use this method to quantify the increase in the pointing accuracy by the addition of gadolinium to water, which tags neutrons from the inverse beta decay background. We also study the dependence of the pointing accuracy on neutrino mixing scenarios and initial spectra. We find that in the ``worst case'' scenario the pointing accuracy is $8^\circ$ at 95% C.L. in the absence of tagging, which improves to $3^\circ$ with a tagging efficiency of 95%. At a megaton detector, this accuracy can be as good as $0.6^\circ$. A TeV-neutrino burst is also expected to be emitted contemporaneously with the SN optical explosion, which may locate the SN to within a few tenths of a degree at a future km$^2$ high-energy neutrino telescope. If the SN is not seen in the electromagnetic spectrum, locating it in the sky through neutrinos is crucial for identifying the Earth matter effects on SN neutrino oscillations.

We have reconsidered the atmospheric neutrino anomaly in light of the laetst data from Super-Kamiokande contained events and from Super-Kamiokande and MACRO up-going muons. We have reanalysed the proposed solution to the atmospheric neutrino anomaly in terms of non-standard neutrino-matter interactions (NSI) as well as the standard nu_mu -> nu_tau oscillations (OSC). Our statistical analysis shows that a pure NSI mechanism is now ruled out at 99%, while the standard nu_mu -> nu_tau OSC mechanism provides a quite remarkably good description of the anomaly. We therefore study an extended mechanism of neutrino propagation which combines both oscillation and non-standard neutrino-matter interactions, in order to derive limits on flavour-changing (FC) and non-universal (NU) neutrino interactions. We obtain that the off-diagonal flavour-changing neutrino parameter epsilon and the diagonal non-universality neutrino parameter epsilon' are confined to -0.03 < epsilon < 0.02 and |epsilon'| < 0.05 at 99.73% CL. These limits are model independent and they are obtained from pure neutrino-physics processes. The stability of the neutrino oscillation solution to the atmospheric neutrino anomaly against the presence of non-standard neutrino interactions establishes the robustness of the near-maximal atmospheric mixing and massive-neutrino hypothesis. The best agreement with the data is obtained for Delta_m^2 = 2.3*10^{-3} eV^2, sin^2(2*theta) = 1, epsilon = 6.7*10^{-3} and epsilon' = 1.1*10^{-3}, although the chi^2 function is quite flat in the epsilon and epsilon' directions for epsilon, epsilon' -> 0.

We investigate the possibility to use the neutrinos coming from a future galactic supernova explosion to perform neutrino oscillation tomography of the Earth's core. We propose to use existing or planned detectors, resulting in an additional payoff. Provided that all of the discussed uncertainties can be reduced as expected, we find that the average matter densities of the Earth's inner and outer cores could be measured with a precision competitive with geophysics. However, since seismic waves are more sensitive to matter density jumps than average matter densities, neutrino physics would give partly complementary information.

We point out that solar neutrino oscillations with large mixing angle as evidenced in current solar neutrino data have a strong impact on strategies for diagnosing collapse-driven supernova (SN) through neutrino observations. Such oscillations induce a significant deformation of the energy spectra of neutrinos, thereby allowing us to obtain otherwise inaccessible features of SN neutrino spectra. We demonstrate that one can determine temperatures and luminosities of non-electron flavor neutrinos by observing bar{nu}_{e} from galactic SN in massive water Cherenkov detectors by the charged current reactions on protons.

We study neutrino oscillations and the level-crossing probability PLZ in power-law potential profiles A(r)\propto r^n. We give local and global adiabaticity conditions valid for all mixing angles theta and discuss different representations for PLZ. For the 1/r^3 profile typical of supernova envelopes we compare our analytical to numerical results and to earlier approximations used in the literature. We then perform a combined likelihood analysis of the observed SN1987A neutrino signal and of the latest solar neutrino data, including the recent SNO CC measurement. We find that, unless all relevant supernova parameters (released binding energy, \bar\nu_e and \bar\nu_{\mu,\tau} temperatures) are near their lowest values found in simulations, the status of large mixing type solutions deteriorates considerably compared to fits using only solar data. This is sufficient to rule out the vacuum-type solutions for most reasonable choices of astrophysics parameters. The LOW solution may still be acceptable, but becomes worse than the SMA-MSW solution which may, in some cases, be the best combined solution. On the other hand the LMA-MSW solution can easily survive as the best overall solution, although its size is generally reduced when compared to fits to the solar data only.

We study the possible impact of the neutrino oscillation which could be induced by a tiny violation of equivalence principle (VEP) for neutrinos emitted from supernova driven by gravitational collapse. Due to the absence of any significant indication of neutrino oscillation in the SN1987A data, we obtain sever bounds on relevant VEP parameters \delta \gamma\lsim O(10^{-31}) for massless or degenerated neutrinos and \delta \gamma \lsim O(10^{-16})\times [\Delta m^2/10^{-5} eV^2] for massive neutrinos.

The present solar and atmospheric neutrino data together with the LSND results and the presence of hot dark matter (HDM) suggest the existence of a sterile neutrino at the eV scale. We have reanalysed the effect of resonant {\sl sterile} neutrino conversions induced by neutrino magnetic moments in a type-II supernova. We analyse the implications of $\nu_e-\nu_s$ and $\bar{\nu}_e-\bar{\nu}_s$ ($\nu_s$ denotes sterile neutrino) conversions for the supernova shock re-heating, the detected $\bar\nu_e$ signal from SN1987A and the $r$-process nucleosynthesis hypothesis. Using reasonable magnetic field profiles we determine the sensitivity of these three arguments to the relevant neutrino parameters, i.e. the value of the transition magnetic moment and the $\nu_e-\nu_s$ mass difference $\Delta m^2 \equiv \Delta m^2_{LSND/HDM}$