We study coupled axial and polar axisymmetric oscillations of a neutron star endowed with a strong magnetic field, having both poloidal and toroidal components. The toroidal component of the magnetic field is driving the coupling between the polar and axial oscillations. The star is composed of a fluid core as well as a solid crust. Using a two dimensional general relativistic simulation and a magnetic field B = 10^16 G, we study the change in the polar and axial spectrum caused by the coupling. We find that the axial spectrum suffers a dramatic change in its nature, losing its continuum character. In fact, we find that only the 'edges' of the continua survive, generating a discrete spectrum. As a consequence the crustal frequencies, that in our previous simulation could be absorbed by the continua, if they were embedded inside it, are now long living oscillations. They may lose their energy only in the very special case that they are in resonance with the 'edges' of the continua.

We study axisymmetric perturbations of neutron star endowed with a strong magnetic field (magnetars), considering the coupled oscillations of the fluid core with the solid crust. We recover discrete oscillations based mainly in the crust and a continuum in the core. We also confirm the presence of "discrete Alfv\'en modes" in the gap between two contiguous continua (see Van Hooven& Levin (2011)) and, in addition, we can resolve some of them also inside the continua. Our results can explain both the lower and the higher observed quasi periodical oscillations (QPOs) in SGR 1806-20 and SGR 1900+14 and put constrains on the mass, radius and crust thickness of the two magnetars.

We study torsional Alfv\'en oscillations of magnetars, i.e., neutron stars with a strong magnetic field. We consider the poloidal and toroidal components of the magnetic field and a wide range of equilibrium stellar models. We use a new coordinate system (X,Y), where $X=\sqrt{a_1} \sin \theta$, $Y=\sqrt{a_1}\cos \theta$ and $a_1$ is the radial component of the magnetic field. In this coordinate system, the 1+2-dimensional evolution equation describing the quasi-periodic oscillations, QPOs, see Sotani et al. (2007), is reduced to a 1+1-dimensional equation, where the perturbations propagate only along the Y-axis. We solve the 1+1-dimensional equation for different boundary conditions and open magnetic field lines, i.e., magnetic field lines that reach the surface and there match up with the exterior dipole magnetic field, as well as closed magnetic lines, i.e., magnetic lines that never reach the stellar surface. For the open field lines, we find two families of QPOs frequencies; a family of "lower" QPOs frequencies which is located near the X-axis and a family of "upper" frequencies located near the Y-axis. According to Levin (2007), the fundamental frequencies of these two families can be interpreted as the turning points of a continuous spectrum. We find that the upper frequencies are constant multiples of the lower frequencies with a constant equaling 2n+1. For the closed lines, the corresponding factor is n+1 . By these relations, we can explain both the lower and the higher observed frequencies in SGR 1806-20 and SGR 1900+14.

We study the effect of the magnetic field geometry on the oscillation spectra of strongly magnetized stars. The magnetic field distributions include both toroidal and poloidal contributions. We observe that the toroidal contribution does not influence significantly the torsional oscillations of the crust. Moreover, in the case that the core is a type I supercontactor and the magnetic fields are confined in the crust, the torsional oscillation spectrum is drastically affected by the presence of the strong magnetic fields. Comparison with results and estimations for the magnetic field strength, from observations, exclude the possibility that magnetars will have a magnetic field solely confined in the crust i.e. our results suggest that the magnetic field in whatever geometry has to permeate the whole star.

We find numerical solutions of the coupled system of Einstein-Maxwell's equations with a linear approach, in which the magnetic field acts as a perturbation of a spherical neutron star. In our study, magnetic fields having both poloidal and toroidal components are considered, and higher order multipoles are also included. We evaluate the deformations induced by different field configurations, paying special attention to those for which the star has a prolate shape. We also explore the dependence of the stellar deformation on the particular choice of the equation of state and on the mass of the star. Our results show that, for neutron stars with mass M = 1.4 Msun and surface magnetic fields of the order of 10^15 G, a quadrupole ellipticity of the order of 10^(-6) - 10^(-5) should be expected. Low mass neutron stars are in principle subject to larger deformations (quadrupole ellipticities up to 10^(-3) in the most extreme case). The effect of quadrupolar magnetic fields is comparable to that of dipolar components. A magnetic field permeating the whole star is normally needed to obtain negative quadrupole ellipticities, while fields confined to the crust typically produce positive quadrupole ellipticities.

The Z-burst model and the direct propagation of UHE proton in negligible extragalactic magnetic fields produce gamma-rays afterglows and precursors halos, respectively at GeVs and TeV energy band a few degree around the UHECR arrival direction. The possible correlation of UHECR clusters (doublet, triplet) with nearby BL Lac sources at E_p = 4*10^19 eV offer a test for this necessary Gamma-UHECR trace. We estimate the secondary gamma energy and spectra and we suggest how to disentangle between the different scenarios. We show why Z-Burst model is still the most realistic model to explain UHECR behaviour and their correlation to known BL Lac sources.