New experiments, exploring the ultra-high energy tail of the cosmic ray spectrum with unprecedented detail, are exerting a severe pressure on extensive air hower modeling. Detailed fast codes are in need in order to extract and understand the richness of information now available. Some hybrid simulation codes have been proposed recently to this effect (e.g., the combination of the traditional Monte Carlo scheme and system of cascade equations or pre-simulated air showers). In this context, we explore the potential of SENECA, an efficient hybrid tridimensional simulation code, as a valid practical alternative to full Monte Carlo simulations of extensive air showers generated by ultra-high energy cosmic rays. We extensively compare hybrid method with the traditional, but time consuming, full Monte Carlo code CORSIKA which is the de facto standard in the field. The hybrid scheme of the SENECA code is based on the simulation of each particle with the traditional Monte Carlo method at two steps of the shower development: the first step predicts the large fluctuations in the very first particle interactions at high energies while the second step provides a well detailed lateral distribution simulation of the final stages of the air shower. Both Monte Carlo simulation steps are connected by a cascade equation system which reproduces correctly the hadronic and electromagnetic longitudinal profile. We study the influence of this approach on the main longitudinal characteristics of proton-induced air showers and compare the predictions of the well known CORSIKA code using the QGSJET hadronic interaction model.

The output parameters from the ground array of the Auger South observatory, were simulated for the typical instrumental and environmental conditions at its Malarg\"ue site using the code sample-sim. Extensive air showers started by photons, protons and iron nuclei at the top of the atmosphere were used as triggers. The study utilized the air shower simulation code Aires with both QGSJet and Sibyll hadronic interaction models. A total of 1850 showers were used to produce more than 35,000 different ground events. We report here on the results of a multivariate analysis approach, including principal component analysis and neural networks, to the development of new primary composition diagnostics.

Introducing a simple Galactic wind model patterned after the solar wind we show that back-tracing the orbits of the highest energy cosmic events suggests that they may all come from the Virgo cluster, and so probably from the active radio galaxy M87. This confirms a long standing expectation. Those powerful radio galaxies that have their relativistic jets stuck in the interstellar medium of the host galaxy, such as 3C147, will then enable us to derive limits on the production of any new kind of particle, expected in some extensions of the standard model in particle physics. New data from HIRES will be crucial in testing the model proposed here.

Intensive monitoring of low-extinction windows towards the galactic bulge has provided in the last years valuable information for studies about the dynamics, kinematics and formation history of this part of the galaxy, mainly by characterizing the bulge stellar populations (Paczy\'nski, 1996). Since 1997, we have been conducting an intensive photometric-astrometric survey of the galactic bulge, with the monitoring of about 120000 stars in 12 windows uniformly distributed in galactic latitude and longitude (Blanco & Terndrup, 1989 e Blanco, 1988) never before submitted to this kind of survey. For this purpose, we have used the IAG/USP CCD Meridian Circle of the Abrah\~ao de Moraes Observatory. The main objective of this work is the identification and classification of variable objects. In this work we present the set up and development of the necessary tools for a project like this and the posterior analysis of our data. We briefly describe the construction of a program to organize and detect variables among the observed stars, including real time alerts (for variations greater than 0.3 magnitudes). The preliminary analysis after the processing of 76 nights of observation yielded 479 variable stars, from which 96.7 % of them are new. We discuss the preliminary classification of this variables, based on: a) the observed amplitude of variation; b) the shape of light curve; c) the expected variable classes among our data and d) the calculated periods, whenever possible. Finally, we discuss the future perspectives for the project and for the applications and analysis of the discovered variable stars.

The Ultra-high energy cosmic ray energy spectrum summarized by the AGASA collaboration indicates clearly that the cosmic ray spectrum extends well beyond the Greisen-Zatsepin-Kuzmin (GZK) cut-off at \sim 5 x 10^19 eV. Furthermore, despite the small number statistics involved, some structure in the spectrum may be emerging. Using numerical simulations, it is demonstrated in the present work that these features are consistent with a spatial distribution of sources that follows the distribution of luminous matter in the local Universe. Therefore, from this point of view, there is no need for a second high-energy component of cosmic rays dominating the spectrum beyond the GZK cut-off.

Several explanations for the existence of Ultra High Energy Cosmic Rays (UHECR) invoke the idea that they originate from the decay of massive particles created in the reheating following inflation. It has been suggested that the decay products can explain the observed isotropic flux of UHECR. We have calculated the anisotropy expected for various models of the dark matter distribution and find that present data are too sparse above 4 x 10^19 eV to discriminate between different models. However, after three years of operation of the southern part of the Pierre Auger Observatory great progress in testing the proposals is expected.

We present in this work secondary catalogs up to $m_{Val} \sim 13$ based on the Tycho reference frame (ESA, 1997) for 12 selected low-extinction fields towards the galactic bulge. The observations have been performed with the Askania-Zeiss Meridian Circle equiped with a CCD camera, located at the Abrah\~ao de Moraes Observatory (Valinhos, Brazil) and operated by the Institute of Astronomy and Geophysics, S\~ao Paulo University. The presented catalog, though not complete, has been designed to help in intensive search programmes (e.g. microlensing and variable searches) and therefore the selected standards have a high astrometric and photometric ($V$ band, approximately) quality. The mean precisions obtained were $0.0018^{s}$ in $\alpha$, 0.013'' in $\delta$, 0.030 for the standard deviation in magnitude and 0.0042 for the magnitude when weighted with the error bars in each night (in the mean, 42 stars for the catalog of each window). Tables B.1 to B.12 are also available in eletronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/Abstract.html.

Several explanations for the existence of Ultra High Energy Cosmic Rays invoke the idea that they originate from the decay of massive particles created in the reheating following inflation. It has been suggested that the decay products can explain the observed isotropic flux of cosmic rays. We have calculated the anisotropy expected for various models of the dark matter distribution and find that at present data are too sparse above $4 \times 10^{19}$ eV to discriminate between different models. However we show that with data from three years of operation of the southern section of the Pierre Auger Observatory significant progress in testing the proposals will be made.

The existence of cosmic ray particles up to the ultra-high energy limit (> 10^20 eV) is now beyond any doubt. The detection of cosmic particles with such energies imposes a challenge for the comprehension of their sources and nature. On one side, particles with such high energies are difficult to be produced by any astrophysical source. On the other side, the interactions of these particles with photons of the cosmic microwave background cause substantial losses of energy which constraint the maximum distances that the particles are able to travel from the sources to the detectors. Aiming to help to elucidate the problem of UHECR source identification, we have performed 3-D simulations of particle trajectories propagated through the stochastic intergalactic and an extended Galactic halo magnetic fields. Going further, we have also performed simulations of proton and Fe nuclei through the spiral Galactic magnetic field (GMF) and built full-sky maps of their arrival direction distribution in both the detector (after deflection in the GMF) and just outside the Galaxy. In this work we summarize the main results of these investigations.

We here investigate the possibility that the observed collimated outflows in luminous infrared galaxies (LIGs) and some Seyfert galaxies can be produced in a starburst (SB) environment. A nuclear disk can be quickly produced by gas infall during star formation in a rotating, stellar cluster. We find that massive nuclear SBs with core disk masses M_d \sim 10^8 - 10^9 M_{\odot}, and supernova rates \nu_{SN} \simeq 5 \times 10^{-3} - 2 yr^{-1} (which are consistent with the \nu_{SN} values inferred from the observed non-thermal radio power in source candidates) may inject kinetic energies which are high enough to blow out directed flows from the accreting disk surface, within the SB lifetimes. In our models, the acceleration and collimation of the nuclear outflow are provided by magnetic fields anchored into the rotating SB-disk. The emerging outflow carries a kinetic power that is only a small fraction (a few percent) of the supernovae energy rate produced in the SB. Based on conditions determined from observed outflows and disks, we find that moderate disk magnetic fields (\gtrsim 8 \times 10^{-4} G) are able to accelerate the outflows up to the observed terminal velocities (\lesssim few 100 km s^{-1} in the case of the Seyfert galaxies, and \sim 400 - 950 km s^{-1} in the case of the LIGs). The outflow is produced within a wind zone in the disk of radius \lesssim 100 pc in the LIGs, and \lesssim 10 pc in the Seyferts, with wind mass loss to disk accretion rate ratios \dot M_w /\dot M_d \gtrsim 0.1 (where \dot M_d \sim 100 M_{\odot} yr^{-1}). The observation of rotating nuclear disks of gas within few 100 pc scales in source candidates like the LIG Arp 220, and magnetized outflows provide observational support for the picture drawn here.

Seven and a half years of continuous monitoring of giant air showers triggered by ultra high-energy cosmic rays have been recently summarized by the AGASA collaboration. The resulting energy spectrum indicates clearly that the cosmic ray spectrum extends well beyond the Greisen-Zatsepin-Kuzmin (GZK) cut-off at $\sim 5 \times 10^{19}$ eV. Furthermore, despite the small number statistics involved, some structure in the spectrum may be emerging. Using numerical simulations, it is demonstrated in the present work that these features are consistent with a spatial distribution of sources that follows the distribution of luminous matter in the local Universe. Therefore, from this point of view, there is no need for a second high-energy component of cosmic rays dominating the spectrum beyond the GZK cut-off.

The interesting possibility of measuring the masses of high energy cosmic ray particles by observing pairs of extensive air showers arriving at the earth nearly simultaneously was proposed some years ago by Gerasimova and Zatsepin (1960). Such showers would be created by the nuclear fragments originating as a result of the photodisintegration of massive nuclei interacting with the solar radiation field. In this paper we re-visit this possibility in the context of existing and proposed detectors of high and ultra-high energy cosmic rays considering a simple, yet realistic, model of the interplanetary magnetic field. The possibilty of observing the mass fragmentation of cosmic rays directly, however, remains challenging.

A general view is presented on the problem of propagation of ultra high-energy cosmic rays through the intergalactic and galactic magnetic fields. Especial emphasis is given to the possibility of correlating the present events with potential sources face to the uncertainties in the intervening magnetic fields. Given high enough statistics, the latter problem can be transformed into a powerful tool for the study of cosmic magnetic fields. Finally, the three pairs of events pointed out by the AGASA collaboration as a possible evidence for clustering are analyzed.

The possibility that the magnetic field is strongly correlated with the large-scale structure of the universe has been recently considered in the literature. In this scenario the intergalactic magnetic field has a strong ($\mu$G) regular component spanning tens of Mpc but localized in sheets and filaments, while the vast voids in between are almost free of magnetic field. If true, this could have important consequences on the propagation of ultra-high energy cosmic rays, and severely affect our capacity of doing astronomy with charged particles. A quantitative discussion of these effects is given in the present work.

The possibility that the magnetic field is strongly correlated with the large-scale structure of the universe has been recently proposed in the literature. In this scenario the intergalactic magnetic field has a strong (m Gauss) regular component spanning tens of Mpc. This could have severe consequences on the propagation of ultra-high energy cosmic rays, and the observed spectra, isotropy and composition. A quantitative discussion of these effects is given in the present work.

Several structures have been discovered embedded in regions of recent or ongoing star formation, which point to the importance of the interaction between fast moving wind-blowing stars and their environment. Using hydrodynamic simulations, we investigate the passage through the interstellar medium of a supersonic stellar wind source, and show how it can naturally lead to the formation of tubes, channels and swamps of globules as interfaces are crossed. The results are in excellent agreement with observation of 30 Doradus.

Three pairs of possibly correlated ultra-high energy cosmic ray events were reported by Hayashida et al (1996). Here we calculate the propagation of the corresponding particles through both the intergalactic and galactic magnetic fields. The large scale disc and halo magnetic components are approximated by the models of Stanev (1997). The intergalactic magnetic field intensity is modulated by the actual density of luminous matter along the corresponding lines of sight, calculated from the CfA redshift catalogue (Huchra et al, 1995). The results indicate that, if the events of each pair had a common source and were simultaneously produced, they either originated inside the galactic halo or otherwise very unlikely events were observed. On the other hand, an estimate of the arrival probability of ultra-high energy cosmic rays, under the assumption that the distribution of luminous matter in the nearby universe traces the distribution of the sources of the particles and intensity of the intergalactic magnetic field, suggests that the pairs are chance clusterings.

The deflection of extragalactic ultra-high energy protons (E > 4E19 eV) travelling to our galaxy is studied assuming that visible matter traces both, the sources of the particles and the intergalactic magnetic field. The reversal scale and the magnitude of the magnetic field are determined by the local density of matter. The CfA Redshift Catalog is used to determine the distribution of galaxies inside 50 Mpc, where the sources are believed to reside. It is demonstrated that the arrival directions of ultra high energy cosmic rays are consistent with the distribution of galaxies inside 50 Mpc and that the proposed clusters of events point to regions of high density of sources instead of individual ones.

We report the results of 3D simulations of the trajectories of ultra-high energy protons and Fe nuclei (with energies $E = 4 \times 10^{19}$ and $2.5 \times 10^{20} eV$) propagating through the galactic magnetic field from the sources to the detector. A uniform distribution of anti-particles is backtracked from the detector, at the Earth, to the halo of the Galaxy. We assume an axisymmetric, large scale spiral magnetic field permeating both the disc and the halo. A normal field component to the galactic plane ($B_z$) is also included in part of the simulations. We find that the presence of a large scale galactic magnetic field does not generally affect the arrival directions of the protons, although the inclusion of a $B_z$ component may cause significant deflection of the lower energy protons ($E = 4 \times 10^{19} $ eV). Error boxes larger than or equal to $\sim 5^{\circ}$ are most expected in this case. On the other hand, in the case of heavy nuclei, the arrival direction of the particles is strongly dependent on the coordinates of the particle source. The deflection may be high enough ($> 20^{\circ}$) as to make extremely difficult any identification of the sources unless the real magnetic field configuration is accurately determined. Moreover, not every incoming particle direction is allowed between a given source and the detector. This generates sky patches which are virtually unobservable from the Earth. In the particular case of the UHE events of Yakutsk, Fly's Eye, and Akeno, they come from locations for which the deflection caused by the assumed magnetic field is not significant.

Using the cosmic ray (CR) data available in the energy interval $(10 - 2 \times 10^{7})$ GeV/particle, we have calculated the profile of the primary $\gamma$-ray spectrum produced by the interaction of these CR with thermal nuclei of the ISM. Normalized to the EGRET measurements, this allows an estimate of the galactic diffuse $\gamma$-ray background due to intermediate and high energy CR at TeV energies. On the other hand, over the last few years, several particles with energies above $10^{20}$ eV (beyond the Greisen-Zatsepin-Kuzmin cut-off) have been detected. These particles are very likely extragalactic protons originated at distances not greater than $30 - 50$ Mpc [e.g., 1]. The propagation of these ultra-high energy protons (UHEP) through the intergalactic medium leads to the development of $\gamma$-ray cascades and an ultimate signature at TeV energies. To assess the statistical significance of this $\gamma$-ray signature by the UHEP, we have also simulated the development of electromagnetic cascades triggered by the decay of a $10^{19}$ eV $\pi^{o}$ in the intergalactic medium after an UHEP collision with a cosmic microwave background photon.

The Ha peak intensity, velocity shift and velocity dispersion maps of the giant HII region NGC 604 in M 33, obtained by two dimensional high spatial resolution Fabry-Perot observations with TAURUS II at the 4.2 m William Herschel Telescope in Spain (Sabalisck, 1995), are analyzed via two point correlation functions. The whole system seems to rotate as a rigid body on scales from 50 to 80 pc (the largest studied scale), with a period of $\sim$ 85 Myr. We demonstrate that the cloud seems to be comprised of eddies with varying characteristic scale lengths which range from 10 pc to the largest observed scales. The calculated kinetic energy spectrum may be interpreted as either a manifestation of a double cascading spectrum of forced two-dimensional turbulence, or as a Kolmogorov three-dimensional turbulence (although this last possibility seems unlikely). According to the first interpretation, turbulence is being forced at scales of $\sim$ 10 pc, while an enstrophy (mean-square vorticity) cascade has developed down to the smallest scales resolved and an inverse kinetic energy cascade extends up to scales of $\sim$ 70 pc where a low wave number turn over is observed; if true, this would be the first time that such a phenomenon has been observed outside the Solar System. As for the second interpretation, energy should be injected at the largest scales, $\sim$ 70 pc. In both cases the average intrinsic optical depth consistent with the results is $\sim$ 20 pc.

We report the results of 3-D simulations of non-diffusive propagation of Ultra-High Energy Cosmic Rays (UHECR) (E > 10^{20} eV) through the intergalactic and extended halo media. We quantify the expected angular and temporal correlations between the events and the sources, and the temporal delay between protons and gamma-ray counterparts with a common origin for both halo and extragalactic origins. It is shown that the proposed UHECR-supergalactic plane source associations require either extremely high values of the halo magnetic field over as much as 100 kpc length scale or a very large correlation length for the IGM, even for the largest possible values of the intergalactic magnetic field. It can be stated that the UHECR seem to point to the sources even more strongly than previously believed. The simulations also show that the calculated time delays between UHE protons and gamma-ray counterparts do not match the claimed GRB-UHECR associations for either cosmological or extended halo distance scales.

A sample of 45 symbiotic stars observed in the direction of the galactic bulge are classified in spectral type using TiO $\lambda \lambda 6180, 7100 \AA$ and VO $\lambda \lambda 7865 \AA$ absorption feature indices and visual comparison. The effective temperatures obtained are used together with K data from the literature to estimate the bolometric luminosity of the giant components if located at d=8.5 kpc. The ($T_{eff}$, $L/L_{\odot}$) pairs are plotted on the HR diagram together with evolutionary tracks for metal rich intermediate mass stars. The mass of the giant components seem to be scattered in the interval ($\sim 1.5 - \sim 4 M_{\odot}$). The giant components of D-type symbiotics appear as if they were more evolved objects than the corresponding giants in S-type systems. Nevertheless, it is unlikely that they form an evolutionary sequence. The results point also to the restriction of the symbiotic phenomenon to the thermally pulsating phase of the AGB, strengthening an accelerated mass loss scenario for the latest evolutionary stages of intermediate mass stars.