We calculate the spectral energy distribution (SED) of electromagnetic radiation and the spectrum of high energy neutrinos from BL Lac objects in the context of the Synchrotron Proton Blazar Model. In this model, the high energy hump of the SED is due to accelerated protons, while most of the low energy hump is due to synchrotron radiation by co-accelerated electrons. To accelerate protons to sufficiently high energies to produce the high energy hump, rather high magnetic fields are required. Assuming reasonable emission region volumes and Doppler factors, we then find that in low-frequency peaked BL Lacs (LBLs), which have higher luminosities than high-frequency peaked BL Lacs (HBLs), there is a significant contribution to the high frequency hump of the SED from pion photoproduction and subsequent cascading, including synchrotron radiation by muons. In contrast, in HBLs we find that the high frequency hump of the SED is dominated by proton synchrotron radiation. We are able to model the SED of typical LBLs and HBLs, and to model the famous 1997 flare of Markarian 501. We also calculate the expected neutrino output of typical BL Lac objects, and estimate the diffuse neutrino intensity due to all BL Lacs. Because pion photoproduction is inefficient in HBLs, as protons lose energy predominantly by synchrotron radiation, the contribution of LBLs dominates the diffuse neutrino intensity. We suggest that nearby LBLs may well be observable with future high-sensitivity TeV gamma-ray telescopes.

Three cometary-shaped objects in the giant HII region NGC 3603, originally found and identified as proto-planetary disks (ProPlyDs) by Brandner et al. 2000 using HST+VLT in the optical and near-infrared, have been detected with the Australia Telescope Compact Array (ATCA) in the radio continuum at 3 and 6cm. All three ProPlyD-like objects are clearly resolved with an extent of a few arcseconds. The integrated 6cm fluxes are up to 1.3 times higher than the 3cm fluxes with spectral indices averaged over the whole clump between $\alpha=-0.1$ and -0.5 ($S_\nu \propto \nu^{\alpha}$), indicating the likely presence of non-thermal emission in at least some of the sources. All measured fluxes are at least one order of magnitude higher than those predicted by Brandner et al. 2000. Due to the unexpectedly large radio luminosities of the ProPlyD-like features and because the radio emission is extended a (proto-)stellar origin of the non-thermal emission from a dust enshrouded star appears unlikely. Instead we propose that magnetized regions within the envelope of the ProPlyD-like nebulae exist. (abridged)

The Synchrotron Proton Blazar model is a promising model to explain high energy emission from gamma-ray loud BL Lac objects like Mkn 421. In contrast to leptonic models, the hadronic explanation of gamma-ray emission predicts ultrahigh energy neutrinos. The predicted neutrino spectra from a typical High-energy cutoff BL Lac Object (HBL) and a Low-energy cutoff BL Lac Object (LBL) are presented. We find that cooling due to muon synchrotron radiation causes a cutoff of the neutrino spectrum at $\sim 10^{18}$ eV, with the exception of $\nu_\mu$ from kaon decay which may extend to higher energies if meson production takes place in the secondary resonance region of the cross section. The impact of the neutrino output from both source populations to the diffuse neutrino background is discussed.

We apply the synchrotron proton blazar (SPB) model to the April 1997 flare of Markarian 501 and find we are able to fit the observed spectral energy distribution. We explore the effect of target photon density on the high energy part of the spectral energy distribution (SED) for fixed assumed magnetic field, emission region size and Doppler factor and find that the luminosity and peak frequency of the high energy part of the SED may depend on the luminosity of the low energy part of the SED in high-frequency peaked BL Lac objects (HBL).

We present a new calculation of the propagation of protons with energies above $10^{19}$ eV over distances of up to several hundred Mpc. The calculation is based on a Monte Carlo approach using the event generator SOPHIA for the simulation of hadronic nucleon-photon interactions and a realistic integration of the particle trajectories in a random extragalactic magnetic field. Accounting for the proton scattering in the magnetic field affects noticeably the nucleon energy as a function of the distance to their source and allows us to give realistic predictions on arrival energy, time delay, and arrival angle distributions and correlations as well as secondary particle production spectra.

(abr.) The spectral energy distribution (SED) of blazars typically has a double-humped appearance usually interpreted in terms of synchrotron self-Compton models. In proton blazar models, the SED is instead explained in terms of acceleration of protons and subsequent cascading. We discuss a variation of the Synchrotron Proton Blazar model, first proposed by M\"ucke & Protheroe (1999), in which the low energy part of the SED is mainly synchrotron radiation by electrons co-accelerated with protons which produce the high energy part of the SED mainly asproton synchrotron radiation. Using a Monte Carlo/numerical technique to simulate the interactions and subsequent cascading of the accelerated protons, we are able to fit the observed SED of Markarian 501 during the April 1997 flare. We find that the emerging cascade spectra initiated by gamma-rays from $\pi^0$ decay and by $e^\pm$ from $\mu^\pm$ decay turn out to be relatively featureless. Synchrotron radiation produced by $\mu^\pm$ from $\pi^\pm$ decay, and even more importantly by protons, and subsequent synchrotron-pair cascading, is able to reproduce well the high energy part of the SED. For this fit we find that synchrotron radiation by protons dominates the TeV emission, pion photoproduction being less important with the consequence that we predict a lower neutrino flux than in other proton blazar models.

The April 1997 giant flare of Mkn 501 is modelled by a stationary Synchrotron-Proton-Blazar model. Our derived model parameters are consistent with X-ray-to-TeV-data in the flare state and diffusive shock acceleration of electron and proton in a Kolmogorov/Kraichnan turbulence sprectrum. While the emerging pair-synchrotron cascade spectra initiated by photons from neutral pi-decay and electrons from pi+- --> mu+- --> e+- decay turn out to be relatively featureless, muon and proton synchrotron radiation and their cascade radiation produce a double-humped spectral energy distribution. For the present model we find proton synchrotron radiation to dominate the TeV emission, while the contribution from the synchrotron radiation of the pairs, produced by the high energy hump, is only minor.

Photomeson production is the main energy loss for relativistic nucleons in dense radiation fields like the cosmic microwave background and the radiation fields in Gamma Ray Bursts (GRB) and jets of Active Galactic Nuclei (AGN). In this paper we study photomeson production in typical GRB and AGN jet radiation fields by using the recently developed Monte Carlo event generator SOPHIA (see these proceedings). We discuss processes that are relevant for the physics of cosmic ray acceleration and the production of neutrinos and gamma rays. We compare our results with widely used approximations, and find significant deviations, particularly for GRBs. The photoproduction of antibaryons as a so far not considered effect in astrophysics is briefly discussed.

A new Monte Carlo program for photohadronic interactions of relativistic nucleons with an ambient photon radiation field is presented. The event generator is designed to fulfil typical astrophysical requirements, but can also be used for radiation and background studies at high energy colliders such as LEP2 and HERA, as well as for simulations of photon induced air showers. We consider the full photopion production cross section from the pion production threshold up to high energies. It includes resonance excitation and decay, direct single pion production and diffractive and non-diffractive multiparticle production. The cross section of each individual process is calculated by fitting experimental data, while the kinematics is determined by the underlying particle production process. We demonstrate that our model is capable of reproducing known accelerator data over a wide energy range.

We discuss the first applications of our newly developed Monte Carlo event generator SOPHIA to multiparticle photoproduction of relativistic protons with thermal and power law radiation fields. The measured total cross section is reproduced in terms of excitation and decay of baryon resonances, direct pion production, diffractive scattering, and non-diffractive multiparticle production. Non--diffractive multiparticle production is described using a string fragmentation model. We demonstrate that the widely used `$\Delta$--approximation' for the photoproduction cross section is reasonable only for a restricted set of astrophysical applications. The relevance of this result for cosmic ray propagation through the microwave background and hadronic models of active galactic nuclei and gamma-ray bursts is briefly discussed.

We present a model for estimating the blazar contribution to the extragalactic gamma-ray background (EGRB) in the EGRET energy range. It is based on inverse-Compton scattering as the dominant gamma-ray production process in the jets of flat spectrum radio quasars (FSRQs) and BL Lac objects, and on the unification scheme of radio-loud AGN. Depending on the evolutionary behaviour of radio-loud AGN, our estimate for the point source contribution to the extragalactic background emission ranges from about ~20%-40% for a redshift drop off at Z_{max}=3 to around ~35%-85% for a cutoff at Z_{max}=5. Roughly 70%-90% of the AGN contribution to the EGRB would be due to BL~Lacs/FR I objects. However, in case of BL Lacs our model parameters are not entirely in agreement with the results from studies at radio wavelengths indicating systematical uncertainties whereas our estimate for the FSRQ contribution seems to be consistent. Thus we end up with two possibilities depending on whether we underpredict or overpredict the BL Lac contribution: either unresolved AGN can not account for the entire EGRB or unresolved BL Lacs produce the observed background.

Starting in 1991 we have performed a regular monitoring of EGRET sources with the Effelsberg 100-m telescope and the IRAM 30-m at Pico Veleta. In comparison with data on a sample of flat-spectrum quasars which have not been seen by EGRET we search for correlations of any kind in the behaviour of FSRQ's in the radio and gamma-ray regime. While there is no radio-to-gamma-ray luminosity correlation for these sources, it seems that gamma-ray high states coincide with increased activity in the radio regime with a strong tendency that gamma-ray outbursts precede radio outbursts. The gamma-ray spectra seem to harden with increasing flux level.