Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Pal 5, Orphan and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a St\"ackel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1 and Pal 5 streams whose combined data set yields $M(< 20\ \mathrm{kpc}) = 3.47^{+0.95}_{-1.44} \times 10^{11} \ M_{\odot}$. When including the Orphan and Helmi stream in the data set, the mass uncertainty increases to $M(< 20\ \mathrm{kpc}) = 3.12^{+5.69}_{-1.07} \times 10^{11} \ M_{\odot}$, indicative of the hidden systematic errors in fits to individual streams. These systematics are likely due to insufficiencies in our St\"ackel model of the potential and to the limited phase space explored by the data set of 4 streams, so the larger uncertainty of the combined fit is the more robust measure of the actual uncertainty in the Milky Way potential.

Several stars detected moving at velocities near to or exceeding the Galactic escape speed likely originated in the Milky Way disc. We quantitatively explore the `binary supernova scenario' hypothesis, wherein these stars are ejected at large peculiar velocities when their close, massive binary companions undergo a core-collapse supernova and the binary is disrupted. We perform an extensive suite of binary population synthesis simulations evolving massive systems to determine the assumptions and parameters which most impact the ejection rate of fast stars. In a simulation tailored to eject fast stars, we find hyper-runaway star progenitor binaries composed of a massive ($\sim$30 M$_{\odot}$) primary and a $\sim$3-4 M$_{\odot}$ companion on an orbital period that shrinks to $\lesssim$1 day prior to the core collapse following a common envelope phase. The black hole remnant formed from the primary must receive a natal kick $\gtrsim$1000 km s$^{-1}$ to disrupt the binary and eject the companion at a large velocity. We compare the fast stars produced in these simulations to a contemporary census of early-type Milky Way hyper-runaway star candidates. We find that these rare objects may be produced in sufficient number only when poorly-constrained binary evolution parameters related to the strength of post-core collapse black hole natal kicks and common envelope efficiency are adjusted to values currently unsupported -- but not excluded -- by the literature. We discuss observational implications that may constrain the existence of these putative progenitor systems.

Tidal disruption events (TDEs) are among the brightest transients in the optical, ultraviolet, and X-ray sky. These flares are set into motion when a star is torn apart by the tidal field of a massive black hole, triggering a chain of events which is -- so far -- incompletely understood. However, the disruption process has been studied extensively for almost half a century, and unlike the later stages of a TDE, our understanding of the disruption itself is reasonably well converged. In this Chapter, we review both analytical and numerical models for stellar tidal disruption. Starting with relatively simple, order-of-magnitude physics, we review models of increasing sophistication, the semi-analytic ``affine formalism,'' hydrodynamic simulations of the disruption of polytropic stars, and the most recent hydrodynamic results concerning the disruption of realistic stellar models. Our review surveys the immediate aftermath of disruption in both typical and more unusual TDEs, exploring how the fate of the tidal debris changes if one considers non-main sequence stars, deeply penetrating tidal encounters, binary star systems, and sub-parabolic orbits. The stellar tidal disruption process provides the initial conditions needed to model the formation of accretion flows around quiescent massive black holes, and in some cases may also lead to directly observable emission, for example via shock breakout, gravitational waves or runaway nuclear fusion in deeply plunging TDEs.

The Laser Interferometer Space Antenna (LISA) will enable Galactic gravitational wave (GW) astronomy by individually resolving $ > 10^4$ signals from double white dwarf (DWD) binaries throughout the Milky Way. In this work we assess for the first time the potential of LISA data to map the Galactic stellar bar and spiral arms, since GWs are unaffected by stellar crowding and dust extinction unlike optical observations of the bulge region. To achieve this goal we combine a realistic population of Galactic DWDs with a high-resolution N-Body simulation a galaxy in good agreement with the Milky Way. We then model GW signals from our synthetic DWD population and reconstruct the structure of the simulated Galaxy from mock LISA observations. Our results show that while the low signal contrast between the background disc and the spiral arms hampers our ability to characterise the spiral structure, the stellar bar will instead clearly appear in the GW map of the bulge. The bar length and bar width derived from these synthetic observations are underestimated, respectively within $1\sigma$ and at a level greater than $2\sigma$, but the resulting axis ratio agrees to well within $1\sigma$, while the viewing angle is recovered to within one degree. These are competitive constraints compared to those from electromagnetic tracers, and they are obtained with a completely independent method. We therefore foresee that the synergistic use of GWs and electromagnetic tracers will be a powerful strategy to map the bar and the bulge of the Milky Way.

Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations.

Milky Way dwarf satellites are unique objects that encode the early structure formation and therefore represent a window into the high redshift Universe. So far, their study was conducted using electromagnetic waves only. The future Laser Interpreter Space Antenna (LISA) has the potential to reveal Milky Way satellites in gravitational waves emitted by double white dwarf (DWD) binaries. We investigate gravitational wave (GW) signals detectable by LISA as a possible tool for the identification and characterisation of the Milky Way satellites. We use the binary population synthesis technique to model the population of DWDs in dwarf satellites and we assess the impact on the number of LISA detections when making changes to the total stellar mass, distance, star formation history and metallicity of satellites. We calibrate predictions for the known Milky Way satellites on their observed properties. We find that DWDs emitting at frequencies $\gtrsim 3\,$mHz can be detected in Milky Way satellites at large galactocentric distances. The number of these high frequency DWDs per satellite primarily depends on its mass, distance, age and star formation history, and only mildly depends on the other assumptions regarding their evolution such as metallicity. We find that dwarf galaxies with $M_\star>10^6\,$M$_{\odot}$ can host detectable LISA sources with a number of detections that scales linearly with the satellite's mass. We forecast that out of the known satellites, Sagittarius, Fornax, Sculptor and the Magellanic Clouds can be detected with LISA. As an all-sky survey that does not suffer from contamination and dust extinction, LISA will provide observations of the Milky Way and dwarf satellites galaxies valuable for Galactic archaeology and near-field cosmology.

We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$).

We explore the prospects for the detection of giant circumbinary exoplanets and brown dwarfs (BDs) orbiting Galactic double white dwarfs binaries (DWDs) with the Laser Interferometer Space Antenna (LISA). By assuming an occurrence rate of 50%, motivated by white dwarf pollution observations, we built a Galactic synthetic population of P-type giant exoplanets and BDs orbiting DWDs. We carried this out by injecting different sub-stellar populations, with various mass and orbital separation characteristics, into the DWD population used in the LISA mission proposal. We then performed a Fisher matrix analysis to measure how many of these three-body systems show a periodic Doppler-shifted gravitational wave perturbation detectable by LISA. We report the number of circumbinary planets (CBPs) and (BDs) that can be detected by LISA for various combinations of mass and semi-major axis distributions. We identify pessimistic and optimistic scenarios corresponding, respectively, to 3 and 83 (14 and 2218) detections of CBPs (BDs), observed during the length of the nominal LISA mission. These detections are distributed all over the Galaxy following the underlying DWD distribution, and they are biased towards DWDs with higher LISA signal-to-noise ratio and shorter orbital period. Finally, we show that if LISA were to be extended for four more years, the number of systems detected will be more than doubled in both the optimistic and pessimistic scenarios. Our results present promising prospects for the detection of post-main sequence exoplanets and BDs, showing that gravitational waves can prove the existence of these populations over the totality of the Milky Way. Detections by LISA will deepen our knowledge on the life of exoplanets subsequent to the most extreme evolution phases of their hosts, clarifying whether new phases of planetary formation take place later in the life of the stars.

Hypervelocity stars (HVSs) travel from the Galactic Centre across the dark matter halo of the Milky Way, where they are observed with velocities in excess of the Galactic escape speed. Because of their quasi-radial trajectories, they represent a unique probe of the still poorly constrained dark matter component of the Galactic potential. In this paper, we present a new method to produce such constraints. Our likelihood is based on the local HVS density obtained by back-propagating the observed phase space position and quantifies the ejection probability along the orbit. To showcase our method, we apply it to simulated Gaia samples of $\sim200$ stars in three realistic Galactic potentials with dark matter components parametrized by spheroidal NFW profiles. We find that individual HVSs exhibit a degeneracy in the scale mass-scale radius plane ($M_s-r_s$) and are able to measure only the combination $\alpha = M_s/r_s^2$. Likewise, a degeneracy is also present between $\alpha$ and the spheroidal axis-ratio $q$. In the absence of observational errors, we show the whole sample can nail down both parameters with {\it sub-per cent} precision (about $1\%$ and $0.1\%$ for $\alpha$ and $q$ respectively) with no systematic bias. This remarkable power to constrain deviations from a symmetric halo is a consequence of the Galactocentric origin of HVSs. To compare our results with other probes, we break the degeneracy in the scale parameters and impose a mass-concentration relation. The result is a competitive precision on the virial mass $M_{200}$ of about $10\%$.

In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various parameters of the IMBH-SMBH binary: IMBH mass, semimajor axis and eccentricity. We also calculate the rates of change of the BH binary orbital elements due to those stellar ejections. One of our new findings is that the ejection rate depends (although mildly) on the rotation of the stellar nucleus (its total angular momentum). We also compare the ejection velocity distribution with that produced by the Hills mechanism (stellar binary disruption) and find that the latter produces faster stars on average. Also, the IMBH mechanism produces an ejection velocity distribution which is flattened towards the BH binary plane while the Hills mechanism produces a spherically symmetric one. The results of this paper will allow us in the future to model the ejection of stars by an evolving BH binary and compare both models with \textit{Gaia} observations, for a wide variety of environments (galactic nuclei, globular clusters, the Large Magellanic Clouds, etc.).

In this work we include black hole (BH) seeding, growth and feedback into our semi-analytic galaxy formation model, Delphi. Our model now fully tracks the, accretion- and merger-driven, hierarchical assembly of the dark matter halo, baryonic and BH masses of high-redshift ($z>5$) galaxies. We use a minimal set of mass- and $z$-independent free parameters associated with star formation and BH growth (and feedback) and include suppressed BH growth in low-mass galaxies to explore a number of physical scenarios including: (i) two types of BH seeds (stellar and those from Direct Collapse BH; DCBH); (ii) the impact of reionization feedback; and (iii) the impact of instantaneous versus delayed galaxy mergers on the baryonic growth. While both reionization feedback and delayed galaxy mergers have no sensible impact on the evolving ultra-violet luminosity function, the latter limits the maximum BH masses achieved at these high-$z$. We then use this model, baselined against all available high-$z$ galaxy and BH data-sets, to predict the LISA detectability of merger events at $z > 5$. As expected, the merger rate is dominated by stellar BH mergers for all scenarios and our model predicts an expected upper limit of about 20 mergers in the case of instantaneous merging and no reionization feedback over the 4-year mission duration. Including the impact of delayed mergers and reionization feedback reduces this to about 12 events over the same observational time-scale.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

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.

The upcoming LISA mission offers the unique opportunity to study the Milky Way through gravitational wave radiation from Galactic binaries. Among the variety of Galactic gravitational wave sources, LISA is expected to individually resolve signals from $\sim 10^5$ ultra-compact double white dwarf (DWD) binaries. DWDs detected by LISA will be distributed across the Galaxy, including regions that are hardly accessible to electromagnetic observations such as the inner part of the Galactic disc, the bulge and beyond. We quantitatively show that the large number of DWD detections will allow us to use these systems as tracers of the Milky Way potential. We demonstrate that density profiles of DWDs detected by LISA may provide constraints on the scale length parameters of the baryonic components that are both accurate and precise, with statistical errors of a few percent to $10$ percent level. Furthermore, the LISA sample is found to be sufficient to disentangle between different (commonly used) disc profiles, by well covering the disc out to sufficiently large radii. Finally, up to $\sim 80$ DWDs can be detected through both electromagnetic and gravitational wave radiation. This enables multi-messenger astronomy with DWD binaries and allows one to extract their physical properties using both probes. We show that fitting the Galactic rotation curve constructed using distances inferred from gravitational waves {\it and} proper motions from optical observations yield a unique and competitive estimate of the bulge mass. Instead robust results for the stellar disc mass are contingent upon knowledge of the Dark Matter content.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

The upcoming LISA mission is the only experiment that will allow us to study the Milky Way's structure using gravitational wave signals from Galactic double white dwarfs (DWDs). The total number of expected detections exceeds $10^5$. Furthermore, up to a hundred DWDs can be simultaneously detected in both gravitational and optical radiation (e.g. with Gaia and LSST as eclipsing), making DWDs ideal sources for performing a multi-messenger tomography of the Galaxy. We show that LISA will detect DWDs everywhere, mapping also the opposite side of the Galaxy. This complete coverage will: (1) provide precise and unbiased constraints on the scale radii of the Milky Way's bulge and disc, and (2) allow us to compute the rotation curve and derive competitive estimates for the bulge and disc masses, when combining gravitational wave and optical observations.

We search for the fastest stars in the subset of stars with radial velocity measurements of the second data release (DR2) of the European Space Agency mission Gaia. Starting from the observed positions, parallaxes, proper motions, and radial velocities, we construct the distance and total velocity distribution of more than $7$ million stars in our Milky Way, deriving the full 6D phase space information in Galactocentric coordinates. These information are shared in a catalogue, publicly available at http://home.strw.leidenuniv.nl/~marchetti/research.html. To search for unbound stars, we then focus on stars with a probability greater than $50 \%$ of being unbound from the Milky Way. This cut results in a clean sample of $125$ sources with reliable astrometric parameters and radial velocities. Of these, $20$ stars have probabilities greater than 80 $\%$ of being unbound from the Galaxy. On this latter sub-sample, we perform orbit integration to characterize the stars' orbital parameter distributions. As expected given the relatively small sample size of bright stars, we find no hypervelocity star candidates, stars that are moving on orbits consistent with coming from the Galactic Centre. Instead, we find $7$ hyper-runaway star candidates, coming from the Galactic disk. Surprisingly, the remaining $13$ unbound stars cannot be traced back to the Galaxy, including two of the fastest stars (around $700$ km/s). If conformed, these may constitute the tip of the iceberg of a large extragalactic population or the extreme velocity tail of stellar streams.

Detached double white dwarf (DWD) binaries are one of the main science case for the Laser Interferometer Space Antenna (LISA). As the most numerous LISA sources, they will provide important contributions towards understanding binary evolution, Supernovae Type Ia (SNIa) formation channels and the structure of the Milky Way. So far only detection prospects for the Milky Way have been computed. In this letter we show that LISA has the potential to detect DWDs in neighboring galaxies up to the border of the Local Group. In particular, we compute quantitative estimates for the number of detections in M31. We expect between a dozen to several tens of DWDs above the nominal detection threshold, for a mission duration between 4 and 10$\,$yr. We show that extra-galactic DWDs detectable by LISA are those with the shortest orbital periods and with the highest chirp masses, that are candidates SNIa progenitor, virtually undetectable at those distances in optical. This implies that LISA will be the only instrument able to provide SNIa merger rates across the Local Group.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

Ultracompact binaries with orbital periods less than a few hours will dominate the gravitational wave signal in the mHz regime. Until recently, 10 systems were expected have a predicted gravitational wave signal strong enough to be detectable by the Laser Interferometer Space Antenna (LISA), the so-called `verification binaries'. System parameters, including distances, are needed to provide an accurate prediction of the expected gravitational wave strength to be measured by LISA. Using parallaxes from {\sl Gaia} Data Release 2 we calculate signal-to-noise ratios (SNR) for $\approx$50 verification binary candidates. We find that 11 binaries reach a SNR$\geq$20, two further binaries reaching a SNR$\geq$5 and three more systems are expected to have a SNR$\approx$5 after four years integration with LISA. For these 16 systems we present predictions of the gravitational wave amplitude ($\mathcal{A}$) and parameter uncertainties from Fisher information matrix on the amplitude ($\mathcal{A}$) and inclination ($\iota$).

The rate of tidal disruption events (TDEs) can vary by orders of magnitude depending on the environment and the mechanism that launches the stars towards the black hole's vicinity. For the largest rates, two disruptions can take place shortly one after the other in a double TDE. In this case, the two debris streams may collide with each other before falling back to the black hole resulting in an electromagnetic emission that is absent from single TDEs. We analytically evaluate the conditions for this streams collision to occur. It requires that the difference in pericenter location between the two disruptions makes up for the time delay between them. In addition, the width of the streams must compensate for the vertical offset induced by the inclination of their orbital planes. If the double TDE happens following the tidal separation of a binary, we find that the streams can collide with a probability as high as $44\%$. We validate our analytical conditions for streams collision through hydrodynamical simulations and find that the associated shocks heat the gas significantly. If photons are able to rapidly escape, a burst of radiation ensues lasting a few days with a luminosity $\sim 10^{43}\, \rm erg\, s^{-1}$, most likely in the optical band. This signal represents a precursor to the main flare of TDEs that could in particular be exploited to determine the efficiency of disc formation from the stellar debris.

Binary stars that are on close orbits around massive black holes (MBH) such as Sgr A* in the center of the Milky Way are liable to undergo tidal disruption and eject a hypervelocity star. We study the interaction between such a MBH and circular binaries for general binary orientations and penetration depths (i.e. binaries penetrate into the tidal radius around the BH). We show that for very deep penetrators, all binaries are disrupted when the binary rotation axis is roughly oriented toward the BH or it is in the opposite direction. The surviving chance becomes significant when the angle between the binary rotation axis and the BH direction is between \pi /4 and 3 \pi /4. The surviving chance is as high as $\sim$ 20$\%$ when the binary rotation axis is perpendicular to the BH direction. The angular dependence is opposite for very shallow penetrators where coplanar prograde orbits have the lowest surviving chance (or equivalently most vulnerable). We provide numerical fits to the disruption probability and energy gain at the the BH encounter as a function of the penetration depth. The latter can be simply rescaled in terms of binary masses, their initial separation and the binary-to-BH mass ratio to evaluate the ejection velocity of a binary members in various systems. We also investigate the disruption of coplanar, eccentric binaries by a MBH. It is shown that for highly eccentric binaries retrograde orbits have a significantly increased disruption probability and ejection velocities compared to the circular binaries.

Hypervelocity stars (HVSs) are amongst the fastest objects in our Milky Way. These stars are predicted to come from the Galactic center (GC) and travel along unbound orbits across the Galaxy. In the coming years, the ESA satellite Gaia will provide the most complete and accurate catalogue of the Milky Way, with full astrometric parameters for more than $1$ billion stars. In this paper, we present the expected sample size and properties (mass, magnitude, spatial, velocity distributions) of HVSs in the Gaia stellar catalogue. We build three Gaia mock catalogues of HVSs anchored to current observations, exploring different ejection mechanisms and GC stellar population properties. In all cases, we predict hundreds to thousands of HVSs with precise proper motion measurements within a few tens of kpc from us. For stars with a relative error in total proper motion below $10 \%$, the mass range extends to ~$10 M_{\odot}$ but peaks at ~$1$ $M_\odot$. The majority of Gaia HVSs will therefore probe a different mass and distance range compared to the current non-Gaia sample. In addition, a subset of a few hundreds to a few thousands of HVSs with $M$ ~ $3$ $M_\odot$ will be bright enough to have a precise measurement of the three-dimensional velocity from Gaia alone. Finally, we show that Gaia will provide more precise proper motion measurements for the current sample of HVS candidates. This will help identifying their birthplace narrowing down their ejection location, and confirming or rejecting their nature as HVSs. Overall, our forecasts are extremely encouraging in terms of quantity and quality of HVS data that can be exploited to constrain both the Milky Way potential and the GC properties.

Whilst astronomy as a science is historically founded on observations at optical wavelengths, studying the Universe in other bands has yielded remarkable discoveries, from pulsars in the radio, signatures of the Big Bang at submm wavelengths, through to high energy emission from accreting, gravitationally-compact objects and the discovery of gamma-ray bursts. Unsurprisingly, the result of combining multiple wavebands leads to an enormous increase in diagnostic power, but powerful insights can be lost when the sources studied vary on timescales shorter than the temporal separation between observations in different bands. In July 2015, the workshop "Paving the way to simultaneous multi-wavelength astronomy" was held as a concerted effort to address this at the Lorentz Center, Leiden. It was attended by 50 astronomers from diverse fields as well as the directors and staff of observatories and spaced-based missions. This community white paper has been written with the goal of disseminating the findings of that workshop by providing a concise review of the field of multi-wavelength astronomy covering a wide range of important source classes, the problems associated with their study and the solutions we believe need to be implemented for the future of observational astronomy. We hope that this paper will both stimulate further discussion and raise overall awareness within the community of the issues faced in a developing, important field.

The paucity of hypervelocity stars (HVSs) known to date has severely hampered their potential to investigate the stellar population of the Galactic Centre and the Galactic Potential. The first Gaia data release gives an opportunity to increase the current sample. The challenge is the disparity between the expected number of hypervelocity stars and that of bound background stars. We have applied a novel data mining algorithm based on machine learning techniques, an artificial neural network, to the Tycho-Gaia astrometric solution (TGAS) catalogue. With no pre-selection of data, we could exclude immediately $\sim 99 \%$ of the stars in the catalogue and find 80 candidates with more than $90\%$ predicted probability to be HVSs, based only on their position, proper motions, and parallax. We have cross-checked our findings with other spectroscopic surveys, determining radial velocities for 30 and spectroscopic distances for 5 candidates. In addition, follow-up observations have been carried out at the Isaac Newton Telescope for 22 stars, for which we obtained radial velocities and distance estimates. We discover 14 stars with a total velocity in the Galactic rest frame > 400 km/s, and 5 of these have a probability $>50\%$ of being unbound from the Milky Way. Tracing back their orbits in different Galactic potential models we find one possible unbound HVS with velocity $\sim$ 520 km/s, 5 bound HVSs, and, notably, 5 runaway stars with median velocity between 400 and 780 km/s. At the moment, uncertainties in the distance estimates and ages are too large to confirm the nature of our candidates by narrowing down their ejection location, and we wait for future Gaia releases to validate the quality of our sample. This test successfully demonstrates the feasibility of our new data mining routine.

When a star gets tidally disrupted by a supermassive black hole, its magnetic field is expected to pervade its debris. In this paper, we study this process via smoothed particle magnetohydrodynamical simulations of the disruption and early debris evolution including the stellar magnetic field. As the gas stretches into a stream, we show that the magnetic field evolution is strongly dependent on its orientation with respect to the stretching direction. In particular, an alignment of the field lines with the direction of stretching induces an increase of the magnetic energy. For disruptions happening well within the tidal radius, the star compression causes the magnetic field strength to sharply increase by an order of magnitude at the time of pericentre passage. If the disruption is partial, we find evidence for a dynamo process occurring inside the surviving core due to the formation of vortices. This causes an amplification of the magnetic field strength by a factor of $\sim 10$. However, this value represents a lower limit since it increases with numerical resolution. For an initial field strength of 1 G, the magnetic field never becomes dynamically important. Instead, the disruption of a star with a strong 1 MG magnetic field produces a debris stream within which magnetic pressure becomes similar to gas pressure a few tens of hours after disruption. If the remnant of one or multiple partial disruptions is eventually fully disrupted, its magnetic field could be large enough to magnetically power the relativistic jet detected from Swift J1644+57. Magnetized streams could also be significantly thickened by magnetic pressure when it overcomes the confining effect of self-gravity.

Double white dwarf (DWD) binaries are expected to be very common in the Milky Way, but their intrinsic faintness challenges the detection of these systems. Currently, only a few tens of detached DWDs are know. Such systems offer the best chance of extracting the physical properties that would allow us to address a wealth of outstanding questions ranging from the nature of white dwarfs, over stellar and binary evolution to mapping the Galaxy. In this paper we explore the prospects for detections of ultra-compact (with binary separations of a few solar radii or less) detached DWDs in: 1) optical radiation with Gaia and the LSST and 2) gravitational wave radiation with LISA. We show that Gaia, LSST and LISA have the potential to detect respectively around a few hundreds, a thousand, and 25 thousand DWD systems. Moreover, Gaia and LSST data will extend by respectively a factor of two and seven the guaranteed sample of binaries detected in electromagnetic and gravitational wave radiation, opening the era of multi-messenger astronomy for these sources.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

The mass assembly history of the Milky Way can inform both theory of galaxy formation and the underlying cosmological model. Thus, observational constraints on the properties of both its baryonic and dark matter contents are sought. Here we show that hypervelocity stars (HVSs) can in principle provide such constraints. We model the observed velocity distribution of HVSs, produced by tidal break-up of stellar binaries caused by Sgr A*. Considering a Galactic Centre (GC) binary population consistent with that inferred in more observationally accessible regions, a fit to current HVS data with significance level > 5\% can only be obtained if the escape velocity from the GC to 50 kpc is $V_G < 850$ km/s, regardless of the enclosed mass distribution. When a NFW matter density profile for the dark matter halo is assumed, haloes with $V_G < 850$ km/s are in agreement with predictions in the $\Lambda$CDM model and that a subset of models around $M_{200} \sim 0.5-1.5 \times 10^{12}$ solar masses and $r_s < 35$ kpc can also reproduce Galactic circular velocity data. HVS data alone cannot currently exclude potentials with $V_G > 850$ km/s. Finally, specific constraints on the halo mass from HVS data are highly dependent on the assumed baryonic mass potentials. This first attempt to simultaneously constrain GC and dark halo properties is primarily hampered by the paucity and quality of data. It nevertheless demonstrates the potential of our method, that may be fully realised with the ESA Gaia mission.

Most massive galaxies are thought to contain a supermassive black hole in their centre surrounded by a tenuous gas environment, leading to no significant emission. In these quiescent galaxies, tidal disruption events represent a powerful detection method for the central black hole. Following the disruption, the stellar debris evolves into an elongated gas stream, which partly falls back towards the disruption site and accretes onto the black hole producing a luminous flare. Using an analytical treatment, we investigate the interaction between the debris stream and the gas environment of quiescent galaxies. Although we find dynamical effects to be negligible, we demonstrate that Kelvin-Helmholtz instability can lead to the dissolution of the stream into the ambient medium before it reaches the black hole, likely dimming the associated flare. This result is robust against the presence of a typical stellar magnetic field and fast cooling within the stream. Furthermore, we find this effect to be enhanced for disruptions involving more massive black holes and/or giant stars. Consequently, although disruptions of evolved stars have been proposed as a useful probe of black holes with masses $\gtrsim 10^8 \, M_{\odot}$, we argue that the associated flares are likely less luminous than expected.

A large number of tidal disruption event (TDE) candidates have been observed recently, often differing in their observational features. Two classes appear to stand out: X-ray and optical TDEs, the latter featuring lower effective temperatures and luminosities. These differences can be explained if the radiation detected from the two categories of events originates from different locations. In practice, this location is set by the evolution of the debris stream around the black hole and by the energy dissipation associated with it. In this paper, we build an analytical model for the stream evolution, whose dynamics is determined by both magnetic stresses and shocks. Without magnetic stresses, the stream always circularizes. The ratio of the circularization timescale to the initial stream period is $t_{\rm ev}/t_{\rm min} = 8.3 (M_{\rm h}/10^6 M_{\odot})^{-5/3} \beta^{-3}$, where $M_{\rm h}$ is the black hole mass and $\beta$ is the penetration factor. If magnetic stresses are strong, they can lead to the stream ballistic accretion. The boundary between circularization and ballistic accretion corresponds to a critical magnetic stresses efficiency $v_{\rm A}/v_{\rm c} \approx 10^{-1}$, largely independent of $M_{\rm h}$ and $\beta$. However, the main effect of magnetic stresses is to accelerate the stream evolution by strengthening self-crossing shocks. Ballistic accretion therefore necessarily occurs on the stream dynamical timescale. The shock luminosity associated to energy dissipation is sub-Eddington but decays as $t^{-5/3}$ only for a slow stream evolution. Finally, we find that the stream thickness rapidly increases if the stream is unable to cool completely efficiently. A likely outcome is its fast evolution into a thick torus, or even an envelope completely surrounding the black hole.

Supermassive black holes are a key ingredient of galaxy evolution. However, their origin is still highly debated. In one of the leading formation scenarios, a black hole of $\sim100$ M$_{\odot}$ results from the collapse of the inner core of a supermassive star ($\gtrsim 10^{4-5}$ M$_{\odot}$), created by the rapid accumulation ($\gtrsim 0.1 $ M$_{\odot}$ yr$^{-1}$) of pristine gas at the centre of newly formed galaxies at $z\sim 15$. The subsequent evolution is still speculative: the remaining gas in the supermassive star can either directly plunge into the nascent black hole, or part of it can form a central accretion disc, whose luminosity sustains a surrounding, massive, and nearly hydrostatic envelope (a system called a "quasi-star"). To address this point, we consider the effect of rotation on a quasi-star, as angular momentum is inevitably transported towards the galactic nucleus by the accumulating gas. Using a model for the internal redistribution of angular momentum that qualitative matches results from simulations of rotating convective stellar envelopes, we show that quasi-stars with an envelope mass greater than a few $10^{5}$ M$_{\odot} \times (\rm black~hole~mass/100 M_{\odot})^{0.82}$ have highly sub-keplerian gas motion in their core, preventing gas circularisation outside the black hole's horizon. Less massive quasi-stars could form but last for only $\lesssim 10^4$ years before the accretion luminosity unbinds the envelope, suppressing the black hole growth. We speculate that this might eventually lead to a dual black hole seed population: (i) massive ($>10^{4}$ M$_{\odot}$) seeds formed in the most massive ($> 10^{8}$ M$_{\odot}$) and rare haloes; (ii) lighter ($\sim 10^{2}$ M$_{\odot}$) seeds to be found in less massive and therefore more common haloes.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

Appreciable star formation, and, therefore, numerous massive stars, are frequently found near supermassive black holes (SMBHs). As a result, core-collapse supernovae in these regions should also be expected. In this paper, we consider the observational consequences of predicting the fate of supernova remnants (SNRs) in the sphere of influence of quiescent SMBHs. We present these results in the context of `autarkic' nuclei, a model that describes quiescent nuclei as steady-state and self-sufficient environments where the SMBH accretes stellar winds with no appreciable inflow of material from beyond the sphere of influence. These regions have properties such as gas density that scale with the mass of the SMBH. Using predictions of the X-ray lifetimes of SNRs originating in the sphere of influence, we make estimates of the number of core collapse SNRs present at a given time. With the knowledge of lifetimes of SNRs and their association with young stars, we predict a number of core-collapse SNRs that grows from ~1 around Milky Way-like (4.3 x 10^6 Msun) SMBHs to ~100 around the highest mass (10^10 Msun) SMBHs. The presence of young SNRs will amplify the X-ray emission near quiescent SMBHs, and we show that the total core-collapse SNR emission has the potential to influence soft X-ray searches for very low-luminosity SMBHs. Our SNR lifetime estimates also allow us to predict star formation rates in these regions. Assuming a steady-state replenishment of massive stars, we estimate a star formation rate density of 2 x 10^-4 Msun/yr/pc^2 around the Milky Way SMBH, and a similar value around other SMBHs due to a weak dependence on SMBH mass. This value is consistent with currently available observations.

The potential of tidal disruption of stars to probe otherwise quiescent supermassive black holes cannot be exploited, if their dynamics is not fully understood. So far, the observational appearance of these events has been derived from analytical extrapolations of the debris dynamical properties just after disruption. By means of hydrodynamical simulations, we investigate the subsequent fallback of the stream of debris towards the black hole for stars already bound to the black hole on eccentric orbits. We demonstrate that the debris circularize due to relativistic apsidal precession which causes the stream to self-cross. The circularization timescale varies between 1 and 10 times the period of the star, being shorter for more eccentric and/or deeper encounters. This self-crossing leads to the formation of shocks that increase the thermal energy of the debris. If this thermal energy is efficiently radiated away, the debris settle in a narrow ring at the circularization radius with shock-induced luminosities of $\sim 10-10^3 \, L_{\rm Edd}$. If instead cooling is impeded, the debris form an extended torus located between the circularization radius and the semi-major axis of the star with heating rates $\sim 1-10^2 \, L_{\rm Edd}$. Extrapolating our results to parabolic orbits, we infer that circularization would occur via the same mechanism in $\sim 1$ period of the most bound debris for deeply penetrating encounters to $\sim 10$ for grazing ones. We also anticipate the same effect of the cooling efficiency on the structure of the disc with associated luminosities of $\sim 1-10 \, L_{\rm Edd}$ and heating rates of $\sim 0.1-1 \, L_{\rm Edd}$. In the latter case of inefficient cooling, we deduce a viscous timescale generally shorter than the circularization timescale. This suggests an accretion rate through the disc tracing the fallback rate, if viscosity starts acting promptly.

We perform smoothed-particle hydrodynamical simulations of the explosion of a helium star in a close binary system, and study the effects of the explosion on the companion star as well as the effect of the presence of the companion on the supernova remnant. By simulating the mechanism of the supernova from just after core bounce until the remnant shell passes the stellar companion, we are able to separate the various effects leading to the final system parameters. In the final system, we measure the mass stripping and ablation from, and the velocity kick imparted to, the companion star, as well as the structure of the supernova shell. The presence of the companion star produces a conical cavity in the expanding supernova remnant, and loss of material from the companion causes the supernova remnant to be more metal-rich on one side and more hydrogen-rich (from the companion material) around the cavity. Following the removal of mass from the companion, we study its subsequent evolution and compare it with a single star not subjected to a supernova impact.

The formation of supermassive black holes is still an outstanding question. In the quasi-star scenario, black hole seeds experience an initial super-Eddington growth, that in less than a million years may leave a $10^4-10^5$ M$_{\odot}$ black hole at the centre of a protogalaxy at $z \sim 20-10$. Super-Eddington accretion, however, may be accompanied by vigorous mass loss that can limit the amount of mass that reaches the black hole. In this paper, we critically assess the impact of radiative driven winds, launched from the surface of the massive envelopes from which the black hole accretes. Solving the full wind equations coupled with the hydrostatic structure of the envelope, we find mass outflows with rates between a few tens and $10^4$ M$_{\odot}$ yr$^{-1}$, mainly powered by advection luminosity within the outflow. We therefore confirm the claim by Dotan, Rossi & Shaviv (2011) that mass losses can severely affect the black hole seed early growth within a quasi-star. In particular, seeds with mass $>10^4$ M$_{\odot}$ can only form within mass reservoirs $ \gtrsim 10^7$ M$_{\odot}$, unless they are refilled at huge rates ($ \gtrsim 100$ M$_{\odot}$ yr$^{-1}$). This may imply that only very massive halos ($>10^9$ M$_{\odot}$) at those redshifts can harbour massive seeds. Contrary to previous claims, these winds are expected to be relatively bright ($10^{44}-10^{47}$ erg s$^{-1}$), blue ($T_{\rm eff} \sim 8000$ K) objects, that while eluding the Hubble Space Telescope, could be observed by the James Webb Space Telescope.

Observational consequences of tidal disruption of stars (TDEs) by supermassive black holes (SMBHs) can enable us to discover quiescent SMBHs, constrain their mass function, study formation and evolution of transient accretion disks and jet formation. A couple of jetted TDEs have been recently claimed in hard X-rays, challenging jet models, previously applied to $\gamma$-ray bursts and active galactic nuclei. It is therefore of paramount importance to increase the current sample. In this paper, we find that the best strategy is not to use up-coming X-ray instruments alone, which will yield between several (e-Rosita) and a couple of hundreds (Einstein Probe) events per year below redshift one. We rather claim that a more efficient TDE hunter will be the Square Kilometer Array (SKA) operating {\it in survey mode} at 1.4 GHz. It may detect up to several hundreds of events per year below $z \sim 2.5$ with a peak rate of a few tens per year at $z\approx 0.5$. Therefore, even if the jet production efficiency is {\it not } $100\%$ as assumed here, the predicted rates should be large enough to allow for statistical studies. The characteristic TDE decay of $t^{-5/3}$, however, is not seen in radio, whose flux is quite featureless. {\it Identification} therefore requires localization and prompt repointing by higher energy instruments. If radio candidates would be repointed within a day by future X-ray observatories (e.g. Athena and LOFT-like missions), it will be possible to detect up to $\approx 400$ X-ray counterparts, almost up to redshift $2$. The shortcome is that only for redshift below $\approx 0.4$ the trigger times will be less than 10 days from the explosion. In this regard the X-ray surveys are better suited to probe the beginning of the flare, and are therefore complementary to SKA.

Observational consequences of the tidal disruption of stars by supermassive black holes (SMBHs) can enable us to discover quiescent SMBHs and constrain their mass function. Moreover, observing jetted TDEs (from previously non-active galaxies) provides us with a new means of studying the early phases of jet formation and evolution in an otherwise "pristine" environment. Although several (tens) TDEs have been discovered since 1999, only two jetted TDEs have been recently discovered in hard X-rays, and only one, Swift J1644+57, has a precise localization which further supports the TDE interpretation. These events alone are not sufficient to address those science issues, which require a substantial increase of the current sample. Despite the way they were discovered, the highest discovery potential for {\em jetted} TDEs is not held by current and up-coming X-ray instruments, which will yield only a few to a few tens events per year. In fact, the best strategy is to use the Square Kilometer Array to detect TDEs and trigger multi-wavelength follow-ups, yielding hundreds candidates per year, up to $z \sim 2$. Radio and X-ray synergy, however, can in principle constrain important quantities such as the absolute rate of jetted TDEs, their jet power, bulk Lorentz factor, the black hole mass function, and perhaps discover massive black holes (MBH) with $<10^{5} M_{\odot}$. Finally, when comparing SKA results with information from optical surveys like LSST, one can more directly constrain the efficiency of jet production.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of jetted tidal disruption events. For a summary, we refer to the paper.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of gamma-ray bursts. For a summary, we refer to the paper.

There is mounting observational evidence that most galactic nuclei host both supermassive black holes (SMBHs) and young populations of stars. With an abundance of massive stars, core-collapse supernovae are expected in SMBH spheres of influence. We develop a novel numerical method, based on the Kompaneets approximation, to trace supernova remnant (SNR) evolution in these hostile environments, where radial gas gradients and SMBH tides are present. We trace the adiabatic evolution of the SNR shock until 50% of the remnant is either in the radiative phase or is slowed down below the SMBH Keplerian velocity and is sheared apart. In this way, we obtain shapes and lifetimes of SNRs as a function of the explosion distance from the SMBH, the gas density profile and the SMBH mass. As an application, we focus here exclusively on quiescent SMBHs, because their light may not hamper detections of SNRs and because we can take advantage of the unsurpassed detailed observations of our Galactic Centre. Assuming that properties such as gas and stellar content scale appropriately with the SMBH mass, we study SNR evolution around other quiescent SMBHs. We find that, for SMBH masses over ~10^7 Msun, tidal disruption of SNRs can occur at less than 10^4 yr, leading to a shortened X-ray emitting adiabatic phase, and to no radiative phase. On the other hand, only modest disruption is expected in our Galactic Centre for SNRs in their X-ray stage. This is in accordance with estimates of the lifetime of the Sgr A East SNR, which leads us to expect one supernova per 10^4 yr in the sphere of influence of Sgr A*.

We consider the process of stellar binaries tidally disrupted by a supermassive black hole. For highly eccentric orbits, as one star is ejected from the three-body system, the companion remains bound to the black hole. Hypervelocity stars (HVSs) observed in the Galactic halo and S-stars observed orbiting the central black hole may originate from such mechanism. In this paper, we predict the velocity distribution of the ejected stars of a given mass, after they have travelled out of the Galactic potential. We use both analytical methods and Monte Carlo simulations. We find that each part of the velocity distribution encodes different information. At low velocities < 800 km/s, the Galactic Potential shapes universally the observed distribution, which rises towards a peak, related to the Galactic escape velocity. Beyond the peak, the velocity distribution depends on binary mass and separation distributions. Finally, the finite star life introduces a break related to their mass. A qualitative comparison of our models with current observations shows the great potential of HVSs to constrain bulge and Galactic properties. Standard choices for parameter distributions predict velocities below and above ~800 km/s with equal probability, while none are observed beyond ~700 km/s and the current detections are more clustered at low velocities 300-400 km/s. These features may indicate that the separation distribution of binaries that reach the tidal sphere is not flat in logarithmic space, as observed in more local massive binaries, but has more power towards larger separations, enhancing smaller velocities. In addition, the binary formation/evolution process or the injection mechanism might also induce a cut-off at ~ 10 solar radii in the separation distribution.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We present the results of a semianalytical model that evolves the masses and spins of massive black holes together with the properties of their host galaxies along the cosmic history. As a consistency check, our model broadly reproduces a number of observations, e.g. the cosmic star formation history, the black hole mass and luminosity function and the galaxy mass function at low redshift, the black hole to bulge mass relation, and the morphological distribution at low redshift. For the first time in a semianalytical investigation, we relax the simplifying assumptions of perfect coherency or perfect isotropy of the gas fueling the black holes. The dynamics of gas is instead linked to the morphological properties of the host galaxies, resulting in different spin distributions for black holes hosted in different galaxy types. We compare our results with the observed sample of spin measurements obtained through broad K$\alpha$ iron line fitting. The observational data disfavor both accretion along a fixed direction and isotropic fueling. Conversely, when the properties of the accretion flow are anchored to the kinematics of the host galaxy, we obtain a good match between theoretical expectations and observations. A mixture of coherent accretion and phases of activity in which the gas dynamics is similar to that of the stars in bulges (i.e., with a significant velocity dispersion superimposed to a net rotation) best describes the data, adding further evidence in support to the coevolution of massive black holes and their hosts.

We present a variability study of a sample of bright gamma-ray (30 MeV -- 50 GeV) sources. This sample is an extension of the first AGILE catalogue of gamma-ray sources (1AGL), obtained using the complete set of AGILE observations in pointing mode performed during a 2.3 year period from July 9, 2007 until October 30, 2009. The dataset of AGILE pointed observations covers a long time interval and its gamma-ray data archive is useful for monitoring studies of medium-to-high brightness gamma-ray sources. In the analysis reported here, we used data obtained with an improved event filter that covers a wider field of view, on a much larger (about 27.5 months) dataset, integrating data on observation block time scales, which mostly range from a few days to thirty days. The data processing resulted in a better characterized source list than 1AGL was, and includes 54 sources, 7 of which are new high galactic latitude (|BII| >= 5) sources, 8 are new sources on the galactic plane, and 20 sources from the previous catalogue with revised positions. Eight 1AGL sources (2 high-latitude and 6 on the galactic plane) were not detected in the final processing either because of low OB exposure and/or due to their position in complex galactic regions. We report the results in a catalogue of all the detections obtained in each single OB, including the variability results for each of these sources. In particular, we found that 12 sources out of 42 or 11 out of 53 are variable, depending on the variability index used, where 42 and 53 are the number of sources for which these indices could be calculated. Seven of the 11 variable sources are blazars, the others are Crab pulsar+nebula, LS I +61{\deg}303, Cyg X-3, and 1AGLR J2021+4030.

Although the internal structure of white dwarfs is considered to be generally well understood, the source and entity of viscosity is still very uncertain. We propose here to study white dwarf viscous properties using short period (< 1 hr), detached white dwarf binaries, such as the newly discovered ~12.8 min system. These binaries are wide enough that mass transfer has not yet started but close enough that the least massive component is subject to a measurable tidal deformation. The associated tidal torque transfers orbital energy, which is partially converted into heat by the action of viscosity within the deformed star. As a consequence, its outer non-degenerate layers expand, and the star puffs up. We self-consistently calculate the fractional change in radius, and the degree of asynchronism (ratio of stellar to orbital spin) as a function of the viscous time. Specializing our calculations to J0651, we find that the discrepancy between the measured radius of the secondary star and He white dwarf model predictions can be interpreted as tidal inflation if the viscous timescale is either ~2 10^5 yr or ~10^4 yr. Such values point to a non-microscopic viscosity, possibly given by tidally induced turbulence, or by magnetic field stresses with a magnetic field strength of 10-100 Gauss. Fortunately, these two timescales produce very different degree of asynchronism, with the shortest one, bringing the system much closer to synchronisation. A measurement of the stellar spin can thus univocally determined the mean viscosity. Extrapolating the secondary's radial expansion, we predict that the star will fill is Roche lobe at a separation which is 1.2-1.3 smaller than the current one. Applying this method to a future sample of systems can allow us to learn whether viscosity changes with mass and/or nuclear composition.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

As we observe in the moon-earth system, tidal interactions in binary systems can lead to angular momentum exchange. The presence of viscosity is generally regarded as the condition for such transfer to happen. In this paper, we show how the orbital evolution can cause a persistent torque between the binary components, even for inviscid bodies. This preferentially occurs at the final stage of coalescence of compact binaries, when the orbit shrinks successively by gravitational waves and plunging on a timescale shorter than the viscous timescale. The total orbital energy transferred to the secondary by this torque is ~0.01 of its binding energy. We further show that this persistent torque induces a differentially rotating quadrupole perturbation. Specializing to the case of a secondary neutron star, we find that this non equilibrium state has an associated free energy of 10^47-10^48 erg, just prior to coalescence. This energy is likely stored in internal fluid motions, with a sizable amount of differential rotation. By tapping this free energy reservoir, a preexisting weak magnetic field could be amplified up to a strength of ~10^15 Gauss. Such a dynamically driven tidal torque can thus recycle an old neutron star into a highly magnetized neutron star, with possible observational consequences at merger.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

We study the tidal disruption of binaries by a massive point mass (e.g. the black hole at the Galactic center), and we discuss how the ejection and capture preference between unequal-mass binary members depends on which orbit they approach the massive object. We show that the restricted three-body approximation provides a simple and clear description of the dynamics. The orbit of a binary with mass m around a massive object M should be almost parabolic with an eccentricity |1-e| < (m/M)^{1/3} << 1 for a member to be captured, while the other is ejected. Indeed, the energy change of the members obtained for a parabolic orbit can be used to describe non-parabolic cases. If a binary has an encounter velocity much larger than (M/m)^{1/3} times the binary rotation velocity, it would be abruptly disrupted, and the energy change at the encounter can be evaluated in a simple disruption model. We evaluate the probability distributions for the ejection and capture of circular binary members and for the final energies. In principle, for any hyperbolic (elliptic) orbit, the heavier member has more chance to be ejected (captured), because it carries a larger fraction of the orbital energy. However, if the orbital energy is close to zero, the difference between the two members becomes small, and there is practically no ejection and capture preference. The preference becomes significant when the orbital energy is comparable to the typical energy change at the encounter. We discuss its implications to hypervelocity stars and irregular satellites around giant planets.

We present results of simulations aimed at tracing the formation of nuclear star clusters (NCs) and black hole (BH) seeds, in a cosmological context. We focus on two mechanisms for the formation of BHs at high redshifts: as end-products of (1) Population III stars in metal free halos, and of (2) runaway stellar collisions in metal poor NCs. Our model tracks the chemical, radiative and mechanical feedback of stars on the baryonic component of the evolving halos. This procedure allows us to evaluate when and where the conditions for BH formation are met, and to trace the emergence of BH seeds arising from the dynamical channel, in a cosmological context. BHs start to appear already at z~30 as remnants of Population III stars. The efficiency of this mechanism begins decreasing once feedbacks become increasingly important. Around redshift z~15, BHs mostly form in the centre of mildly metal enriched halos inside dense NCs. The seed BHs that form along the two pathways have at birth a mass around 100-1000M\odot. The occupation fraction of BHs is a function of both halo mass and mass growth rate: at a given z, heavier and faster growing halos have a higher chance to form a native BH, or to acquire an inherited BH via merging of another system. With decreasing z, the probability of finding a BH shifts toward progressively higher mass halo intervals. This is due to the fact that, at later cosmic times, low mass systems rarely form a seed, and already formed BHs are deposited into larger mass systems due to hierarchical mergers. Our model predict that at z=0, all halos above 10^11M\odot should host a BH (in agreement with observational results), most probably inherited during their lifetime. Halos less massive then 10^9M\odot have a higher probability to host a native BH, but their occupation fraction decreases below 10%.

The tidal disruption of a solar mass star around a supermassive black hole has been extensively studied analytically and numerically. In these events the star develops into an elongated banana-shaped structure. After completing an eccentric orbit, the bound debris fall onto the black hole, forming an accretion disk and emitting radiation. The same process may occur on planetary scales, if a minor body orbits too close to its star. In the Solar System, comets fall directly onto our Sun or onto planets. If the star is a compact object, the minor body can become tidally disrupted. Indeed, one of the first mechanisms invoked to produce strong gamma-ray emission involved accretion of comets onto neutron stars in our Galaxy. Here we report that the peculiarities of the `Christmas' gamma-ray burst (GRB 101225A) can be explained by a tidal disruption event of a minor body around a Galactic isolated neutron star. This result would indicate either that minor bodies can be captured by compact stellar remnants more frequently than it occurs in our Solar System or that minor body formation is relatively easy around millisecond radio pulsars. A peculiar Supernova associated to a GRB may provide an alternative explanation.

Unambiguous detection of the tidal disruption of a star would allow an assessment of the presence and masses of supermassive black holes in quiescent galaxies. It would also provide invaluable information on bulge scale stellar processes (such as two-body relaxation) via the rate at which stars are injected into the tidal sphere of influence of the black holes. This rate, in turn, is essential to predict gravitational radiation emission by compact object inspirals. The signature of a tidal disruption event is thought to be a fallback rate for the stellar debris onto the black hole that decreases as $t^{-5/3}$. This mass flux is often assumed to yield a luminous signal that decreases in time at the same rate. In this paper, we calculate the monochromatic lightcurves arising from such an accretion event. Differently from previous studies, we adopt a more realistic description of the fallback rate and of the super-Eddigton accretion physics. We also provide simultaneous lightcurves in optical, UV and X-rays. We show that, after a few months, optical and UV lightcurves scale as $t^{-5/12}$, and are thus substantially flatter than the $t^{-5/3}$ behaviour, which is a prerogative of the bolometric lightcurve, only. At earlier times and for black hole masses $< 10^7 M_{\sun}$, the wind emission dominates: after reaching a peak of $10^{41}-10^{43}$ erg/s at roughly a month, the lightcurve decreases steeply as $\sim t^{-2.6}$, until the disc contribution takes over. The X-ray band, instead, is the best place to detect the $t^{-5/3}$ "smoking gun" behaviour, although it is displayed only for roughly a year, before the emission steepens exponentially.

We consider a scenario where supermassive black holes form through direct accumulation of gas at the centre of proto-galaxies. In the first stage, the accumulated gas forms a super-massive star whose core collapses when the nuclear fuel is exhausted, forming a black hole of $M_{\rm BH} \approx 100 M_{\sun}$. As the black hole starts accreting, it inflates the surrounding dense gas into an almost hydrostatic self-gravitating envelope, with at least 10-100 times the mass of the hole. We find that these "quasistars" suffer extremely high rates of mass loss through winds from their envelopes, in analogy to very massive stars such as eta-Carinae. Only for envelope masses greater than 2.8 \times 10^{5} (M_{\rm BH}/100 M_{\sun})^{9/11} is the envelope evaporation time-scale longer than the accretion time-scale of the black hole. This relation thus constitutes a "threshold growth line" above which quasistars can grow their internal black holes. Accretion rates can be 10 to 100 times the Eddington rate. The quasistars born in this "growth region" with 10^7-10^8 M_{\sun} can grow black holes with masses between 10^4$ to $10^5 M_{\sun}, before crossing the threshold growth line and dispersing their envelopes in less than $10^4$ yr. This scenario therefore predicts that massive black hole seeds can be found only in dark matter halos with total masses larger than about 10^9 M_{\sun}, which can provide sufficiently high accretion rates to form such massive quasistars.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

Pulsars are known to power winds of relativistic particles that can produce bright nebulae by interacting with the surrounding medium. These pulsar wind nebulae (PWNe) are observed in the radio, optical, x-rays and, in some cases, also at TeV energies, but the lack of information in the gamma-ray band prevents from drawing a comprehensive multiwavelength picture of their phenomenology and emission mechanisms. Using data from the AGILE satellite, we detected the Vela pulsar wind nebula in the energy range from 100 MeV to 3 GeV. This result constrains the particle population responsible for the GeV emission, probing multivavelength PWN models, and establishes a class of gamma-ray emitters that could account for a fraction of the unidentified Galactic gamma-ray sources.

The study of relativistic particle acceleration is a major topic of high-energy astrophysics. It is well known that massive black holes in active galaxies can release a substantial fraction of their accretion power into energetic particles, producing gamma-rays and relativistic jets. Galactic microquasars (hosting a compact star of 1-10 solar masses which accretes matter from a binary companion) also produce relativistic jets. However, no direct evidence of particle acceleration above GeV energies has ever been obtained in microquasar ejections, leaving open the issue of the occurrence and timing of extreme matter energization during jet formation. Here we report the detection of transient gamma-ray emission above 100 MeV from the microquasar Cygnus X-3, an exceptional X-ray binary which sporadically produces powerful radio jets. Four gamma-ray flares (each lasting 1-2 days) were detected by the AGILE satellite simultaneously with special spectral states of Cygnus X-3 during the period mid-2007/mid-2009. Our observations show that very efficient particle acceleration and gamma-ray propagation out of the inner disk of a microquasar usually occur a few days before major relativistic jet ejections. Flaring particle energies can be thousands of times larger than previously detected maximum values (with Lorentz factors of 105 and 102 for electrons and protons, respectively). We show that the transitional nature of gamma-ray flares and particle acceleration above GeV energies in Cygnus X-3 is clearly linked to special radio/X-ray states preceding strong radio flares. Thus gamma-rays provide unique insight into the nature of physical processes in microquasars.

Motivated by detections of hypervelocity stars that may originate from the Galactic Center, we revist the problem of a binary disruption by a passage near a much more massive point mass. The six order of magnitude mass ratio between the Galactic Center black hole and the binary stars allows us to formulate the problem in the restricted parabolic three-body approximation. In this framework, results can be simply rescaled in terms of binary masses, its initial separation and binary-to-black hole mass ratio. Consequently, an advantage over the full three-body calculation is that a much smaller set of simulations is needed to explore the relevant parameter space. Contrary to previous claims, we show that, upon binary disruption, the lighter star does not remain preferentially bound to the black hole. In fact, it is ejected exactly in 50% of the cases. Nonetheless, lighter objects have higher ejection velocities, since the energy distribution is independent of mass. Focusing on the planar case, we provide the probability distributions for disruption of circular binaries and for the ejection energy. We show that even binaries that penetrate deeply into the tidal sphere of the black hole are not doomed to disruption, but survive in 20% of the cases. Nor do these deep encounters produce the highest ejection energies, which are instead obtained for binaries arriving to 0.1-0.5 of the tidal radius in a prograde orbit. Interestingly, such deep-reaching binaries separate widely after penetrating the tidal radius, but always approach each other again on their way out from the black hole.[shortened]

SuperAGILE is the hard X-ray monitor of the AGILE gamma ray mission, in orbit since 23$^{rd}$ April 2007. It is an imaging experiment based on a set of four independent silicon strip detectors, equipped with one-dimensional coded masks, operating in the nominal energy range 18-60 keV. The main goal of SuperAGILE is the observation of cosmic sources simultaneously with the main gamma-ray AGILE experiment, the Gamma Ray Imaging Detector (GRID). Given its $\sim$steradian-wide field of view and its $\sim$15 mCrab day-sensitivity, SuperAGILE is also well suited for the long-term monitoring of Galactic compact objects and the detection of bright transients. The SuperAGILE detector properties and design allow for a 6 arcmin angular resolution in each of the two independent orthogonal projections of the celestial coordinates. Photon by photon data are continuously available by the experiment telemetry, and are used to derive images and fluxes of individual sources, with integration times depending on the source intensity and position in the field of view. In this paper we report on the main scientific results achieved by SuperAGILE over its first two years in orbit, until April 2009.

The coalescence of supermassive black hole binaries occurs via the emission of gravitational waves, that can impart a substantial recoil to the merged black hole. We consider the energy dissipation, that results if the recoiling black hole is surrounded by a thin circumbinary disc. Our results differ significantly from those of previous investigations. We show analytically that the dominant source of energy is often potential energy, released as gas in the outer disc attempts to circularize at smaller radii. Thus, dimensional estimates, that include only the kinetic energy gained by the disc gas, underestimate the real energy loss. This underestimate can exceed an order of magnitude, if the recoil is directed close to the disc plane. We use three dimensional Smooth Particle Hydrodynamics (SPH) simulations and two dimensional finite difference simulations to verify our analytic estimates. We also compute the bolometric light curve, which is found to vary strongly depending upon the kick angle. A prompt emission signature due to this mechanism may be observable for low mass (10^6 Solar mass) black holes whose recoil velocities exceed about 1000 km/s. Emission at earlier times can mainly result from the response of the disc to the loss of mass, as the black holes merge. We derive analytically the condition for this to happen.

We present the first catalog of high-confidence gamma-ray sources detected by the AGILE satellite during observations performed from July 9, 2007 to June 30, 2008. Catalogued sources are detected by merging all the available data over the entire time period. AGILE, launched in April 2007, is an ASI mission devoted to gamma-ray observations in the 30 MeV - 50 GeV energy range, with simultaneous X-ray imaging capability in the 18-60 keV band. This catalog is based on Gamma-Ray Imaging Detector (GRID) data for energies greater than 100 MeV. For the first AGILE catalog we adopted a conservative analysis, with a high-quality event filter optimized to select gamma-ray events within the central zone of the instrument Field of View (radius of 40 degrees). This is a significance-limited (4 sigma) catalog, and it is not a complete flux-limited sample due to the non-uniform first year AGILE sky coverage. The catalog includes 47 sources, 21 of which are associated with confirmed or candidate pulsars, 13 with Blazars (7 FSRQ, 4 BL Lacs, 2 unknown type), 2 with HMXRBs, 2 with SNRs, 1 with a colliding-wind binary system, 8 with unidentified sources.

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

When double neutron star or neutron star-black hole binaries merge, the final remnant may comprise a central solar-mass black hole surrounded by a 0.01-0.1 solar masses torus. The subsequent evolution of this disc may be responsible for short gamma-ray bursts (SGRBs). A comparable amount of mass is ejected into eccentric orbits and will eventually fall back to the merger site after approximately 0.01 seconds. In this Paper, we investigate analytically the fate of the fallback matter, which may provide a luminous signal long after the disc is exhausted. We find that matter in the eccentric tail returns at a super-Eddington rate and is eventually (0.1 sec) unable to cool via neutrino emission and accrete all the way to the black hole. Therefore, contrary to previous claims, our analysis suggests that fallback matter is not an efficient source of late time accretion power and is unlikely to cause the late flaring activity observed in SGRB afterglows. The fallback matter rather forms a radiation-driven wind or a bound atmosphere. In both cases, the emitting plasma is very opaque and photons are released with a degraded energy in the X-ray band. We therefore suggest that compact binary mergers could be followed by an "X-ray renaissance", as late as several days to weeks after the merger. This might be observed by the next generation of X-ray detectors.

The Mini-Calorimeter (MCAL) instrument on-board the AGILE satellite is a non-imaging gamma-ray scintillation detector sensitive in the 300keV-100MeV energy range with a total on-axis geometrical area of 1400cm^2. Gamma-Ray Bursts (GRBs) are one of the main scientific targets of the AGILE mission and the MCAL design as an independent self-triggering detector makes it a valuable all-sky monitor for GRBs. Furthermore MCAL is one of the very few operative instruments with microsecond timing capabilities in the MeV range. In this paper the results of GRB detections with MCAL after one year of operation in space are presented and discussed. A flexible trigger logic implemented in the AGILE payload data-handling unit allows the on-board detection of GRBs. For triggered events, energy and timing information are sent to telemetry on a photon-by-photon basis, so that energy and time binning are limited by counting statistics only. When the trigger logic is not active, GRBs can be detected offline in ratemeter data, although with worse energy and time resolution. Between the end of June 2007 and June 2008 MCAL detected 51 GRBs, with a detection rate of about 1 GRB/week, plus several other events at a few milliseconds timescales. Since February 2008 the on-board trigger logic has been fully active. Comparison of MCAL detected events and data provided by other space instruments confirms the sensitivity and effective area estimations. MCAL also joined the 3rd Inter-Planetary Network, to contribute to GRB localization by means of triangulation.

Hyper-accreting discs occur in compact-object mergers and collapsars, and may power gamma-ray bursts (GRBs). We calculate the microscopic viscosity and resistivity of plasma in these discs, and discuss the implications for their global structure and evolution. In the neutrino-cooled innermost regions, the viscosity is provided mainly by mildly degenerate electrons, while the resistivity is modified from the Spitzer value due to the effects of both relativity and degeneracy. The plasma behaves as an almost ideal MHD fluid. Among the non-ideal MHD effects the Hall term is relatively the most important, while the magnetic Prandtl number, Pr (the ratio of viscosity to resistivity), is typically larger than unity: 10 < Pr < 6000. The outer radiatively inefficient regions also display high Pr. Numerical simulations of the magneto-rotational instability indicate that the saturation level and angular momentum transport efficiency may be greatly enhanced at high Pr. If this behaviour persists in the presence of a strong Hall effect we would expect that hyper-accreting discs should be strongly magnetised and highly variable. The expulsion of magnetic field that cannot be dissipated at small scales may also favour a magnetic outflow. We note the similaries between the Prandtl number in hyper-accreting discs and X-ray binary discs, which suggests that a comparison between late-time activity in GRBs and X-ray binary accretion states may be fruitful. Our results imply that the behavior of high Prandtl number MHD flows needs to be considered in studies of hyper-accreting discs.

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.

We present and discuss the results of an extensive observational campaign devoted to GRB071010A, a long-duration gamma-ray burst detected by the Swift satellite. This event was followed for almost a month in the optical/near-infrared (NIR) with various telescopes starting from about 2min after the high-energy event. Swift-XRT observations started only later at about 0.4d. The light-curve evolution allows us to single out an initial rising phase with a maximum at about 7min, possibly the afterglow onset in the context of the standard fireball model, which is then followed by a smooth decay interrupted by a sharp rebrightening at about 0.6d. The rebrightening was visible in both the optical/NIR and X-rays and can be interpreted as an episode of discrete energy injection, although various alternatives are possible. A steepening of the afterglow light curve is recorded at about 1d. The entire evolution of the optical/NIR afterglow is consistent with being achromatic. This could be one of the few identified GRB afterglows with an achromatic break in the X-ray through the optical/NIR bands. Polarimetry was also obtained at about 1d, just after the rebrightening and almost coincident with the steepening. This provided a fairly tight upper limit of 0.9% for the polarized-flux fraction.

We study the structure and evolution of "quasistars", accreting black holes embedded within massive hydrostatic gaseous envelopes. These configurations may model the early growth of supermassive black hole seeds. The accretion rate onto the black hole adjusts so that the luminosity carried by the convective envelope equals the Eddington limit for the total mass. This greatly exceeds the Eddington limit for the black hole mass alone, leading to rapid growth of the black hole. We use analytic models and numerical stellar structure calculations to study the structure and evolution of quasistars. We derive analytically the scaling of the photospheric temperature with the black hole mass and envelope mass, and show that it decreases with time as the black hole mass increases. Once the photospheric temperature becomes lower than 10,000 K, the photospheric opacity drops precipitously and the photospheric temperature hits a limiting value, analogous to the Hayashi track for red giants and protostars, below which no hydrostatic solution for the convective envelope exists. For metal-free (Population III) opacities this limiting temperature is approximately 4000 K. After a quasistar reaches this limiting temperature, the envelope is rapidly dispersed by radiation pressure. We find that black hole seeds with masses between 1000 and 10000 solar masses could form via this mechanism in less than a few Myr.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

Measurements at 100 TeV and above are an important goal for the next generation of high energy gamma-ray astronomy experiments to solve the still open problem of the origin of galactic cosmic rays. The most natural experimental solution to detect very low radiation fluxes is provided by the Extensive Air Shower (EAS) arrays. They benefit from a close to 90% duty cycle and a very large field of view (about 2 sr), but the sensitivity is limited by their angular resolution and their poor cosmic ray background discrimination. Above 10 TeV the standard technique for rejecting the hadronic background consists in looking for "muon-poor" showers. In this paper we discuss the capability of a large muon detector (A=2500 m2) operated with an EAS array at very high altitude (>4000 m a.s.l.) to detect gamma-ray fluxes around 100 TeV. Simulation-based estimates of energy ranges and sensitivities are presented.

The ARGO-YBJ experiment is a full coverage EAS-array installed at the YangBaJing Cosmic Ray Laboratory (4300 m a.s.l., Tibet, P.R. China). We present the results on the angular resolution measured with different methods with the full central carpet. The comparison of experimental results with MC simulations is discussed.

We study the spectrum of the prompt emission and the X-ray and optical afterglow fluxes of 54 X-Ray Rich Gamma Ray Burst (XRRs) and X-Ray Flashes (XRFs), observed by BeppoSax and HETE-2. A comparison is then performed with classical Gamma Ray Bursts (GRBs). The goal of this paper is to investigate the nature of XRRs/XRFs, as high redshift GRBs or off-axis GRBs, analyzing both their prompt and afterglow properties. We find that the XRR/XRF spectral indexes of the Band function are similar to those of classical GRBs, whereas the peak energy is lower by a factor of 4. We study the optical and X-Ray afterglow properties of the XRRs/XRFs; in particular we analyze the XRR 011030 afterglow. We find that the X-ray and optical flux distributions and the lightcurves of the XRRs/XRFs sample are consistent with those of classical GRBs; in particular, they show evidence of a break and no rising temporal slope. We compare these results with the afterglow predictions of the high redshift scenario, where XRFs are GRBs at higher redshift and of the off-axis scenario, where the observed differences are due to viewing angle effects. In this last framework, we consider jets with a homogeneous, a -2 power-law shaped and a Gaussian luminosity angular distribution. We find that the high redshift scenario can explain some events but not the total sample of XRRs/XRFs. The off-axis model may be consistent with our findings when a homogeneous jet is considered. However, given the uncertainties on the selection effects in our sample, a Gaussian jet viewed at small angles from the Gaussian core and a power-law shaped cannot be ruled out.

We predict the redshift distribution of long Gamma-Ray Bursts (GRBs) with Monte Carlo simulations. Our improved analysis constrains free parameters with three kinds of observation: (i) the log(N)-log(P) diagram of BATSE bursts; (ii) the peak energy distribution of bright BATSE bursts; (iii) the HETE2 fraction of X-ray rich GRBs and X-ray flashes. The statistical analysis of the Monte Carlo simulation results allow us to carefully study the impact of the uncertainties in the GRB intrinsic properties on the redshift distribution. The comparison with SWIFT data then leads to the following conclusions. The Amati relation should be intrinsic, if observationally confirmed by SWIFT. The progenitor and/or the GRB properties have to evolve to reproduce the high mean redshift of SWIFT bursts. Our results favor an evolution of the efficiency of GRB production by massive stars, that would be ~6-7 times higher at z~7 than at z~2. We finally predict around 10 GRBs detected by SWIFT at redshift z>6 for a three year mission. These may be sufficient to open a new observational window over the high redshift Universe.

Relativistic neutron-loaded outflows in gamma-ray bursts are studied at their early stages, before deceleration by a surrounding medium. The outflow has four components: radiation, electrons, protons and neutrons. The components interact with each other and exchange energy as the outflow expands. The presence of neutrons significantly changes the outflow evolution. Before neutrons decouple from protons, friction between the two components increases their temperatures by many orders of magnitude. After the decoupling, the gradual neutron decay inside the outflow has a drag effect on the protons and reduces their final Lorentz factor.

We review some open problems in the physics of afterglows, namely their polarization properties and the existence of dark/faint bursts. Polarization studies yield precious insights in the physical structure and dynamical evolution of GRB jets, revealing their magnetization properties and their energy profile. Polarimetric observations of GRB 020813 already allowed to exclude a homogeneous jet for this event. We then present observations of faint/dark bursts, showing that some of them may be obscured by dust, while others are possibly just intrinsically dim.

Polarimetry of Gamma-Ray Burst (GRB) afterglows in the last few years has been considered one of the most effective tool to probe the geometry, energetic, dynamics and the environment of GRBs. We report some of the most recent results and discuss their implications and future perspectives.

We numerically compute light and polarisation curves of gamma-ray burst afterglows for various configurations of the jet luminosity structure and for different dynamical evolutions. We especially consider the standard homogeneous ``top hat'' jet and the ``universal structured jet'' with power-law wings. We also investigate a possible more physical variation of the ``top hat'' model: the ``Gaussian jet''. The polarisation curves for the last two jet types are shown here for the first time together with the computation of X-ray and radio polarised fluxes. We show that the lightcurves of the total flux from these configurations are very similar to each other, and therefore only very high quality data could allow us to pin down the underlying jet structure. We demonstrate instead that polarisation curves are a powerful means to solve the jet structure, since the predicted behaviour of polarisation and its position angle at times around the jet break are very different if not opposite. We conclude that the afterglow polarisation measurements provide clear footprints of any outflow energy distribution (unlike the lightcurves of the total flux) and the joint analysis of the total and polarised flux should reveal GRBs jet structure.

Fireballs with high energy per baryon rest mass (> 400) contain a relatively slow neutron component. We show here that in this situation the thermal history of fireballs is very different from the standard adiabatic cooling.

The polarization curve of GRB 020813 is discussed and compared to different models for the structure, evolution and magnetisation properties of the jet and the interstellar medium onto which the fireball impacts. GRB 020813 is best suited for this kind of analysis for the smoothness of its afterglow light curve, ensuring the applicability of current models. The polarization dataset allows us to rule out the standard GRB jet, in which the energy and Lorentz factor have a well defined value inside the jet opening angle and the magnetic field is generated at the shock front. We explore alternative models finding that a structured jet or a jet with a toroidal component of the magnetic field can fit equally well the polarization curve. Stronger conclusions cannot be drawn due to the incomplete sampling of the polarization curve. A more dense sampling, especially at early times, is required to pin down the structure of the jet and the geometry of its magnetic field.

Thanks to its extraordinary brightness, the X-ray afterglow of GRB030329 could be studied by XMM-Newton up to two months after the prompt Gamma-ray emission. We present the results of two XMM-Newton observations performed on May 5 and 29, as well as an analysis of the Rossi-XTE data of the early part of the afterglow, discussing in particular the stability of the X-ray spectrum and presenting upper limits on the presence of X-ray emission lines.

The recent claim by Coburn & Boggs to have detected a very high degree of linear polarization in the prompt emission of GRB 021206 has stimulated interest in how much polarization could arise in gamma-ray bursts from synchrotron emission. Alternatively, as Shaviv & Dar have shown, GRB polarization could be produced by inverse Compton scattering in the point-source limit. We discuss polarization from a fireball that upscatters a soft radiation field. We show that, after the proper angular integration, the residual polarization can be large, in some cases approaching the point-source limit. We discuss the probability of realizing the geometrical conditions in which a large polarization is obtained showing that, for a particularly bright burst as GRB 021206, the detection of polarization at the first attempt in the Compton drag scenario is not unlikely.

PICsIT (Pixellated Imaging CaeSium Iodide Telescope) is the high energy detector of the IBIS telescope on-board the INTEGRAL satellite. PICsIT operates in the gamma-ray energy range between 175 keV and 10 MeV, with a typical energy resolution of 10% at 1 MeV, and an angular resolution of 12 arcmin within a \~100 square degree field of view, with the possibility to locate intense point sources in the MeV region at the few arcmin level. PICsIT is based upon a modular array of 4096 independent CsI(Tl) pixels, ~0.70 cm^2 in cross-section and 3 cm thick. In this work, the PICsIT on-board data handling and science operative modes are described. This work presents the in-flight performances in terms of background count spectra, sensitivity limit, and imaging capabilities.

The determination of the arrival direction of gamma-induced air showers in the ARGO-YBJ experiment has been investigated using different algorithms in the reconstruction procedure. The calculation has been performed, as a function of pad multiplicity, for different primary energies. The performance of a conical correction to the shower front with a suitable fixed cone slope is discussed.

In this paper we present a procedure able to identify and reject showers with the core outside the ARGO-YBJ carpet boundaries. The efficiency of this method is investigated for different primary energies and fiducial areas. A comparison of the results for gamma and proton induced showers is also presented.

We report on XMM-Newton and Rossi-XTE observations of the bright (fluence $\sim$ 10$^{-4}$ erg cm$^{-2}$) and nearby (z=0.1685) Gamma-Ray Burst GRB030329 associated to SN2003dh. The first Rossi-XTE observation, 5 hours after the burst, shows a flux decreasing with time as a power law with index 0.9$\pm$0.3. Such a decay law is only marginally consistent with a further Rossi-XTE measurement (at t-t$_{GRB}\sim$30 hr). Late time observations of this bright afterglow at X-ray wavelengths have the advantage, compared to optical observations, of not being affected by contributions from the supernova and host galaxy. A first XMM-Newton observation, at t-t$_{GRB}\sim$37 days, shows a flux of 4$\times10^{-14}$ erg cm$^{-2}$ s$^{-1}$ (0.2-10 keV). The spectrum is a power law with photon index $\Gamma$=1.9 and absorption $<2.5\times10^{20}$ cm$^{-2}$, consistent with the Galactic value. A further XMM-Newton pointing at t-t$_{GRB}\sim$61 days shows a flux fainter by a factor $\sim$2. The combined Rossi-XTE and XMM-Newton measurements require a break at t$\sim$0.5 days in the afterglow decay, with a power law index increasing from 0.9 to 1.9, similar to what is observed in the early part of the optical afterglow. The extrapolation of the XMM-Newton spectra to optical frequencies lies a factor of $\sim10$ below simultaneous measurements. This is likely due to the presence of SN2003dh.

In any EAS array, the rejection of events with shower cores outside the detector boundaries is of great importance. A large difference between the true and the reconstructed shower core positions may lead to a systematic miscalculation of some shower characteristics. Moreover, an accurate determination of the shower core position for selected internal events is important to reconstruct the primary direction using conical fits to the shower front, improving the detector angular resolution, or to performe an efficient gamma/hadron discrimination. In this paper we present a procedure able to identify and reject showers with cores outside the ARGO-YBJ carpet boundaries. A comparison of the results for gamma and proton induced showers is reported.

We present a model for the early optical afterglow of GRB 021004. This burst had one of the earliest detected optical afterglows, allowing a dense optical sampling. The lightcurve was peculiar, since bright bumps were superimposed to the regular power-law decay observed in many other events. We show that, given their time scale and shape, the bumps are due to the interaction of the fireball with moderate enhancements in the ambient medium density. The enhancements have a density contrast of order 10, modifying only slightly the dynamics of the fireball, which therefore surfs on them rather than colliding onto them. A relativistic reverse shock does not develop. Since the interaction takes place when the fireball is still hyper-relativistic it is not possible to understand if the overdensities are localized in clumps or are spherically symmetric around the GRB progenitor. The monotonic decrease of the contrast of successive rebrightenings suggests however the presence of clumps in an average uniform environment. Such an interpretation, complemented by the detection of several high velocity absorption systems in the optical spectrum, strongly suggests that GRB 021004 exploded within the remnant of a previous explosion.

In models for gamma ray burst afterglows, it is normally assumed that the external shock strongly amplifies the magnetic field and that this field maintains a steady value throughout the shocked region. We discuss the effects of modifying this (probably simplistic) assumption by allowing for a short decay time. The observations are incompatible with a post-shock field that decays too rapidly. However if the field pervades only a few percent of the total thickness of the shocked shell (and the electrons undergo only inverse Compton losses in the remainder) the model could be compatible with data. This would suggest a strong dependence of model parameters on the uncertainties of shock physics, therefore calling for independent external density estimates. We claim that afterglow emission should be instead seen as a laboratory where we can understand relativistic shock physics. The model we propose here together with afterglow observations in all wavebands could help to pin down the field structure.

It has been recently proposed that current burst and afterglow observations can be accounted for by a top hat jet structure as well as by a more universal structured outflow, in which the energy carried by the flow is maximum at the jet spine, decreasing outwards as a power-law with the off-axis angle. To pin down the GRBs jet structure is important in order to quantify the total energy released by the unknown engine and put constrains on its nature. In paper we will review the most popular GRBs jet structures, focusing on the predicted observed properties that could discriminate between them. In particular the behaviour of linearly polarised emission during the afterglow phase seems very promising, since the degree and temporal dependence are predicted to be very different from model to model. In addition, polarization measurements and detections are accumulating fast enouqgh to lead soon to a comprehension of magnetic field and GRBs jet configurations.

The emission features observed in the X-ray afterglow of Gamma Ray Bursts are extremely powerful. Since they last at least for several hours, they imply energies of the order of 1e49 ergs. This in turn implies that the energy contained in the illuminating continuum thought to be responsible of the line production must exceed 1e51 ergs. This is a strong lower limit to the energy reservoir of Gamma Ray Bursts, which is independent of collimation and beaming, and bears important consequences on the possible collimation of the fireball radiation and the density of the medium surrounding the burst.

Gamma ray bursts are often modelled as jet-like outflows directed towards the observer; the cone angle of the jet is then commonly inferred from the time at which there is a steepening in the power-law decay of the afterglow. We consider an alternative model in which the jet has a beam pattern where the luminosity per unit solid angle (and perhaps also the initial Lorentz factor) decreases smoothly away from the axis, rather than having a well-defined cone angle within which the flow is uniform. We show that the break in the afterglow light curve then occurs at a time that depends on the viewing angle. Instead of implying a range of intrinsically different jets - some very narrow, and others with similar power spread over a wider cone - the data on afterglow breaks could be consistent with a standardized jet, viewed from different angles. We discuss the implication of this model for the luminosity function.

The imager IBIS on board INTEGRAL is composed of two detector planes: ISGRI (15 keV - 1 MeV) and PICsIT (170 keV - 10 MeV). Here we describe the data structure of PICsIT and the instrument specific software, including the standard scientific analysis. Some examples of shadowgrams and spectra are shown, by using data collected during the End-to-End Tests of January 2002.

The VST (VLT Survey Telescope) is a 2.6 m class Alt-Az telescope to be installed at Mount Paranal in the Atacama desert, Chile, in the European Southern Observatory (ESO) site. The VST is a wide-field imaging facility planned to supply databases for the ESO Very Large Telescope (VLT) science and carry out stand-alone observations in the UV to I spectral range. This paper will focus on the distributed control system of active optics based on CAN bus and PIC microcontrollers. Both axial and radial pads of the primary mirror will be equipped by astatic lever supports controlled by microcontroller units. The same CAN bus + microcontroller boards approach will be used for the temperature acquisition modules.