We present a detailed analysis of the properties of twisted, force-free magnetospheres of non-rotating neutron stars, which are of interest in the modelling of magnetar properties and evolution. In our models the magnetic field smoothly matches to a current-free (vacuum) solution at some large external radius, and they are specifically built to avoid pathological surface currents at any of the interfaces. By exploring a large range of parameters, we find a few remarkable general trends. We find that the total dipolar moment can be increased by up to $40\%$ with respect to a vacuum model with the same surface magnetic field, due to the contribution of magnetospheric currents to the global magnetic field. Thus, estimates of the surface magnetic field based on the large-scale dipolar braking torque are slightly overestimating the surface value by the same amount. Consistently, there is a moderate increase in the total energy of the model with respect to the vacuum solution of up to $25\%$, which would be the available energy budget in the event of a fast, global magnetospheric reorganization commonly associated with magnetar flares. We have also found the interesting result of the existence of a critical twist ($\varphi_{\rm max} \lesssim 1.5$ rad), beyond which we cannot find any more numerical solutions. Combining the models considered in this paper with the evolution of the interior of neutron stars will allow us to study the influence of the magnetosphere on the long-term magnetic, thermal, and rotational evolution.

We study the long-term quasi-steady evolution of the force-free magnetosphere of a magnetar coupled to its internal magnetic field. We find that magnetospheric currents can be maintained on long timescales of the order of thousands of years. Meanwhile, the energy, helicity and twist stored in the magnetosphere all gradually increase over the course of this evolution, until a critical point is reached, beyond which a force-free magnetosphere cannot be constructed. At this point, some large-scale magnetospheric rearrangement, possibly resulting in an outburst or a flare, must occur, releasing a large fraction of the stored energy, helicity and twist. After that, the quasi-steady evolution should continue in a similar manner from the new initial conditions. The timescale for reaching this critical point depends on the overall magnetic field strength and on the relative fraction of the toroidal field. The energy stored in the force-free magnetosphere is found to be up to $\sim 30\%$ larger than the corresponding vacuum energy. This implies that for a $10^{14}$ G field at the pole, the energy budget available for fast magnetospheric events is of the order of a few $10^{44}$ erg. The spindown rate is estimated to increase by up to $\sim 60\%$, since the dipole content in the magnetosphere is enhanced by the currents present there. A rough estimate of the braking index $n$ reveals that it is systematically $n < 3$ for the most part of the evolution, consistent with actual measurements for pulsars and early estimates for several magnetars.

We extend our previous study of equilibrium solutions of non-rotating force-free magnetospheres of neutron stars. We show that multiple solutions exist for the same sets of parameters, implying that the solutions are degenerate. We are able to obtain configurations with disconnected field lines, however, in nearly all cases these correspond to degenerate higher energy solutions. We carry out a wide parametric search in order to understand the properties of the solutions. We confirm our previous results that the lower energy solutions have up to $\sim 25\%$ more energy than the vacuum case, helicity of the order of $\sim 5$ (in some defined units), maximum twist of $\sim 1.5$ rad, and a dipole strength that is up to $\sim 40\%$ larger than the vacuum dipole. Including the degenerate higher energy solutions allows for larger theoretical limits of up to $\sim 80\%$ more energy with respect to the vacuum case, helicity of the order of $\sim 8$, and a dipole strength that can now be up to four times that of the vacuum dipole, while the twist can be significantly larger and even diverge for configurations with disconnected domains. The higher energy solutions are probably unstable, therefore, it is unlikely that such magnetospheres exist under normal conditions in magnetars and high magnetic field pulsars.

We study the coupling of the force-free magnetosphere to the long-term internal evolution of a magnetar. We allow the relation between the poloidal and toroidal stream functions - that characterizes the magnetosphere - to evolve freely without constraining its particular form. We find that, on time-scales of the order of kyr, the energy stored in the magnetosphere gradually increases, as the toroidal region grows and the field lines expand outwards. This continues until a critical point is reached beyond which force-free solutions for the magnetosphere can no longer be constructed, likely leading to some large-scale magnetospheric reorganization. The energy budget available for such events can be as high as several $10^{45}\,$erg for fields of $10^{14}\,$G. Subsequently, starting from the new initial conditions, the evolution proceeds in a similar manner. The time-scale to reach the critical point scales inversely with the magnetic field amplitude. Allowing currents to pass through the last few meters below the surface, where the magnetic diffusivity is orders of magnitude larger than in the crust, should give rise to a considerable amount of energy deposition through Joule heating. We estimate that the effective surface temperature could increase locally from $\sim 0.1\,$keV to $\sim 0.3 - 0.6\,$keV, in good agreement with observations. Similarly, the power input from the interior into the magnetosphere could be as high as $10^{35} - 10^{36}\,$erg/s, which is consistent with peak luminosities observed during magnetar outbursts. Therefore, a detailed treatment of currents flowing through the envelope may be needed to explain the thermal properties of magnetars.

We revisit the various approximations employed to study the long-term evolution of the magnetic field in neutron star cores and discuss their limitations and possible improvements. A recent controversy on the correct form of the induction equation and the relevant evolution timescale in superconducting neutron star cores is addressed and clarified. We show that this ambiguity in the estimation of timescales arises as a consequence of nominally large terms that appear in the induction equation, but which are, in fact, mostly irrotational. This subtlety leads to a discrepancy by many orders of magnitude when velocity fields are absent or ignored. Even when internal velocity fields are accounted for, only the solenoidal part of the electric field contributes to the induction equation, which can be substantially smaller than the irrotational part. We also argue that stationary velocity fields must be incorporated in the slow evolution of the magnetic field as the next level of approximation.

We study ambipolar diffusion in strongly magnetised neutron stars, with special focus on the effects of neutrino reaction rates and the impact of a superfluid/superconducting transition in the neutron star core. For axisymmetric magnetic field configurations, we determine the deviation from $\beta-$equilibrium induced by the magnetic force and calculate the velocity of the slow, quasi-stationary, ambipolar drift. We study the temperature dependence of the velocity pattern and clearly identify the transition to a predominantly solenoidal flow. For stars without superconducting/superfluid constituents and with a mixed poloidal-toroidal magnetic field of typical magnetar strength, we find that ambipolar diffusion proceeds fast enough to have a significant impact on the magnetic field evolution only at low core temperatures, $T \lesssim 1-2\times10^8$ K. The ambipolar diffusion timescale becomes appreciably shorter when fast neutrino reactions are present, because the possibility to balance part of the magnetic force with pressure gradients is reduced. We also find short ambipolar diffusion timescales in the case of superconducting cores for $T \lesssim 10^9$ K, due to the reduced interaction between protons and neutrons. In the most favourable scenario, with fast neutrino reactions and superconducting cores, ambipolar diffusion results in advection velocities of several km/kyr. This velocity can substantially reorganize magnetic fields in magnetar cores, in a way that can only be confirmed by dynamical simulations.

A large fraction of Gamma Ray Bursts (GRBs) displays an X-ray plateau phase within <10^{5} s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, between January 2005 and August 2014. From the derived initial magnetic field distribution for the possible magnetars left behind by the GRBs, we study the evolution and properties of a simulated GRB-magnetar population using numerical simulations of magnetic field evolution, coupled with Monte Carlo simulations of Pulsar Population Synthesis in our Galaxy. We find that if the GRB X-ray plateaus are powered by the rotational energy of a newly formed magnetar, the current observational properties of the Galactic magnetar population are not compatible with being formed within the GRB scenario (regardless of the GRB type or rate at z=0). Direct consequences would be that we should allow the existence of magnetars and "super-magnetars" having different progenitors, and that Type Ib/c SNe related to Long GRBs form systematically neutron stars with higher initial magnetic fields. We put an upper limit of <16 "super-magnetars" formed by a GRB in our Galaxy in the past Myr (at 99% c.l.). This limit is somewhat smaller than what roughly expected from Long GRB rates, although the very large uncertainties do not allow us to draw strong conclusion in this respect.

Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include $X$-ray thermal emitting pulsars within the same evolutionary model as radio-pulsars. We find that the cumulative distribution of the number of X-ray pulsars can be well reproduced by several models that, simultaneously, reproduce the characteristics of the radio-pulsar distribution. However, even considering the most favourable magneto-thermal evolution models with fast field decay, log-normal distributions of the initial magnetic field over-predict the number of visible sources with periods longer than 12 s. We then show that the problem can be solved with different distributions of magnetic field, such as a truncated log-normal distribution, or a binormal distribution with two distinct populations. We use the observational lack of isolated NSs with spin periods P>12 s to establish an upper limit to the fraction of magnetars born with B > 10^{15} G (less than 1\%). As future detections keep increasing the magnetar and high-B pulsar statistics, our approach can be used to establish a severe constraint on the maximum magnetic field at birth of NSs.

Isolated neutron stars show a diversity in timing and spectral properties, which has historically led to a classification in different sub-classes. The magnetic field plays a key role in many aspects of the neutron star phenomenology: it regulates the braking torque responsible for their timing properties and, for magnetars, it provides the energy budget for the outburst activity and high quiescent luminosities (usually well above the rotational energy budget). We aim at unifying this observational variety by linking the results of the state-of-the-art 2D magneto-thermal simulations with observational data. The comparison between theory and observations allows to place two strong constraints on the physical properties of the inner crust. First, strong electrical currents must circulate in the crust, rather than in the star core. Second, the innermost part of the crust must be highly resistive, which is in principle in agreement with the presence of a novel phase of matter so-called nuclear pasta phase.

We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in parameter space that are more consistent with the observational data. In particular, we find that any broad distribution of birth spin periods with $P_0\lesssim 0.5$ s can fit the data, and that if the alignment angle is allowed to vary consistently with the torque model, realistic magnetospheric models are favoured compared to models with classical magneto-dipolar radiation losses. Assuming that the initial magnetic field is given by a lognormal distribution, our optimal model has mean strength $\langle\log B_0{\rm [G]}\rangle \approx 13.0-13.2$ with width $\sigma (\log B_0) = 0.6-0.7$. However, there are strong correlations between parameters. This degeneracy in the parameter space can be broken by an independent estimate of the pulsar birth rate or by future studies correlating this information with the population in other observational bands (X-rays and $\gamma$-rays).

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radio-pulsars, thus confirming the idea that magnetic fields are involved in their X-ray emission. Here we present the results of 2D simulations of the fully-coupled evolution of temperature and magnetic field in neutron stars, including the state-of-the-art kinetic coefficients and, for the first time, the important effect of the Hall term. After gathering and thoroughly re-analysing in a consistent way all the best available data on isolated, thermally emitting neutron stars, we compare our theoretical models to a data sample of 40 sources. We find that our evolutionary models can explain the phenomenological diversity of magnetars, high-B radio-pulsars, and isolated nearby neutron stars by only varying their initial magnetic field, mass and envelope composition. Nearly all sources appear to follow the expectations of the standard theoretical models. Finally, we discuss the expected outburst rates and the evolutionary links between different classes. Our results constitute a major step towards the grand unification of the isolated neutron star zoo.

The characteristic formulation of the relativistic hydrodynamic equations (Donat et al 1998 J. Comput. Phys. 146 58), which has been implemented in many relativistic hydro-codes that make use of Godunov-type methods, has to be slightly modified in the case of evolving barotropic flows. For a barotropic equation of state, a removable singularity appears in one of the eigenvectors. The singularity can be avoided by means of a simple renormalization which makes the system of eigenvectors well defined and complete. An alternative strategy for the particular case of barotropic flows is discussed.

Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.

Nowadays, the analysis of the X-ray spectra of magnetically powered neutron stars or magnetars is one of the most valuable tools to gain insight into the physical processes occurring in their interiors and magnetospheres. In particular, the magnetospheric plasma leaves a strong imprint on the observed X-ray spectrum by means of Compton up-scattering of the thermal radiation coming from the star surface. Motivated by the increased quality of the observational data, much theoretical work has been devoted to develop Monte Carlo (MC) codes that incorporate the effects of resonant Compton scattering in the modeling of radiative transfer of photons through the magnetosphere. The two key ingredients in this simulations are the kinetic plasma properties and the magnetic field (MF) configuration. The MF geometry is expected to be complex, but up to now only mathematically simple solutions (self-similar solutions) have been employed. In this work, we discuss the effects of new, more realistic, MF geometries on synthetic spectra. We use new force-free solutions in a previously developed MC code to assess the influence of MF geometry on the emerging spectra. Our main result is that the shape of the final spectrum is mostly sensitive to uncertain parameters of the magnetospheric plasma, but the MF geometry plays an important role on the angle-dependence of the spectra.

The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.

In a recent paper (Ant\'on et al. 2010) we have derived sets of right and left eigenvectors of the Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. We present a summary of the main steps followed in the above derivation and the numerical experiments carried out with the linearized (Roe-type) Riemann solver we have developed, and some note on the (non-)convex character of the relativistic MHD equations.

We perform population synthesis studies of different types of neutron stars taking into account the magnetic field decay. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G]) > \sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible. We apply the obtained field distribution and the model of decay to study long-term evolution of neuton stars till the stage of accretion from the interstellar medium. It is shown that though the subsonic propeller stage can be relatively long, initially highly magnetized neutron stars ($B_0 > \sim 10^{13}$ G) reach the accretion regime within the Galactic lifetime if their kick velocities are not too large. The fact that in previous studies made $>$10 years ago, such objects were not considered results in a slight increase of the Accretor fraction in comparison with earlier conclusions. Most of the neutron stars similar to the Magnificent seven are expected to become accreting from the interstellar medium after few billion years of their evolution. They are the main predecestors of accreting isolated neutron stars.

Here we introduce PHAT, the PHoto-z Accuracy Testing programme, an international initiative to test and compare different methods of photo-z estimation. Two different test environments are set up, one (PHAT0) based on simulations to test the basic functionality of the different photo-z codes, and another one (PHAT1) based on data from the GOODS survey. The accuracy of the different methods is expressed and ranked by the global photo-z bias, scatter, and outlier rates. Most methods agree well on PHAT0 but produce photo-z scatters that can differ by up to a factor of two even in this idealised case. A larger spread in accuracy is found for PHAT1. Few methods benefit from the addition of mid-IR photometry. Remaining biases and systematic effects can be explained by shortcomings in the different template sets and the use of priors on the one hand and an insufficient training set on the other hand. Scatters of 4-8% in Delta_z/(1+z) were obtained, consistent with other studies. However, somewhat larger outlier rates (>7.5% with Delta_z/(1+z)>0.15; >4.5% after cleaning) are found for all codes. There is a general trend that empirical codes produce smaller biases than template-based codes. The systematic, quantitative comparison of different photo-z codes presented here is a snapshot of the current state-of-the-art of photo-z estimation and sets a standard for the assessment of photo-z accuracy in the future. The rather large outlier rates reported here for PHAT1 on real data should be investigated further since they are most probably also present (and possibly hidden) in many other studies. The test data sets are publicly available and can be used to compare new methods to established ones and help in guiding future photo-z method development. (abridged)

We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wavefront in the fluid rest frame. Proper expressions of the renormalized eigenvectors in conserved variables are obtained through the corresponding matrix transformations. Our work completes previous analysis that present different sets of right eigenvectors for non-degenerate and degenerate states, and can be seen as a relativistic generalization of earlier work performed in classical MHD. Based on the full wave decomposition (FWD) provided by the the renormalized set of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) Riemann solver. Extensive testing against one- and two-dimensional standard numerical problems allows us to conclude that our solver is very robust. When compared with a family of simpler solvers that avoid the knowledge of the full characteristic structure of the equations in the computation of the numerical fluxes, our solver turns out to be less diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The amount of operations needed by the FWD solver makes it less efficient computationally than those of the HLL family in one-dimensional problems. However its relative efficiency increases in multidimensional simulations.

We perform population synthesis studies of different types of neutron stars (thermally emitting isolated neutron stars, normal radio pulsars, magnetars) taking into account the magnetic field decay and using results from the most recent advances in neutron star cooling theory. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. For this purpose, we fix a baseline neutron star model (all microphysics input), and other relevant parameters to standard values (velocity distribution, mass spectrum, birth rates ...), allowing to vary the initial magnetic field strength. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G])>\sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life (only about a factor of 2 for low field neutron stars but more than an order of magnitude for magnetars), and a mass distribution function dominated by low mass objects. This model explains satisfactorily all known populations. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible.

(Abridged) We present numerical simulations of the spectral evolution and emission of radio components in relativistic jets. We have developed an algorithm (SPEV) for the transport of a population of non-thermal electrons including radiative losses. For large values of the ratio of gas pressure to magnetic field energy density, \ab \sim 6\times 10^4, quiescent jet models show substantial spectral evolution, with observational consequences only above radio frequencies. Larger values of the magnetic field (\ab \sim 6\times 10^2), such that synchrotron losses are moderately important at radio frequencies, present a larger ratio of shocked-to-unshocked regions brightness than the models without radiative losses, despite the fact that they correspond to the same underlying hydrodynamic structure. We also show that jets with a positive photon spectral index result if the lower limit \gamma_min of the non-thermal particle energy distribution is large enough. A temporary increase of the Lorentz factor at the jet inlet produces a traveling perturbation that appears in the synthetic maps as a superluminal component. We show that trailing components can be originated not only in pressure matched jets, but also in over-pressured ones, where the existence of recollimation shocks does not allow for a direct identification of such features as Kelvin-Helmholtz modes, and its observational imprint depends on the observing frequency. If the magnetic field is large (\ab \sim 6\times 10^2), the spectral index in the rarefaction trailing the traveling perturbation does not change much with respect to the same model without any hydrodynamic perturbation. If the synchrotron losses are considered the spectral index displays a smaller value than in the corresponding region of the quiescent jet model.

We study the mutual influence of thermal and magnetic evolution in a neutron star's crust in axial symmetry. Taking into account realistic microphysical inputs, we find the heat released by Joule effect consistent with the circulation of currents in the crust, and we incorporate its effects in 2D cooling calculations. We solve the induction equation numerically using a hybrid method (spectral in angles, but a finite--differences scheme in the radial direction), coupled to the thermal diffusion equation. We present the first long term 2D simulations of the coupled magneto-thermal evolution of neutron stars. This substantially improves previous works in which a very crude approximation in at least one of the parts (thermal or magnetic diffusion) has been adopted. Our results show that the feedback between Joule heating and magnetic diffusion is strong, resulting in a faster dissipation of the stronger fields during the first million years of a NS's life. As a consequence, all neutron stars born with fields larger than a critical value (about 5 10^13 G) reach similar field strengths (approximately 2-3 10^{13} G) at late times. Irrespectively of the initial magnetic field strength, after $10^6$ years the temperature becomes so low that the magnetic diffusion timescale becomes longer than the typical ages of radio--pulsars, thus resulting in apparently no dissipation of the field in old NS. We also confirm the strong correlation between the magnetic field and the surface temperature of relatively young NSs discussed in preliminary works. The effective temperature of models with strong internal toroidal components are systematically higher than those of models with purely poloidal fields, due to the additional energy reservoir stored in the toroidal field that is gradually released as the field dissipates.

Context: Many thermally emitting isolated neutron stars have magnetic fields larger than 10^13 G. A realistic cooling model that includes the presence of high magnetic fields should be reconsidered. Aims: We investigate the effects of anisotropic temperature distribution and Joule heating on the cooling of magnetized neutron stars. Methods: The 2D heat transfer equation with anisotropic thermal conductivity tensor and including all relevant neutrino emission processes is solved for realistic models of the neutron star interior and crust. Results: The presence of the magnetic field affects significantly the thermal surface distribution and the cooling history during both, the early neutrino cooling era and the late photon cooling era. Conclusions: There is a large effect of the Joule heating on the thermal evolution of strongly magnetized neutron stars. Both magnetic fields and Joule heating play a key role in keeping magnetars warm for a long time. Moreover, this effect is important for intermediate field neutron stars and should be considered in radio-quiet isolated neutron stars or high magnetic field radio-pulsars.

We present two-dimensional simulations for the cooling of neutron stars with strong magnetic fields (B > 1e13 Gauss). We study how the cooling curves are influenced by magnetic field decay. We show that the Joule heating effects are very large and in some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field and predict the surface temperature distribution for the early and late stages of the evolution of isolated neutron stars, comparing our results with available observational data of isolated neutron stars.

The impact of strong magnetic fields B>10e13 G on the thermal evolution of neutron stars is investigated, including crustal heating by magnetic field decay. For this purpose, we perform 2D cooling simulations with anisotropic thermal conductivity considering all relevant neutrino emission processes for realistic neutron stars. The standard cooling models of neutron stars are called into question by showing that the magnetic field has relevant (and in many cases dominant) effects on the thermal evolution. The presence of the magnetic field significantly affects the thermal surface distribution and the cooling history of these objects during both, the early neutrino cooling era and the late photon cooling era. The minimal cooling scenario is thus more complex than generally assumed. A consistent magneto-thermal evolution of magnetized neutron stars is needed to explain the observations.

We present 2D simulations of the cooling of neutron stars with strong magnetic fields (B \geq 10^{13} G). We solve the diffusion equation in axial symmetry including the state of the art microphysics that controls the cooling such as slow/fast neutrino processes, superfluidity, as well as possible heating mechanisms. We study how the cooling curves depend on the the magnetic field strength and geometry. Special attention is given to discuss the influence of magnetic field decay. We show that Joule heating effects are very large and in some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field for the early and late stages of the evolution of isolated neutron stars.

This is the first paper of a series describing our measurement of weak lensing by large-scale structure using archival observations from the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST). In this work we present results from a pilot study testing the capabilities of the ACS for cosmic shear measurements with early parallel observations and presenting a re-analysis of HST/ACS data from the GEMS survey and the GOODS observations of the Chandra Deep Field South (CDFS). We describe our new correction scheme for the time-dependent ACS PSF based on observations of stellar fields. This is currently the only technique which takes the full time variation of the PSF between individual ACS exposures into account. We estimate that our PSF correction scheme reduces the systematic contribution to the shear correlation functions due to PSF distortions to < 2*10^{-6} for galaxy fields containing at least 10 stars. We perform a number of diagnostic tests indicating that the remaining level of systematics is consistent with zero for the GEMS and GOODS data confirming the success of our PSF correction scheme. For the parallel data we detect a low level of remaining systematics which we interpret to be caused by a lack of sufficient dithering of the data. Combining the shear estimate of the GEMS and GOODS observations using 96 galaxies arcmin^{-2} with the photometric redshift catalogue of the GOODS-MUSIC sample, we determine a local single field estimate for the mass power spectrum normalisation sigma_{8,CDFS}=0.52^{+0.11}_{-0.15} (stat) +/- 0.07 (sys) (68% confidence assuming Gaussian cosmic variance) at fixed Omega_m=0.3 for a LambdaCDM cosmology. We interpret this exceptionally low estimate to be due to a local under-density of the foreground structures in the CDFS.

Qualitatively new aspects of the (linear and non-linear) stability of sheared relativistic (slab) jets are analyzed. The linear problem has been solved for a wide range of jet models well inside the ultrarelativistic domain (flow Lorentz factors up to 20; specific internal energies $\approx 60c^2$). As a distinct feature of our work, we have combined the analytical linear approach with high-resolution relativistic hydrodynamical simulations, which has allowed us i) to identify, in the linear regime, resonant modes specific to the relativistic shear layer ii) to confirm the result of the linear analysis with numerical simulations and, iii) more interestingly, to follow the instability development through the non-linear regime. We find that very high-order reflection modes with dominant growth rates can modify the global, long-term stability of the relativistic flow. We discuss the dependence of these resonant modes on the jet flow Lorentz factor and specific internal energy, and on the shear layer thickness. The results could have potential applications in the field of extragalactic relativistic jets.

We derive the relativistic equations for stellar perturbations, including in a consistent way shear viscosity in the stress-energy tensor, and we numerically integrate our equations in the case of large viscosity. We consider the slow rotation approximation, and we neglect the coupling between polar and axial perturbations. In our approach, the frequency and damping time of the emitted gravitational radiation are directly obtained. We find that, approaching the inviscid limit from the finite viscosity case, the continuous spectrum is regularized. Constant density stars, polytropic stars, and stars with realistic equations of state are considered. In the case of constant density stars and polytropic stars, our results for the viscous damping times agree, within a factor two, with the usual estimates obtained by using the eigenfunctions of the inviscid limit. For realistic neutron stars, our numerical results give viscous damping times with the same dependence on mass and radius as previously estimated, but systematically larger of about 60%.

We show the existence of a strong trend between neutron star surface temperature and the dipolar component of the magnetic field extending through three orders of field magnitude, a range that includes magnetars, radio-quiet isolated neutron stars, and many ordinary radio pulsars. We suggest that this trend can be explained by the decay of currents in the crust over a time scale of few Myr. We estimate the minimum temperature that a NS with a given magnetic field can reach in this interpretation.

This paper presents new five passbands (UBVRI) optical wide-field imaging data accumulated as part of the DEEP Public Survey (DPS) carried out as a public survey by the ESO Imaging Survey (EIS) project. Out of the 3 square degrees originally proposed, the survey covers 2.75 square degrees, in at least one band (normally R), and 1.00 square degrees in five passbands. The median seeing, as measured in the final stacked images, is 0.97", ranging from 0.75" to 2.0". The median limiting magnitudes (AB system, 2" aperture, 5 sigma detection limit) are U_(AB)=25.65, B_(AB)=25.54, V_(AB)=25.18, R_(AB) = 24.8 and I_(AB)=24.12 mag, consistent with those proposed in the original survey design. The paper describes the observations and data reduction using the EIS Data Reduction System and its associated EIS/MVM library. The quality of the individual images were inspected, bad images discarded and the remaining used to produce final image stacks in each passband, from which sources have been extracted. Finally, the scientific quality of these final images and associated catalogs was assessed qualitatively by visual inspection and quantitatively by comparison of statistical measures derived from these data with those of other authors as well as model predictions, and from direct comparison with the results obtained from the reduction of the same dataset using an independent (hands-on) software system. Finally to illustrate one application of this survey, the results of a preliminary effort to identify sub-mJy radio sources are reported. To the limiting magnitude reached in the R and I passbands the success rate ranges from 66 to 81% (depending on the fields). These data are publicly available at CDS.

This paper presents new J and Ks data obtained from observations conducted at the ESO 3.5m New Technology Telescope using the SOFI camera. These data were taken as part of the ESO Imaging Survey Deep Public Survey (DPS) and significantly extend the earlier optical/infrared EIS-DEEP survey presented in a previous paper. The DPS-IR survey comprises two observing strategies: shallow Ks observations providing nearly full coverage of pointings with complementary multi-band optical data and deeper J and Ks observations of the central parts of these fields. The DPS-IR survey provides a coverage of roughly 2.1 square degrees in Ks with 0.63 square degrees to fainter magnitudes and also covered in J, over three independent regions of the sky. The goal of the present paper is to describe the observations, the data reduction procedures, and to present the final survey products. The astrometric solution with an estimated accuracy of <0.15" is based on the USNO catalog. The final stacked images presented here number 89 and 272, in J and Ks, respectively, the latter reflecting the larger surveyed area. The J and Ks images were taken with a median seeing of 0.77" and 0.8". The images reach a median 5sigma limiting magnitude of J_AB~23.06 in an aperture of 2", while the corresponding limiting magnitude in Ks_AB is ~21.41 and ~22.16 mag for the shallow and deep strategies. Overall, the observed limiting magnitudes are consistent with those originally proposed. The quality of the data has been assessed by comparing the measured magnitude of sources at the bright end directly with those reported by the 2MASS survey and at the faint end by comparing the counts of galaxies and stars with those of other surveys to comparable depth and to model predictions. The final science-grade catalogs and images are available at CDS.

We present a unified explanation for the observed properties of the isolated neutron star RX J0720.4-3125 by obtaining, for the first time, a self-consistent model that accounts simultaneously for the observed X-ray spectrum and optical excess, the pulsed fraction, the observed spectral feature around 0.3 keV, and the long--term spectral evolution. By fitting the parameters of our realistic self--consistent models to all archival XMM--Newton observations and available optical data, we show that all observed properties are consistent with a normal neutron star with a proper radius of about 12 km, a temperature at the magnetic pole of about 100 eV and a magnetic field strength of 2-3 x 10^{13} G, with no need to invoke additional emission and absorption components nor exotic internal composition. The observed variability of the effective temperature, strength of the spectral feature, and pulsed fraction are in good agreement with the predictions of our model in which the star is subject to free precession, producing changes in the angle between the magnetic field and the rotation axis of tens of degrees with a periodicity of 7 years, as pointed out by other authors on the basis of phenomenological models. In addition to the evidence of internal toroidal components, we also find strong evidence of non-dipolar magnetic fields, since all spectral properties are better reproduced with models with strong quadrupolar components.

The thermal emission from isolated neutron stars is not well understood. The X-ray spectrum is very close to a blackbody but there is a systematic optical excess flux with respect to the extrapolation to low energy of the best blackbody fit. This fact, in combination with the observed pulsations in the X-ray flux, can be explained by anisotropies in the surface temperature distribution.We study the thermal emission from neutron stars with strong magnetic fields in order to explain the origin of the anisotropy. We find (numerically) stationary solutions in axial symmetry of the heat transportequations in the neutron star crust and the condensed envelope. The anisotropy in the conductivity tensor is included consistently. The presence of magnetic fields of the expected strength leads to anisotropy in the surface temperature. Models with toroidal components similar to or larger than the poloidal field reproduce qualitatively the observed spectral properties and variability of isolated neutron stars. Our models also predict spectral features at energies between 0.2 and 0.6 keV.

This paper presents infrared data obtained from observations carried out at the ESO 3.5m New Technology Telescope (NTT) of the Hubble Deep Field South (HDF-S) and the Chandra Deep Field South (CDF-S). These data were taken as part of the ESO Imaging Survey (EIS) program, a public survey conducted by ESO to promote follow-up observations with the VLT. In the HDF-S field the infrared observations cover an area of ~53 square arcmin, encompassing the HST WFPC2 and STIS fields, in the JHKs passbands. The seeing measured in the final stacked images ranges from 0.79" to 1.22" and the median limiting magnitudes (AB system, 2" aperture, 5sigma detection limit) are J_AB~23.0, H_AB~22.8 and K_AB~23.0 mag. Less complete data are also available in JKs for the adjacent HST NICMOS field. For CDF-S, the infrared observations cover a total area of \~100 square arcmin, reaching median limiting magnitudes (as defined above) of J_AB~23.6 and K_AB~22.7 mag. For one CDF-S field H-band data are also available. This paper describes the observations and presents the results of new reductions carried out entirely through the un-supervised, high-throughput EIS Data Reduction System and its associated EIS/MVM C++-based image processing library developed, over the past 5 years, by the EIS project and now publicly available. The paper also presents source catalogs extracted from the final co-added images which are used to evaluate the scientific quality of the survey products, and hence the performance of the software. This is done comparing the results obtained in the present work with those obtained by other authors from independent data and/or reductions carried out with different software packages and techniques. The final science-grade catalogs and co-added images are available at CDS.

(Abridged) This paper presents the data recently released for the XMM-Newton/WFI survey carried out as part of the ESO Imaging Survey (EIS) project. The aim of this survey is to provide optical imaging follow-up data in BVRI for identification of serendipitously detected X-ray sources in selected XMM-Newton fields. In this paper, fully calibrated individual and stacked images of 12 fields as well as science-grade catalogs for the 8 fields located at high-galactic latitude are presented. The data covers an area of \sim 3 square degrees for each of the four passbands. The median limiting magnitudes (AB system, 2" aperture, 5\sigma detection limit) are 25.20, 24.92, 24.66, and 24.39 mag for B-, V-, R-, and I-band, respectively. These survey products, together with their logs, are available to the community for science exploitation in conjunction with their X-ray counterparts. Preliminary results from the X-ray/optical cross-correlation analysis show that about 61% of the detected X-ray point sources in deep XMM-Newton exposures have at least one optical counterpart within 2" radius down to R \simeq 25 mag, 50% of which are so faint as to require VLT observations thereby meeting one of the top requirements of the survey, namely to produce large samples for spectroscopic follow-up with the VLT, whereas only 15% of the objects have counterparts down to the DSS limiting magnitude.

We have extended the procedure to find the exact solution of the Riemann problem in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics in which the magnetic field of the initial states is tangential to the discontinuity and orthogonal to the flow velocity. The wave pattern produced after the break up of the initial discontinuity is analogous to the non--magnetic case and we show that the problem can be understood as a purely relativistic hydrodynamical problem with a modified equation of state. The new degree of freedom introduced by the non-zero component of the magnetic field results in interesting effects consisting in the change of the wave patterns for given initial thermodynamical states, in a similar way to the effects arising from the introduction of tangential velocities. Secondly, when the magnetic field dominates the thermodynamical pressure and energy, the wave speeds approach the speed of light leading to fast shocks and fast and arbitrarily thin rarefaction waves. Our approach is the first non-trivial exact solution of a Riemann problem in relativistic magnetohydrodynamics and it can also be of great interest to test numerical codes against known analytical or exact solutions.

We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e. Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic magnetohydrodynamics equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic magnetohydrodynamics in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick accretion disks around a black hole prone to the magnetorotational instability.

We have shown that the cluster-mass reconstruction method which combines strong and weak gravitational lensing data, developed in the first paper in the series, successfully reconstructs the mass distribution of a simulated cluster. In this paper we apply the method to the ground-based high-quality multi-colour data of RX J1347.5-1145, the most X-ray luminous cluster to date. A new analysis of the cluster core on very deep, multi-colour data analysis of VLT/FORS data reveals many more arc candidates than previously known for this cluster. The combined strong and weak lensing reconstruction confirms that the cluster is indeed very massive. If the redshift and identification of the multiple-image system as well as the redshift estimates of the source galaxies used for weak lensing are correct, we determine the enclosed cluster mass in a cylinder to M(<360 h^-1 kpc)= (1.2 +/- 0.3) 10^15 Msun. In addition the reconstructed mass distribution follows the distribution found with independent methods (X-ray measurements, SZ). With higher resolution (e.g. HST imaging data) more reliable multiple imaging information can be obtained and the reconstruction can be improved to accuracies greater than what is currently possible with weak and strong lensing techniques.

Within a fully relativistic framework, we derive and solve numerically the perturbation equations of relativistic stars, including the stresses produced by a non-vanishing shear viscosity in the stress-energy tensor. With this approach, the real and imaginary parts of the frequency of the modes are consistently obtained. We find that, approaching the inviscid limit from the finite viscosity case, the continuous spectrum is regularized and we can calculate the quasi-normal modes for stellar models that do not admit solutions at first order in perturbation theory when the coupling between the polar and axial perturbations is neglected. The viscous damping time is found to agree within factor 2 with the usual estimate obtained by using the eigenfunctions of the inviscid limit and some approximation for the energy dissipation integrals. We find that the frequencies and viscous damping times for relativistic $r-$modes lie between the Newtonian and Cowling results. We compare the results obtained with homogeneous, polytropic and realistic equations of state and find that the frequencies depend only on the rotation rate and on the compactness parameter (M/R), being almost independent of the equation of state. Our numerical results for realistic neutron stars give viscous damping times with the same dependence on mass and radius as previously estimated, but systematically larger of about 60%.

We investigate the thermal emission from magnetic neutron star surfaces in which the cohesive effects of the magnetic field have produced the condensation of the atmosphere and the external layers. This may happen for sufficiently cool atmospheres with moderately intense magnetic fields. The thermal emission from an isothermal bare surface of a neutron star shows no remarkable spectral features, but it is significantly depressed at energies below some threshold energy. However, since the thermal conductivity is very different in the normal and parallel directions to the magnetic field lines, the presence of the magnetic field is expected to produce a highly anisotropic temperature distribution, depending on the magnetic field geometry. In this case, the observed flux of such an object looks very similar to a BB spectrum, but depressed in a nearly constant factor at all energies. This results in a systematic underestimation of the area of the emitter (and therefore its size) by a factor 5-10 (2-3).

We present our image processing system for the reduction of optical imaging data from multi-chip cameras. In the framework of the Garching Bonn Deep Survey (GaBoDS; Schirmer et al. 2003) consisting of about 20 square degrees of high-quality data from WFI@MPG/ESO 2.2m, our group developed an imaging pipeline for the homogeneous and efficient processing of this large data set. Having weak gravitational lensing as the main science driver, our algorithms are optimised to produce deep co-added mosaics from individual exposures obtained from empty field observations. However, the modular design of our pipeline allows an easy adaption to different scientific applications. Our system has already been ported to a large variety of optical instruments and its products have been used in various scientific contexts. In this paper we give a thorough description of the algorithms used and a careful evaluation of the accuracies reached. This concerns the removal of the instrumental signature, the astrometric alignment, photometric calibration and the characterisation of final co-added mosaics. In addition we give a more general overview on the image reduction process and comment on observing strategies where they have significant influence on the data quality.

We study the conditions for convective instability in rotating, non-magnetic proto--neutron stars. The criteria that determine stability of nascent neutron stars are analogous to the Solberg--Hoiland conditions but including the presence of lepton gradients. Our results show that, for standard angular velocity profiles, convectively unstable modes with wave-vectors parallel to the rotation axis are suppressed by a stable angular momentum profile, while unstable modes with wave-vectors perpendicular to the axis remain unaltered. Since the wave-vector is perpendicular to the velocity perturbation, the directional selection of the unstable modes may result in fluid motions along the direction of the rotation axis. This occurs in rigidly rotating stars as well as in the inner core of differentially rotating stars. Our results provide a natural source of asymmetry for proto--neutron stars with the only requirement that angular velocities be of the order of the convective characteristic frequency.

Following the detection of a cosmic shear signal at the 30'' scale using archival parallel data from the STIS CCD camera onboard HST in Haemmerle et al. (2002), we analyzed a larger data set obtained from an HST GO pure parallel program. Although this data set is considerably larger than the one analyzed previously, we do not obtain a significant detection of the cosmic shear signal. The potential causes of this null result are the multiple systematics that plague the STIS CCD data, and in particular the degradation of the CCD charge transfer efficiency after 4 years in space.

We perform a weak lensing analysis on a previously reported dark-lens candidate on STIS Parallel data (Miralles et al. 2002). New VLT-data indicate that the reported signal originates from a small number of galaxies tangentially aligned towards the center of the STIS field but no signature for an extended mass distribution is found. We argue that we should be able to detect a massive cluster ($M\geq 3.2\times 10^{14}M_{\odot}$) through its lensing signal up to a redshift of $z\approx 0.6$ with our data. Also the double image hypothesis of two galaxies with very similar morphologies in the STIS data is ruled out with colour information.

We consider the transverse structure and stability properties of relativistic jets formed in the course of the collapse of a massive progenitor. Our numerical simulations show the presence of a strong shear in the bulk velocity of such jets. This shear can be responsible for a very rapid shear--driven instability that arises for any velocity profile. This conclusion has been confirmed both by numerical simulations and theoretical analysis. The instability leads to rapid fluctuations of the main hydrodynamical parameters (density, pressure, Lorentz factor, etc.). However, the perturbations of the density are effectively decoupled from those of the pressure because the beam of the jet is radiation--dominated. The characteristic growth time of instability is much shorter than the life time of the jet and, therefore, may lead to a complete turbulent beam. In the course of the non-linear evolution, these fluctuations may yield to internal shocks which can be randomly distributed in the jet. In the case that internal shocks in a ultrarelativistic outflow are responsible for the observed phenomenology of gamma-ray bursts, the proposed instability can well account for the short-term variability of gamma-ray light curves down to milliseconds.

We report on the marginal detection of cosmic shear on sub-arcmin scales with archive data from the STIS camera on board HST. For the measurement 121 galaxy fields with a field of view of 51'' x 51'' are used to obtain an rms cosmic shear of ~ 4% with 1.5sigma significance. This value is consistent with groundbased results obtained by other groups on larger scales, and with theoretical predictions for a standard LambdaCDM cosmology. To show the suitability of STIS for weak shear measurements we carefully investigated the stability of the PSF. We demonstrate that small temporal changes do not affect the cosmic shear measurement by more than ~ 10%. We also discuss the influence of various weighting and selection schemes for the galaxy ellipticities.

The stability properties of newly born neutron stars, or proto--neutron stars, are considered. We take into account dissipative processes, such as neutrino transport and viscosity, in the presence of a magnetic field. In order to find the regions of the star subject to different sorts of instability, we derive the general instability criteria and apply it to evolutionary models of PNSs. The influence of the magnetic field on instabilities is analyzed and the critical magnetic field stabilizing the star is obtained. In the light of our results, we estimate of the maximum poloidal magnetic field that might be present in young pulsars or magnetars.

We report the serendipituous discovery of a conspicuous alignment of galaxies in a field obtained through the STIS Parallel Shear Survey. This project collects randomly distributed 50'' x 50'' fields to investigate the cosmic shear effect on this scale. Analyzing the parallel observations having the Seyfert galaxy NGC625 as primary target, we recognized over the whole field of view a strong apparent tangential alignment of galaxy ellipticities towards the image center. The field shows several arclet-like features typical for images of massive galaxy clusters, but no obvious over-density of bright foreground galaxies. We also find a multiple image candidate. On the basis of the possible strong and weak lensing effect within the data, we discuss whether this could be compatible with a massive halo with no clear optical counterpart.

We present the X-ray energy spectrum of the lensed BAL QSO H1413+117 (the Cloverleaf) at z=2.56 observed with the Chandra X-ray observatory. We detected 293 photons in a 40 ks Advanced CCD Imaging Spectrometer (ACIS-S) observation. The X-ray image consists of four lensed image components, thus the photons are from the lensed QSO itself. The overall spectrum can be described with a power-law function heavily absorbed by neutral matter at a redshift consistent with the QSO redshift. This supports the idea that intrinsic absorption is significant for BAL QSOs. The spectral fit significantly (99% confidence) improves when we include an emission line. The centroid energy and intrinsic width (Gaussian \sigma) of the line are 6.21 \pm 0.16 keV and 220 ^{+270}_{-130} eV (90% errors), respectively, in the QSO rest frame, assuming the absorbed power-law as the continuum. The equivalent width of the line in the QSO rest frame is 960 ^{+1400}_{-480} eV. We suggest that the large equivalent width, the centroid energy, and the line broadness can be explained by iron K emission arising from X-ray reprocessing in the BAL flow, assuming it has a conical thin-sheet structure.

Since the Universe is inhomogeneous on scales well below the Hubble radius, light bundles from distant galaxies are deflected and distorted by the tidal gravitational field of the large-scale matter distribution as they propagate through the Universe. Two-point statistical measures of the observed ellipticities, like the dispersion within a finite aperture or the ellipticity cross-correlation, can be related to the power spectrum of the large-scale structure. The measurement of cosmic shear (especially on small angular scales) can thus be used to constrain cosmological parameters and to test cosmological structure formation in the non-linear regime, without any assumptions about the relation between luminous and dark matter. In this paper we will present preliminary cosmic shear measurements on sub-arcminute scales, obtained from archival STIS parallel data. The high angular resolution of HST, together with the sensitivity and PSF-stability of STIS, allows us to measure cosmic shear along many independent lines-of-sight. Ongoing STIS parallel observations, currently being carried out in the frame of a big GO program (8562+9248), will greatly increase the number of available useful fields and will enable us to measure cosmic shear with higher accuracy on sub-arcminute scales.

The measurement of cosmic shear requires deep imaging with high image quality on many lines of sight to sample the statistics of large-scale structure. The expected distortion of galaxy images by cosmic shear on the STIS angular scale is a few percent, therefore the PSF anisotropy has to be understood and controlled to an accuracy better than 1%. In this poster we present the analysis of the PSF of STIS and a preliminary cosmic shear measurement using archival data from the STIS pure parallel program to show that the STIS camera on-board HST is well suited for our project. The data reduction and catalog production are described in an accompanying paper (astro-ph/0102330).

In June 1997, parallel observations using the Space Telescope Imaging Spectrograph (STIS) on the HST started to be taken in substantial numbers along many different lines-of-sight. We are using the imaging data to investigate the distortion of background galaxies by the gravitational field of the large scale matter distribution, also known as Cosmic Shear. This poster presents the data and the catalog production that leads to the cosmic shear result presented in the poster "First Cosmic Shear results from STIS parallel program archive data" (Haemmerle et al., this conference). The data is publicly available also at http://www.stecf.org/projects/shear .

Following the second HST servicing mission in 1997 when the STIS instrument was installed and the capability for parallel observations was enhanced, a substantial archive of non-proprietary parallel data has been accumulating. In this paper, we discuss the use of unfiltered STIS imaging data for a project that requires deep observations along as many independent lines-of-sight as possible. We have developed a technique to determine which datasets in the archive can safely be co-added together and have developed an iterative co-addition technique which enabled us to produce 498 high-quality, deep images. The principal motivation for this work is to measure the Cosmic Shear on small angular scales and a value derived from these data will be presented in a subsequent paper. A valuable by-product of this work is a set of high quality combined fields which can be used for other projects. The data are publicly available at http://www.stecf.org/projects/shear/

In this paper we study the accuracy of photometric redshifts computed through a standard SED fitting procedure, where SEDs are obtained from broad-band photometry. We present our public code hyperz, which is presently available on the web. We introduce the method and we discuss the expected influence of the different observational conditions and theoretical assumptions. In particular, the set of templates used in the minimization procedure (age, metallicity, reddening, absorption in the Lyman forest, ...) is studied in detail, through both real and simulated data. The expected accuracy of photometric redshifts, as well as the fraction of catastrophic identifications and wrong detections, is given as a function of the redshift range, the set of filters considered, and the photometric accuracy. Special attention is paid to the results expected from real data.

We present simulations of the evolution of a proto-neutron star in which kaon-condensed matter might exist, including the effects of finite temperature and trapped neutrinos. The phase transition from pure nucleonic matter to the kaon condensate phase is described using Gibbs' rules for phase equilibrium, which permit the existence of a mixed phase. A general property of neutron stars containing kaon condensates, as well as other forms of strangeness, is that the maximum mass for cold, neutrino-free matter can be less than the maximum mass for matter containing trapped neutrinos or which has a finite entropy. A proto-neutron star formed with a baryon mass exceeding that of the maximum mass of cold, neutrino-free matter is therefore metastable, that is, it will collapse to a black hole at some time during the Kelvin-Helmholtz cooling stage. The effects of kaon condensation on metastable stars are dramatic. In these cases, the neutrino signal from a hypothetical galactic supernova (distance $\sim8.5$ kpc) will stop suddenly, generally at a level above the background in the SuperK and SNO detectors, which have low energy thresholds and backgrounds. This is in contrast to the case of a stable star, for which the signal exponentially decays, eventually disappearing into the background. We find the lifetimes of kaon-condensed metastable stars to be restricted to the range 40--70 s and weakly dependent on the proto-neutron star mass, in sharp contrast to the significantly larger mass dependence and range (1--100 s) of hyperon-rich metastable stars.

The linear hydrodynamic stability of proto-neutron stars (PNSs) is considered taking into account dissipative processes such as neutrino transport and viscosity. We obtain the general instability criteria which differ essentially from the well-known Ledoux criterion used in previous studies. We apply the criteria to evolutive models of PNSs that, in general, can be subject to the various known regimes such as neutron fingers and convective instabilities. Our results indicate that the fingers instability arises in a more extended region of the stellar volume and lasts a longer time than expected.

We study the mathematical character of the angular moment equations of radiative transfer in spherical symmetry and conclude that the system is hyperbolic for general forms of the closure relation found in the literature. Hyperbolicity and causality preservation lead to mathematical conditions allowing to establish a useful characterization of the closure relations. We apply numerical methods specifically designed to solve hyperbolic systems of conservation laws (the so-called Godunov-type methods), to calculate numerical solutions of the radiation transport equations in a static background. The feasibility of the method in any kind of regime, from diffusion to free-streaming, is demonstrated by a number of numerical tests and the effect of the choice of the closure relation on the results is discussed.

Our contribution concerns with the numerical solution of the 3D general relativistic hydrodynamical system of equations within the framework of the 3+1 formalism. We summarize the theoretical ingredients which are necessary in order to build up a numerical scheme based on the solution of local Riemann problems. Hence, the full spectral decomposition of the Jacobian matrices of the system, i.e., the eigenvalues and the right and left eigenvectors, is explicitly shown. An alternative approach consists in using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits an accurate description of numerical general relativistic hydrodynamics.

We discuss the first results obtained on the study of a sample of high-z galaxies (2 < z < 7), using the gravitational amplification effect in the core of lensing clusters. Sources are located close to the critical lines in clusters with well constrained mass distributions, and selected through photometric redshifts, computed on a large wavelength domain, and lens inversion techniques.

We discuss a method aiming to use photometric redshifts in lensing clusters to access the population of distant background sources. The amplification provided by gravitational lensing allows to calibrate photometric redshifts 1 to 3 magnitudes deeper with respect to field studies. We present and discuss some results obtained with this procedure on extensively studied clusters, as well as the future issues of this project.

We present the first results on the identification and study of very distant field galaxies in the core of cluster-lenses, using a selection criterium based on both lens modelling and photometric redshifts. We concentrate on two multiple-imaged sources at z=4.05 in the cluster A2390. The 2 objects presented in this paper, namely H3 and H5, were identified through lens modelling as multiple images of high-redshift sources at z>3.5. We confirm the excellent agreement between this identification and both their photometric redshifts and morphologies. Our CFHT/WHT program for a systematic redshift survey of arcs in clusters has allowed to obtain a set of spectra on 3 different images at z~4: the brightest image of H3, which redshift was already confirmed by Frye & Broadhurst (1998), and the two brightest images of H5. The later is then confirmed spectroscopically as a multiple image, giving a strong support to the lens model. The main feature in each of these spectra is a strong emission line, identified as Ly-alpha, leading to z=4.05 for both H3 and H5. The spectrophotometric properties of these galaxies are studied, in particular the degeneracy in the parameter-space defined by the SFR type, age, metallicity and reddening. H3 and H5 are intrinsically bright and clumpy sources located ~100 kpc part on the source plane, with mean metallicities compatible with a fraction of solar or even solar values. All these results seem to favour a hierarchical merging scenario, where we are actually seeing a relatively advanced step for these 2 z~4 objects, with stars forming locally and efficiently from a preenriched gas.

We discuss a method to build up and study a sample of distant galaxies, with 2<z<5, using the gravitational amplification effect in cluster-lenses for which the mass distribution is well known. The candidates are selected close to the critical lines at high-z, through photometric redshift estimates on a large wavelength range. With respect to other methods, this procedure allows to reduce the selection biases in luminosity, thanks to the cluster magnification of ~1 magnitude, and towards active SFRs when the near-IR domain is included. Some new results obtained with this technique are presented and discussed.

We present a general and practical procedure to solve the general relativistic hydrodynamic equations by using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits accurate numerical calculations of general relativistic hydrodynamics problems using the numerical tools developed for the special relativistic case with negligible computational cost. The feasibility of the method has been confirmed by a number of numerical experiments.

We study the thermal and chemical evolution during the Kelvin-Helmholtz phase of the birth of a neutron star, employing neutrino opacities that are consistently calculated with the underlying equation of state (EOS). Expressions for the diffusion coefficients appropriate for general relativistic neutrino transport in the equilibrium diffusion approximation are derived. The diffusion coefficients are evaluated using a field-theoretical finite temperature EOS that includes the possible presence of hyperons. The variation of the diffusion coefficients is studied as a function of EOS and compositional parameters. We present results from numerical simulations of protoneutron star cooling for internal stellar properties as well as emitted neutrino energies and luminosities. We discuss the influence of the initial stellar model, the total mass, the underlying EOS, and the addition of hyperons on the evolution of the protoneutron star and upon the expected signal in terrestrial detectors.

We consider Joule heating caused by dissipation of the magnetic field in the neutron star crust. This mechanism may be efficient in maintaining a relatively high surface temperature in very old neutron stars. Calculations of the thermal evolution show that, at the late evolutionary stage ($t \geq 10$ Myr), the luminosity of the neutron star is approximately equal to the energy released due to the field dissipation and is practically independent of the atmosphere models. At this stage, the surface temperature can be of the order of $3 \times 10^{4} - 10^{5}$K. Joule heating can maintain this high temperature during extremely long time ($\geq 100$ Myr), comparable with the decay time of the magnetic field.

We calculate the Legendre expansion of the rate of the process $\nu + \bar{\nu} \leftrightarrow e^+ + e^-$ up to 3rd order extending previous results of other authors which only consider the 0th and 1st order terms. Using different closure relations for the moment equations of the radiative transfer equation we discuss the physical implications of taking into account quadratic and cubic terms on the energy deposition outside the neutrinosphere in a simplified model. The main conclusion is that 2nd order is necessary in the semi-transparent region and gives good results if an appropriate closure relation is used.

This paper addresses the problem of detecting high-redshift clustering in deep photometric surveys. We have used photometric redshifts to select different samples of galaxies in the HDF, in order to study their clustering properties within the 2.5<z<4.5 domain, in redshift bins of 0.5, using different techniques. A strong clustering signal is found for redshift bins including z~3.4 and leading to a present correlation length of r_0 = 4.1 +- 0.8 h^-1 Mpc (q_0=0.1) assuming linear evolution. An excess appears in the correlation function of chip 2 with respect to the fit. This excess could be associated to the main structure detected in this field, which contains 20% of the objects identified at 3.4<z<3.9. Its dimensions are 3 h^-1 Mpc x 0.5 h^-1 Mpc. The galaxies at 3.4<z<3.9 exhibit a SFR of a few solar masses per year, but their comoving density is a factor of ~50 higher than the population of star-forming galaxies reported by Steidel et al. (1996b). The resulting star formation rate density is at least 1.1 10^-2 h Msun/yr/Mpc^3 (1.8 10^-2 h Msun/yr/Mpc^3) with q_0=0.1(0.5), slightly higher than the results by Madau et al. (1996) at 2.5<z<3.5, and then incompatible with a global decrease of the star formation in this redshift domain. These results on star formation and clustering are consistent with a hierarchical scenario for galaxy formation.

We present an optical and near-infrared survey of galaxies in nearby clusters aimed at determining fundamental quantities of galaxies, such as multivariate luminosity function and color distribution for each Hubble type. The main characteristics of our survey are completeness in absolute magnitude, wide wavelength coverage and faint limiting magnitudes.

This paper describes new results on the identification of the complex gravita- tional lens responsible for the double quasar Q2345+007. A gravitational shear field was detected recently about 45'' away from the QSO, centered on an excess of faint blue galaxies. We present deep photometric data in the near-IR (J & K'), which are used together with optical BRI photometry to build spectral ener- gy distributions for all the objects, and to derive a photometric redshift esti- mate by comparison with synthetic spectrophotometric data. We propose a statis- tical method to analyse the distribution in z. An excess of galaxies at z about 0.75 is detected in this field, with a 2D distribution showing a maximum located at the center of the weak-shear field. The z inferred for this overdensity (a distant cluster of galaxies) corresponds to an absorption system in the spectrum of the quasar B. Most cluster-member candidates at z=0.75 are undergoing a star- formation process or are burst systems 1 to 3 Gyr old. We also discuss the exis- tence of other possible excesses of galaxies at different z planes.

In this paper we present a full general relativistic one-dimensional hydro-code which incorporates a modern high-resolution shock-capturing algorithm, with an approximate Riemann solver, for the correct modelling of formation and propagation of strong shocks. The efficiency of this code in treating strong shocks is demonstrated by some numerical experiments. The interest of this technique in several astrophysical scenarios is discussed.

The equation of state of an ideal Fermi gas is expressed in terms of Fermi-Dirac integrals. We give formulae for evaluation the Fermi-Dirac integrals of orders 1/2, 3/2, and 5/2 and their derivatives in various limits of non- and extreme degeneracy and relativity. We provide tables and a Fortran subroutine for numerical evaluation of the integrals and derivatives when a limit does not apply. The functions can be evaluated to better than 1% accuracy for any temperature and density using these methods.