As the sensitivity of gravitational wave (GW) instruments improves and new networks start operating, hundreds of merging stellar-mass black holes (SBHs) and intermediate-mass black holes (IMBHs) are expected to be observed in the next few years. The origin and distribution of SBH and IMBH binaries in various dynamical environments is a fundamental scientific question in GW astronomy. In this paper we discuss the electromagnetic window provided by tidal disruption events (TDEs) of stars, into the assembly and merger of binary SBHs and IMBHs in nuclear star clusters (NSCs). We discuss how the host NSC mass and density and the slope of the black-hole mass function set the orbital properties and the masses of the binaries that undergo a TDE. For typical NSC properties, we predict a merger rate of $\sim 10^{-6}$--$10^{-7}\ {\rm yr}^{-1}$ per galaxy. The lightcurve of TDEs in NSCs could be interrupted and modulated by the companion black hole on the orbital period of the binary. The ejected mass associated with TDEs can produce optical transients of luminosity $\sim 10^{42}$--$10^{44}$ erg s$^{-1}$ over timescales of days to months. These should be readily detectable by optical transient surveys such as the Zwicky Transient Facility and LSST.

In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analytic- and population synthesis-based calculations performed using the \texttt{SeBa} code to constrain the expected frequency of such systems, given a realistic initial population of MS triples, and provide the predicted distributions of orbital periods. We identify the parameter space for triples that can accommodate a CBD, to inform future numerical simulations of suitable initial conditions. We find that $\lesssim$ 10\% of all MS triples should be able to accommodate a CBD around the inner binary, and compute lower limits for the production rates. This scenario broadly predicts mergers of near equal-mass binaries, producing blue stragglers (BSs), Type Ia supernovae, gamma ray bursts and gravitational wave-induced mergers, along with the presence of an outer WD tertiary companion. We compare our predicted distributions to a sample of field BS binaries, and argue that our proposed mechanism explains the observed range of orbital periods. Finally, the mechanism considered here could produce hypervelocity MS stars, WDs and even millisecond pulsars with masses close to the Chandrasekhar mass limit, and be used to constrain the maximum remnant masses at the time of any supernova explosion.

GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15 to 30 degrees away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativistic baryon wind from the merger remnant, resulting in a structured outflow with a narrow core and broad wings. In this paper, we couple the observational constraints on the structured outflow with a model for the jet-wind interaction to constrain the intrinsic properties with which the jet was launched by the central engine, including its time delay from the merger event. Using wind prescriptions inspired by magnetized BNS merger simulations, we find that the jet was launched within about 0.4 s from the merger, implying that the 1.74 s observed delay was dominated by the fireball propagation up to the photospheric radius. We also constrain, for the first time for any gamma-ray burst, the jet opening angle at injection and set a lower limit to its asymptotic Lorentz factor. These findings suggest an initially Poynting-flux dominated jet, launched via electromagnetic processes. If the jet was powered by an accreting black hole, they also provide a significant constraint on the survival time of the metastable neutron star remnant.

The magnetar Swift ,J1818.0-1607 was discovered in March 2020 when Swift detected a 9 ms hard X-ray burst and a long-lived outburst. Prompt X-ray observations revealed a spin period of 1.36 s, soon confirmed by the discovery of radio pulsations. We report here on the analysis of the Swift burst and follow-up X-ray and radio observations. The burst average luminosity was $L_{\rm burst} \sim2\times 10^{39}$ erg/s (at 4.8 kpc). Simultaneous observations with XMM-Newton and NuSTAR three days after the burst provided a source spectrum well fit by an absorbed blackbody ($N_{\rm H} = (1.13\pm0.03) \times 10^{23}$ cm$^{-2}$ and $kT = 1.16\pm0.03$ keV) plus a power-law ($\Gamma=0.0\pm1.3$) in the 1-20 keV band, with a luminosity of $\sim$$8\times10^{34}$ erg/s, dominated by the blackbody emission. From our timing analysis, we derive a dipolar magnetic field $B \sim 7\times10^{14}$ G, spin-down luminosity $\dot{E}_{\rm rot} \sim 1.4\times10^{36}$ erg/s and characteristic age of 240 yr, the shortest currently known. Archival observations led to an upper limit on the quiescent luminosity $<$$5.5\times10^{33}$ erg/s, lower than the value expected from magnetar cooling models at the source characteristic age. A 1 hr radio observation with the Sardinia Radio Telescope taken about 1 week after the X-ray burst detected a number of strong and short radio pulses at 1.5 GHz, in addition to regular pulsed emission; they were emitted at an average rate 0.9 min$^{-1}$ and accounted for $\sim$50% of the total pulsed radio fluence. We conclude that Swift ,J1818.0-1607 is a peculiar magnetar belonging to the small, diverse group of young neutron stars with properties straddling those of rotationally and magnetically powered pulsars. Future observations will make a better estimation of the age possible by measuring the spin-down rate in quiescence.

After the identification of a candidate $\gamma$-ray transient in the error region of the binary black hole (BBH) merger GW150914 by the \textit{Fermi} satellite, the question of whether BBH mergers can be associated to electromagnetic counterparts remains highly debated. Here, we present radio follow-up observations of GW170608, a BBH merger that occurred during the second observing run (O2) of the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO). Our radio follow up focused on a specific field contained in the GW170608 sky localization area, where a candidate high-energy transient was detected by the \textit{Fermi} Large Area Telescope (LAT). We make use of data collected at 1.4\,GHz with the Karl G. Jansky Very Large Array (VLA), as well as with the VLA Low-band Ionosphere and Transient Experiment (VLITE). Our analysis is sensitive to potential radio afterglows with luminosity densities $L_{\rm 1.4\,GHz}\gtrsim 6\times10^{28}$\,erg\,s$^{-1}$\,Hz$^{-1}$. In the most optimistic theoretical models, $\approx 20\%$ of BBH events occurring in massive hosts could be associated with outflows as radio luminous as this. Although we find no evidence for the presence of a radio counterpart associated with the LAT candidate in the GW170608 error region, our analysis demonstrates the feasibility of future radio follow-up observations of well localized BBHs. Comparing our radio upper-limits with theoretical expectations for the radio afterglows potentially associated with jets launched in BBH mergers, we find that for jets of energy $\approx 10^{49}$\,erg seen on-axis, only jet angles $\theta_{jet}\gtrsim 40^{\circ}$ are compatible with the observations.

Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.

The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has opened interesting debates about the number, age and nature of pulsars expected in the Galactic center region. In this work, we present extensive X-ray monitoring of the outburst of SGR J1745-2900 using the Chandra X-ray Observatory, the only instrument with the spatial resolution to distinguish the magnetar from the supermassive black hole (2.4" angular distance). It was monitored from its outburst onset in April 2013 until August 2019, collecting more than fifty Chandra observations for a total of more than 2.3 Ms of data. Soon after the outburst onset, the magnetar emission settled onto a purely thermal emission state that cooled from a temperature of about 0.9 to 0.6 keV over 6 years. The pulsar timing properties showed at least two changes in the period derivative, increasing by a factor of about 4 during the outburst decay. We find that the long-term properties of this outburst challenge current models for the magnetar outbursts.

The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not predicted by the core accretion theory of planet formation. Here we suggest a novel mechanism, in which the giant planet is born around a more typical Sun-like star ($M_{*,1}$), but is subsequently exchanged during a dynamical interaction with a flyby low-mass star ($M_{*,2}$). We perform state-of-the-art $N$-body simulations with $M_{*,1}=1M_\odot$ and $M_{*,2}=0.1M_\odot$ to study the statistical outcomes of this interaction, and show that exchanges result in high eccentricities for the new orbit around the low-mass star, while about half of the outcomes result in tighter orbits than the planet had around its birth star. We numerically compute the cross section for planet exchange, and show that an upper limit for the probability per planetary system to have undergone such an event is $\Gamma\sim 4.4(M_{\rm c}/100M_\odot)^{-2}(a_{\rm p}/{\rm AU}) (\sigma/1\,{\rm km}\,{\rm s}^{-1})^{5}$Gyr$^{-1}$, where $a_{\rm p}$ is the planet semi-major axis around the birth star, $\sigma$ the velocity dispersion of the star cluster, and $M_{\rm c}$ the total mass of the star cluster. Hence these planet exchanges could be relatively common for stars born in open clusters and groups, should already be observed in the exoplanet database, and provide new avenues to create unexpected planetary architectures.

The discovery of Exoplanetary Systems has challenged some of the theories of planet formation, which assume unperturbed evolution of the host star and its planets. However, in star clusters the interactions with flyby stars and binaries may be relatively common during the lifetime of a planetary system. Here, via high-resolution $N$-body simulations of star-planet systems perturbed by interlopers (stars and binaries), we explore the reconfiguration to the planetary system due to the encounters. In particular, via an exploration focused on the strong scattering regime, we derive the fraction of encounters which result in planet ejections, planet transfers and collisions by the interloper star/binary, as a function of the characteristics of the environment (density, velocity dispersion), and for different masses of the flyby star/binary. We find that binary interlopers can significantly increase the cross section of planet ejections and collisions, while they only slightly change the cross section for planet transfers. Therefore, in environments with high binary fractions, floating planets are expected to be relatively common, while in environments with low binary fractions, where the cross sections of planet ejection and transfer are comparable, the rate of planet exchanges between two stars will be comparable to the rate of production of free-floating planets.

Stars passing too close to a black hole can produce tidal disruption events (TDEs), when the tidal force across the star exceeds the gravitational force that binds it. TDEs have usually been discussed in relation to massive black holes that reside in the centres of galaxies or lurk in star clusters. In this paper, we investigate the possibility that triple stars hosting a stellar black hole (SBH) may be sources of TDEs. We start from a triple system made up of three main sequence (MS) stars and model the supernova (SN) kick event that led to the production of an inner binary comprised of a SBH. We evolve these triples in isolation with a high precision $N$-body code and study their TDEs as a result of Kozai-Lidov oscillations. We explore a variety of distributions of natal kicks imparted during the SN event, various maximum initial separations for the triples, and different distributions of eccentricities. We show that the main parameter that governs the properties of the SBH-MS binaries which produce a TDE in triples is the mean velocity of the natal kick distribution. Smaller $\sigma$'s lead to larger inner and outer semi-major axes of the systems that undergo a TDE, smaller SBH masses, and longer timescales. We find that the fraction of systems that produce a TDE is roughly independent of the initial conditions, while estimate a TDE rate of $7.3\times 10^{-5}-4.1 \ \mathrm{yr}^{-1}$, depending on the prescriptions adopted for the SBH natal kicks. This rate is almost comparable to the expected TDE rate for massive black holes.

The discovery of GW170817, the merger of a binary neutron star (NS) triggered by a gravitational wave detection by LIGO and Virgo, has opened a new window of exploration in the physics of NSs and their cosmological role. Among the important quantities to measure are the mass and velocity of the ejecta produced by the tidally disrupted NSs and the delay - if any - between the merger and the launching of a relativistic jet. These encode information on the equation of state of the NS, the nature of the merger remnant, and the jet launching mechanism, as well as yielding an estimate of the mass available for r-process nucleosynthesis. Here we derive analytic estimates for the structure of jets expanding in environments with different density, velocity, and radial extent. We compute the jet-cocoon structure and the properties of the broadband afterglow emission as a function of the ejecta mass, velocity, and time delay between merger and launch of the jet. We show that modeling of the afterglow light curve can constrain the ejecta properties and, in turn, the physics of neutron density matter. Our results increase the interpretative power of electromagnetic observations by allowing for a direct connection with the merger physics.

Nuclear star clusters that surround supermassive black holes (SMBHs) in galactic nuclei are thought to contain large numbers of black holes (BHs) and neutron stars (NSs), a fraction of which form binaries and could merge by Kozai-Lidov oscillations (KL). Triple compact objects are likely to be present, given what is known about the multiplicity of massive stars, whose life ends either as a NS or a BH. In this paper, we present a new possible scenario for merging BHs and NSs in galactic nuclei. We study the evolution of a triple black hole (BH) or neutron star (NS) system orbiting an SMBH in a galactic nucleus by means of direct high-precision $N$-body simulations, including Post-Newtonian terms. We find that the four-body dynamical interactions can increase the KL angle window for mergers compared to the binary case and make BH and NS binaries merge on shorter timescales. We show that the merger fraction can be up to $\sim 5$--$8$ times higher for triples than for binaries. Therefore, even if the triple fraction is only $\sim 10\%$--$20\%$ of the binary fraction, they could contribute to the merger events observed by LIGO/VIRGO in comparable numbers.

We combine sophisticated high precision scattering experiments, together with results from the Millenium-II simulation, to compute the cosmic merger rate of bound compact object (CO) binaries dynamically interacting with supermassive black hole binaries (SMBHBs). We consider binaries composed of white dwarfs (WDs), neutron stars (NSs) and black holes (BHs). The overall merger rates for WD-WD, NS-NS, BH-BH, BH-NS binaries and EBBH (eccentric binaries of black holes) from redshift $\sim 5$ are found to be $4.32\times 10^3\,{\rm yr}^{-1}$($5.93\times10^2\,{\rm yr}^{-1}$ for Type Ia SNe), $82.7\,{\rm yr}^{-1}$, $96.3\,{\rm yr}^{-1}$, $13.1\,{\rm yr}^{-1}$ and $148\,{\rm yr}^{-1}$ , respectively, for a nominal CO binary fraction in the Galactic centre of 0.1. We calculate the distance ($R$) distribution of the merger sites with respect to the host galaxies of the binaries. The distribution shows a wide range of distances up to $\sim$Mpc; this tail is produced by escaped hyper-velocity CO binaries. Due to the differences in the matter density of the surrounding environment, merger events with different $R$ are expected to display significantly different signatures in their EM counterparts. In particular, merger events (and especially NS-NS) producing a relativistic jet but occurring in the intergalactic medium will have very weak afterglow radiation relative to their prompt emission. These events, which we call 'off-center', can only be produced from a close encounter between CO binaries and SMBHBs; hence the detection of such merger events would indicate the existence of nearby SMBHBs, and in particular with high mass ratio, produced in the aftermath of a major galaxy merger.

Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies, but for the first time perform an extensive statistical study of BH-BH, NS-NS and BH-NS binary mergers by means of direct high-precision regularized $N$-body simulations, including Post-Newtonian (PN) terms up to order PN2.5. We consider different SMBH masses, slopes for the BH mass function, binary semi-major axis and eccentricity distributions, and different spatial distributions for the binaries. We find that the merger rates are a decreasing function of the SMBH mass and are in the ranges $\sim 0.17$-$0.52 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, $\sim 0.06$-$0.10 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ and $\sim 0.04$-$0.16 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ for BH-BH, BH-NS and NS-NS binaries, respectively. However, the rate estimate from this channel remains highly uncertain and depends on the specific assumptions regarding the star-formation history in galactic nuclei and the supply rate of compact objects. We find that $\sim 10\%$--$20\%$ of the mergers enter the LIGO band with eccentricities $\gtrsim 0.1$. We also compare our results to the secular approximation, and show that $N$-body simulations generally predict a larger number of mergers. Finally, these events can also be observable via their electromagnetic counterparts, thus making these compact object mergers especially valuable for cosmological and astrophysical purposes.

Prior to the detection of black holes (BHs) via the gravitational waves (GWs) they generate at merger, the presence of BHs was inferred in X-ray binaries, mostly via dynamical measurements, with masses in the range between $\sim 5-20~M_\odot$. The LIGO discovery of the first BHs via GWs was surprising in that the two BHs that merged had masses of $35.6^{+4.8}_{-3.0}$ and $30.6^{+3.0}_{-4.4}\,M_\odot$, which are both above the range inferred from X-ray binaries. With 20 BH detections from the O1/O2 runs, the distribution of masses remains generally higher than the X-ray inferred one, while the effective spins are generally lower, suggesting that, at least in part, the GW-detected population might be of dynamical origin rather than produced by the common evolution of field binaries. Here we perform high-resolution N-body simulations of a cluster of isolated BHs with a range of initial mass spectra and upper mass cut-offs, and study the resulting binary mass spectrum resulting from the dynamical interactions. Our clusters have properties similar to those of the massive remnants in an OB association $\sim 10 \, \mathrm{Myr}$ after formation. We perform a likelihood analysis for each of our dynamically-formed binary population against the data from the O1 and O2 LIGO/Virgo runs. We find that an initial mass spectrum $M_{\rm BH}\propto M^{-2.35}$ with an upper mass cutoff $M_{\rm max}\sim 50M_\odot$ is favored by the data, together with a slight preference for a merger rate that increases with redshift.

The conventional view of stellar-mass binary black hole (sBBH) mergers is that there should not be enough matter present to produce a detectable electromagnetic transient. However, there ARE a number of mechanisms for producing such a counterpart through accretion of matter from various reserves, charged black holes, or interactions with magnetic or exotic fields. After reviewing these mechanisms, we describe what we can learn from multi-messenger observations of sBBH mergers in the areas of stellar evolution and compact binary formation, cosmological measurements, searches for charged black holes, and tests of general relativity and fundamental physics. We end with some remarks about the need to support both observing capabilities and modeling in the next decade and beyond.

The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}

Unlike mergers of two compact objects containing a neutron star (NS), binary black hole (BBH) mergers are not accompanied by the production of tidally disrupted material, and hence lack the most direct source of accretion to power a jet and generate electromagnetic (EM) radiation. However, following a tentative detection by the Fermi GBM of a $\gamma$-ray counterpart to GW150914, several alternative ideas were proposed for driving a jet and producing EM radiation. {\em If} such jets were in fact produced, they would however lack the bright cocoon emission that makes jets from binary NSs bright also at large viewing angles. Here, via Monte Carlo simulations of a population of BBH mergers with properties consistent with those inferred from the existing LIGO/Virgo observations, and the angular emission characteristic of jets propagating into the interstellar medium, we derive limits on the allowed energetics and Lorentz factors of such jets from EM follow ups to GW-detected BBH merger events to date, and we make predictions which will help tighten these limits with broadband EM follow ups to events in future LIGO/Virgo runs. The condition that $\lesssim 1$ event out of 10 GW-detected BBH mergers be above the Fermi/GBM threshold imposes that any currently allowed emission model has to satisfy the condition $E_{\rm iso}/(10^{49} {\rm erg})(\theta_{jet}/20^\circ)\lesssim 1$.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

After two major outbursts in 2006 and 2011, on 2017 May 16 the magnetar CXOU J164710.2-455216, hosted within the massive star cluster Westerlund I, emitted a short (20 ms) burst, which marked the onset of a new active phase. We started a long-term monitoring campaign with Swift (45 observations), Chandra (5 observations) and NuSTAR (4 observations) from the activation until 2018 April. During the campaign, Swift BAT registered the occurrence of multiple bursts, accompanied by two other enhancements of the X-ray persistent flux. The long time span covered by our observations allowed us to study the spectral as well as the timing evolution of the source. After 11 months since the 2017 May outburst onset, the observed flux was about 15 times higher than its historical minimum level and a factor of 3 higher than the level reached after the 2006 outburst. This suggests that the crust has not fully relaxed to the quiescent level, or that the source quiescent level has changed following the multiple outburst activities in the past 10 years or so. This is another case of multiple outbursts from the same source on a yearly time scale, a somehow recently discovered behaviour in magnetars.

The degeneracy between disc and halo contributions in spiral galaxy rotation curves makes it difficult to obtain a full understanding of the distribution of baryons and dark matter in disc galaxies like our own Milky Way. Using mock data, we study how constraints on dark matter profiles obtained from kinematics, strong lensing, or a combination of the two are affected by assumptions about the halo model. We compare four different models: spherical isothermal and Navarro-Frenk-White halos, along with spherical and elliptical Burkert halos. For both kinematics and lensing we find examples where different models fit the data well but give enclosed masses that are inconsistent with the true (i.e., input) values. This is especially notable when the input and fit models differ in having cored or cuspy profiles (such as fitting an NFW model when the underlying dark matter distribution follows a different profile). We find that mass biases are more pronounced with lensing than with kinematics, and using both methods can help reduce the bias and provide stronger constraints on the dark matter distributions.

Hypervelocity binaries (HVBs) can be produced by a close encounter between stellar binaries and supermassive black hole (SMBH)-intermediate massive black hole (IMBH) binaries. We perform scattering experiments between tight stellar binaries and an SMBH-IMBH binary (within the observationally allowed parameters), using high-precision N-body integrations. Simulation results show that only a few HVBs can be produced in the Galactic center if it harbors a $4\times10^3M_\odot$ IMBH. However, we find that tight HVBs can be efficiently produced in prograde coplanar orbits with respect to the SMBH-IMBH binary. Therefore, the trajectories of the HVBs provide a promising way to constrain the IMBH orbital plane, if one exists. More importantly, we find that the ejected tight binaries are eccentric. For other galaxies hosting equal-mass SMBH binaries for which the HVB ejection rate is much higher relative to SMBH-IMBH binaries, these ejected compact binaries constitute promising extragalactic merger sources that could be detected by gravitational wave detectors.

We have performed high resolution $3-$dimensional compressible hydrodynamics simulations to investigate the effects of shocks and turbulence on energy transport into hot Jupiter atmospheres, under a variety of shear gradients. We focus on a local atmospheric region to accurately follow the small-scale structures of turbulence and shocks. We find that the effects of turbulence above and below a shear layer are different in scale and magnitude: below the shear layer, the effects of turbulence on the vertical energy transfer are local, generally $\lesssim2~\times(\text{scale height})$. However, turbulence can have a spatially and thermally-large influence on almost the entire region above the shear layer. We also find that shock formation is local and transient. Once the atmosphere becomes steady, the time-averaged heat energy flux at $P\sim 1$ bar is insignificant, on the order of 0.001\% of the incoming stellar flux with a shear motion at $P\sim 1$ mbar, and 0.1\% with a deeper shear layer at $P\sim 100$ mbar. Accordingly, the diffusion coefficient is higher for the deeper shear layer. Therefore, our results suggest that turbulence near less dense regions ($P\sim 1$ mbar) does not cause a sufficiently deep and large penetration of thermal energy to account for radius inflation in hot Jupiters, regardless of how violent the turbulence is. However, as the shear layer gets deeper, heat energy transfer becomes more effective throughout the atmosphere (upwards and downwards) due to a larger kinetic energy budget. Therefore, it is more important how deep turbulence occurs in the atmosphere, than how unstable the atmosphere is for effective energy transfer.

The binary neutron star (BNS) merger GW170817 was the first astrophysical source detected in gravitational waves and multi-wavelength electromagnetic radiation. The almost simultaneous observation of a pulse of gamma-rays proved that BNS mergers are associated with at least some short gamma-ray bursts (GRBs). However, the gamma-ray pulse was faint, casting doubts on the association of BNS mergers with the luminous, highly relativistic outflows of canonical short GRBs. Here we show that structured jets with a relativistic, energetic core surrounded by slower and less energetic wings produce afterglow emission that brightens characteristically with time, as recently seen in the afterglow of GW170817. Initially, we only see the relatively slow material moving towards us. As time passes, larger and larger sections of the outflow become visible, increasing the luminosity of the afterglow. The late appearance and increasing brightness of the multi-wavelength afterglow of GW170817 allow us to constrain the geometry of its ejecta and thus reveal the presence of an off-axis jet pointing about 30 degrees away from Earth. Our results confirm a single origin for BNS mergers and short GRBs: GW170817 produced a structured outflow with a highly relativistic core and a canonical short GRB. We did not see the bright burst because it was beamed away from Earth. However, approximately one in 20 mergers detected in gravitational waves will be accompanied by a bright, canonical short GRB.

We report on a new NuSTAR observation and on the ongoing Swift XRT monitoring campaign of the peculiar source 1E 161348-5055, located at the centre of the supernova remnant RCW 103, which is recovering from its last outburst in June 2016. The X-ray spectrum at the epoch of the NuSTAR observation can be described by either two absorbed blackbodies ($kT_{BB_1}$ ~ 0.5 keV, $kT_{BB_2}$ ~ 1.2 keV) or an absorbed blackbody plus a power law ($kT_{BB_1}$ ~ 0.6 keV, $\Gamma$ ~ 3.9). The observed flux was ~ 9 $\times$ 10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, ~ 3 times lower than what observed at the outburst onset, but about one order of magnitude higher than the historical quiescent level. A periodic modulation was detected at the known 6.67 hr periodicity. The spectral decomposition and evolution along the outburst decay are consistent with 1E 161348-5055 being a magnetar, the slowest ever detected.

Ultra-long gamma ray bursts (ULGRBs) are a distinct class of GRBs characterized by durations of several thousands of seconds, about two orders of magnitude longer than those of standard long GRBs (LGRBs). The driving engine of these events has not been uncovered yet, and ideas range from magnetars, to tidal disruption events, to extended massive stars, such as blue super giants (BSG). BSGs, a possible endpoint of stellar evolution, are attractive for the relatively long free-fall times of their envelopes, allowing accretion to power a long-lasting central engine. At the same time, their large radial extension poses a challenge to the emergence of a jet. Here we perform an end-to-end simulation aimed at assessing the viability of BSGs as ULGRB progenitors. The evolution to core collapse of a BSG star model is calculated with the MESA code. We then compute the accretion rate for the fraction of envelope material with enough angular momentum to circularize and form an accretion disk, and input the corresponding power into a jet which we evolve through the star envelope with the FLASH code. Our simulation shows that the jet can emerge, and the resulting light curves resemble those observed in ULGRBs, with durations $T_{90}$ ranging from $\approx$4000~s to $\approx10^4$ s depending on the viewing angle.

Binary black holes (BBHs) are one of the endpoints of isolated binary evolution, and their mergers a leading channel for gravitational wave events. Here, using the evolutionary code \textsc{StarTrack}, we study the statistical properties of the BBH population from isolated binary evolution for a range of progenitor star metallicities and BH natal kicks. We compute the mass function and the distribution of the primary BH spin $a$ as a result of mass accretion during the binary evolution, and find that this is not an efficient process to spin up BHs, producing an increase by at most $a\sim$~0.2--0.3 for very low natal BH spins. We further compute the distribution of merger sites within the host galaxy, after tracking the motion of the binaries in the potentials of a massive spiral, a massive elliptical, and a dwarf galaxy. We find that a fraction of 70-90\% of mergers in massive galaxies and of 40-60\% in dwarfs (range mostly sensitive to the natal kicks) is expected to occur inside of their hosts. The number density distribution at the merger sites further allows us to estimate the broadband luminosity distribution that BBH mergers would produce, \textit{if} associated with a kinetic energy release in an outflow, {which, as a reference, we assume at the level inferred for the \textit{Fermi} GBM counterpart to GW150914, with the understanding that current limits from the O1 and O2 runs would require such emission to be produced within a jet of angular size within $\lesssim 50^\circ$.}

We present a half-life formalism for describing the disruption of gravitationally-bound few-body systems, with a focus on binary-binary scattering. For negative total encounter energies, the four-body problem has three possible decay products in the point particle limit. For each decay product and a given set of initial conditions, we obtain directly from numerical scattering simulations the half-life for the distribution of disruption times. As in radioactive decay, the half-lives should provide a direct prediction for the relative fractions of each decay product. We test this prediction with simulated data and find good agreement with our hypothesis. We briefly discuss applications of this feature of the gravitational four-body problem to populations of black holes in globular clusters. This paper, the second in the series, builds on extending the remarkable similarity between gravitational chaos at the macroscopic scale and radioactive decay at the microscopic scale to larger-N systems.

The mass of the stellar precursor of supernova (SN) 1987A and the burst of neutrinos observed at the moment of the explosion are consistent with the core-collapse formation of a neutron star. However, no compelling evidence for the presence in SN 1987A of a compact object of any kind has been found yet in any band of the electromagnetic spectrum, prompting questions on whether the neutron star survived and, if it did, on its properties. Starting from the analysis of recent Chandra observations, here we appraise the current observational situation. We derived limits on the X-ray luminosity of a compact object with a nonthermal, Crab-pulsar-like spectrum of the order of $\approx$$(1$-$5)\times10^{35}$ erg s$^{-1}$, corresponding to limits on the rotational energy loss of a possible X-ray pulsar in SN 1987A of $\approx$$(0.5$-$1.5)\times10^{38}$ erg s$^{-1}$. However, a much brighter X-ray source cannot be excluded if, as is likely, it is enshrouded in a cloud of absorbing matter with metallicity similar to that expected in the outer layers of a massive star towards the end of its life. We found that other limits obtained from various arguments and observations in other energy ranges either are unbinding or allow a similar maximum luminosity of the order of $\approx$$10^{35}$ erg s$^{-1}$. We conclude that while a pulsar alike the one in the Crab Nebula in both luminosity and spectrum is hardly compatible with the observations, there is ample space for an `ordinary' X-ray-emitting young neutron star, born with normal initial spin period, temperature and magnetic field, to be hiding inside the evolving remnant of SN 1987A.

The evolution of main sequence binaries that reside in the Galactic Centre can be heavily influenced by the central supermassive black hole (SMBH). Due to these perturbative effects, the stellar binaries in dense environments are likely to experience mergers, collisions or ejections through secular and/or non-secular interactions. More direct interactions with the central SMBH are thought to produce hypervelocity stars (HVSs) and tidal disruption events (TDEs). In this paper, we use N-body simulations to study the dynamics of stellar binaries orbiting a central SMBH primary with an outer SMBH secondary orbiting this inner triple. The effects of the secondary SMBH on the event rates of HVSs, TDEs and stellar mergers are investigated, as a function of the SMBH-SMBH binary mass ratio. Our numerical experiments reveal that, relative to the isolated SMBH case, the TDE and HVS rates are enhanced for, respectively, the smallest and largest mass ratio SMBH-SMBH binaries. This suggests that the observed event rates of TDEs and HVSs have the potential to serve as a diagnostic of the mass ratio of a central SMBH-SMBH binary. The presence of a secondary SMBH also allows for the creation of hypervelocity binaries. Observations of these systems could thus constrain the presence of a secondary SMBH in the Galactic Centre.

Understanding of the role of X-rays for driving the thermal evolution of the intergalactic medium (IGM) at high redshifts is one of important questions in astrophysics. High-mass X-ray binaries (HMXBs) in early stellar populations are prime X-ray source; however, their formation efficiency is not well understood. Using $N$-body simulations, we estimate the HMXB formation rate via mutual gravitational interactions of nascent, small groups of the Population~III stars. We find that HMXBs form at a rate of one per $\gtrsim 10^{4}M_{\odot}$ in newly born stars, and that they emit with a power of $\sim 10^{41} {\rm erg}~{\rm s}^{-1}$ in the $2-10$ keV band per star formation rate (SFR). This value is a factor $\sim 10^{2}$ larger than what is observed in star forming galaxies at lower redshifts; the X-ray production from early HMXBs would have been even more copious, if they also formed \textit{in situ} or via migration in protostellar disks. Combining our results with earlier studies suggests that early HMXBs were highly effective at heating the IGM and leaving a strong 21 cm signature. We discuss broader implications of our results, such as the rate of long gamma-ray bursts from Population~III stars and the direct collapse channel for massive black hole formation.

We present the results of numerical simulations of the prompt emission of short-duration gamma-ray bursts. We consider emission from the relativistic jet, the mildly relativistic cocoon, and the non-relativistic shocked ambient material. We find that the cocoon material is confined between off-axis angles 15<theta<45 degrees and gives origin to X-ray transients with a duration of a few to ~10 seconds, delayed by a few seconds from the time of the merger. We also discuss the distance at which such transients can be detected, finding that it depends sensitively on the assumptions that are made about the radiation spectrum. Purely thermal cocoon transients are detectable only out to a few Mpc, Comptonized transients can instead be detected by the FERMI GBM out to several tens of Mpc.

Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evolution, and allow us to separately explore the role of dynamical friction and that of multi-body BH interactions on BH mergers. Using the computed merger rates, we infer the stochastic gravitational wave background (GWB). Our dynamical approach is a first attempt to study the dynamical evolution of multiple SMBHs in the host galaxies undergoing mergers with various mass ratios ($10^{-4} < q_{\star} < 1$). Our main result demonstrates that SMBH binaries are able to merge in both scenarios. In the empty loss cone case, we find that BHs merge via multi-body interactions, avoiding the `final parsec' problem, and entering the PTA band with substantial orbital eccentricity. Our full loss cone treatment, albeit more approximate, suggests that the eccentricity becomes even higher when GWs become dominant, leading to rapid coalescences (binary lifetime $\lesssim1 {\rm ~Gyr}$). Despite the lower merger rates in the empty loss cone case, due to their higher mass ratios and lower redshifts, the GWB in the full/empty loss cone models are comparable ($0.70\times10^{-15}$ and $0.53\times10^{-15}$ at a frequency of $1~{\rm yr}^{-1}$, respectively). Finally, we compute the effects of high eccentricities on the GWB spectrum.

We study how the viscosity of neutron star (NS) matter affects the distribution of tilt angles ($\chi$) between the spin and magnetic axes in young pulsars. Under the hypothesis that the NS shape is determined by the magnetically-induced deformation, and that the toroidal component of the internal magnetic field exceeds the poloidal one, we show that the dissipation of precessional motions by bulk viscosity can naturally produce a bi-modal distribution of tilt angles, as observed in radio/$\gamma$-ray pulsars, with a low probability of achieving $\chi \sim (20^\circ - 70^\circ)$ if the interior B-field is $\sim (10^{11} - 10^{15})$~G and the birth spin period is $\sim 10 - 300$~ms. As a corollary of the model, the idea that the NS shape is solely determined by the poloidal magnetic field, or by the centrifugal deformation of the crust, is found to be inconsistent with the tilt angle distribution in young pulsars. When applied to the Crab pulsar, with $\chi \sim 45^\circ - 70^\circ$ and birth spin $\sim$ 20 ms, our model implies that: (i) the magnetically-induced ellipticity is $\epsilon_B \gtrsim 3 \times 10^{-6}$; (ii) the measured positive $\dot{\chi} \sim 3.6 \times 10^{-12}$ rad s$^{-1}$ requires an additional viscous process, acting on a timescale $\lesssim 10^4$ yrs. We interpret the latter as crust-core coupling via mutual friction in the superfluid NS interior. One critical implication of our model is a GW signal at (twice) the spin frequency of the NS, due to $\epsilon_B \sim 10^{-6}$. This could be detectable by Advanced LIGO/Virgo operating at design sensitivity.

We study the formation of runaway stars due to binary-binary (2+2) interactions in young star-forming clusters and/or associations. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We focus on interactions that produce two single stars and a binary, and study the outcomes as a function of the depth of the background potential, within a range typical of cluster cores. As reference parameters for the observational properties, we use those observed for the system of runaway stars AE Aur and $\mu$ Col and binary $\iota$ Ori. We find that the outcome fractions have no appreciable dependence on the depth of the potential, and neither do the velocities of the ejected single stars. However, as the potential gets deeper and a larger fraction of binaries remain trapped, two binary populations emerge, with the escaped component having higher speeds and shorter semi-major axes than the trapped one. Additionally, we find that the relative angles between the ejected products are generally large. In particular, the angle between the ejected fastest star and the escaped binary is typically $\gtrsim 120-135^{\circ}$, with a peak at around $160^{\circ}$. However, as the potential gets deeper, the angle distribution becomes broader. Finally, we discuss the implications of our results for the interpretation of the properties of the runaway stars AE Aur and $\mu$ Col.

In this paper, we study the 3-body products (two single stars and a binary) of binary-binary (2+2) scattering interactions. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We derive analytically a simple formula relating the angle between the velocity vectors of the two ejected single stars and the orbital separation of the remaining binary. We compare our analytic formulation to numerical scattering simulations, and illustrate that the agreement is excellent, both in isolation and in a homogeneous background potential. Our results are ideally suited for application to the GAIA database, which is expected to identify many hundred runaway stars. The analytic relation presented here has the potential to identify runaway stars formed dynamically with high confidence. Finally, by applying our method to the runaways AE Aur and $\mu$ Col, we illustrate that it can be used to constrain the history of the background potential, which was denser than the presently observed density in the case of the Trapezium cluster.

Merging binary neutron stars (BNSs) represent the ultimate targets for multimessenger astronomy, being among the most promising sources of gravitational waves (GWs), and, at the same time, likely accompanied by a variety of electromagnetic counterparts across the entire spectrum, possibly including short gamma-ray bursts (SGRBs) and kilonova/macronova transients. Numerical relativity simulations play a central role in the study of these events. In particular, given the importance of magnetic fields, various aspects of this investigation require general relativistic magnetohydrodynamics (GRMHD). So far, most GRMHD simulations focused the attention on BNS mergers leading to the formation of a hypermassive NS, which, in turn, collapses within few tens of ms into a black hole surrounded by an accretion disk. However, recent observations suggest that a significant fraction of these systems could form a long-lived NS remnant, which will either collapse on much longer timescales or remain indefinitely stable. Despite the profound implications for the evolution and the emission properties of the system, a detailed investigation of this alternative evolution channel is still missing. Here, we follow this direction and present a first detailed GRMHD study of BNS mergers forming a long-lived NS. We consider magnetized binaries with different mass ratios and equations of state and analyze the structure of the NS remnants, the rotation profiles, the accretion disks, the evolution and amplification of magnetic fields, and the ejection of matter. Moreover, we discuss the connection with the central engine of SGRBs and provide order-of-magnitude estimates for the kilonova/macronova signal. Finally, we study the GW emission, with particular attention to the post-merger phase.

Following merger, a neutron star (NS) binary can produce roughly one of three different outcomes: (1) a stable NS, (2) a black hole (BH), or (3) a supra-massive, rotationally-supported NS, which then collapses to a BH following angular momentum losses. Which of these fates occur and in what proportion has important implications for the electromagnetic transient associated with the mergers and the expected gravitational wave signatures, which in turn depend on the high density equation of state (EOS). Here we combine relativistic calculations of NS masses using realistic EOSs with Monte Carlo population synthesis based on the mass distribution of NS binaries in our Galaxy to predict the distribution of fates expected. For many EOSs, a significant fraction of the remnants are NSs or supra-massive NSs. This lends support to scenarios where a quickly spinning, highly magnetized NS may be powering an electromagnetic transient. This also indicates that it will be important for future gravitational wave (GW) observatories to focus on high frequencies to study the post merger GW emission. Even in cases where individual GW events are too low in signal to noise to study the post merger signature in detail, the statistics of how many mergers produce NSs versus BHs can be compared with our work to constrain the EOS. To match short gamma-ray burst (SGRB) X-ray afterglow statistics, we find that the stiffest EOSs are ruled out. Furthermore, many popular EOSs require ~60-70% of SGRBs to be from NS-BH mergers rather than just binary NSs.

We perform a large suite of $N=4$ numerical scattering experiments between two identical binaries consisting of identical point particles in a (continuous) background potential. For investigative purposes, albeit without loss of generality, we assume that the potential corresponds to a uniform (natal or star-forming) gas medium. We explore a range of constant gas densities, from $n=10~ {\rm cm}^{-3}$ to $10^{5}~ {\rm cm}^{-3}$. These densities are relevant for various astrophysical environments, including molecular clouds (i.e., star-forming regions) and denser, fragmented cores within these clouds. Our primary goal is to characterize the effects of the background potential on the subsequent stellar dynamics. We consider the outcome probabilities as well as the properties of any binaries formed during the binary-binary encounters, such as the distributions of binary binding energies and eccentricities. We also present the final velocity distributions of the ejected single stars. The background potential has two important effects on the stellar dynamics: 1) The potential acts to reset the zero-point of the total system energy, which in turn affects the types and properties of the products of the encounter; 2) For higher $n$ and weakly bound systems (i.e., large semimajor axes), the stellar dynamics are significantly affected when stars become trapped in the potential, oscillating around the system centre of mass (CM). This, in turn, increases the number of scattering events between stars (single, binary or triple) near the CM and makes it harder for single stars to escape to infinity. This ultimately leads to the preferential ionization of triples and wide binaries and the survival of compact binaries, with the single stars escaping at very high ejection velocities.

A peculiar infrared ring-like structure was discovered by {\em Spitzer} around the strongly magnetised neutron star SGR 1900$+$14. This infrared structure was suggested to be due to a dust-free cavity, produced by the SGR Giant Flare occurred in 1998, and kept illuminated by surrounding stars. Using a 3D dust radiative transfer code, we aimed at reproducing the emission morphology and the integrated emission flux of this structure assuming different spatial distributions and densities for the dust, and different positions for the illuminating stars. We found that a dust-free ellipsoidal cavity can reproduce the shape, flux, and spectrum of the ring-like infrared emission, provided that the illuminating stars are inside the cavity and that the interstellar medium has high gas density ($n_H\sim$1000 cm$^{-3}$). We further constrain the emitting region to have a sharp inner boundary and to be significantly extended in the radial direction, possibly even just a cavity in a smooth molecular cloud. We discuss possible scenarios for the formation of the dustless cavity and the particular geometry that allows it to be IR-bright.

We study the contribution of binary black hole (BH-BH) mergers from the first, metal-free stars in the Universe (Pop III) to gravitational wave detection rates. Our study combines initial conditions for the formation of Pop III stars based on N-body simulations of binary formation (including rates, binary fraction, initial mass function, orbital separation and eccentricity distributions) with an updated model of stellar evolution specific for Pop III stars. We find that the merger rate of these Pop III BH-BH systems is relatively small (< 0.1 Gpc^-3 yr^-1) at low redshifts (z<2), where it can be compared with the LIGO empirical estimate of 9-240 Gpc^-3 yr^-1 (Abbott et al. 2016). The predicted rates are even smaller for Pop III double neutron star and black hole neutron star mergers. Our rates are compatible with those of Hartwig et al. (2016), but significantly smaller than those found in previous work (Bond & Carr 1984; Belczynski et al. 2004; Kinugawa et al. 2014, 2016). We explain the reasons for this discrepancy by means of detailed model comparisons and point out that (i) identification of Pop III BH-BH mergers may not be possible by advanced LIGO, and (ii) the level of stochastic gravitational wave background from Pop III mergers may be lower than recently estimated (Kowalska et al. 2012; Inayoshi et al. 2016; Dvorkin et al. 2016). We further estimate gravitational wave detection rates for third-generation interferometric detectors. Our calculations are relevant for low to moderately rotating Pop III stars. We can now exclude significant (> 1 per cent) contribution of these stars to low-redshift BH-BH mergers. However, it remains to be tested whether (and at what level) rapidly spinning Pop III stars (homogeneous evolution) can contribute to BH-BH mergers in the local Universe.

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

Flaring activity following gamma-ray bursts (GRBs), observed in both long and short GRBs, signals a long-term activity of the central engine. However, its production mechanism has remained elusive. Here we develop a quantitative model of the idea proposed by Perna et al. of a disc whose outer regions fragment due to the onset of gravitational instability. The self-gravitating clumps migrate through the disc and begin to evolve viscously when tidal and shearing torques break them apart. Our model consists of two ingredients: theoretical bolometric flare lightcurves whose shape (width, skewness) is largely insensitive to the model parameters, and a spectral correction to match the bandpass of the available observations, that is calibrated using the observed spectra of the flares. This simple model reproduces, with excellent agreement, the empirical statistical properties of the flares as measured by their width-to-arrival time ratio and skewness (ratio between decay and rise time). We present model fits to the observed lightcurves of two well-monitored flares, of GRB 060418 and of GRB 060904B. To the best of our knowledge, this is the first quantitative model able to reproduce the flare lightcurves and explain their global statistical properties.

We present the detection of persistent soft X-ray radiation with L_x ~ 10^41-10^42 erg/s at the location of the extremely luminous, double-humped transient ASASSN-15lh as revealed by Chandra and Swift. We interpret this finding in the context of observations from our multiwavelength campaign, which revealed the presence of weak narrow nebular emission features from the host-galaxy nucleus and clear differences with respect to superluminous supernova optical spectra. Significant UV flux variability on short time-scales detected at the time of the re-brightening disfavors the shock interaction scenario as the source of energy powering the long-lived UV emission, while deep radio limits exclude the presence of relativistic jets propagating into a low-density environment. We propose a model where the extreme luminosity and double-peaked temporal structure of ASASSN-15lh is powered by a central source of ionizing radiation that produces a sudden change of the ejecta opacity at later times. As a result, UV radiation can more easily escape, producing the second bump in the light-curve. We discuss different interpretations for the intrinsic nature of the ionizing source. We conclude that, IF the X-ray source is physically associated with the optical-UV transient, ASASSN-15lh most likely represents the tidal disruption of a main-sequence star by the most massive spinning black hole detected to date. In this case, ASASSN-15lh and similar events discovered in the future would constitute the most direct probes of very massive, dormant, spinning, supermassive black holes in galaxies. Future monitoring of the X-rays may allow us to distinguish between the supernova and TDE hypothesis.

We report the formation of intermediate-mass black holes (IMBHs) in suites of numerical $N$-body simulations of Population III remnant black holes (BHs) embedded in gas-rich protogalaxies at redshifts $z\gtrsim10$. We model the effects of gas drag on the BHs' orbits, and allow BHs to grow via gas accretion, including a mode of hyper-Eddington accretion in which photon trapping and rapid gas inflow suppress any negative radiative feedback. Most initial BH configurations lead to the formation of one (but never more than one) IMBH in the center of the protogalaxy, reaching a mass of $10^{3-5}\mathrm{M}_{\odot}$ through hyper-Eddington growth. Our results suggest a viable pathway to forming the earliest massive BHs in the centers of early galaxies. We also find that the nuclear IMBH typically captures a stellar-mass BH companion, making these systems observable in gravitational waves as extreme mass-ratio inspirals (EMRIs) with \textit{eLISA}.

Short Gamma-Ray Bursts (GRBs) are explosions of cosmic origin believed to be associated with the merger of two compact objects, either two neutron stars, or a neutron star and a black hole. The presence of at least one neutron star has long been thought to be an essential element of the model: its tidal disruption provides the needed baryonic material whose rapid accretion onto the post-merger black hole powers the burst. The recent tentative detection by the Fermi satellite of a short GRB in association with the gravitational wave signal GW150914 produced by the merger of two black holes has challenged this standard paradigm. Here we show that the evolution of two high-mass, low-metallicity stars with main sequence rotational speeds a few tens of percent of the critical speed eventually undergoing a weak supernova explosion {\em can} produce a short gamma-ray burst. The outer layers of the envelope of the last exploding star remain bound and circularize at large radii. With time, the disk cools and becomes neutral, suppressing the magneto-rotational instability, and hence the viscosity. The disk remains 'long-lived dead' until tidal torques and shocks during the pre-merger phase heat it up and re-ignite accretion, rapidly consuming the disk and powering the short gamma-ray burst.

We analyze archival Chandra HRC observations of the ultra luminous accreting pulsar M82-X2 (NuSTAR J095551+6940.8), and determine an upper limit of $< 1.7\times 10^{38}$~erg/s to its luminosity at an epoch at which it was undetected. Combined with other recent measurements, this confirms that the source X-ray emission has been highly variable during the last 15 years, ranging from a maximum of $10^{40}$ erg/s through intermediate values $\sim$ a few $\times 10^{39}$ erg/s, and down to a minimum that must be below the current detection threshold $\sim (2-3) \times 10^{38}$ erg/s . We interpret these results by means of a magnetically-threaded disk model: when at peak luminosity, the neutron star (NS) is close to spin equilibrium, its inner disk edge r_m ~ 10^8 cm is approximately half the corotation radius r_{co}, and radiation pressure dominates the disk out to r_{tr} ~ 10^9 cm. In the radiation pressure-dominated regime, r_m grows very slowly as the mass inflow rate drops: as a result, r_m < r_{co} remains valid until the mass accretion rate becomes ~ the Eddington accretion rate, allowing a wide range of accretion luminosities to the NS. Once the mass accretion rate is below Eddington, accretion onto the NS is inhibited because r_m > r_{co}, and the source luminosity is expected to drop by a large factor. We conclude that a magnetically threaded, radiation pressure-dominated disk, around a highly magnetized NS (B~10^{13} G) offers the best intepretation for all the currently observed properties of NuSTAR J095551+6940.8. This source offers an unprecedented opportunity to study the disk-magnetosphere interaction in a new regime of supercritical accretion, and across the transition between-radiation pressure and gas-pressure dominance inside the disk.

We study the timing and spectral properties of the low-magnetic field, transient magnetar SWIFT J1822.3-1606 as it approached quiescence. We coherently phase-connect the observations over a time-span of ~500 days since the discovery of SWIFT J1822.3-1606 following the Swift-BAT trigger on 2011 July 14, and carried out a detailed pulse phase spectroscopy along the outburst decay. We follow the spectral evolution of different pulse phase intervals and find a phase and energy-variable spectral feature, which we interpret as proton cyclotron resonant scattering of soft photon from currents circulating in a strong (10^14 G) small-scale component of the magnetic field near the neutron star surface, superimposed to the much weaker (~3 x 10^13 G) magnetic dipole field. We discuss also the implications of the pulse-resolved spectral analysis for the emission regions on the surface of the cooling magnetar.

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.

The first astrophysical objects shaped the cosmic environment by reionizing and heating the intergalactic medium (IGM). In particular, X-rays are very efficient at heating the IGM before it became completely ionized, an effect that can be measured through the 21 cm line of neutral hydrogen. High-mass X-ray binaries (HMXBs), known to be prolific X-ray sources in star-forming galaxies at lower redshifts, are prime candidates for driving the thermal evolution of the IGM at redshifts $z > 20$. Despite their importance, the formation efficiency of HMXBs from the first stellar populations is not well understood---as such, their collective X-ray emission and the subsequent imprint on the 21 cm signature are usually evaluated using free parameters. Using $N$-body simulations, we estimate the rate of HMXB formation via mutual gravitational interactions of nascent, small groups of the first stars (Pop III stars). We run two sets of calculations: one in which stars form in small groups of five in nearly Keplerian initial orbits, and another in which two such groups collide (expected from mergers of host protogalaxies). We find that HMXBs form at a rate of one per $~10^{4}~{\rm M}_{\odot}$ in newly born stars, and that they emit with a power of $\sim 10^{41} ~{\rm erg}~{\rm s} ^{-1}$ in the $2-10$ keV band per solar mass per year of star formation. This value is a factor $\sim 10^{2}$ larger than what is observed in star forming galaxies at lower redshifts; the X-ray production from early HMXBs would have been even more copious, if they also formed $in$ $situ$ or via migration in protostellar disks. Combining our results with earlier studies suggests that early HMXBs were highly effective at heating the IGM and leaving a strong 21 cm signature. We discuss broader implications of our results, such as the rate of long GRBs from PopIII stars and the direct collapse black hole formation.

The redshifted 21 cm transition line of hydrogen tracks the thermal evolution of the neutral intergalactic medium (IGM) at "cosmic dawn," during the emergence of the first luminous astrophysical objects (~100 Myr after the Big Bang) but before these objects ionized the IGM (~400-800 Myr after the Big Bang). Because X-rays, in particular, are likely to be the chief energy courier for heating the IGM, measurements of the 21 cm signature can be used to infer knowledge about the first astrophysical X-ray sources. Using analytic arguments and a numerical population synthesis algorithm, we argue that the progenitors of supermassive black holes (SMBHs) should be the dominant source of hard astrophysical X-rays---and thus the primary driver of IGM heating and the 21 cm signature---at redshifts $z < 20$, if (i) they grow readily from the remnants of Population III stars and (ii) produce X-rays in quantities comparable to what is observed from active galactic nuclei and high-mass X-ray binaries. We show that models satisfying these assumptions dominate over contributions to IGM heating from stellar populations, and cause the 21 cm brightness temperature to rise at $z > 20$. An absence of such a signature in the forthcoming observational data would imply that SMBH formation occurred later (e.g. via so-called direct collapse scenarios), that it was not a common occurrence in early galaxies and protogalaxies, or that it produced far fewer X-rays than empirical trends at lower redshifts, either due to intrinsic dimness (radiative inefficiency) or Compton-thick obscuration close to the source.

We present the first uniform treatment of long duration gamma-ray burst (GRB) host galaxy detections and upper limits over the redshift range 3<z<5, a key epoch for observational and theoretical efforts to understand the processes, environments, and consequences of early cosmic star formation. We contribute deep imaging observations of 13 GRB positions yielding the discovery of eight new host galaxies. We use this dataset in tandem with previously published observations of 31 further GRB positions to estimate or constrain the host galaxy rest-frame ultraviolet (UV; 1600 A) absolute magnitudes M_UV. We then use the combined set of 44 M_UV estimates and limits to construct the M_UV luminosity function (LF) for GRB host galaxies over 3<z<5 and compare it to expectations from Lyman break galaxy (LBG) photometric surveys with the Hubble Space Telescope. Adopting standard prescriptions for the luminosity dependence of galaxy dust obscuration (and hence, total star formation rate), we find that our LF is compatible with LBG observations over a factor of 600x in host luminosity, from M_UV = -22.5 mag to >-15.6 mag, and with extrapolations of the assumed Schechter-type LF well beyond this range. We review proposed astrophysical and observational biases for our sample, and find they are for the most part minimal. We therefore conclude, as the simplest interpretation of our results, that GRBs successfully trace UV metrics of cosmic star formation over the range 3<z<5. Our findings suggest GRBs are providing an accurate picture of star formation processes from z ~3 out to the highest redshifts.

The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale dynamo activity. In order to properly resolve such instabilities and obtain the correct magnetic field amplification, one would need to employ resolutions that are currently unfeasible in global general relativistic magnetohydrodynamic (GRMHD) simulations of BNS mergers. Here, we present a subgrid model that allows global simulations to take into account the small-scale amplification of the magnetic field which is caused by the development of turbulence during BNS mergers. Assuming dynamo saturation, we show that magnetar-level fields ($\sim 10^{16}\,{\rm G}$) can be easily reached, and should therefore be expected from the merger of magnetized BNSs. The total magnetic energy can reach values up to $\sim 10^{51}\,{\rm erg}$ and the post-merger remnant can therefore emit strong electromagnetic signals and possibly produce short gamma-ray bursts.

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.

We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6x10^12 G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. Nonetheless, the modelling allows us to place constraints on the system geometry (i.e. the angles $\psi$ and $\xi$ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte-Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to $65\deg < \psi+\xi < 125\deg$ or $235\deg < \psi+\xi <295\deg$; (b) the temperature contrast between the poles and the equator must be at least a factor of $\sim6$, and (c) the size of the hottest region ranges between 0.2-0.7 km (including uncertainties on the source distance). Last, we interpret our findings within the context of internal and external heating models.

In 2013 April a new magnetar, SGR 1745-2900, was discovered as it entered an outburst, at only 2.4 arcsec angular distance from the supermassive black hole at the Centre of the Milky Way, Sagittarius A*. SGR 1745-2900 has a surface dipolar magnetic field of ~ 2x10^{14} G, and it is the neutron star closest to a black hole ever observed. The new source was detected both in the radio and X-ray bands, with a peak X-ray luminosity L_X ~ 5x10^{35} erg s^{-1}. Here we report on the long-term Chandra (25 observations) and XMM-Newton (8 observations) X-ray monitoring campaign of SGR 1745-2900, from the onset of the outburst in April 2013 until September 2014. This unprecedented dataset allows us to refine the timing properties of the source, as well as to study the outburst spectral evolution as a function of time and rotational phase. Our timing analysis confirms the increase in the spin period derivative by a factor of ~2 around June 2013, and reveals that a further increase occurred between 2013 Oct 30 and 2014 Feb 21. We find that the period derivative changed from 6.6x10^{-12} s s^{-1} to 3.3x10^{-11} s s^{-1} in 1.5 yr. On the other hand, this magnetar shows a slow flux decay compared to other magnetars and a rather inefficient surface cooling. In particular, starquake-induced crustal cooling models alone have difficulty in explaining the high luminosity of the source for the first ~200 days of its outburst, and additional heating of the star surface from currents flowing in a twisted magnetic bundle is probably playing an important role in the outburst evolution.

Long-Duration Gamma-Ray Bursts (GRBs) are powerful probes of the Universe star formation history, but correlation between the two depends on the highly debated presence/strength of a metallicity bias. To investigate this correlation, we use a phenomenological model that successfully describes star formation rates, luminosities and stellar masses of star forming galaxies, applying it to GRB production. We predict luminosities, stellar masses, and metallicities of host galaxies depending on the metallicity bias. Our best-fitting model includes a moderate metallicity bias, broadly consistent with the large majority of long-duration GRBs in metal-poor environments originating from collapsars (probability ~83%), but with a secondary contribution (~17%) from metal-independent production channels, such as binary evolution. Because of the mass-metallicity relation of galaxies, the maximum likelihood model predicts that the metal-independent channel becomes dominant at z<2, where hosts have higher metallicities and collapsars are suppressed. This possibly explains why some studies find no clear evidence of a metal-bias based on low-z samples. However, while metallicity predictions match observations well at high redshift, there is tension with low redshift observations, since a significant fraction of GRB hosts are predicted to have (near-)solar metallicity. This is in contrast to observations, unless obscured, metal-rich hosts are preferentially missed in current datasets, and suggests that lower efficiencies of the metal-independent GRB channel might be preferred following a comprehensive fit from complete samples. Overall, we are able to establish the presence of a metallicity bias for GRB production, but continued characterization of GRB host galaxies is needed to quantify its strength.

The identification of the Ultraluminous X-ray source (ULX) X-2 in M82 as an accreting pulsar has shed new light on the nature of a subset of ULXs, while rising new questions on the nature of the super-Eddington accretion. Here, by numerically solving the torque equation of the accreting pulsar within the framework of the magnetically threaded-disk scenario, we show that three classes of solutions, corresponding to different values of the magnetic field, are mathematically allowed. We argue that the highest magnetic field one, corresponding to B $\sim 10^{13}$ G, is favoured based on physical considerations and the observed properties of the source. In particular, that is the only solution which can account for the observed variations in $\dot{P}$ (over four time intervals) without requiring major changes in $\dot{M}$, which would be at odds with the approximately constant X-ray emission of the source during the same time. For this solution, we find that the source can only accomodate a moderate amount of beaming, 0.5 $\lesssim b < 1$. Last, we show that the upper limit on the luminosity, L$_X < 2.5 \times 10^{38}$ erg s$^{-1}$ from archival observations, is consistent with a highly-magnetized neutron star being in the propeller phase at that time.

We present high resolution transmission spectra of giant planet atmospheres from a coupled 3-D atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9 to 55 day orbital periods around solar-type stars. The results of our 3-D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple 1-D models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blue shifts of up to 3 km s$^{-1}$, whereas less irradiated planets show almost no net Doppler shifts. Compared to 1-D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3-D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1-D atmospheric models may be inadequate, as 3-D atmospheric motions can produce a noticeable effect on the absorption signatures.

Binary neutron star (NS) mergers are among the most promising sources of gravitational waves (GWs), as well as candidate progenitors for short Gamma-Ray Bursts (SGRBs). Depending on the total initial mass of the system, and the NS equation of state (EOS), the post-merger phase can see a prompt collapse to a black hole, or the formation of a supramassive NS, or even a stable NS. In the case of post-merger NS (PMNS) formation, magnetic field amplification during the merger will produce a magnetar with a large induced mass quadrupole moment, and millisecond spin. If the timescale for orthogonalization of the magnetic symmetry axis with the spin axis is sufficiently short the NS will radiate its spin down energy primarily via GWs. Here we study this scenario for various outcomes of NS formation: we generalise the set of equilibrium states for a twisted torus magnetic configuration to include solutions that, at a fixed exterior dipole field, carry a larger magnetic energy reservoir; we hence compute their magnetic ellipticity and the strength of the expected GW signal as a function of the magnitude of the dipole and toroidal field. The relative number of GW detections from PMNSs and from binary NSs is a strong function of the NS equation of state (EOS), being higher (~ 1%) for the stiffest EOSs and negligibly small for the softest ones. For intermediate-stiffness EOSs, such as the n=4/7 polytrope recently used by Giacomazzo \& Perna or the GM1 used by Lasky et al., the relative fraction is ~0.3%; correspondingly we estimate a GW detection rate from stable PMNSs of ~ (0.1-1) yr$^{-1}$ with Advanced detectors, and of ~ (100-1000) yr$^{-1}$ with third generation detectors such as the Einstein Telescope. Measurement of such GW signal would provide strong constraints on the NS EOS and on the nature of the binary progenitors giving rise to SGRBs.

The thermal X-ray spectra of several isolated neutron stars display deviations from a pure blackbody. The accurate physical interpretation of these spectral features bears profound implications for our understanding of the atmospheric composition, magnetic field strength and topology, and equation of state of dense matter. With specific details varying from source to source, common explanations for the features have ranged from atomic transitions in the magnetized atmospheres or condensed surface, to cyclotron lines generated in a hot ionized layer near the surface. Here we quantitatively evaluate the X-ray spectral distortions induced by inhomogeneous temperature distributions of the neutron star surface. To this aim, we explore several surface temperature distributions, we simulate their corresponding general relativistic X-ray spectra (assuming an isotropic, blackbody emission), and fit the latter with a single blackbody model. We find that, in some cases, the presence of a spurious 'spectral line' is required at a high significance level in order to obtain statistically acceptable fits, with central energy and equivalent width similar to the values typically observed. We also perform a fit to a specific object, RX J0806.4-4123, finding several surface temperature distributions able to model the observed spectrum. The explored effect is unlikely to work in all sources with detected lines, but in some cases it can indeed be responsible for the appearance of such lines. Our results enforce the idea that surface temperature anisotropy can be an important factor that should be considered and explored also in combination with more sophisticated emission models like atmospheres.

The presence of fallback disks around young neutron stars has been invoked over the years to explain a large variety of phenomena. Here we perform a numerical investigation of the formation of such disks during a supernova explosion, considering both neutron star (NS) and black hole (BH) remnants. Using the public code MESA, we compute the angular momentum distribution of the pre-supernova material, for stars with initial masses M in the range 13 - 40 Msun, initial surface rotational velocities vsurf between 25% and 75% of the critical velocity, and for metallicities Z of 1%, 10% and 100% of the solar value. These pre SN models are exploded with energies E varying between 10^{50} - 3x10^{52} ergs, and the amount of fallback material is computed. We find that, if magnetic torques play an important role in angular momentum transport, then fallback disks around NSs, even for low-metallicity main sequence stars, are not an outcome of SN explosions. Formation of such disks around young NSs can only happen under the condition of negligible magnetic torques and a fine-tuned explosion energy. For those stars which leave behind BH remnants, disk formation is ubiquitous if magnetic fields do not play a strong role; however, unlike the NS case, even with strong magnetic coupling in the interior, a disk can form in a large region of the {Z,M,vsurf,E} parameter space. Together with the compact, hyperaccreting fallback disks widely discussed in the literature, we identify regions in the above parameter space which lead to extended, long-lived disks around BHs. We find that the physical conditions in these disks may be conducive to planet formation, hence leading to the possible existence of planets orbiting black holes.

The duration of the prompt emission of long gamma-ray bursts is generally considered to be fairly similar to the duration of the activity of the engine in the center of the progenitor star. Here, we investigate the relation between the duration of the engine activity and that of the observed light curve, using inputs from both numerical simulations and observations. We find that the observed burst duration is a good proxy for the engine duration after the time necessary for the jet to break out the star's surface is subtracted. However, the observed duration is a function of the viewing angle and can be significantly shorter than the duration of the engine activity. We also show that the observed, redshift-corrected burst duration evolves only moderately with redshift for both observations and synthetic light curves. We conclude that the broad distribution of the observed duration of long BATSE gamma-ray bursts is mostly accounted for by an engine lasting ~20 s, the dispersion being due to viewing and redshift effects. Our results do not rule out the existence of engines with very long duration. However, we find that they are constrained to be a small minority of the BATSE detected bursts.

The center of our Galaxy hosts a supermassive black hole, Sagittarius (Sgr) A*. Young, massive stars within 0.5 pc of SgrA* are evidence of an episode of intense star formation near the black hole a few Myr ago, which might have left behind a young neutron star traveling deep into SgrA*'s gravitational potential. On 2013 April 25, a short X-ray burst was observed from the direction of the Galactic center. Thanks to a series of observations with the Chandra and the Swift satellites, we pinpoint the associated magnetar at an angular distance of 2.4+/-0.3 arcsec from SgrA*, and refine the source spin period and its derivative (P=3.7635537(2) s and \dot{P} = 6.61(4)x10^{-12} s/s), confirmed by quasi simultaneous radio observations performed with the Green Bank (GBT) and Parkes antennas, which also constrain a Dispersion Measure of DM=1750+/-50 pc cm^{-3}, the highest ever observed for a radio pulsar. We have found that this X-ray source is a young magnetar at ~0.07-2 pc from SgrA*. Simulations of its possible motion around SgrA* show that it is likely (~90% probability) in a bound orbit around the black hole. The radiation front produced by the past activity from the magnetar passing through the molecular clouds surrounding the Galactic center region, might be responsible for a large fraction of the light echoes observed in the Fe fluorescence features.

Gamma Ray Bursts (GRBs) and galaxies at high redshift represent complementary probes of the star formation history of the Universe. In fact, both the GRB rate and the galaxy luminosity density are connected to the underlying star formation. Here, we combine a star formation model for the evolution of the galaxy luminosity function from z=0 to z=10 with a metallicity-dependent efficiency for GRB formation to simultaneously predict the comoving GRB rate. Our model sheds light on the physical origin of the empirical relation often assumed between GRB rate and luminosity density-derived star formation rate: Rgrb(z) = \epsilon(z)*SFR_{obs}(z), with \epsilon(z) (1+z)^{1.2}. At z<4, \epsilon(z) is dominated by the effects of metallicity evolution in the GRB efficiency. Our best-fitting model only requires a moderate preference for low-metallicity, that is a GRB rate per unit stellar mass about four times higher for log(Z/Zsun)<-3 compared to log(Z/Zsun)>0. Models with total suppression of GRB formation at log(Z/Zsun)>0 are disfavored. At z>4, most of the star formation happens in low-metallicity hosts with nearly saturated efficiency of GRB production per unit stellar mass. However at the same epoch, galaxy surveys miss an increasing fraction of the predicted luminosity density because of flux limits, driving an accelerated evolution of \epsilon(z) compared to the empirical power-law fit from lower z. Our findings are consistent with the non-detections of GRB hosts in ultradeep imaging at z>5, and point toward current galaxy surveys at z>8 only observing the top 15-20 % of the total luminosity density.

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50-60% level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated "observed" spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ~ 1-2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages, and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements.

By performing fully general relativistic magnetohydrodynamic simulations of binary neutron star mergers, we investigate the possibility that the end result of the merger is a stable magnetar. In particular, we show that, for a binary composed of two equal-mass neutron stars (NSs) of gravitational mass M~1.2 Msun and equation of state similar to Shen et al. at high densities, the merger product is a stable NS. Such NS is found to be differentially rotating and ultraspinning with spin parameter J/M^2~0.86, where J is its total angular momentum, and it is surrounded by a disk of ~0.1 Msun. While in our global simulations the magnetic field is amplified by about two orders of magnitude, local simulations have shown that hydrodynamic instabilities and the onset of the magnetorotational instability could further increase the magnetic field strength up to magnetar levels. This leads to the interesting possibility that, for some NS mergers, a stable and magnetized NS surrounded by an accretion disk could be formed. We discuss the impact of these new results for the emission of electromagnetic counterparts of gravitational wave signals and for the central engine of short gamma-ray bursts.

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.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

We report on the long term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729, using all the available X-ray data obtained with RXTE, SWIFT, Chandra, and XMM-Newton observations, from the discovery of the source in June 2009, up to August 2012. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: \dot{P}=4(1)x10^{-15} s/s, significant at ~3.5 sigma confidence level. This leads to a surface dipolar magnetic field of B_dip ~6x 10^{12} G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ~1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ~0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about 3 orders of magnitude: from about 1.4x10^{-11} to 1.2x10^{-14} erg/s/cm^2 . Deep radio, millimeter, optical and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ~550 kyr, and a dipolar magnetic field at birth of ~10^{14} G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such result in several other astrophysical contexts.

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^\alpha$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $\alpha$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$.

Central Compact Objects are peculiar young neutron stars, with very low external magnetic fields combined with high fluxes in the X-ray band and surface temperature anisotropies. However, in their crust the magnetic field can be strong, result of its burial during a short post-supernova hypercritical accretion episode. The implications of this latter scenario for the temperature anisotropy, pulsed fraction and luminosity are discussed.

In recent years, detailed observations and accurate numerical simulations have provided support to the idea that mergers of compact binaries containing either two neutron stars (NSs) or an NS and a black hole (BH) may constitute the central engine of short gamma-ray bursts (SGRBs). The merger of such compact binaries is expected to lead to the production of a spinning BH surrounded by an accreting torus. Several mechanisms can extract energy from this system and power the SGRBs. Here we connect observations and numerical simulations of compact binary mergers, and use the current sample of SGRBs with measured energies to constrain the mass of their powering tori. By comparing the masses of the tori with the results of fully general-relativistic simulations, we are able to infer the properties of the binary progenitors which yield SGRBs. By assuming a constant efficiency in converting torus mass into jet energy, epsilon_{jet}=10%, we find that most of the tori have masses smaller than 0.01M_{sun}, favoring "high-mass" binary NSs mergers, i.e., binaries with total masses >~1.5 the maximum mass of an isolated NS. This has important consequences for the gravitational-wave signals that may be detected in association with SGRBs, since "high-mass" systems do not form a long-lived hypermassive NS after the merger. While NS-BH systems cannot be excluded to be the engine of at least some of the SGRBs, the BH would need to have an initial spin of ~0.9, or higher.

We present a long-term study of the 2011 outburst of the magnetar Swift J1834.9-0846 carried out using new Chandra observations, as well as all the available Swift, RXTE, and XMM-Newton data. The last observation was performed on 2011 November 12, about 100 days after the onset of the bursting activity that had led to the discovery of the source on 2011 August 07. This long time span enabled us to refine the rotational ephemeris and observe a downturn in the decay of the X-ray flux. Assuming a broken power law for the long-term light curve, the break was at ~46 d after the outburst onset, when the decay index changed from alpha ~ 0.4 to ~4.5. The flux decreased by a factor ~2 in the first ~50 d and then by a factor ~40 until November 2011 (overall, by a factor ~70 in ~100 d). At the same time, the spectrum, which was well described by an absorbed blackbody all along the outburst, softened, the temperature dropping from ~1 to ~0.6 keV. Diffuse X-ray emission extending up to 20" from the source was clearly detected in all Chandra observations. Its spatial and spectral properties, as well as its time evolution, are consistent with a dust-scattering halo due to a single cloud located at a distance of $\approx$200 pc from Swift J1834.9-0846, which should be in turn located at a distance of ~5 kpc. Considering the time delay of the scattered photons, the same dust cloud might also be responsible for the more extended emission detected in XMM-Newton data taken in September 2011. We searched for the radio signature of Swift J1834.9-0846 at radio frequencies using the Green Bank Radio Telescope and in archival data collected at Parkes from 1998 to 2003. No evidence for radio emission was found, down to a flux density of 0.05 mJy (at 2 GHz) during the outburst and ~0.2-0.3 mJy (at 1.4 GHz) in the older data.

We present a general formalism for computing self-consistent, numerical solutions to the time-dependent radiative transfer equation in low velocity, multi-level ions undergoing radiative interactions. Recent studies of time-dependent radiative transfer have focused on radiation hydrodynamic and magnetohydrodynamic effects without lines, or have solved time-independent equations for the radiation field simultaneously with time-dependent equations for the state of the medium. In this paper, we provide a fully time-dependent numerical solution to the radiative transfer and atomic rate equations for a medium irradiated by an external source of photons. We use Accelerated Lambda Iteration to achieve convergence of the radiation field and atomic states. We perform calculations for a three-level atomic model that illustrates important time-dependent effects. We demonstrate that our method provides an efficient, accurate solution to the time-dependent radiative transfer problem. Finally, we characterize astrophysical scenarios in which we expect our solutions to be important.

Neutron stars (NSs) in the astrophysical universe are often surrounded by accretion disks. Accretion of matter onto an NS may increase its mass above the maximum value allowed by its equation of state, inducing its collapse to a black hole (BH). Here we study this process for the first time, in three-dimensions, and in full general relativity. By considering three initial NS configurations, each with and without a surrounding disk (of mass ~7% M_{NS}), we investigate the effect of the accretion disk on the dynamics of the collapse and its imprint on both the gravitational wave (GW) and electromagnetic (EM) signals that can be emitted by these sources. We show in particular that, even if the GW signal is similar for the accretion induced collapse (AIC) and the collapse of an NS in vacuum (and detectable only for Galactic sources), the EM counterpart could allow us to discriminate between these two types of events. In fact, our simulations show that, while the collapse of an NS in vacuum leaves no appreciable baryonic matter outside the event horizon, an AIC is followed by a phase of rapid accretion of the surviving disk onto the newly formed BH. The post-collapse accretion rates, on the order of ~10^{-2} M_{sun} s^{-1}, make these events tantalizing candidates as engines of short gamma-ray bursts.

Observations of high-redshift quasars at z>6 imply that supermassive black holes (SMBHs) with masses over a billion solar masses were in place less than 1 Gyr after the Big Bang. If these SMBHs assembled from "seed" BHs left behind by the first stars, then they must have accreted gas at close to the Eddington limit during a large fraction (>50%) of the time. A generic problem with this scenario, however, is that the mass density in million-solar-mass SMBHs at z=6 already exceeds the locally observed SMBH mass density by several orders of magnitude; in order to avoid this overproduction, BH seed formation and growth must become significantly less efficient in less massive protogalaxies, while proceeding uninterrupted in the most massive galaxies that formed first. Using Monte-Carlo realizations of the merger and growth history of BHs, we show that X-rays from the earliest accreting BHs can provide such a feedback mechanism. Our calculations paint a self-consistent picture of black-hole-made climate change, in which the first miniquasars---among them the ancestors of the z>6 quasar SMBHs---globally warm the IGM and suppress the formation and growth of subsequent generations of BHs. We present two specific models with global miniquasar feedback that provide excellent agreement with recent estimates of the z=6 SMBH mass function. For each of these models, we estimate the rate of BH mergers at z>6 that could be detected by the proposed gravitational-wave observatory eLISA/NGO.

We report on the long term X-ray monitoring with Swift, RXTE, Suzaku, Chandra and XMM-Newton of the outburst of the newly discovered magnetar Swift J1822.3-1606 (SGR 1822-1606), from the first observations soon after the detection of the short X-ray bursts which led to its discovery, through the first stages of its outburst decay (covering the time-span from July 2011, until end of April 2012). We also report on archival ROSAT observations which witnessed the source during its likely quiescent state, and on upper limits on Swift J1822.3-1606's radio-pulsed and optical emission during outburst, with the Green Bank Telescope (GBT) and the Gran Telescopio Canarias (GTC), respectively. Our X-ray timing analysis finds the source rotating with a period of P=8.43772016(2) s and a period derivative \dot{P}=8.3(2)x10^{-14} s s^{-1} , which entails an inferred dipolar surface magnetic field of B~2.7x10^{13} G at the equator. This measurement makes Swift J1822.3-1606 the second lowest magnetic field magnetar (after SGR 0418+5729; Rea et al. 2010). Following the flux and spectral evolution from the beginning of the outburst, we find that the flux decreased by about an order of magnitude, with a subtle softening of the spectrum, both typical of the outburst decay of magnetars. By modeling the secular thermal evolution of Swift J1822.3-1606, we find that the observed timing properties of the source, as well as its quiescent X-ray luminosity, can be reproduced if it was born with a poloidal and crustal toroidal fields of B_{p}~1.5x10^{14} G and B_{tor}~7x10^{14} G, respectively, and if its current age is ~550 kyr.

Fallback discs around neutron stars (NSs) are believed to be an expected outcome of supernova explosions. Here we investigate the consequences of such a common outcome for the timing and spectral properties of the associated NS population, using Monte Carlo population synthesis models. We find that the long-term torque exerted by the fallback disc can substantially influence the late-time period distribution, but with quantitative differences which depend on whether the initial spin distribution is dominated by slow or fast pulsars. For the latter, a single-peaked initial spin distribution becomes bimodal at later times. Timing ages tend to underestimate the real age of older pulsars, and overestimate the age of younger ones. Braking indices cluster in the range 1.5 <~ n <~ 3 for slow-born pulsars, and -0.5 <~ n <~ 5 for fast-born pulsars, with the younger objects found predominantly below n <~ 3. Large values of n, while not common, are possible, and associated with torque transitions in the NS+disc system. The 0.1-10 keV thermal luminosity of the NS+disc system is found to be generally dominated by the disc emission at early times, t <~ 10^3 yr, but this declines faster than the thermal surface emission of the NS. Depending on the initial parameters, there can be occasional periods in which some NSs switch from the propeller to the accretion phase, increasing their luminosity up to the Eddington limit for ~ 10^3-10^4 years.

Hot Jupiters, due to the proximity to their parent stars, are subjected to a strong irradiating flux which governs their radiative and dynamical properties. We compute a suite of 3D circulation models with dual-band radiative transfer, exploring a relevant range of irradiation temperatures, both with and without temperature inversions. We find that, for irradiation temperatures T \lesssim 2000 K, heat redistribution is very efficient, producing comparable day- and night-side fluxes. For Tirr \approx 2200-2400 K, the redistribution starts to break down, resulting in a high day-night flux contrast. Our simulations indicate that the efficiency of redistribution is primarily governed by the ratio of advective to radiative timescales. Models with temperature inversions display a higher day-night contrast due to the deposition of starlight at higher altitudes, but we find this opacity-driven effect to be secondary compared to the effects of irradiation. The hotspot offset from the substellar point is large when insolation is weak and redistribution is efficient, and decreases as redistribution breaks down. The atmospheric flow can be potentially subjected to the Kelvin-Helmholtz instability (as indicated by the Richardson number) only in the uppermost layers, with a depth that penetrates down to pressures of a few millibars at most. Shocks penetrate deeper, down to several bars in the hottest model. Ohmic dissipation generally occurs down to deeper levels than shock dissipation (to tens of bars), but the penetration depth varies with the atmospheric opacity. The total dissipated Ohmic power increases steeply with the strength of the irradiating flux and the dissipation depth recedes into the atmosphere, favoring radius inflation in the most irradiated objects.

Measuring the star formation rate (SFR) at high redshift is crucial for understanding cosmic reionization and galaxy formation. Two common complementary approaches are Lyman-Break-Galaxy (LBG) surveys for large samples and Gamma-Ray-Burst (GRB) observations for sensitivity to SFR in small galaxies. The z>4 GRB-inferred SFR is higher than the LBG rate, but this difference is difficult to understand, as both methods rely on several modeling assumptions. Using a physically motivated galaxy luminosity function model, with star formation in dark-matter halos with virial temperature Tvir>2e4 K (M_DM>2e8 M_sun), we show that GRB and LBG-derived SFRs are consistent if GRBs extend to faint galaxies (M_AB<-11). To test star formation below the detection limit L_lim~0.05L^*_{z=3} of LBG surveys, we propose to measure the fraction f_det(L>L_lim,z) of GRB hosts with L>L_lim. This fraction quantifies the missing star formation fraction in LBG surveys, constraining the mass-suppression scale for galaxy formation, with weak dependence on modeling assumptions. Because f_det(L>L_lim,z) corresponds to the ratio of star formation rates derived from LBG and GRB surveys, if these estimators are unbiased, measuring f_det(L>L_lim,z) also constrains the redshift evolution of the GRB production rate per unit mass of star formation. Our analysis predicts significant success for GRB host detections at z~5 with f_det(L>L_lim,z)~0.4, but rarer detections at z>6. By analyzing the upper limits on host-galaxy luminosities of six z>5 GRBs from literature data, we infer that galaxies with M_AB>-15 were present at z>5 at 95% confidence, demonstrating the key role played by very faint galaxies during reionization.

Conversion of photons into axions under the presence of a strong magnetic field can dim the radiation from magnetized astrophysical objects. Here we perform a detailed calculation aimed at quantifying the signatures of photon-axion conversion in the spectra, light curves, and polarization of neutron stars (NSs). We take into account the energy and angle-dependence of the conversion probability and the surface thermal emission from NSs. The latter is computed from magnetized atmosphere models that include the effect of photon polarization mode conversion due to vacuum polarization. The resulting spectral models, inclusive of the general-relativistic effects of gravitational redshift and light deflection, allow us to make realistic predictions for the effects of photon to axion conversion on observed NS spectra, light curves, and polarization signals. We identify unique signatures of the conversion, such as an increase of the effective area of a hot spot as it rotates away from the observer line of sight. For a star emitting from the entire surface, the conversion produces apparent radii that are either larger or smaller (depending on axion mass and coupling strength) than the limits set by NS equations of state. For an emission region that is observed phase-on, photon-axion conversion results in an inversion of the plane of polarization with respect to the no-conversion case. While the quantitative details of the features that we identify depend on NS properties (magnetic field strength, temperature) and axion parameters, the spectral and polarization signatures induced by photon-axion conversion are distinctive enough to make NSs very interesting and promising probes of axion physics.

Highly magnetized neutron stars (NSs) are characterized by a bewildering range of astrophysical manifestations. Here, building on our simulations of the evolution of magnetic stresses in the NS crust and its ensuing fractures (Perna & Pons 2011), we explore in detail, for the middle-age and old NSs, the dependence of starquake frequency and energetics on the relative strength of the poloidal (B_p) and toroidal (B_tor) components. We find that, for B_p >~10^{14}G, since a strong crustal toroidal field B_tor B_p is quickly formed on a Hall timescale, the initial toroidal field needs to be B_tor >> B_p to have a clear influence on the outbursting behaviour. For initial fields B_p <~ 10^{14}G, it is very unlikely that a middle-age (t~10^5 years) NS shows any bursting activity. This study allows us to solve the apparent puzzle of how NSs with similar dipolar magnetic fields can behave in a remarkably different way: an outbursting 'magnetar' with a high X-ray luminosity, or a quiet, low-luminosity, "high-$B$" radio pulsar. As an example, we consider the specific cases of the magnetar 1E2259+586 and the radio pulsar PSRJ1814-1744, which at present have a similar dipolar field ~6x10^{13}G. We determine for each object an initial magnetic field configuration that reproduces the observed timing parameters at their current age. The same two configurations also account for the differences in quiescent X-ray luminosity and for the 'magnetar/outbursting' behaviour of 1E2259+586 but not of PSRJ1814-1744. We further use the theoretically predicted surface temperature distribution to compute the light-curve for these objects. In the case of 1E2259+586, for which data are available, our predicted temperature distribution gives rise to a pulse profile whose double-peaked nature and modulation level is consistent with the observations.

The return to the quiescent state of the Anomalous X-ray pulsar XTE J1810-197 following its 2003 outburst represents a unique opportunity to probe the surface emission properties of a magnetar. The quiescent emission of XTE J1810-197 is composed of two thermal components, one arising from the whole star surface, and the other from a small warm spot on it. By modeling the magnitude and shape of the pulse profile in narrow spectral bands, we have been able to constrain the physical characteristics and geometrical parameters of the system: the two angles that the line of sight and the spin axis make with respect to the warm spot axis (\psi\ and \xi\ respectively), the angular size of the spot, and the overall surface temperature distribution. Our modeling accounts for the general relativistic effects of gravitational redshift and light bending near the stellar surface, and allows for local anisotropic emission. We found that the surface temperature distribution on the neutron star is consistent with the expectations of a dipole magnetic field configuration; the local radiation requires a pencil-beamed emission pattern, suggesting the presence of a magnetized atmosphere. For a typical value of the radius, R=13 km, the viewing parameters (symmetric for an interchange between \psi\ and \xi), range from \psi=\xi=38 deg to (\psi,\xi)=(52 deg, 29 deg). These angles are consistent with those obtained by modeling the AXP in outburst, with uncertainty contours reduced by a factor of 2.5.

Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) are young neutron stars (NSs) characterized by high X-ray quiescent luminosities, outbursts, and, in the case of SGRs, sporadic giant flares. They are believed to be powered by ultra-strong magnetic fields (hence dubbed magnetars). The diversity of their observed behaviours is however not understood, and made even more puzzling by the discovery of magnetar-like bursts from "low-field" pulsars. Here we perform long-term 2D simulations that follow the evolution of magnetic stresses in the crust; these, together with recent calculations of the breaking stress of the neutron star crust, allow us to establish when starquakes occur. For the first time, we provide a quantitative estimate of the burst energetics, event rate, and location on the neutron star surface, which bear a direct relevance for the interpretation of the overall magnetar phenomenology. Typically, an "SGR-like" object tends to be more active than an "AXP-like" object or a "high-$B$ radio pulsar", but there is no fundamental separation among what constitutes the apparent different classes. Among the key elements that create the variety of observed phenomena, age is more important than a small variation in magnetic field strength. We find that outbursts can also be produced in old, lower-field pulsars (B ~ a few x 10^{12} G), but those events are much less frequent than in young, high-field magnetars.

Both types of long and short gamma ray bursts involve a stage of a hyper-Eddington accretion of hot and dense plasma torus onto a newly born black hole. The prompt gamma ray emission originates in jets at some distance from this 'central engine' and in most events is rapidly variable, having a form of spikes and subpulses. This indicates at the variable nature of the engine itself, for which a plausible mechanism is an internal instability in the accreting flow. We solve numerically the structure and evolution of the neutrino-cooled torus. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. The models are calculated for both Schwarzschild and Kerr black holes. We find that for sufficiently large accretion rates (> 10 Msun/s for non-rotating black hole, and >1 Msun/s for rotating black hole, depending on its spin), the inner regions of the disk become opaque, while the helium nuclei are being photodissociated. The sudden change of pressure in this region leads to the development of a viscous and thermal instability, and the neutrino pressure acts similarly to the radiation pressure in sub-Eddington disks. In the case of rapidly rotating black holes, the instability is enhanced and appears for much lower accretion rates. We also find the important and possibly further destabilizing role of the energy transfer from the rotating black hole to the torus via the magnetic coupling.

We show that a pair of thermal, antipodal hot-spots on the neutron star surface is able to fully account for the pulsar's double blackbody spectrum and energy-dependent pulse profile, including the observed 180 degree phase reversal at approximately 1.2 keV. By comparing the observed pulse modulation and phase to the model predictions, we strongly constrain the hot-spot pole (xi) and the line-of-sight (psi) angles with respect to the spin axis. For a nominal radius of R = 12 km and distance D = 2.2 kpc, we find (xi,psi) = (86d,6d), with 1-sigma error ellipse of (2d,1d); this solution is degenerate in the two angles. The best-fit spectral model for this geometry requires that the temperatures of the two emission spots differ by a factor of 2 and their areas by a factor of ~ 20. Including a cosine-beamed pattern for the emitted intensity modifies the result, decreasing the angles to (84d,3d); however this model is not statistically distinguishable from the isotropic emission case. We also present a new upper limit on the period derivative of Pdot < 3.5E-16 (2-sigma), which limits the global dipole magnetic field to B_s < 2.0E11 G, confirming PSR J0821-4300 as an "anti-magnetar." We discuss the results in the context of observations and theories of nonuniform surface temperature on isolated NSs of both weak and strong magnetic field. To explain the nonuniform temperature of PSR J0821-4300 may require a crustal field that is much stronger than the external, global dipole field.

Hot Jupiter atmospheres exhibit fast, weakly-ionized winds. The interaction of these winds with the planetary magnetic field generates drag on the winds and leads to ohmic dissipation of the induced electric currents. We study the magnitude of ohmic dissipation in representative, three-dimensional atmospheric circulation models of the hot Jupiter HD 209458b. We find that ohmic dissipation can reach or exceed 1% of the stellar insolation power in the deepest atmospheric layers, in models with and without dragged winds. Such power, dissipated in the deep atmosphere, appears sufficient to slow down planetary contraction and explain the typically inflated radii of hot Jupiters. This atmospheric scenario does not require a top insulating layer or radial currents that penetrate deep in the planetary interior. Circulation in the deepest atmospheric layers may actually be driven by spatially non-uniform ohmic dissipation. A consistent treatment of magnetic drag and ohmic dissipation is required to further elucidate the consequences of magnetic effects for the atmospheres and the contracting interiors of hot Jupiters.

We report on the detailed study of the 2008 October outburst from the anomalous X-ray pulsar (AXP) 1E 1547.0-5408 discovered through the Swift/Burst Alert Telescope (BAT) detection of SGR-like short X-ray bursts on 2008 October 3. The Swift/X-ray Telescope (XRT) started observing the source after less than 100 s since the BAT trigger, when the flux (about 6E-11 erg/cm^2/s in the 2-10 keV range) was >50 times higher than its quiescent level. Swift monitored the outbursting activity of 1E 1547.0-5408 on a daily basis for approximately three weeks. This strategy allowed us to find a phase-coherent solution for the source pulsations after the burst, which, besides period and period derivative, requires a positive Period second derivative term (spin-down increase). The time evolution of the pulse shape is complex and variable, with the pulsed fraction increasing from 20% to 50% within the Swift observational window. The XRT spectra can be fitted well by means of a single component, either a power-law (PL) or a blackbody (BB). During the very initial phases of the outburst the spectrum is hard, with a PL photon index about 2 (or kT about 1.4 keV) which steepens to about 4 (or kT about 0.8 keV) within one day from the BAT trigger, though the two components are likely present simultaneously during the first day spectra. An INTEGRAL observation carried out five days after the trigger provided an upper limit of about 2E-11 erg/cm^2/s to the emission of 1E 1547.0-5408 in the 18-60 keV band.

We report on submillimetre bolometer observations of the isolated neutron star RX J1856.5--3754 using the LABOCA bolometer array on the Atacama Pathfinder Experiment (APEX) Telescope. No cold dust continuum emission peak at the position of RX J1856.5--3754 was detected. The 3 sigma flux density upper limit of 5 mJy translates into a cold dust mass limit of a few earth masses. We use the new submillimetre limit, together with a previously obtained H-band limit, to constrain the presence of a gaseous, circumpulsar disc. Adopting a simple irradiated-disc model, we obtain a mass accretion limit of dM/dt less than 10^{14} g/s, and a maximum outer disc radius of around 10^{14} cm. By examining the projected proper motion of RX J1856.5--3754, we speculate about a possible encounter of the neutron star with a dense fragment of the CrA molecular cloud a few thousand years ago.

The properties of a massive star prior to its final explosion are imprinted in the circumstellar medium (CSM) created by its wind and termination shock. We perform a detailed, comprehensive calculation of the time-variable and angle-dependent transmission spectra of an average-luminosity Gamma-Ray Burst (GRB) which explodes in the CSM structure produced by the collapse of a 20 Msun, rapidly rotating, Z=0.001 progenitor star. We study both the case in which metals are initially in the gaseous phase, as well as the situation in which they are heavily depleted into dust. We find that high-velocity lines from low-ionization states of silicon, carbon, and iron are initially present in the spectrum only if the metals are heavily depleted into dust prior to the GRB explosion. However, such lines disappear on timescales of a fraction of a second for a burst observed on-axis, and of a few seconds for a burst seen at high-latitude, making their observation virtually impossible. Rest-frame lines produced in the termination shock are instead clearly visible in all conditions. We conclude that time-resolved, early-time spectroscopy is not a promising way in which the properties of the GRB progenitor wind can be routinely studied. Previous detections of high velocity features in GRB UV spectra must have been due either due to a superposition of a physically unrelated absorber or to a progenitor star with very unusual properties.

Hot Jupiters, with atmospheric temperatures T ~ 1000 K, have residual thermal ionization levels sufficient for the interaction of the ions with the planetary magnetic field to result in a sizable magnetic drag on the (neutral) atmospheric winds. We evaluate the magnitude of magnetic drag in a representative three-dimensional atmospheric model of the hot Jupiter HD 209458b and find that it is a plausible mechanism to limit wind speeds in this class of atmospheres. Magnetic drag has a strong geometrical dependence, both meridionally and from the day to the night side (in the upper atmosphere), which could have interesting consequences for the atmospheric flow pattern. By extension, close-in eccentric planets with transiently heated atmospheres will experience time-variable levels of magnetic drag. A robust treatment of magnetic drag in circulation models for hot atmospheres may require iterated solutions to the magnetic induction and Saha equations as the hydrodynamical flow is evolved.

While there is mounting evidence that long Gamma-Ray Bursts (GRBs) are associated with the collapse of massive stars, the detailed structure of their pre-supernova stage is still debatable. Particularly uncertain is the degree of mixing among shells of different composition, and hence the role of magnetic torques and convection in transporting angular momentum. Here we show that early-time afterglow observations with the Swift satellite place constraints on the allowed GRB pre-supernova models. In particular, they argue against pre-supernova models in which different elemental shells are unmixed. These types of models would produce energy injections into the GRB engine on timescales between several hundreds of seconds to a few hours. Flaring activity has {\em not} been observed in a large fraction of well-monitored long GRBs. Therefore, if the progenitors of long GRBs have common properties, then the lack of flares indicates that the massive stars which produce GRBs are mostly well mixed, as expected in low-metallicity, rapidly rotating massive stars.

This contribution summarizes, as of early 2008, the observational and theoretical understanding of the origin, physics, and emission properties of short gamma-ray bursts in both electromagnetic and gravitational waves.

Aims. With an observed and rest-frame duration of < 2s and < 0.5s, respectively, GRB090426 could be classified as a short GRB. The prompt detection, both from space and ground-based telescopes, of a bright optical counterpart to this GRB offered a unique opportunity to complete a detailed study. Methods. Based on an extensive ground-based observational campaign, we obtained the spectrum of the optical afterglow of GRB090426, measuring its redshift and obtaining information about the medium in which the event took place. We completed follow-up observation of the afterglow optical light curve down to the brightness level of the host galaxy that we firmly identified and studied. We also retrieved and analyzed all the available high-energy data of this event, and compared the results with our findings in the optical. This represents one of the most detailed studies of a short-duration event presented so far. Results. The time properties qualify GRB090426 as a short burst. In this case, its redshift of z = 2.61 would be the highest yet found for a GRB of this class. On the other hand, the spectral and energy properties are more similar to those of long bursts. LBT late-time deep imaging identifies a star-forming galaxy at a redshift consistent with that of the GRB. The afterglow lies within the light of its host and shows evidence of local absorption.

Emission spectra from magnetars in the soft X-ray band likely contain a thermal component emerging directly from the neutron star surface. However, the lack of observed absorption-like features in quiescent spectra makes it difficult to directly constrain physical properties of the atmosphere. We argue that future X-ray polarization measurements represent a promising technique for directly constraining the magnetar magnetic field strength and geometry. We construct models of the observed polarization signal from a finite surface hotspot, using the latest NS atmosphere models for magnetic fields B = 4 x 10^13--5 x 10^14 G. Our calculations are strongly dependent on the NS magnetic field strength and geometry, and are more weakly dependent on the NS equation of state and atmosphere composition. We discuss how the complementary dependencies of phase-resolved spectroscopy and polarimetry might resolve degeneracies that currently hamper the determination of magnetar physical parameters using thermal models.

We study the Gamma Ray Burst (GRB) environment and intervening absorbers by analyzing the optical absorption features produced by gas surrounding the GRB or along its line of sight. We analyzed high resolution spectroscopic observations (R=40000, S/N=3 - 6) of the optical afterglow of GRB080330, taken with UVES at the VLT ~ 1.5 hours after the GRB trigger. The spectrum illustrates the complexity of the ISM of the GRB host galaxy at z = 1.51 which has at least four components in the main absorption system. We detect strong FeII, SiII, and NiII excited absorption lines associated with the bluemost component only. In addition to the host galaxy, at least two more absorbers lying along the line of sight to the afterglow have been detected in the redshift range 0.8 < z < 1.1, each exhibiting MgII absorption. For the bluemost component in the host galaxy, we derive information about its distance from the site of the GRB explosion. We do so by assuming that the excited absorption lines are produced by indirect UV pumping, and compare the data with a time dependent photo-excitation code. The distance of this component is found to be 280+40-50 pc, which is lower than found for other GRBs (1 - 6 kpc). We identify two additional MgII absorbers, one of them with a rest frame equivalent width larger than 1A. The distance between the GRB and the absorber measured in this paper confirms that the power of the GRB radiation can influence the conditions of the interstellar medium up to a distance of at least several hundred pc. For the intervening absorbers, we confirm the trend that on average one strong intervening system is found per afterglow, as has been noted in studies exhibiting an excess of strong MgII absorbers along GRB sightlines compared to quasars.

In strong gravitational lensing, the multiple images we see correspond to light rays that leave the source in slightly different directions. If the source emission is anisotropic, the images may differ from conventional lensing predictions (which assume isotropy). To identify scales on which source anisotropy may be important, we study the angle delta between the light rays emerging from the source, for different lensing configurations. If the lens has a power law profile M propto R^gamma, the angle delta initially increases with lens redshift and then either diverges (for a steep profile gamma<1), remains constant (for an isothermal profile gamma=1), or vanishes (for a shallow profile gamma>1) as zl approaches zs. The scaling with lens mass is roughly delta propto M^(1/(2-gamma)). The results for an NFW profile are qualitatively similar to those for a shallow power law, with delta peaking at about half the redshift of the source (not half the distance). In practice, beaming could modify the statistics of beamed sources lensed by massive clusters: for an opening angle theta_jet, there is a probability as high as P ~ 0.02-0.07 (theta_jet/{0.5} deg)^-1 that one of the lensed images may be missed (for 2 < zs < 6). Differential absorption within Active Galactic Nuclei could modify the flux ratios of AGNs lensed by clusters; a sample of AGNs lensed by clusters could provide further constraints on the sizes of absorbing regions. Source anisotropy is not likely to be a significant effect in galaxy-scale strong lensing.

We report here on the outburst onset and evolution of the new Soft Gamma Repeater SGR 0501+4516. We monitored the new SGR with XMM-Newton starting on 2008 August 23, one day after the source became burst-active, and continuing with 4 more observations, with the last one on 2008 September 30. Combining the data with the Swift-XRT and Suzaku data, we modelled the outburst decay over 160 days, and we found that the source flux decreased exponentially with a timescale of t_c=23.8 days. In the first XMM-Newton observation a large number of short X-ray bursts were observed, the rate of which decayed drastically in the following observations. We found large changes in the spectral and timing behavior of the source during the outburst, with softening emission as the flux decayed, and the non-thermal soft X-ray spectral component fading faster than the thermal one. Almost simultaneously to our XMM-Newton observations (on 2008 August 29 and September 2), we observed the source in the hard X-ray range with INTEGRAL, which clearly detected the source up to ~100keV in the first pointing, while giving only upper limits during the second pointing, discovering a variable hard X-ray component fading in less than 10 days after the bursting activation. We performed a phase-coherent X-ray timing analysis over about 160 days starting with the burst activation and found evidence of a strong second derivative period component (\ddot{P} = -1.6(4)x10^{-19} s/s^{-2}). Thanks to the phase-connection, we were able to study the the phase-resolved spectral evolution of SGR 0501+4516 in great detail. We also report on the ROSAT quiescent source data, taken back in 1992 when the source exhibits a flux ~80 times lower than that measured during the outburst, and a rather soft, thermal spectrum.

Short GRBs are commonly thought to originate from the merging of double compact object binaries but direct evidence for this scenario is still missing. Optical observations of short GRBs allow us to measure redshifts, firmly identify host galaxies, characterize their properties, and accurately localize GRBs within them. Multiwavelength observations of GRB afterglows provide useful information on the emission mechanisms at work. These are all key issues that allow one to discriminate among different models of these elusive events. We carried out photometric observations of the short/hard GRB 051227, GRB 061006, and GRB 071227 with the ESO-VLT starting from several hours after the explosion down to the host galaxy level several days later. For GRB 061006 and GRB 071227 we also obtained spectroscopic observations of the host galaxy. We compared the results obtained from our optical observations with the available X-ray data of these bursts. For all the three above bursts, we discovered optical afterglows and firmly identified their host galaxies. About half a day after the burst, the optical afterglows of GRB 051227 and GRB 061006 present a decay significatly steeper than in the X-rays. In the case of GRB 051227, the optical decay is so steep that it likely indicates different emission mechanisms in the two wavelengths ranges. The three hosts are blue, star forming galaxies at moderate redshifts and with metallicities comparable to the Solar one. The projected offsets of the optical afterglows from their host galaxies centers span a wide range, but all afterglows lie within the light of their hosts and present evidence for local absorption in their X-ray spectra. We discuss our findings in light of the current models of short GRB progenitors.

Despite about a decade of observations, very little is known about the optical and infrared (IR) emission properties of the Soft Gamma-ray Repeaters (SGRs) and of the Anomalous X-ray Pulsars (AXPs), the magnetar candidates, and about the physical processes which drive their emission at these wavelengths. This is mainly due to the limited number of identifications achieved so far, five in total, and to the sparse spectral coverage obtained from multi-band optical/IR photometry. Aim of this work is to search for a likely candidate counterpart to the recently discovered transient radio AXP 1E 1547.0-5408. We performed the first deep near-IR (NIR) observations (Ks band) of 1E 1547.0-5408 with the VLT on three nights (July 8th, 12th, and August 17th), after the X-ray source rebrightening and during the subsequent decay reported around June 2007. We detected four objects within, or close to, the 3sigma radio position of 1E 1547.0-5408. The faintest of them (object 1) has a brightness K = 20.27 +/- 0.05, which would yield an unabsorbed X-ray to NIR flux ratio FX/FKs ~ 800 for 1E 1547.0-5408, i.e. on average lower than those derived for other magnetars. The non-detection of object 1 on the nights of July 8th and August 17th only allowed us to set an upper limit of Delta Ks ~ 0.2 on its NIR variability, which prevented us to search for correlations with the radio or X-ray flux. We detected no other object at the radio position down to a limit of Ks ~ 21.7 (at 5sigma), computed in our deepest VLT image (July 12th). From our observations we can not confidently propose a NIR counterpart to 1E 1547.0-5408. More NIR observations of object 1, e.g. to determine its colors and to monitor variability, would be conclusive to determine whether or not it can be considered a plausible candidate.

XTE J1810-197 is the first transient Anomalous X-ray Pulsar ever discovered. Its highly variable X-ray flux allowed us to study the timing and spectral emission properties of a magnetar candidate over a flux range of about two orders of magnitude. We analyzed nine XMM-Newton observations of XTE J1810-197 collected over a four years baseline (September 2003 - September 2007). EPIC PN and MOS data were reduced and used for detailed timing and spectral analysis. Pulse phase spectroscopic studies were also carried out for observations with sufficiently high signal to noise. We find that: (i) a three blackbodies model reproduces the spectral properties of the source over the entire outburst statistically better than the two blackbodies model previously used in the literature, (ii) the coldest blackbody is consistent with the thermal emission from the whole surface, and has temperature and radius similar to those inferred from ROSAT observations before the outburst onset, (iii) there is a spectral feature around 1.1 keV during six consecutive observations (since March 2005); if due to proton resonant cyclotron scattering, it would imply a magnetic field of around 2E14 G. This is in a very good agreement with the value from the spin period measurements.

We have studied the X-ray properties of ageing historical core-collapse supernovae in nearby galaxies, using archival data from Chandra, XMM-Newton and Swift. We found possible evidence of a young X-ray pulsar in SN 1968D and in few other sources, but none more luminous than ~ a few 10^{37} erg/s. We compared the observational limits to the X-ray pulsar luminosity distribution with the results of Monte Carlo simulations for a range of birth parameters. We conclude that a pulsar population dominated by periods <~ 40 ms at birth is ruled out by the data.

The aim of this paper is to investigate the environment of gamma ray bursts (GRBs) and the interstellar matter of their host galaxies. We use to this purpose high resolution spectroscopic observations of the afterglow of GRB050922C, obtained with UVES/VLT 3.5 hours after the GRB event. We found that, as for most high resolution spectra of GRBs, the spectrum of the afterglow of GRB050922C is complex. At least seven components contribute to the main absorption system at z=2.1992. The detection of lines of neutral elements like MgI and the detection of fine-structure levels of the ions FeII, SiII and CII allows us to separate components in the GRB ISM along the line of sight. Moreover, in addition to the main system, we have analyzed the five intervening systems between z = 2.077 and z = 1.5664 identified along the GRB line of sight. GRB afterglow spectra are very complex, but full of information. This can be used to disentangle the contribution of the different parts of the GRB host galaxy and to study their properties. Our metallicity estimates agree with the scenario of GRBs exploding in low metallicity galaxies

GRB080319B reached 5th optical magnitude during the burst prompt emission. Thanks to the VLT/UVES rapid response mode, we observed its afterglow just 8m:30s after the GRB onset when the magnitude was R ~ 12. This allowed us to obtain the best signal-to-noise, high resolution spectrum of a GRB afterglow ever (S/N per resolution element ~ 50). The spectrum is rich of absorption features belonging to the main system at z=0.937, divided in at least six components spanning a total velocity range of 100 km/s. The VLT/UVES observations caught the absorbing gas in a highly excited state, producing the strongest Fe II fine structure lines ever observed in a GRB. A few hours later the optical depth of these lines was reduced by a factor of 4-20, and the optical/UV flux by a factor of ~ 60. This proves that the excitation of the observed fine structure lines is due to "pumping" by the GRB UV photons. A comparison of the observed ratio between the number of photons absorbed by the excited state and those in the Fe II ground state suggests that the six absorbers are ~ 2-6 kpc from the GRB site, with component I ~ 3 times closer to the GRB site than components III to VI. Component I is characterized also by the lack of Mg I absorption, unlike all other components. This may be due both to a closer distance and a lower density, suggesting a structured ISM in this galaxy complex.

We present the results from an ESO/VLT campaign aimed at studying the afterglow properties of the short/hard gamma ray burst GRB 070707. Observations were carried out at ten different epochs from ~0.5 to ~80 days after the event. The optical flux decayed steeply with a power-law decay index greater than 3, later levelling off at R~27.3 mag; this is likely the emission level of the host galaxy, the faintest yet detected for a short GRB. Spectroscopic observations did not reveal any line features/edges that could unambiguously pinpoint the GRB redshift, but set a limit z < 3.6. In the range of allowed redshifts, the host has a low luminosity, comparable to that of long-duration GRBs. The existence of such faint host galaxies suggests caution when associating short GRBs with bright, offset galaxies, where the true host might just be too dim for detection. The steepness of the decay of the optical afterglow of GRB 070707 challenges external shock models for the optical afterglow of short/hard GRBs. We argue that this behaviour might results from prolonged activity of the central engine or require alternative scenarios.

The paucity of reliable achromatic breaks in Gamma-Ray Burst afterglow light curves motivates independent measurements of the jet aperture. Serendipitous searches of afterglows, especially at radio wavelengths, have long been the classic alternative. These survey data have been interpreted assuming a uniformly emitting jet with sharp edges (``top-hat'' jet), in which case the ratio of weakly relativistically beamed afterglows to GRBs scales with the jet solid angle. In this paper, we consider, instead, a very wide outflow with a luminosity that decreases across the emitting surface. In particular, we adopt the universal structured jet (USJ) model, that is an alternative to the top-hat model for the structure of the jet. However, the interpretation of the survey data is very different: in the USJ model we only observe the emission within the jet aperture and the observed ratio of prompt emission rate to afterglow rate should solely depend on selection effects. We compute the number and rate of afterglows expected in all-sky snapshot observations as a function of the survey sensitivity. We find that the current (negative) results for OA searches are in agreement with our expectations. In radio and X-ray bands this was mainly due to the low sensitivity of the surveys, while in the optical band the sky-coverage was not sufficient. In general we find that X-ray surveys are poor tools for OA searches, if the jet is structured. On the other hand, the FIRST radio survey and future instruments like the Allen Telescope Array (in the radio band) and especially GAIA, Pan-Starrs and LSST (in the optical band) will have chances to detect afterglows.

It is now recognized that long-duration Gamma-Ray Bursts (GRBs) are linked to the collapse of massive stars, based on the association between (low-redshift) GRBs and (type Ic) core-collapse supernovae (SNe). The census of massive stars and GRBs reveals, however, that not all massive stars do produce a GRB. Only ~1% of core collapse SNe are able to produce a highly relativistic collimated outflow, and hence a GRB. The extra crucial parameter has long been suspected to be metallicity and/or rotation. We find observational evidence strongly supporting that both ingredients are necessary in order to make a GRB out of a core-collapsing star. A detailed study of the absorption pattern in the X-ray spectrum of GRB060218 reveals evidence of material highly enriched in low atomic number metals ejected before the SN/GRB explosion. We find that, within the current scenarios of stellar evolution, only a progenitor star characterized by a fast stellar rotation and sub-solar initial metallicity could produce such a metal enrichment in its close surrounding.

We have studied the X-ray emission from four historical Type-II supernovae (the newly-discovered 1941C in NGC 4631 and 1959D in NGC 7331; and 1968D, 1980K in NGC 6946), using Chandra ACIS-S imaging. In particular, the first three are the oldest ever found in the X-ray band, and provide constraints on the properties of the stellar wind and circumstellar matter encountered by the expanding shock at more advanced stages in the transition towards the remnant phase. We estimate emitted luminosities ~ 5 x 10^{37} erg/s for SN 1941C, ~ a few x 10^{37} erg/s for SN 1959D, ~ 2 x 10^{38} erg/s for SN 1968D, and ~ 4 x 10^{37} erg/s for SN 1980K, in the 0.3-8 keV band. X-ray spectral fits to SN 1968D suggest the presence of a harder component, possibly a power law with photon index ~ 2, contributing ~ 10^{37} erg/s in the 2-10 keV band. We speculate that it may be evidence of non-thermal emission from a Crab-like young pulsar.

We compute the average luminosity of X-ray flares as a function of time, for a sample of 10 long-duration gamma-ray burst afterglows. The mean luminosity, averaged over a timescale longer than the duration of the individual flares, declines as a power-law in time with index ~-1.5. We elaborate on the properties of the central engine that can produce such a decline. Assuming that the engine is an accreting compact object, and for a standard conversion factor between accretion rate and jet luminosity, the switch between a neutrino-cooled thin disk and a non-cooled thick disk takes place at the transition from the prompt to the flaring phase. We discuss the implications of this coincidence under different scenarios for the powering of the GRB outflow. We also show that the interaction of the outflow with the envelope of the progenitor star cannot produce flares out of a continuous relativistic flow, and conclude that it is the dynamics of the disk or the jet-launching mechanism that generates an intrinsically unsteady outflow on timescales much longer than the dynamical timescale of the system. This is consistent with the fact that X-ray flares are observed in short-duration GRBs as well as in long-duration ones.

Flares overlaid on the smooth power-law decay of Swift X-ray afterglows are rather common, appearing in roughly half the observed light curves. They are a manifestation of the late time activity of the inner engine, since their temporal evolution is too fast to be linked to activity taking place in the external shock blastwave. In this paper we show that the energy emitted in the form of flares decreases with time as a power-law. We discuss several possibilities in which the flares can be powered and the source of the observed variability. We show that late time accretion from a disk can provide the necessary energy input in both classes of short duration and long duration gamma-ray bursts.

X-ray observations unveiled the existence of enigmatic point-like sources at the centre of young supernova remnants (SNRs). These sources, dubbed Central Compact Objects (CCOs), are thought to be neutron stars formed by the supernova explosion. However, their multi-wavelength phenomenology is surprisingly different from that of most young neutron stars.The aim of this work is to understand the nature of the CCO 1WGA J1713.4-3949 in the G347.3-0.5 SNR, through deep optical and IR observations, the first ever performed for this source.By exploiting its derived Chandra X-ray position we carried out optical (BVI) observations with the NTT and Adaptive Optics IR (JHKs) observations with the VLT. We detected two faint (I~23.5, I~24.3.) patchy objects in the NTT images, close to the Chandra error circle. They were clearly resolved in our VLT images which unveiled a total of six candidate counterparts (17.8<H<20.3) with quite red colours (H-Ks~0.6). If they are stars, none of them can be associated with 1WGA J1713.4-3949 for the most likely values of distance and hydrogen column density. The identification of the faintest candidate with the neutron star itself can not be firmly excluded, while the identification with a fallback disk is ruled out by its non-detection in the I band. No other candidates are detected down to B~26, V~26.2, I~24.7, H~21.3 and K~20.5. Our high-resolution IR imaging of unveiled a few objects close/within the Chandra X-ray position of 1WGA J1713.4-3949. However, at present none of them can be firmly identified as its likely counterpart.

We present a proof-of-concept study that dust extinction curves can be extracted from the infrared (IR), optical, ultraviolet (UV) and X-ray afterglow observations of GRBs without assuming known extinction laws. We focus on GRB 050525A (z = 0.606), for which we also present IR observations from the Spitzer Space Telescope at about 2.3 days post-burst. We construct the spectral energy distribution (SED) of the afterglow and use it to derive the dust extinction curve of the host galaxy in 7 optical/UV wavebands. By comparing our derived extinction curve to known templates, we see that the Galactic or Milky Way extinction laws are disfavored versus the Small and Large Magellanic Cloud (SMC and LMC) ones, but that we cannot rule out the presence of a LMC-like 2175 angstrom bump in our extinction curve. The dust-to-gas ratio present within the host galaxy of GRB 050525A is similar to that found in the LMC, while about 10 to 40% more dust is required if the SMC template is assumed. Our method is useful to observatories that are capable of simultaneously observing GRB afterglows in multiple wavebands from the IR to the X-ray.

The temporal decay of the flux components of Transient Anomalous X-ray Pulsar XTE J1810-197 following its 2002 outburst presents a unique opportunity to probe the emission geometry of a magnetar. Toward this goal, we model the magnitude of the pulsar's modulation in narrow spectral bands over time. Following previous work, we assume that the post-outburst flux is produced in two distinct thermal components arising from a hot spot and a warm concentric ring. We include general relativistic effects on the blackbody spectra due to gravitational redshift and light bending near the stellar surface, which strongly depend on radius. This affects the model fits for the temperature and size of the emission regions. For the hot spot, the observed temporal and energy-dependent pulse modulation is found to require an anisotropic, pencil-beamed radiation pattern. We are able to constrain an allowed range for the angles that the line-of-sight (psi) and the hot spot pole (xi) make with respect to the spin-axis. Within errors, this is defined by the locus of points in the xi-psi plane that lie along the line (xi+beta(R))(psi+beta(R)) ~ constant, where beta(R) is a function of the radius R of the star. For a canonical value of R=12 km, the viewing parameters range from psi=xi=37 deg to (psi,xi)=(85 deg,15 deg). We discuss our results in the context of magnetar emission models.

<Context>. We report on near-infrared (IR) observations of the three anomalous X-ray pulsars XTE J1810-197, 1RXS J1708-4009, 1E 1841-045 and the soft gamma-ray repeater SGR 1900+14, taken with the ESO-VLT, the Gemini, and the CFHT telescopes. <Aims>. This work is aimed at identifying and/or confirming the IR counterparts of these magnetars, as well as at measuring their possible IR variability. <Methods>. In order to perform photometry of objects as faint as Ks~20, we have used data taken with the largest telescopes, equipped with the most advanced IR detectors and in most of the cases with Adaptive Optics devices. The latter are critical to achieve the sharp spatial accuracy required to pinpoint faint objects in crowded fields. <Results>. We confirm with high confidence the identification of the IR counterpart to XTE J1810-197, and its IR variability. For 1E 1841-045 and SGR 1900+14 we propose two candidate IR counterparts based on the detection of IR variability. For 1RXS J1708-4009 we show that none of the potential counterparts within the source X-ray error circle can be yet convincingly associated with this AXP. <Conclusions>. The IR variability of the AXP XTE J1810-197 does not follow the same monotonic decrease of its post-outburst X-ray emission. Instead, the IR variability appears more similar to the one observed in radio band, although simultaneous IR and radio observations are crucial to draw any conclusion in this respect. For 1E 1841-045 and SGR 1900+14, follow-up observations are needed to confirm our proposed candidates with higher confidence.

Traditionally, studies aimed at inferring the distribution of birth periods of neutron stars are based on radio surveys. Here we propose an independent method to constrain the pulsar spin periods at birth based on their X-ray luminosities. In particular, the observed luminosity distribution of supernovae poses a constraint on the initial rotational energy of the embedded pulsars, via the L_X-dot{E}_{rot} correlation found for radio pulsars, and under the assumption that this relation continues to hold beyond the observed range. We have extracted X-ray luminosities (or limits) for a large sample of historical SNe observed with Chandra, XMM and Swift, that have been firmly classified as core-collapse supernovae. We have then compared these observational limits with the results of Monte Carlo simulations of the pulsar X-ray luminosity distribution, for a range of values of the birth parameters. We find that a pulsar population dominated by millisecond periods at birth is ruled out by the data.

4U 2129+47 was discovered in the early 80's and classified as an accretion disk corona source due to its broad and partial X-ray eclipses. The 5.24 hr binary orbital period was inferred from the X-ray and optical light curve modulation, implying a late K or M spectral type companion star. The source entered a low state in 1983, during which the optical modulation disappeared and an F8 IV star was revealed, suggesting that 4U 2129+47 might be part of a triple system. The nature of 4U 2129+47 has since been investigated, but no definitive conclusion has been reached. Here, we present timing and spectral analyses of two XMM-Newton observations of this source, carried out in May and June, 2005. We find evidence for a delay between two mid-eclipse epochs measured ~22 days apart, and we show that this delay can be naturally explained as being due to the orbital motion of the binary 4U 2129+47 around the center of mass of a triple system. This result thus provides further support in favor of the triple nature of 4U 2129+47.

We report on two new XMM-Newton observations of the Anomalous X-ray Pulsar (AXP) 4U 0142+614 performed in March and July 2004, collecting the most accurate spectrum for this source to date. Furthermore, we analyse two short archival observations performed in February 2002 and January 2003 (the latter already reported by Gohler et al. 2005) in order to study the long term behaviour of this AXP. 4U 0142+614 appears to be relatively steady in flux between 2002 and 2004, and the phase-averaged spectrum does not show any significant variability between the four epochs. We derive the deepest upper limits to date on the presence of lines in the 4U 0142+614 spectrum as a function of energy: equivalent width in the 1-3 keV energy range < 4 eV and < 8 eV for narrow and broad lines, respectively. A remarkable energy dependence in both the pulse profile and the pulsed fraction is detected, and consequently pulse-phase spectroscopy shows spectral variability as a function of phase. By making use of XMM-Newton and INTEGRAL data, we successfully model the 1-250 keV spectrum of 4U 0142+614 with three models presented in Rea et al. (2007a), namely the canonical absorbed blackbody plus two power-laws, a resonant cyclotron scattering model plus one power-law and two log-parabolic functions.

X-ray observations have unveiled the existence of a family of radio-quiet Isolated Neutron Stars whose X-ray emission is purely thermal, hence dubbed X-ray Dim Isolated Neutron Stars (XDINSs). While optical observations have allowed to relate the thermal emission to the neutron star cooling and to build the neutron star surface thermal map, IR observations are critical to pinpoint a spectral turnover produced by a so far unseen magnetospheric component, or by the presence of a fallback disk. The detection of such a turnover can provide further evidence of a link between this class of isolated neutron stars and the magnetars, which show a distinctive spectral flattening in the IR. Here we present the deepest IR observations ever of five XDINSs, which we use to constrain a spectral turnover in the IR and the presence of a fallback disk. The data are obtained using the ISAAC instrument at the VLT. For none of our targets it was possible to identify the IR counterpart down to limiting magnitudes H = 21.5 - 22.9. Although these limits are the deepest ever obtained for neutron stars of this class, they are not deep enough to rule out the existence and the nature of a possible spectral flattening in the IR. We also derive, by using disk models, the upper limits on the mass inflow rate in a fallback disk. We find the existence of a putative fallback disk consistent (although not confirmed) with our observations.

Recent radio observations have unveiled the existence of a number of radio pulsars with spin-down derived magnetic fields in the magnetar range. However, their observational properties appears to be more similar to classical radio pulsars than to magnetars. To shed light on this puzzle we first have to determine whether the spin-down derived magnetic field values for these radio pulsars are indeed representative of the actual neutron star magnetic field or if they are polluted, e.g. by the effects of a torque from a fallback disk. To investigate this possibility, we have performed deep IR observations of one of these high magnetic field radio pulsars (PSR J1119-6127) with the ESO VLT to search for IR emission which can be associated with a disk. No IR emission is detected from the pulsar position down to J=24, H=23, Ks=22. By comparing our flux upper limits with the predictions of fallback disk models, we have found that we can only exclude the presence of a disk with accretion rate dot M >3x10^16 g/s. This lower limit cannot rule out the presence of a substantial disk torque on the pulsar, which would then lead to overestimate the value of the magnetic field inferred from P and dot P.

We investigate the properties and evolution of accretion tori formed after the coalescence of two compact objects. At these extreme densities and temperatures, the accreting torus is cooled mainly by neutrino emission produced primarily by electron and positron capture on nucleons (beta reactions). We solve for the disc structure and its time evolution by introducing a detailed treatment of the equation of state which includes photodisintegration of helium, the condition of beta-equilibrium, and neutrino opacities. We self-consistently calculate the chemical equilibrium in the gas consisting of helium, free protons, neutrons and electron-positron pairs and compute the chemical potentials of the species, as well as the electron fraction throughout the disc. We find that, for sufficiently large accretion rates (> 10 solar masses per second), the inner regions of the disk become opaque and develop a viscous and thermal instability. The identification of this instability might be relevant for GRB observations.

On September 21, 2006, an intense (~10^39 erg/s) and short (20 ms) burst was detected by Swift BAT at a position consistent with that of the candidate Anomalous X-ray Pulsar, CXOU J164710.2-455216, discovered by Chandra in 2005. Swift follow-up observations began about 13 hours after the event and found the source at a 1-10keV flux level of about 4.5 x 10^-11 erg/s/cm^2, i.e. ~300 times brighter than measured 5 days earlier by XMM. We report the results obtained from Swift BAT observations of the burst and subsequent Swift XRT observations carried out during the first four months after the burst. These data are complemented with those from two XMM observations (carried out just before and after the BAT event) and four archival Chandra observations carried out between 2005 and 2007. We find a phase coherent solution for the source pulsations after the burst. The evolution of the pulse phase comprises an exponential component decaying with timescale of 1.4d which we interpret as the recovery stage following a large glitch (Delta nu / nu about 6 x 10^-5). We also detect a quadratic component corresponding to a spin-down rate of Pdot ~ 9 x 10^-13, implying a magnetic field strength of 10^14 Gauss. During the first Swift XRT observation taken 0.6 days after the burst, the spectrum showed a kT = 0.65keV blackbody (R_BB = 1.5km) plus a Gamma=2.3 power-law accounting for about 60% of the 1-10 keV observed flux. Analysis of Chandra archival data, taken during 2005 when the source was in quiescence, reveal that the modulation in quiescence is 100% pulsed at energies above 4 keV and consistent with the (unusually small-sized) blackbody component being occulted by the neutron star as it rotates.

Long gamma-ray bursts (GRBs) are likely associated with the collapse of massive stars, which produce dust and are born in dusty environments. Absorption and scattering of ultraviolet/X-ray photons from the prompt, optical flash and afterglow emission of the GRB produce dust echoes. We perform time-dependent calculations of these echoes, accounting for the evolution of the dust grain distribution due to selective grain destruction by the GRB radiation, and for off-axis beaming. We explore cloud configurations of differing density and size -- the echo light curve and spectrum depend on the cloud radius, with larger clouds peaking at longer wavelengths. For a region ~ 3 pc in size with hydrogen density ~ 1000 per cubic centimeter, the echo spectrum peaks at ~ 3.6 microns and ~ 8.8 eV for thermal and scattered components, respectively. Dust echoes should be detectable with the Very Large Telescope up to z ~ 0.1, IRAC onboard the Spitzer Space Telescope up to z ~ 0.2, and NICMOS onboard the Hubble Space Telescope up to z ~ 0.3. Furthermore, the shape of the echo light curve allows one to infer: the jet opening angle; the inclination of the jet axis with respect to the line of sight; the size of the dust-emitting region. For sources with symmetric, bipolar jets, dust echoes exhibit two bumps in the light curve, making them easily distinguishable from the rebrightening due to an underlying supernova.

We analyze high resolution spectroscopic observations of the optical afterglow of GRB050730, obtained with UVES@VLT about hours after the GRB trigger. The spectrum shows that the ISM of the GRB host galaxy at z = 3.967 is complex, with at least five components contributing to the main absorption system. We detect strong CII*, SiII*, OI* and FeII* fine structure absorption lines associated to the second and third component. For the first three components we derive information on the relative distance from the site of the GRB explosion. Component 1, which has the highest redshift, does not present any fine structure nor low ionization lines; it only shows very high ionization features, such as CIV and OVI, suggesting that this component is very close to the GRB site. From the analysis of low and high ionization lines and fine structure lines, we find evidences that the distance of component 2 from the site of the GRB explosion is 10-100 times smaller than that of component 3. We evaluated the mean metallicity of the z=3.967 system obtaining values about 0.01 of the solar metallicity or less. However, this should not be taken as representative of the circumburst medium, since the main contribution to the hydrogen column density comes from the outer regions of the galaxy while that of the other elements presumably comes from the ISM closer to the GRB site. Furthermore, difficulties in evaluating dust depletion correction can modify significantly these values. The mean [C/Fe] ratio agrees well with that expected by single star-formation event models. Interestingly the [C/Fe] of component 2 is smaller than that of component 3, in agreement with GRB dust destruction scenarios, if component 2 is closer than component 3 to the GRB site.

We compute the probability that intermediate mass black holes (IMBHs) capture companions due to dynamical interactions and become accreting sources, and explore the possibility that the accreting IMBHs would appear as ultra-luminous X-ray sources (ULXs). We focus on IMBHs originating from low-metallicity Population III stars. Two channels of IMBH formation are considered: from primordial halos in the framework of hierarchical clustering, and from non-mixed, zero-metallicity primeval gas in galactic discs. IMBHs can form binary systems due to tidal captures of single stars and exchange interactions with existing binary systems in galactic discs. We find that neither formation mechanism of the accreting IMBH binary is able to provide enough sources to explain the observed population of ULXs. Even at sub-ULX luminosity, the total number of accreting IMBHs with $L> 10^{36}$ erg s$^{-1}$ with dynamically captured companions is found to be $< 0.01$ per galaxy.

Aims. We derive simple estimates of the maximum efficiency with which matter can be ejected by the propeller mechanism in disk-fed, rotating magnetic neutron stars. Some binary evolution scenarios envisage that this mechanism is responsible for expelling to infinity the mass inflowing at a low rate from the companion star, therefore limiting the total amount of mass that can be accreted by the neutron star. Methods. We demonstrate that, for typical neutron star parameters, a maximum of \eta_{pro} < 5.7 (P_{-3})^{1/3} times more matter than accreted can be expelled through the propeller mechanism at the expenses of the neutron star rotational energy (P_{-3} is the NS spin period in unit of 10E-3 s). Approaching this value, however, would require a great deal of fine tuning in the system parameters and the properties of the interaction of matter and magnetic field at the magnetospheric boundary. Results. We conclude that some other mechanism must be invoked in order to prevent that too much mass accretes onto the neutron stars of some low mass X-ray binaries.

The detection of bright X-ray flares superimposed on the regular afterglow decay in Swift gamma-ray bursts has triggered theoretical speculations on their origin. We study the temporal properties of flares due to internal dissipation and external shock mechanisms. We first show that at least a sizable fraction of the flares cannot be related to external shock mechanisms, since external shock flares evolve on much longer time scales than observed. We then study flares from internal dissipation, showing that the temporal properties allow us to distinguish the emission of slow early shells from that of late faster shells. We show that, due to the rapid evolution of the detected flares, it is most likely that the flares are produced by relatively fast shells ejected by the central engine shortly before they are observed. This implies that the central engine must be active for, in some cases, as long as one day. We finally discuss the constraints and implications that this observation has on the properties and physics of the inner engine, and we elaborate on possible future observational tests on the flare sample to further understand their origin and physics.

GRB050904 is the gamma-ray burst with the highest measured redshift. We performed time resolved X-ray spectroscopy of the late GRB and early afterglow emission. We find robust evidence for a decrease with time of the soft X-ray absorbing column. We model the evolution of the column density due to the flash ionization of the GRB and early afterglow photons. This allows us to constrain the metallicity and geometry of the absorbing cloud. We conclude that the progenitor of GRB050904 was a massive star embedded in a dense metal enriched molecular cloud with Z~0.03 Z_solar. This is the first local measurement of metallicity in the close environment of a GRB and one of the highest redshift metallicity measurements. We also find that the dust associated with the cloud cannot be similar to that of our Galaxy but must be either sizably depleted or dominated by silicate grains. We discuss the implications of these results for GRB progenitors and high redshift star formation.

High resolution spectroscopy of GRB 021004 revealed a wealth of absorption lines from several intermediate ionization species. The velocity structure of the absorber is complex and material with velocity up to >3000 km/s is observed. Since only the blueshifted component is observed, the absorber is very likely to be material closely surrounding the gamma-ray burst. We use a time-dependent photoionization code to track the abundance of the ions over time. Thanks to the presence of absorption from intermediate ionization states at long times, we can estimate the location and mass of the components of the absorber. We interpret those constraints within the hypernova scenario showing that the mass loss rate of the progenitor must have been ~10^{-4} solar masses per year, suggestive of a very massive star. In addition, the wind termination shock must lie at a distance of at least 100 pc, implying a low density environment. The velocity structure of the absorber also requires clumping of the wind at those large distances.

In the last years, optical studies of Isolated Neutron Stars (INSs) have expanded from the more classical rotation-powered ones to other categories, like the Anomalous X-ray Pulsars (AXPs) and the Soft Gamma-ray Repeaters (SGRs), which make up the class of the magnetars, the radio-quiet INSs with X-ray thermal emission and, more recently, the enigmatic Compact Central Objects (CCOs) in supernova remnants. Apart from 10 rotation-powered pulsars, so far optical/IR counterparts have been found for 5 magnetars and for 4 INSs. In this work we present some of the latest observational results obtained from optical/IR observations of different types of INSs.

Two sensitive Chandra X-ray observations of the heavily-reddened galactic starburst cluster Westerlund 1 in May and June 2005 detected a previously unknown X-ray pulsar (CXO J164710.20-455217). Its slow 10.6 s pulsations, moderate X-ray temperature kT $\approx$ 0.5 keV, and apparent lack of a massive companion tentatively suggest that it is an Anomalous X-ray Pulsar (AXP). An isothermal blackbody model yields an acceptable spectral fit but the inferred source radius is much less than that of a neutron star, a result that has also been found for other AXPs. We analyze the X-ray spectra with more complex models including a model that assumes the pulsar is a strongly magnetized neutron star (``magnetar'') with a light element atmosphere. We conclude that the observed X-ray emission cannot be explained as global surface emission arising from the surface of a cooling neutron star or magnetar. The emission likely arises in one or more localized regions (``hot spots'') covering a small fraction of the surface. We discuss these new results in the context of both accretion and magnetar interpretations for the X-ray emission.

We present a detailed study of the spectral and temporal properties of the X-ray and optical emission of GRB050713a up to 0.5 day after the main GRB event. The X-ray light curve exhibits large amplitude variations with several rebrightenings superposed on the underlying three-segment broken powerlaw that is often seen in Swift GRBs. Our time-resolved spectral analysis supports the interpretation of a long-lived central engine, with rebrightenings consistent with energy injection in refreshed shocks as slower shells generated in the central engine prompt phase catch up with the afterglow shock at later times. Our sparsely-sampled light curve of the optical afterglow can be fitted with a single power law without large flares. The optical decay index appears flatter than the X-ray one, especially at later times.

The energy output in the gamma-ray burst (GRB) prompt emission and afterglow phase is expected to photoionize the surrounding medium out to large radii. Cooling of this gas produces line emission, particularly strong in the optical, whose variability is a strong diagnostics of the gas density and geometry in the close environment of the burst. We present the results of a spectral time series analysis of the host galaxy of GRB 990712 observed up to ~6 years after the burst. We analyze the emission line fluxes together with those of the previous observations of the same GRB, in search for photoionization signatures. We find that the emission line fluxes show no variation within the uncertainities up to 6 years after the burst, and we use the measured line intensities to set a limit on the density of the gas within a few parsecs of the burst location. This is the first time that emission from cooling GRB remnants is probed on years time scales.

We present a theoretical study of double compact objects as potential short/hard gamma-ray burst (GRB) progenitors. An updated population synthesis code StarTrack is used to calculate properties of double neutron stars and black-hole neutron star binaries. We obtain their formation rates, estimate merger times and finally predict their most likely merger locations and afterglow properties for different types of host galaxies. Our results serve for a direct comparison with the recent HETE-II and SWIFT observations of several short bursts, for which afterglows and host galaxies were detected. We also discuss the possible constraints these observations put on the evolutionary models of double compact object formation. We emphasize that our double compact object models can successfully reproduce at the same time short GRBs within both young, star-forming galaxies (e.g., GRB 050709 and GRB 051221A), as well as within old, elliptical hosts (e.g., GRB 050724 and probably GRB 050509B)

Understanding the reasons for the faintness of the optical/near-infrared afterglows of the so-called dark bursts is essential to assess whether they form a subclass of GRBs, and hence for the use of GRBs in cosmology. With VLT and other ground-based telescopes, we searched for the afterglows of the INTEGRAL bursts GRB 040223, GRB 040422 and GRB 040624 in the first hours after the triggers. A detection of a faint afterglow and of the host galaxy in the K band was achieved for GRB 040422, while only upper limits were obtained for GRB 040223 and GRB 040624, although in the former case the X-ray afterglow was observed. A comparison with the magnitudes of a sample of afterglows clearly shows the faintness of these bursts, which are good examples of a population that an increasing usage of large diameter telescopes is beginning to unveil.

We present optical observations of the short/hard gamma-ray burst GRB 050709, the first such event with an identified optical counterpart. The object is coincident with a weak X-ray source and is located inside a galaxy at redshift z = 0.1606 +- 0.0002. Multiband photometry allowed us to study the broad-band spectral energy distribution. Late-time monitoring places strong limits on any supernova simultaneous with the GRB. The host galaxy is not of early type. Spectra show that the dominant stellar population is relatively young (~ 1 Gyr), and that ongoing star formation is present at a level of 2-3 (L/L*) Msun/yr. This is at least 2 orders of magnitude larger than that observed in the elliptical hosts of the short GRB 050509B and GRB 050724. This shows that at least some short GRBs originate in a young population. Short/hard GRB models based on the merger of a binary degenerate system are compatible with the host galaxy characteristics, although there is still the possibility of a connection between young stars and at least a fraction of such events.

Early-time X-ray observations of GRBs with the Swift satellite have revealed a more complicated phenomenology than was known before. In particular, the presence of flaring activity on a wide range of time scales probably requires late-time energy production within the GRB engine. Since the flaring activity is observed in both long and short GRBs, its origin must be within what is in common for the two likely progenitors of the two classes of bursts: a hyperaccreting accretion disk around a black hole of a few solar masses. Here, we show that some of the observational properties of the flares, such as the duration-time scale correlation, and the duration-peak luminosity anticorrelation displayed by most flares within a given burst, are qualitatively consistent with viscous disk evolution, provided that the disk at large radii either fragments or otherwise suffers large amplitude variability. We discuss the physical conditions in the outer parts of the disk, and conclude that gravitational instability, possibly followed by fragmentation, is the most likely candidate for this variability.

Accreting X-Ray Binaries display a wide range of behaviours. Some of them are observed to spin up steadily, others to alternate between spin-up and spin-down states, sometimes superimposed on a longer trend of either spin up or spin down. Here we interpret this rich phenomenology within a new model of the disk-magnetosphere interaction. Our model, based on the simplest version of a purely material torque, accounts for the fact that, when a neutron star is in the propeller regime, a fraction of the ejected material does not receive enough energy to completely unbind, and hence falls back into the disk. We show that the presence of this feedback mass component causes the occurrence of multiple states available to the system, for a given, constant value of the mass accretion rate dot{M}_* from the companion star. If the angle chi of the magnetic dipole axis with respect to the perpendicular to the disk is larger than a critical value chi_crit, the system eventually settles in a cycle of spin-up/spin-down transitions for a constant value of dot{M}_* and independent of the initial conditions. No external perturbations are required to induce the torque reversals. The transition from spin up to spin down is often accompanied by a large drop in luminosity. The frequency range spanned in each cycle and the timescale for torque reversals depend on dot{M}_*, the magnetic field of the star, the magnetic colatitude chi, and the degree of elasticity regulating the magnetosphere-disk interaction. The critical angle chi_crit ranges from \~25-30 deg for a completely elastic interaction to ~40-45 deg for a totally anelastic one. For chi ~< chi_crit, cycles are no longer possible and the long-term evolution of the system is a pure spin up. We specifically illustrate our model in the cases of the X-ray binaries GX 1+4 and 4U 1626-67.

At the extreme densities and temperatures typical of the central engine of GRBs, the accreting torus is cooled mainly by advection and by neutrino emission. The latter process is dominated by electron and positron capture onto nucleons ($\beta$ reactions). We calculate the reaction rates and the nuclear composition of matter, assuming that the torus consists of helium, eletron-positron pairs, free neutrons and protons. After determining the equation of state and solving for the disk structure for a given initial accretion rate, we subsequently follow its time evolution. We find that, for accretion rates of the order of $10 M_{\odot}$/s, likely typical for the early stages of the accretion event, the disk becomes unstable, giving rise to variable energy output. This instability may play an important role for producing internal shocks.

We present the spectral analysis of 14 gamma-ray bursts (GRB) X-ray afterglows in order to investigate the properties of interstellar matter (ISM) along the line of sight of GRB. We carried out a simultaneous analysis of the NIR-optical and X-band for those afterglows with an optical counterpart too, in order to evaluate and strongly constrain the absorption effect on the spectral energy distribution due to dust extinction from GRB environment. We evaluated the equivalent hydrogen column density Nh from X-ray spectroscopy and rest frame visual extinction Av by assuming different type of ISM composition and dust grain size distribution. From our analysis we obtained a distribution of the GRB rest frame consistent with the one expected if GRB were embedded in a Galactic-like molecular cloud. Moreover, values of the visual extinction estimated from the simultaneous analysis of NIR-to-X band favour an environment where small dust grain are destroyed by the interaction with the X-ray and UV photons from GRB.

We present an optical-to-X-ray spectral analysis of the afterglow of GRB 020405. The optical spectral energy distribution not corrected for the extragalactic extinction is significantly below the X-ray extrapolation of the single powerlaw spectral model suggested by multiwavelength studies. We investigate whether considerable extinction could explain the observed spectral ``mismatch'' by testing several types of extinction curves. For the first time we test extinction curves computed with time-dependent numerical simulations of dust grains destruction by the burst radiation. We find that an extinction law weakly depen dent on wavelength can reconcile the unabsorbed optical and X-ray data with the expected synchrotron spectrum. A gray extinction law can be provided by a dust grain size distribution biased toward large grains.

The sky region including the Chandra position of SGR 1806-20 was monitored in the IR band during 2004, following its increased high energy bursting activity. Observations were performed using NAOS-CONICA, the adaptive optics IR camera mounted on Yepun VLT, which provided images of unprecedented quality (FWHM better than 0.1"). After the 2004 December 27th giant flare, the source position has been nailed by VLA observations of its radio counterpart, reducing the positional uncertainty to 0.04". Using IR data from our monitoring campaign, we discovered the likely IR counterpart to SGR 1806-20 based on positional coincidence with the Chandra and VLA uncertainty regions and flux variability of a factor of about 2 correlated with that at higher energies. We compare our findings with other isolated neutron star classes thought to be related, at some level, with SGRs.

Cosmological models that include suppression of the power spectrum of density fluctuations on small scales exhibit an exponential reduction of high-redshift, non-linear structures, including a reduction in the rate of gamma ray bursts (GRBs). Here we quantify the constraints that the detection of distant GRBs would place on structure formation models with reduced small-scale power. We compute the number of GRBs that could be detectable by the Swift satellite at high redshifts (z > 6), assuming that the GRBs trace the cosmic star formation history, which itself traces the formation of non-linear structures. We calibrate simple models of the intrinsic luminosity function of the bursts to the number and flux distribution of GRBs observed by the Burst And Transient Source Experiment (BATSE). We find that a discovery of high-z GRBs would imply strong constraints on models with reduced small-scale power. For example, a single GRB at z > 10, or 10 GRBs at z > 5, discovered by Swift during its scheduled two-year mission, would rule out an exponential suppression of the power spectrum on scales below R_c=0.09 Mpc (exemplified by warm dark matter models with a particle mass of m_x=2 keV). Models with a less sharp suppression of small-scale power, such as those with a red tilt or a running scalar index, n_s, are more difficult to constrain, because they are more degenerate with an increase in the power spectrum normalization, sigma_8, and with models in which star-formation is allowed in low-mass minihalos. We find that a tilt of \delta n_s ~ 0.1 is difficult to detect; however, an observed rate of 1 GRB/yr at z > 12 would yield an upper limit on the running of the spectral index, alpha = d(n_s)/d(ln k) > -0.05.

GRB 040422 was detected by the INTEGRAL satellite at an angle of only 3 degrees from the Galactic plane. Analysis of the prompt emission observed with the SPI and IBIS instruments on INTEGRAL are presented. The IBIS spectrum is well fit by the Band model with a break energy of Eo=56+/-2 keV and Epeak=41+/-3 keV. The peak flux is 1.8 10^(-7) erg/cm2/s and fluence 3.4 10^(-7) erg/cm2 in the range 20-200 keV. We then present the observations of the afterglow of GRB 040422, obtained with the ISAAC and FORS 2 instruments at the VLT less than 2 hours after the burst. We report the discovery of its near-infrared afterglow, for which we give here the astrometry and photometry. No detection could have been obtained in the R and I bands, partly due to the large extinction in the Milky Way. We re-imaged the position of the afterglow two months later in the Ks band, and detected a likely bright host galaxy. We compare the magnitude of the afterglow with a those of a compilation of promptly observed counterparts of previous GRBs, and show that the afterglow of GRB 040422 lies at the very faint end of the distribution, brighter only than that of GRB 021211, singled out later and in the optical bands, and GRB 040924 after accounting for Milky Way extinction. This observation suggests that the proportion of dark GRBs can be significantly lowered by a more systematic use of 8-m class telescopes in the infrared in the very early hours after the burst.

X-ray flashes (XRFs) and X-ray rich gamma-ray bursts (XRGRBs) share many observational characteristics with long duration GRBs, but the reason for which their prompt emission peaks at lower photon energies, $E_p$, is still under debate. Although many different models have been invoked in order to explain the lower $E_p$ values, their implications for the afterglow emission were not considered in most cases, mainly because observations of XRF afterglows have become available only recently. Here we examine the predictions of the various XRF models for the afterglow emission, and test them against the observations of XRF 030723 and XRGRB 041006, the events with the best monitored afterglow light curves in their respective class. We show that most existing XRF models are hard to reconcile with the observed afterglow light curves, which are very flat at early times. Such light curves are, however, naturally produced by a roughly uniform jet with relatively sharp edges that is viewed off-axis (i.e. from outside of the jet aperture). This type of model self consistently accommodates both the observed prompt emission and the afterglow light curves of XRGRB 041006 and XRF 030723, implying viewing angles $\theta_{obs}$ from the jet axis of $(\theta_{obs}-\theta_0)\sim 0.15\theta_0$ and $\sim \theta_0$, respectively, where $\theta_0\sim 3$ deg is the jet half-opening angle. This suggests that GRBs, XRGRBs and XRFs are intrinsically similar relativistic jets viewed from different angles, corresponding to $\gamma(\theta_{obs}-\theta_0)$ of less than 1, between 1 and a few, and more than a few, respectively, where $\gamma$ is the Lorentz factor. Future observations with Swift could help test this unification scheme in which GRBs, XRGRBs and XRFs share the same basic physics and differ only by their orientation relative to our line of sight.

We present the first high resolution (R=20000--45000, corresponding to 14 km/s at 4200A to 6.6 km/s at 9000A) observations of the optical afterglow of Gamma Ray Bursts. GRB020813 and GRB021004 were observed by UVES@VLT 22.19 hours and 13.52 hours after the trigger, respectively. These spectra show that the inter--stellar matter of the GRB host galaxies is complex, with many components contributing to each main absorption system, and spanning a total velocity range of up to about 3000 km/s. Several narrow components are resolved down to a width of a few tens of km/s. In the case of GRB021004 we detected both low and high ionization lines. Combined with photoionization results obtained with CLOUDY, the ionization parameters of the various systems are consistent with a remarkably narrow range with no clear trend with system velocity. This can be interpreted as due to density fluctuations on top of a regular R^-2 wind density profile.

We investigate the consequences of a model of the assembly and growth of massive black holes from primordial seeds, remnants of the first generation of stars in a hierarchical structure formation scenario. Our model traces the build-up of MBHs from an early epoch, and follows the merger history of dark matter halos and their associated holes via Monte Carlo realizations of the merger hierarchy from early times to the present time. The sequence of minor and major mergers experienced by galactic halos in their hierarchical growth affects the merger history of MBHs embedded in their nuclei. So, if the formation route for the assembly of SMBHs dates back to the early universe, a large number of BH interactions is inevitable. Binary black holes coalescence timescales can be long enough for a third BH to fall in and interact with the central binary. These BH triple interactions lead typically to the final expulsion of one of the three bodies and to the recoil of the binary. Also, asymmetric emission of gravitational waves in the last stages of the black hole merging can give a recoil velocity to the centre of mass of the coalescing binary. This scenario leads to the prediction of a population of intermediate mass BHs (IMBHs) wandering in galaxy halos at the present epoch. We compute the luminosity distribution produced by these IMBHs accreting from their circumstellar medium. We find that in a Milky Way-sized galaxy they are unable to account for sources with luminosities > 1e39 erg/s unless they carry a baryonic remnant from which they are able to accrete for a long time. We also find that, for typical spiral galaxies, the bright end of the point source distribution correlates with the mass of the galaxy, and the most luminous sources are expected to be found in the disk.

We study the probability that three GRBs (980425, 030329, 031203) are found within z=0.17, given the luminosity functions consistent with the log N-log S relationship for classical cosmological bursts (i.e., those observed by BATSE). We show that, in order for the probability of these three low-z events to be non-negligible (thus making it more likely that they belong to the same class of the classical cosmological bursts), the bursts' luminosity function must be a broken powerlaw. By reasoning in analogy with beamed AGNs, we show that observations are consistent with the expectations if GRB 980425 and GRB 031203 are indeed normal bursts seen sideways. Within this model, no bright burst within z=0.17 should be observed by a HETE--like instrument within the next \sim 20 yr.

Rapid, hyper-Eddington accretion is likely to power the central engines of gamma-ray bursts (GRBs). In the extreme conditions of densities and temperatures the accreting torus is cooled by neutrino emission rather than by radiation. Another important cooling mechanism is the advection of energy into the central black hole. We compute the time evolution of a neutrino-dominated disc that proceeds during the burst and investigate the changes in its density and temperature. The discrimination between short and long bursts is made on the basis of the different rates of material inflow to the outer parts of the disc, thus favoring the binary merger scenario for the short GRBs and the collapsar scenario for the long ones. Within the context of the collapsar model, we also study the evolution of the photon luminosity of the remnant disc up to times of the order of 1 day, and discuss its implications for the production of emission lines in GRB spectra.

We report on observations aimed at searching for flux variations from the proposed IR counterpart of the Anomalous X-ray Pulsar XTE J1810-197. These data, obtained in March 2004 with the adaptive optics camera NAOS-CONICA at the ESO VLT, show that the candidate proposed by Israel et al. (2004) was fainter by Delta H=0.7+/-0.2 and Delta Ks=0.5+/-0.1 with respect to October 2003, confirming it as the IR counterpart of XTE J1810-197. We also report on an XMM-Newton observation carried out the day before the VLT observations. The 0.5-10 keV absorbed flux of the source was 2.2x10^-11 erg/s/cm^2, which is less by a factor of about two compared to the previous XMM-Newton observation on September 2003. Therefore, we conclude that a similar flux decrease took place in the X-ray and IR bands. We briefly discuss these results in the framework of the proposed mechanism(s) responsible for the IR variable emission of Anomalous X-ray Pulsars.

We investigate the possible contribution from the emission of accretion flows around supermassive black holes in early type galaxies to current measurements of the Cosmic Microwave Background (CMB) at radio frequencies. We consider a range of luminosities suggested by targeted radio observations and accretion models and compute the residual contribution of these sources to the spectrum and bispectrum of the observed CMB maps. As for high-resolution CMB experiments, we find that the unresolved component of these sources could make up to ~40-50% of the observed CBI and BIMA power spectrum at l > 2000. As a consequence, the inferred sigma_8^{SZ} value could be biased high by up to 6-7%. As for all sky experiments, we find that the contribution of accretion-flow sources to the WMAP bispectrum is at the 2-3 per cent level at most. At the flux limit that Planck will achieve, however, these sources may contribute up to 15 per cent of the bispectrum in the 60-100 GHz frequency range. Moreover, Planck should detect hundreds of these sources in the 30-300 GHz frequency window. These detections, possibly coupled with galaxy type confirmation from optical surveys, will allow number counts to put tighter constraints on early-type galaxies radio luminosity and accretion flows properties. These sources may also contribute up to the 30 per cent level to the residual radio sources power spectrum in future high-resolution SZ surveys (like ACT or APEX) reaching mJy flux limits.

Recent Chandra observations of nearby galaxies have revealed a number of ultraluminous X-ray sources (ULXs) with super-Eddington luminosities, away from the central regions of non-active galaxies. The nature of these sources is still debated. We argue that a fraction of them could be young, Crab-like pulsars, the X-ray luminosity of which is powered by rotation. We use the pulsar birth parameters estimated from radio pulsar data to compute the steady-state pulsar X-ray luminosity distribution as a function of the star formation rate (SFR) in the galaxy. We find that ~10% of optically dull galaxies are expected to have a source with L_x >~ 10^{39} erg/s, while starbursts galaxies should each have several of these sources. We estimate that the X-ray luminosity of a few percents of galaxies is dominated by a single bright pulsar with L_x >~10^{39} erg/s, roughly independently of its SFR. We discuss observational diagnostics that can help distinguish the young pulsar population in ULXs.

We report observation results of the prompt X-- and $\gamma$--ray emission from GRB 000528. This event was detected with the Gamma-Ray Burst Monitor and one of the two Wide Field Cameras aboard the BeppoSAX satellite. The GRB was promptly followed on with the BeppoSAX Narrow Field Instruments and with ground optical and radio telescopes. The X-ray afterglow source was identified, but no optical or radio counterpart found. We report here results from the prompt and afterglow emission analysis. The main feature derived from spectral evolution of the prompt emission is a high hydrogen-equivalent column density with evidence of its decrease with time. We model this behavior in terms of a time-dependent photoionization of the local circum-burst medium, finding that a compact and dense environment is required by the data. We also find a fading of the late part of the 2--10 keV prompt emission which is consistent with afterglow emission. We discuss this result in the light of the external shock model scenario.

We report the accurate sub-arcsec X-ray position of the new Anomalous X-ray Pulsar (AXP) XTE J1810-197, derived with a Chndra-HRC Target of Opportunity observation carried out in November 2003. We also report the discovery of a likely IR counterpart based on a VLT (IR band) Target of Opportunity observation carried out in October 2003. Our proposed counterpart is the only IR source (Ks=20.8) in the X-ray error circle. Its IR colors as well as the X-ray/IR flux ratio, are consistent with those of the counterparts of all other AXPs (at variance with field star colors). Deep Gunn-i band images obtained at the 3.6m ESO telescope detected no sources down to a limiting magnitude of 24.3. Moreover, we find that the pulsed fraction and count rates of XTE J1810-197 remained nearly unchanged since the previous Chandra and XMM-Newton observations (2003 August 27th and September 8th, respectively). We briefly discuss the implications of these results. In particular, we note that the transient (or at least highly variable) nature of this AXP might imply a relatively large number of hidden members of this class.

Traditional searches for radio pulsars have targeted individual small regions such as supernova remnants or globular clusters, or have covered large contiguous regions of the sky. None of these searches has been specifically directed towards giant supershells, some of which are likely to have been produced by multiple supernova (SN) explosions from an OB association. Here we perform a Montecarlo simulation of the pulsar population associated with supershells powered by multiple SNe. We predict that several tens of radio pulsars could be detected with current instruments associated with the largest Galactic supershells (with kinetic energies >~ 10^{53} ergs), and a few pulsars with the smaller ones. We test these predictions for some of the supershells which lie in regions covered by past pulsar surveys. For the smaller supershells, our results are consistent with the few detected pulsars per bubble. For the giant supershell GSH 242-03+37, we find the multiple SN hypothesis inconsistent with current data at the 95% level. We stress the importance of undertaking deep pulsar surveys in correlation with supershells. Failure to detect any pulsar enhancement in the largest of them would put serious constraints on the multiple SN origin for them. Conversely, the discovery of the pulsar population associated with a supershell would allow a different/independent approach to the study of pulsar properties.

Gamma Ray Bursts (GRBs) are expected to leave behind GRB remnants, similar to how ``standard'' supernovae (SN) leave behind SN remnants. The identification of these remnants in our own and in nearby galaxies would allow a much closer look at GRB birth sites, and possibly lead to the discovery of the compact object left behind. It would also provide independent constraints on GRB rates and energetics. We have initiated an observational program to search for GRB remnants in nearby galaxies. The identification is based on specific line ratios, such as OIII/$H_{\beta}$ and HeII/$H_{\beta}$, which are expected to be unusually high in case of GRB remnants according to the theoretical predictions of Perna et al. (2000). The observing strategies and preliminary studies from a test run at 2.34 m VBT as well as archival data from planetary nebulae surveys of spiral galaxies are discussed.

For many astrophysical sources with jets, there is evidence for the contemporaneous presence of disks. In contrast, pulsars such as the Crab and Vela show jets but have not yet revealed direct evidence for accretion disks. Here we show that for such pulsars, an accretion disk radiating below detectable thresholds may simultaneously account for (1) observed deviations in the braking indices from that of the simple dipole, (2) observed pulsar timing ages, and (3) possibly even the jet morphology via a disk outflow that interacts with the pulsar wind within, collimating and/or redirecting it.

During 2002-2003 the number of IR-identified counterparts to the Anomalous X-ray Pulsars (AXPs) has grown to four (4U0142+614, 2E2259+584, 1E 1048-59 and 1RXS J170849-400910) out of the six assessed objects of this class, plus two candidates. More importantly, some new common observational characteristics have been identified, such as the IR variability, the IR flattening in the broad band energy spectrum, the X-ray spectral variability as a function of pulse phase (which are not predicted by the magnetar model), and the SGR-like bursts (which can not be explained in terms of standard accretion models). We present the results obtained from an extensive multi-wavelength observational campaign carried out collecting data from the NTT, CFHT for the optical/IR bands, and XMM, Chandra (plus BeppoSAX archival data) in the X-rays. Based on these results and those reported in the literature, the IR-to-X-ray band emission of AXPs has been compared and studied.

We report Chandra observations of the dwarf nova WX Hyi in quiescence. The X-ray spectrum displays strong and narrow emission lines of N, O, Mg, Ne, Si, S and Fe. The various ionization states implied by the lines suggest that the emission is produced within a flow spanning a wide temperature range, from T ~ 10^6 K to T >~ 10^8 K. Line diagnostics indicate that most of the radiation originates from a very dense region, with n ~ 10^{13}-10^{14} cm^{-3}. The Chandra data allow the first tests of specific models proposed in the literature for the X-ray emission in quiescent dwarf novae. We have computed the spectra for a set of models ranging from hot boundary layers, to hot settling flows solutions, to X-ray emitting coronae. WX Hyi differs from other dwarf novae observed at minimum in having much stronger low temperature lines, which prove difficult to fit with existing models, and possibly a very strong, broad O VII line, perhaps produced in a wind moving at a few x 10^3 km/s. The accretion rate inferred from the X-rays is lower than the value inferred from the UV. The presence of high-velocity mass ejection could account for this discrepancy while at the same time explaining the presence of the broad O VII line. If this interpretation is correct, it would provide the first detection of a wind from a dwarf nova in quiescence.

We present a pulsar-atmosphere (PA) model for modulated thermal X-ray emission from cooling magnetized neutron stars. The model synthesizes the spectral properties of detailed stellar atmosphere calculations with the non-uniform surface properties anticipated for isolated, aging radio pulsars, and general relativistic effects on photon trajectories. We analyze the archival Chandra observations of the middle-aged radio pulsar PSR 0656+14 with the PA model and find it is an excellent representation of the phase-averaged X-ray spectrum for a broad range of polar effective temperature T_p and column density N_H ~ 10^{20} cm^{-2}. The spectral fits favor a sub-solar neutron star mass M < 1.0 M_{\sun}, a large radius R=15-16 km for a source distance d ~ 190 pc, while retaining consistency with theoretical neutron star equations-of-state. The modulated spectrum constrains the angular displacement of the pulsar spin axis to 30 +/- 2 \arcdeg with respect to the line of sight.

A large number of neutron stars (NSs), ~10^9, populate the Galaxy, but only a tiny fraction of them is observable during the short radio pulsar lifetime. The majority of these isolated NSs, too cold to be detectable by their own thermal emission, should be visible in X-rays as a result of accretion from the interstellar medium. The ROSAT all sky survey has however shown that such accreting isolated NSs are very elusive: only a few tentative candidates have been identified, contrary to theoretical predictions that up to several thousands should be seen. We suggest that the fundamental reason for this discrepancy lies in the use of the standard Bondi formula to estimate the accretion rates. We compute the expected source counts using updated estimates of the pulsar velocity distribution, realistic hydrogen atmosphere spectra, and a modified expression for the Bondi accretion rate as suggested by recent MHD simulations, and supported by direct observations in the case of accretion around supermassive black holes in nearby galaxies and in our own. We find that, whereas the inclusion of atmospheric spectra partly compensates for the reduction in the counts due to the higher mean velocities of the new distribution, the modified Bondi formula dramatically suppresses the source counts. The new predictions are consistent with a null detection at the ROSAT sensitivity.

The two leading interpretations of achromatic breaks that are observed in the light curves of GRBs afterglow are (i) the manifestation of the edge of a jet, which has a roughly uniform energy profile and a sharp edge and (ii) a line of sight effect in jets with a variable energy profile. The first scenario requires the inner engine to produce a jet with a different opening angle each explosion, while the latter requires a standard engine. The physical structure of the jet is a crucial factor in understanding GRB progenitors, and therefore discriminating the two jet scenarios is particularly relevant. In the structured jet case, specific predictions can be made for the distribution of observed break angles $\theta_{\rm break}$, while that distribution is arbitrary in the first scenario. We derive the theoretical distribution for the structured jet model. Specifically, we predict the most common angle to be about 0.12 rad, in rough agreement with the sample. If this agreement would hold as the sample size increases, it would strengthen the case for the standard jet hypothesis. We show that a prediction of this model is that the average viewing angle is an increasing function of the survey sensitivity, and in particular that a mission like {\em Swift} will find the typical viewing angle to be about 0.3 rad. The local event rate predicted by this model is $R_{\rm GRB}(z=0)\sim 0.5 $ Gpc$^{-3}$ yr$^{-1}$.

We report the discovery of the likely IR counterpart to the Anomalous X-ray Pulsar (AXP) 1RXS J1708-4009, based on the combination of Chandra HRC-I X-ray position, and deep optical/IR observations carried out from ESO telescopes and the Canada France Hawaii Telescope (CFHT) during 1999-2002. Within the narrow uncertainty region we found two relatively faint (K'=20.0 and K'=17.53) IR objects. Based on their color and position in the J-K' versus J-H diagram only the brighter object is consistent with the known IR properties of the counterparts to other AXPs. No variability was detected for this source, similarly to what is observed in the case of 4U 0142+614. Like in other AXPs, we found that the IR flux of 1RXSJ1708-4009 is higher then expected for a simple blackbody component extrapolated from the X-ray data. If confirmed, this object would be the fourth IR counterpart to a source of the AXP class, and would make the IR excess a likely new characteristic of AXPs.

We present a broadband modeling of the afterglow of GRB 021004. The optical transient of this burst has been detected very early and densely sampled in several bands. Its light curve shows significant deviations from a simple power law. We use the data from the X-ray to the $I-$band gathered in the first month of observations, and examine three models. Two involve variations in the energy of the shock. The first (energy injection) allows only increases to the shock energy, while the second (patchy shell) allows the energy to increase or decrease. In the final model (clumpy medium) the energy of the shock is constant while the density varies. While all three models reproduce well the optical bands, the variable density model can best account for the X-ray data, and the energy-injection model has the poorest fit. None of the models can account for the modest color variations observed during the first few days of the burst.

We have examined the complete set of X-ray afterglow observations of dark and optically bright GRBs performed by BeppoSAX until February 2001. X-ray afterglows are detected in 90% of the cases. We do not find significant differences in the X-ray spectral shape, in particular no higher X-ray absorption in GRBs without optical transient (dark GRBs) compared to GRBs with optical transient (OTGRBs). Rather, we find that the 1.6-10 keV flux of OTGRBs is on average about 5 times larger than that of the dark GRBs. A K-S test shows that this difference is significant at 99.8% probability. Under the assumption that dark and OTGRB have similar spectra, this could suggest that the first are uncaught in the optical band because they are just faint sources. In order to test this hypothesis, we have determined the optical-to-X ray flux ratios of the sample. OTGRBs show a remarkably narrow distribution of flux ratios, which corresponds to an average optical-to-x spectral index 0.794\pm 0.054. We find that, while 75% of dark GRBs have flux ratio upper limits still consistent with those of OT GRBs, the remaining 25% are 4 - 10 times weaker in optical than in X-rays. The significance of this result is equal to or higher than 2.6 sigma. If this sub-population of dark GRBs were constituted by objects assimilable to OTGRBs, they should have shown optical fluxes higher than upper limits actually found. We discuss the possible causes of their behaviour, including a possible occurrence in high density clouds or origin at very high redshift and a connection with ancient, Population III stars.

We use the time-dependent photoionization and dust destruction code developed by Perna & Lazzati (2002) to study the time evolution of the medium in a dusty gaseous cloud illuminated by a bright central source that sets on at time zero. We study the case of a bright source, which lasts for a time scale much smaller than the recombination and dust creation time scales. For this reason an equilibrium is never reached. We show that the presence of dust and its properties, such as its composition, can have a big effect on the time scale for the evaporation of the soft X-ray absorption, in particular for ionizing sources with hard spectra. We discuss the profile of evaporation of the soft X-ray absorbing column, as well as how the apparent ionization state of the cloud evolves in time. We finally consider the apparent metallicity of the cloud that is left behind as a function of the cloud and ionizing source properties.

The interaction of a radiation field with a dusty medium is a relevant issue in several astrophysical contexts. We use the time-dependent photoionization code in a dusty medium developed by Perna & Lazzati (2002), to study the modifications in the dust distribution and the relative optical opacities when a strong X-ray UV radiation flux propagates into a medium. We find that silicates are preferentially destroyed with respect to graphite, and the extinction curve becomes significantly flatter (hence implying less reddening), with the characteristic bump at lambda 2175 A highly suppressed, due to the destruction of the small graphite grains. This could explain the observational lack of such a feature in GRB afterglow and AGN spectra. For a very intense and highly variable source irradiating a compact and dense region, time variability in the optical opacity resulting from dust destruction can be observed on a relatively short timescale. We show that, under these circumstances, monitoring the time variability of the opacity can yield powerful clues on the properties of dust in the environment of the source. In particular, it allows to break the observational degeneracy as to whether a grey extinction is the result of a low dust-to-gas ratio or of a dust grain distribution that is skewed towards large grains.

Many models of Gamma Ray Bursts invoke a central engine consisting of a black hole of a few solar masses accreting matter from a disk at a rate of a fraction to a few solar masses per second. Popham et al. and Narayan et al. have shown that, for Mdot >~ 0.1 Msun/s, accretion proceeds via neutrino cooling and neutrinos can carry away a significant amount of energy from the inner regions of the disks. We improve on these calculations by including a simple prescription for neutrino transfer and neutrino opacities in such regions. We find that the flows become optically thick to neutrinos inside a radius R~6-40R_s for Mdot in the range of 0.1 -10 Msun/s, where R_s is the black hole Schwarzchild radius. Most of the neutrino emission comes from outside this region and, the neutrino luminosity stays roughly constant at a value L_{\nu} ~ 10^{53} erg/s. We show that, for Mdot > 1 Msun/s, neutrinos are sufficiently trapped that energy advection becomes the dominant cooling mechanism in the flow. These results imply that neutrino annihilation in hyperaccreting black holes is an inefficient mechanism for liberating large amounts of energy. Extraction of rotational energy by magnetic processes remains the most viable mechanism.

We present a time-dependent photoionization code that combines self-consistently metal evolution and dust destruction under an intense X-ray UV radiation field. Firstly, we extend the mathematical formulation of the time-dependent evolution of dust grains under an intense radiation flux with the inclusion of the process of ion field emission (IFE). We determine the relative importance of IFE with respect to X-ray and UV sublimation as a function of grain size, intensity and hardness of the incident spectrum. We then combine the processes of dust destruction with a photoionization code that follows the evolution of the ionization states of the metals and the relative radiative transitions. Our code treats, self-consistently, the gradual recycling of metals into gas as dust is sublimated away; it allows for any initial dust grain distribution and follows its evolution in space and time. In this first paper, we use our code to study the time-dependent behaviour of the X-ray and optical opacities in the nearby environment of a Gamma-ray Burst, and show how the time variability of the low-energy and high-energy opacities can yield powerful clues on the characteristics of the medium in which the bursts occur.

We derive a fully relativistic, self-similar solution to describe the propagation of a shock along an exponentially decreasing atmosphere, in the limit of very large Lorentz factor. We solve the problem in planar symmetry and compute the acceleration of the shock in terms of the density gradient crossed during its evolution. We apply our solution to the acceleration of shocks within the atmosphere of a HyperNova, and show that velocities consistent with the requirements of GRB models can be achieved with exponential atmospheres spanning a wide density range.

Despite great observational and theoretical effort, the burst progenitor is still a mysterious object. It is generally accepted that one of the best ways to unveil its nature is the study of the properties of the close environment in which the explosion takes place. We discuss the potentiality and feasibility of time resolved X-ray spectroscopy, focusing on the prompt gamma-ray phase. We show that the study of absorption features (or continuum absorption) can reveal the radial structure of the close environment, unaccessible with different techniques. We discuss the detection of absorption in the prompt and afterglow spectra of several bursts, showing how these are consistent with gamma-ray bursts taking place in dense regions. In particular, we show that the radius and density of the surrounding cloud can be measured through the evolution of the column density in the prompt burst phase. The derived cloud properties are similar to those of the star forming cocoons and globules within molecular clouds. We conclude that the burst are likely associated with the final evolutionary stages of massive stars.

GRB data accumulated over the years have shown that the distribution of their time duration is bimodal. While there is some evidence that long bursts are associated with star-forming regions, nothing is known regarding the class of short bursts. Their very short timescales are hard to explain with the collapse of a massive star, but would be naturally produced by the merger of two compact objects, such as two neutron stars (NS), or a neutron star and a black hole (BH). As for the case of long bursts, afterglow obervations for short bursts should help reveal their origin. By using updated population synthesis code calculations, we simulate a cosmological population of merging NS-NS and NS-BH, and compute the distribution of their galactic off-sets, the density distribution of their environment, and, if indeed associated with GRBs, their expected afterglow characteristics.

We show that accretion disks in several compact X-ray binaries with hydrogen-depleted donors are likely subject to a thermal ionization instability, unless they are strongly irradiated. These disks are particularly interesting in that their MHD-turbulent properties in the neutral phase may be quite different from those of standard, hydrogen-rich disks.

We investigate through dedicated numerical simulations the evolution of the soft X-ray absorption properties of a cloud surrounding a gamma-ray burst source. We show that the absorption properties of the material are strongly modified by the ionization induced by the intense burst flux. We derive the temporal evolution of the measured column density as a function of the density and size of the absorbing medium. Even if their statistical significance is not extremely compelling, we find that the detection in several bursts of variable absorption during the gamma-ray phase can be accounted for if these bursts are associated to overdense regions in molecular clouds with properties similar to those of star formation globules. We fit our model variable column density to the data of GRB 980329 and GRB 780506, showing that with this method the size, density and density distribution of the material surrounding a burst can be measured.

Absorption or emission against the cosmic microwave background radiation (CMB) may be observed in the redshifted 21cm line if the spin temperature of the neutral intergalactic medium prior to reionization differs from the CMB temperature. This so-called 21cm tomography should reveal important information on the physical state of the intergalactic medium at high redshifts. The fluctuations in the redshifted 21 cm, due to gas density inhomogeneities at early times, should be observed at meter wavelengths by the next generation radio telescopes such as the proposed {\it Square Kilometer Array (SKA)}. Here we show that the extra-galactic radio sources provide a serious contamination to the brightness temperature fluctuations expected in the redshifted 21 cm emission from the IGM at high redshifts. Unless the radio source population cuts off at flux levels above the planned sensitivity of SKA, its clustering noise component will dominate the angular fluctuations in the 21 cm signal. The integrated foreground signal is smooth in frequency space and it should nonetheless be possible to identify the sharp spectral feature arising from the non-uniformities in the neutral hydrogen density during the epoch when the first UV sources reionize the intergalactic medium.

The recent detection of a transient absorption feature in the X-ray prompt emission of GRB 990705 showed the importance of such observations in the understanding of gamma-ray bursts and their progenitors. We investigate the time dependence of photoionization edges during the prompt emission of bursts in different environments. We show that their variability can be used to infer the density and geometry of the surrounding medium, giving important clues to unveil the nature of the burst progenitor.

We show that thin accretion disks made of Carbon or Oxygen are subject to the same thermal ionization instability as Hydrogen and Helium disks. We argue that the instability applies to disks of any metal content. The relevance of the instability to supernova fallback disks probably means that their power-law evolution breaks down when they first become neutral. We construct simple analytical models for the viscous evolution of fallback disks to show that it is possible for these disks to become neutral when they are still young (ages of a few 10^3 to 10^4 years), compact in size (a few 10^9 cm to 10^11 cm) and generally accreting at sub-Eddington rates (Mdot ~ a few 10^14 - 10^18 g/s). Based on recent results on the nature of viscosity in the disks of close binaries, we argue that this time may also correspond to the end of the disk activity period. Indeed, in the absence of a significant source of viscosity in the neutral phase, the entire disk will likely turn to dust and become passive. We discuss various applications of the evolutionary model, including anomalous X-ray pulsars and young radio pulsars. Our analysis indicates that metal-rich fallback disks around newly-born neutron stars and black holes become neutral generally inside the tidal truncation radius (Roche limit) for planets, at \~10^11 cm. Consequently, the efficiency of the planetary formation process in this context will mostly depend on the ability of the resulting disk of rocks to spread via collisions beyond the Roche limit. It appears easier for the merger product of a doubly degenerate binary, whether it is a massive white dwarf or a neutron star, to harbor planets because it can spread beyond the Roche limit before becoming neutral.[Abridged]

We present a model for the spindown of young radio pulsars in which the neutron star loses rotational energy not only by emitting magnetic dipole radiation but also by torquing a surrounding accretion disk produced by supernova fallback. The braking index predicted in our model is in general less than n=3 (the value for pure dipole magnetic radiation), in agreement with the reported values of n<3 for five young radio pulsars. With an accuracy of 30% or better, our model reproduces the age, braking index and third frequency derivative of the Crab pulsar for a disk accretion rate in the range 3 x 10^16 - 10^17 g/s.

By now there is substantial evidence that Gamma-Ray Bursts (GRBs) originate at cosmological distances from very powerful explosions. The interaction between a GRB and its surrounding environment has dramatic consequences on the environment itself. At early times, the strong X-ray UV afterglow flux photoionizes the medium on distance scales on the order of 100 pc or more. Here I discuss the long-term effects resulting from the interaction between a GRB and its environment, and in particular the signatures of the emission spectrum produced while the heated and ionized gas slowly cools and recombines. Besides photoionizing the medium with its afterglow, a GRB explosion drives a blast wave which is expected to have a very long lifetime. I discuss possible candidates for such GRB remnants in our own and in nearby galaxies, and ways to distinguish them from remnants due to other phenomena, such as multiple supernova (SN) explosions.

We present a simple, unified model which accounts for properties of the X-ray emission from the three middle-aged pulsars PSR 1055-52, PSR 0656+14 and PSR 0630+18 (Geminga). The X-ray radiation from these objects is pulsed more strongly at energies above a transition point around 0.5 keV. In addition, the phase of the pulses shifts by about 80-100 degrees around the same point. Geminga also has the peculiarity that its pulsed fraction decreases in the 0.3-0.5 keV energy range, attaining a minimum near 0.5 keV. We show that a two-component hydrogen atmosphere is able to account for these disparate features. In our model, the hotter component is powered by particle bombardment and is restricted to the polar regions, while the cooler one covers the entire stellar surface. The two components also differ in their emission patterns, with the hard and soft contributions coming from areas radiating into fan and pencil beams, respectively.

Most of the baryons in the low-redshift universe reside in a warm/hot component which is difficult to detect with standard absorption/emission line techniques. We propose to use quasar refractive scintillation as an useful, complementary probe of such ionized, intergalactic gas. In particular, an application to the case of the intracluster medium is presented. We show that clusters located at z approx 0.02 should produce a source rms intensity fluctuation at 50-100 GHz of several tens of percent and on time scales ranging from days to months, depending on the projected location of the source on the foreground cluster. However, in order to produce such a signal, the source needs to be very compact. This effect, if observed, can be used as an independent test of the baryonic mass fraction in clusters.

The energy source powering the X-ray emission from anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) is still uncertain. In one scenario, the presence of an ultramagnetized neutron star, or ``magnetar'', with B on the order of 10^{14} - 10^{15} G is invoked. To investigate this hypothesis, we have analyzed archival ASCA data for several known AXPs and SGRs, and fitted them with a model in which all or part of the X-ray flux originates as thermal emission from a magnetar. Our magnetar spectral model includes the effects of the anisotropy of the heat flow through an ultramagnetized neutron star envelope, reprocessing by a light element atmosphere, and general relativistic corrections to the observed spectrum. We obtain good fits to the data with radii for the emitting areas which are generally consistent with those expected for neutron stars, in contrast to blackbody (BB) fits, which imply much smaller radii. Furthermore, the inclusion of atmospheric effects results in inferred temperatures which are lower than those implied by BB fits, but however still too high to be accounted by thermal cooling alone. An extra source of heating (possibly due to magnetic field decay) is needed. Despite the harder tail in the spectrum produced by reprocessing of the outgoing flux through the atmosphere, spectral fits still require a considerable fraction of the flux to be in a power-law component.

Rather than accreting via a disk, some White Dwarfs (WDs) in quiescent Dwarf Novae (DN) could accrete via an Advection-Dominated Accretion Flow (ADAF) possibly responsible for the X-ray Bremsstrahlung emission observed. Such a hot accretion flow is also expected to produce characteristic thermal line emission. Using SS Cyg as a specific example, we show that knowing a priori the inclination and the WD mass in quiescent DN makes X-ray line diagnostics powerful probes of the flow structure in these systems. Current X-ray instruments can discriminate, from their width, between lines emitted from a flow with a Keplerian rotation rate and those emitted at a substantially sub-Keplerian rate. This could be used to observationally test the property of energy advection, which is at the origin of the partial radial pressure support by the hot gas and the sub-Keplerian rotation rate in an ADAF.

We examine the reprocessing of X-ray radiation from compact objects by accretion disks when the X-ray emission from the star is highly beamed. The reprocessed flux for various degrees of beaming and inclinations of the beam axis with respect to the disk is determined. We find that, in the case where the beam is produced by a non-relativistic object, the intensity of the emitted spectrum is highly suppressed if the beam is pointing away from the disk. However, for beams produced by compact objects, general relativistic effects cause only a small reduction in the reradiated flux even for very narrow beams oriented perpendicularly to the disk. This is especially relevant in constraining models for the anomalous X-ray pulsars, whose X-ray emission is highly beamed. We further discuss other factors that can influence the emission from disks around neutron stars.

We examine properties of thermal radiation emitted by strongly magnetized neutron stars (NSs). In particular, we show that the pulsation amplitudes of the energy-integrated flux are an increasing function of the absorption column density to the source. This is especially important for the interpretation of the Anomalous $X$-ray Pulsars (AXPs) as cooling neutron stars with high magnetic fields. We show that the high-pulsation amplitudes observed in these objects are consistent with cooling models, if the large amount of absorption to these sources is taken into account. We also show that cooling models imply inferred radii of the emitting regions on the order of $\sim 5-6$ times smaller than the actual NS radii, again in agreement with observations.

The amount and properties of high-redshift galactic and intergalactic (IG) dust are largely unknown, but could be investigated using multi-wavelength photometry of high-z objects that have a known intrinsic spectrum. Observations of gamma-ray burst (GRB) afterglows appear to support the theoretical model of an adiabatic blast wave expanding into an external medium. In this model, the synchrotron peak flux is independent of frequency, providing a flat spectrum when observed over time, and therefore allowing straightforward measurement of the relative attenuation of afterglow flux in widely separated bands. Applying this method to dust extinction, we show that for a sample of afterglows which have been corrected by galactic extinction, comparison between the number counts of peak fluxes in X-ray versus optical can provide constraints on an intergalactic component of dust. A similar technique can probe the redshift-dependence of extinction in GRB-forming regions without requiring an assumed relation between extinction and reddening by the dust. Probing systematic changes in extinction with redshift - particularly in IG and/or non-reddening dust - is crucial to a proper interpretation of the Type Ia Supernova Hubble diagram and similar observations, and useful in understanding GRB progenitor environments.

We compute X-ray emission lines from thermal plasma in hot accretion flows. We show that line profiles are strong probes of the gas dynamics, and we present line-ratio diagnostics which are sensitive to the distribution of mass with temperature in the flow. We show how these can be used to constrain the run of density with radius, and the size of the hot region. We also present diagnostics which are primarily sensitive to the importance of recombination versus collisional ionization, and which could help discriminate ADAFs from photoionization-dominated accretion disk coronae. We apply our results to the Galactic center source Sagittarius A* and to the nucleus of M87. We find that the brightest predicted lines are within the detection capability of current $X$-ray instruments.

Current estimates of number counts of radio sources in the frequency range where the most sensitive Cosmic Microwave Background (CMB) experiments are carried out significantly under-represent sources with strongly inverted spectra. Hot accretion flows around supermassive black holes in the nuclei of nearby galaxies are expected to produce inverted radio spectra by thermal synchrotron emission. We calculate the temperature fluctuations and power spectra of these sources in the Planck Surveyor 30 GHz energy channel, where their emission is expected to peak. We find that their potential contribution is generally comparable to the instrumental noise, and approaches the CMB anisotropy level at small angular scales. Forthcoming CMB missions, which will provide a large statistical sample of inverted-spectra sources, will be crucial for determining the distribution of hot accretion flows in nearby quiescent galactic nuclei. Detection of these sources in different frequency channels will help constrain their spectral characteristics, hence their physical properties.

The energy source of the giant HI supershells in nearby galaxies and in our own is still an unsettled issue. Proposed scenarios for production of near-complete ringlike supershells are multiple supernovae (SNe) and gamma-ray bursts (GRBs). In the late phase of evolution it is difficult to tell these models apart. We show that, if a supershell has been powered by multiple SNe, peculiar metal abundances should be found in the medium within the bubble. We present line diagnostics that could detect such unusual abundances, such as unusually high ratios of [O]/[Fe] and [Ne]/[Fe]. Among ions of the same element, a higher enhancement is expected in lines from a high-ionization state than in lines from a low-ionization state. Identification of the energy source of HI supershells would set strong constraints on the rates and energetics of GRBs, as well as on their location within a galaxy.

The nature of the energy source powering anomalous X-ray pulsars is uncertain. Proposed scenarios involve either an ultramagnetized neutron star, or accretion onto a neutron star. We consider the accretion model proposed recently by Chatterjee, Hernquist & Narayan, in which a disk is fed by fallback material following a supernova. We compute the optical, infrared, and submillimeter emission expected from such a disk, including both viscous dissipation and reradiation of X-ray flux impinging on the disk from the pulsar. We find that it is possible with current instruments to put serious constraints on this and on other accretion models of AXPs. Fallback disks could also be found around isolated radio pulsars and we compute the corresponding spectra. We show that the excess emission in the R and I bands observed for the pulsar PSR 0656+14 is broadly consistent with emission from a disk.

We study the spectral signatures arising from cooling and recombination of an interstellar medium whose equilibrium state has been altered over \sim 100 pc by the radiation of a Gamma-Ray Burst (GRB) and its afterglow. We identify signatures in the line diagnostics which are indicative of a photo-ionized GRB remnant which is \la 5 x 10^4 years old . We estimate that at least a few such remnants should be detectable in the Virgo cluster of galaxies. If the gamma-ray emission from GRBs is beamed to a fraction f_b of their sky, then the expected number of Virgo remnants is larger by a factor of f_b^{-1}. Virgo remnants can be resolved with arcsecond imaging, and are likely to be center-filled using narrow-band filters of high ionization lines (such as [O III] \lambda5007 or He II \lambda4686), and limb-brightened for low-ionization lines (such as [S II] \lambda6717). The non-relativistic blast wave might be visible separately, since it does not reach the outer edge of these young photo-ionized remnants. The remnants should show evidence for ionization cones if the prompt or afterglow UV emission from GRBs is beamed.

The degree of beaming in Gamma-Ray Bursts (GRBs) is currently unknown. The uncertainty in the gamma-ray beaming angle, theta, leaves the total energy release (proportional to theta^2) and the event rate per galaxy (proportional to 1/theta^2) unknown to within orders of magnitude. Since the delayed radio emission of GRB sources originates from a mildly relativistic shock and receives only weak relativistic beaming, the rate of radio-selected transients with no GRB counterpart can be used to set an upper limit on 1/theta^2. We find that a VLA survey with a sensitivity of 0.1 mJy at 10 GHz could identify more than 2e4/(theta/10degrees)^2 radio afterglows across the sky if each source is sampled at least twice over a period of one month or longer. From the total number of >0.1 mJy sources observed at 8.44 GHz and the fraction of fading sources at 1.44 GHz, we get the crude limit, theta>6 degrees.

Gamma-Ray Bursts (GRBs) are thought to originate at cosmological distances from the most powerful explosions in the Universe. If GRBs are not beamed then the distribution of their number as a function of Gamma-ray flux implies that they occur once per (0.3-40) million years per bright galaxy and that they deposit >10^{53} ergs into their surrounding interstellar medium. The blast wave generated by a GRB explosion would be washed out by interstellar turbulence only after tens of millions of years when it finally slows down to a velocity of 10 km/s. This rather long lifetime implies that there could be up to several tens of active GRB remnants in each galaxy at any given time. For many years, radio observations have revealed the enigmatic presence of expanding neutral-hydrogen (HI) supershells of kpc radius in the Milky Way and in other nearby galaxies. The properties of some supershells cannot be easily explained in terms of conventional sources such as stellar winds or supernova explosions. However, the inferred energy and frequency of the explosions required to produce most of the observed supershells agree with the above GRB parameters. More careful observations and analysis might reveal which fraction of these supershells are GRB remnants. We show that if this link is established, the data on HI supershells can be used to constrain the energy output, the rate per galaxy, the beaming factor, and the environment of GRB sources in the Universe.

The current census of observed baryons in the local Universe is still missing a significant fraction of them according to standard Big-Bang nucleosynthesis. Numerical simulations predict that most of the missing baryons are in a hot intergalactic medium, which is difficult to observe through its X-ray emission or Sunyaev-Zel'dovich effect. We show that the next generation of X-ray satellites will be able to detect this gas through the X-ray absorption lines imprinted by its highly-ionized metals on the spectrum of a background quasar. For the metallicity typically found in intracluster gas, up to 70% of the baryons produce O VIII absorption lines with an equivalent width >0.1eV. The spectrum of any high redshift quasar is expected to show several such lines per unit redshift due to intervening gaseous halos of galaxy groups. These lines will be detectable at a signal-to-noise ratio >5 after a day of integration with the future Constellation-X telescope for any of the several hundred brightest quasars across the sky.

The discovery of Gamma-Ray Burst (GRB) afterglows revealed a new class of variable sources at optical and radio wavelengths. At present, the environment and precise redshift of the detected afterglows are still unknown. We show that if a GRB source resides in a compact (<100pc) gas-rich environment, the afterglow spectrum will show time-dependent absorption features due to the gradual ionization of the surrounding medium by the afterglow radiation. Detection of this time-dependence can be used to constrain the size and density of the surrounding gaseous system. For example, the MgII absorption line detected in GRB970508 should have weakened considerably during the first month if the absorption occurred in a gas cloud of size <100pc around the source. The time-dependent HI or metal absorption features provide a precise determination of the GRB redshift.

The afterglow of a cosmological Gamma-Ray Burst (GRB) should appear on the sky as a narrow emission ring of radius 3x10^{16}cm(t/day)^{5/8} which expands faster than light. After a day, the ring radius is comparable to the Einstein radius of a solar mass lens at a cosmological distance. Thus, microlensing by an intervening star can modify significantly the lightcurve and polarization signal from a GRB afterglow. We show that the achromatic amplification signal of the afterglow flux can be used to determine the impact parameter and expansion rate of the source in units of the Einstein radius of the lens, and probe the superluminal nature of the expansion. If the synchrotron emission from the afterglow photosphere originates from a set of coherent magnetic field patches, microlensing would induce polarization variability due to the transient magnification of the patches behind the lens. The microlensing interpretation of the flux and polarization data can be confirmed by a parallax experiment which would probe the amplification peak at different times. The fraction of microlensed afterglows can be used to calibrate the density parameter of stellar-mass objects in the Universe.

The damped Lyman alpha absorbers (DLAs) in quasar spectra are believed to be the progenitors of present-day disk galaxies. We examine the probability for microlensing of background quasars by stars in their DLAs. Microlensing by an individual star should magnify the continuum but not the broad emission lines of the quasars. Consequently, the equivalent width distribution of microlensed quasars would be distorted. We model a representative spiral galaxy as a closed system composed of a bulge, a disk, and a halo, and evolve the mass fraction of stars in the disk based on the observed metallicity of DLAs at high redshifts. The microlensing signatures are stronger if the halo of the galaxy is made of Massive Compact Halo Objects (MACHOs). In this case, the distortion imprinted by microlensing on the equivalent width distribution of quasar emission lines can be detected with high significance in a sample of about 10 DLAs with HI column densities N>10^21 cm^-2 and absorption redshifts z_abs<1. About a tenth of all quasars with DLAs (N>10^20 cm^-2) might show excess variability on timescales shorter than five years. A search for these signals would complement microlensing searches in local galaxies and calibrate the MACHO mass fraction in galactic halos at high redshifts.

Current microlensing searches calibrate the mass fraction of the Milky Way halo which is in the form of Massive Compact Halo Objects (MACHOs). We show that surveys like the Sloan Digital Sky Survey (SDSS) can probe the same quantity in halos of distant galaxies. Microlensing of background quasars by MACHOs in intervening galaxies would distort the equivalent width distribution of the quasar emission lines by an amplitude that depends on the projected quasar-galaxy separation. For a statistical sample of detectable at the >2sigma level out to a quasar-galaxy impact parameter of several tens of kpc, as long as extragalactic halos are made of MACHOs. Detection of this signal would test whether the MACHO fraction inferred for the Milky-Way halo is typical of other galaxies.

Gravitational lensing of an optical QSO by a spiral galaxy is often counteracted by dust obscuration, since the line-of-sight to the QSO passes close to the center of the galactic disk. The dust in the lens is likely to be correlated with neutral hydrogen, which in turn should leave a Lyman-alpha absorption signature on the QSO spectrum. We use the estimated dust-to-gas ratio of the Milky-Way galaxy as a mean and allow a spread in its values to calculate the effects of dust on lensing by low redshift spiral galaxies. Using a no-evolution model for spirals at z<1 we find (in Lambda=0 cosmologies) that the magnification bias due to lensing is stronger than dust obscuration for QSO samples with a magnitude limit B<16. The density parameter of neutral hydrogen, Omega_HI, is overestimated in such samples and is underestimated for fainter QSOs.

The microlensing monitoring programs have studied large numbers of standard light curves which seem to be due to lensing by a dark point mass. Theory predicts that many microlensing events should display significant deviations from the standard form. Lens binarity in particular is expected to be common. So far, however, only a handful of light curves exhibit evidence that the lens is a binary; all of these display dramatic deviations from the standard light curve, exhibiting pronounced multiple peaks and/or caustic crossings. Binary-lens events in which the light curve is less dramatically perturbed should also exist in the data set. Why, then, have we not detected them? The answer may lie in the fact that the perturbations, though often significant, tend to be less distinctive than those associated with caustic crossings. We present a method to determine whether a light curve is due to lensing by a binary. The method works for both gently and dramatically perturbed binary-lens light curves. Our method identifies all degenerate solutions-- i.e., all possible lensing events that might have given rise to the observed light curve. It also enables us to eliminate from consideration large ranges of possible false positive identifications associated with light curves that might mimic microlensing by a binary. This method, or a generalization of it, can also be applied to the analysis of light curves that deviate from the standard point-mass lens form because of astronomical effects other than lens binarity.