Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface to the expansion of its remnant through the surrounding nebula by 3D MHD simulations. The model considers the radiative cooling, the deviations from equilibrium of ionization, the deviation from temperature-equilibration between electrons and ions, and a plausible configuration of the pre-SN ambient magnetic field. From the simulations, we synthesize both thermal X-ray and non-thermal radio emission and compare the model results with observations. Results. The presence of an ambient magnetic field reduces the erosion and fragmentation of the dense equatorial ring after the impact of the SN blast wave. As a result, the ring survives the passage of the blast, at least, during the time covered by the simulations (40 years). Our model is able to reproduce the morphology and lightcurves of SN 1987A in both X-ray and radio bands. The model reproduces the observed radio emission if the flux originating from the reverse shock is heavily suppressed. In this case, the radio emission originates mostly from the forward shock traveling through the H II region. Possible mechanisms for the suppression of emission from the reverse shock are investigated. We find that synchrotron self-absorption and free-free absorption have negligible effects on the emission during the interaction with the nebula. We suggest that the emission from the reverse shock at radio frequencies might be limited by highly magnetized ejecta.

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH.

About half of the expected total baryon budget in the local Universe is `missing'. Hydrodynamical simulations suggest that most of the missing baryons are located in a mildly overdense, warm-hot intergalactic medium (WHIM), which is difficult to be detected at most wavelengths. In this paper we explore multi-wavelength synthetic observations of a massive galaxy cluster developed in a full Eulerian-AMR cosmological simulation. A novel numerical procedure is applied on the outputs of the simulation, which are post-processed with a full-radiative transfer code that allows to compute the change of the intensity at any frequency along the null-geodesic of photons. We compare the emission from the whole inter-galactic medium (IGM) and from the WHIM component (defined as the gas with a temperature in the range $10^5-10^7$ K) at three observational bands associated to thermal X-rays, thermal and kinematic Sunyaev-Zel'dovich effect, and radio emission. The synthetic maps produced by this procedure could be directly compared with existing observational maps and could be used as a guide for future observations with forthcoming instruments. The analysis of the different emissions associated to a high-resolution galaxy cluster is in broad agreement with previous simulated and observational estimates of both gas components.

Internal shocks occurring in blazars may accelerate both thermal and non-thermal electrons. In this paper we examine the consequences that such a hybrid (thermal/non-thermal) EED has on the spectrum of blazars. Since the thermal component of the EED may extend to very low energies. We replace the standard synchrotron process by the more general magneto-bremsstrahlung (MBS). Significant differences in the energy flux appear at low radio frequencies when considering MBS instead of the standard synchrotron emission. A drop in the spectrum appears in the all the radio band and a prominent valley between the infrared and soft X-rays bands when a hybrid EED is considered, instead of a power-law EED. In the $\gamma$-ray band an EED of mostly thermal particles displays significant differences with respect to the one dominated by non-thermal particles. A thermally-dominated EED produces a synchrotron self-Compton (SSC) peak extending only up to a few MeV, and the valley separating the MBS and the SSC peaks is much deeper than if the EED is dominated by non-thermal particles. The combination of these effects modifies the Compton dominance of a blazar, suggesting that the vertical scatter in the distribution of FSRQs and BL Lac objects in the peak synchrotron frequency - Compton dominance parameter space could be attributed to different proportions of thermal/non-thermal particles in the EED of blazars. Finally, the temperature of the electrons in the shocked plasma is shown to be a degenerated quantity for different magentizations of the ejected material.

Radiatively inefficient accretion flow models have been shown to accurately account for the spectrum and luminosity observed from Sgr A* in the X-ray regime down to mm wavelengths. However, observations at a few GHz cannot be explained by thermal electrons alone but require the presence of an additional non-thermal particle population. Here, we propose a model for the origin of such a population in the accretion flow via means of a pulsar orbiting the supermassive black hole in our Galaxy. Interactions between the relativistic pulsar wind with the disc lead to the formation of a bow shock in the wind. During the pulsar's transit through the accretion disc, relativistic pairs, accelerated at the shock front, are injected into the disc. The radio-emitting particles are long-lived and remain within the disc long after the pulsar's transit. Periodic pulsar transits through the disc result in regular injection episodes of non-thermal particles. We show that for a pulsar with spin-down luminosity $L_{\rm sd} \sim 3\times 10^{35}$ erg s$^{-1}$ and a wind Lorentz factor of $\gamma_{\rm w} \sim 10^4$ a quasi-steady synchrotron emission is established with luminosities in the $1-10$ GHz range comparable to the observed one.

Dozens of stellar tidal disruption events (TDEs) have been identified at optical, UV and X-ray wavelengths. A small fraction of these, most notably Swift J1644+57, produce radio synchrotron emission, consistent with a powerful, relativistic jet shocking the surrounding circumnuclear gas. The dearth of similar non-thermal radio emission in the majority of TDEs may imply that powerful jet formation is intrinsically rare, or that the conditions in galactic nuclei are typically unfavorable for producing a detectable signal. Here we explore the latter possibility by constraining the radial profile of the gas density encountered by a TDE jet using a one-dimensional model for the circumnuclear medium which includes mass and energy input from a stellar population. Near the jet Sedov radius of 10$^{18}$ cm, we find gas densities in the range of $n_{18} \sim$ 0.1$-$1000 cm$^{-3}$ across a wide range of plausible star formation histories. Using one- and two-dimensional relativistic hydrodynamical simulations, we calculate the synchrotron radio light curves of TDE jets (as viewed both on and off-axis) across the allowed range of density profiles. We find that bright radio emission would be produced across the plausible range of nuclear gas densities by jets as powerful as Swift J1644+57, and we quantify the relationship between the radio luminosity and jet energy. We use existing radio detections and upper limits to constrain the energy distribution of TDE jets. Radio follow up observations several months to several years after the TDE candidate will strongly constrain the energetics of any relativistic flow.

Sgr A* is an ideal target to study low-luminosity accreting systems. It has been recently proposed that properties of the accretion flow around Sgr A* can be probed through its interactions with the stellar wind of nearby massive stars belonging to the S-cluster. When a star intercepts the accretion disk, the ram and thermal pressures of the disk terminate the stellar wind leading to the formation of a bow shock structure. Here, a semi-analytical model is constructed which describes the geometry of the termination shock formed in the wind. With the employment of numerical hydrodynamic simulations, this model is both verified and extended to a region prone to Kelvin-Helmholtz instabilities. Because the characteristic wind and stellar velocities are in $\sim10^{8}$ cm s$^{-1}$ range, the shocked wind may produce detectable X-rays via thermal bremsstrahlung emission. The application of this model to the pericenter passage of S2, the brightest member of the S-cluster, shows that the shocked wind produces roughly a month long X-ray flare with a peak luminosity of $L\approx 4 \times 10^{33}$ erg s$^{-1}$ for a stellar mass-loss rate, disk number density, and thermal pressure strength of $\dot{M}_{\rm w}= 10^{-7} M_\odot\, {\rm yr}^{-1}$, $n_{\rm d} = 10^{5}$ cm$^{-3}$, and $\alpha=0.1$, respectively. This peak luminosity is comparable to the quiescent X-ray emission detected from Sgr A* and is within the detection capabilities of current X-ray observatories. Its detection could constrain the density and thickness of the disk at a distance of $\sim 3000$ gravitational radii from the supermassive black hole.

In this work we focus on the technical details of the numerical simulations of the non- thermal transient Swift J1644+57, whose emission is probably produced by a two-component jet powered by a tidal disruption event. In this context we provide details of the coupling between the relativistic hydrodynamic simulations and the radiative transfer code. First, we consider the technical demands of one-dimensional simulations of a fast relativistic jet, and show to what extent (for the same physical parameters of the model) do the computed light curves depend on the numerical parameters of the different codes employed. In the second part we explain the difficulties of computing light curves from axisymmetric two dimensonal simulations and discuss a procedure that yields an acceptable tradeoff between the computational cost and the quality of the results.

High resolution Very Long Baseline Interferometry (VLBI) observations of Active Galactic Nuclei (AGN) revealed traveling and stationary or quasi-stationary radio-components in several blazar jets. The traveling ones are in general interpreted as shock waves generated by pressure perturbations injected at the jet nozzle. The stationary features can be interpreted as recollimation shocks in non-pressure matched jets if they show a quasi-symmetric bump in the spectral index distribution. In some jets there may be interactions between the two kinds of shocks. These shock--shock interactions are observable with VLBI techniques, and their signature should also be imprinted on the single--dish light curves. We performed relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets. To simulate the shock interaction we injected a perturbation at the jet nozzle once a steady-state was reached. We computed the non-thermal emission (including adiabatic and synchotron losses) resulting from the simulation. We show that the injection of perturbations in a jet can produce a bump in emission at GHz frequencies previous to the main flare, which is produced when the perturbation fills the jet in the observer's frame. The detailed analysis of our simulations and the non-thermal emission calculations show that interaction between a recollimation shock and traveling shock produce a typical and clear signature in both the single--dish light curves and in the VLBI observations: the flaring peaks are higher and delayed with respect to the evolution of a perturbation through a conical jet. This fact can allow to detect such interactions for stationary components lying outside of the region in where the losses are dominated by inverse Compton scattering.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

We model the non-thermal transient Swift J1644+57 as resulting from a relativistic jet powered by the accretion of a tidally-disrupted star onto a super-massive black hole. Accompanying synchrotron radio emission is produced by the shock interaction between the jet and the dense circumnuclear medium, similar to a gamma-ray burst afterglow. An open mystery, however, is the origin of the late-time radio rebrightening, which occurred well after the peak of the jetted X-ray emission. Here, we systematically explore several proposed explanations for this behavior by means of multi-dimensional hydrodynamic simulations coupled to a self-consistent radiative transfer calculation of the synchrotron emission. Our main conclusion is that the radio afterglow of Swift J1644+57 is not naturally explained by a jet with a one-dimensional top-hat angular structure. However, a more complex angular structure comprised of an ultra-relativistic core (Lorentz factor $\Gamma \sim 10$) surrounded by a slower ($\Gamma \sim $ 2) sheath provides a reasonable fit to the data. Such a geometry could result from the radial structure of the super-Eddington accretion flow or as the result of jet precession. The total kinetic energy of the ejecta that we infer of $\sim$ few $10^{53}\,$erg requires a highly efficient jet launching mechanism. Our jet model providing the best fit to the light curve of the on-axis event Swift J1644+57 is used to predict the radio light curves for off-axis viewing angles. Implications for the presence of relativistic jets from TDEs detected via their thermal disk emission, as well as the prospects for detecting orphan TDE afterglows with upcoming wide-field radio surveys and resolving the jet structure with long baseline interferometry, are discussed.

The internal shocks scenario in relativistic jets is used to explain the variability of the blazar emission. Recent studies have shown that the magnetic field significantly alters the shell collision dynamics, producing a variety of spectral energy distributions and light-curves patterns. However, the role played by magnetization in such emission processes is still not entirely understood. In this work we numerically solve the magnetohydodynamic evolution of the magnetized shells collision, and determine the influence of the magnetization on the observed radiation. Our procedure consists in systematically varying the shell Lorentz factor, relative velocity, and viewing angle. The calculations needed to produce the whole broadband spectral energy distributions and light-curves are computationally expensive, and are achieved using a high-performance parallel code.

Blackbody-dominated (BBD) gamma-ray bursts (GRBs) are events characterized by the absence of a typical afterglow, long durations and the presence of a significant thermal component following the prompt gamma-ray emission. GRB 101225A (the `Christmas burst') is a prototype of this class. A plausible progenitor system for it, and for the BBD-GRBs, is the merger of a neutron star (NS) and a helium core of an evolved, massive star. Using relativistic hydrodynamic simulations we model the propagation of an ultrarelativistic jet through the enviroment created by such a merger and we compute the whole radiative signature, both thermal and non-thermal, of the jet dynamical evolution. We find that the thermal emission originates from the interaction between the jet and the hydrogen envelope ejected during the NS/He merger.

The observational signature of supernova remnants (SNRs) is very complex, in terms of both their geometrical shape and their spectral properties, dominated by non-thermal synchrotron and inverse-Compton scattering. We propose a post-processing method to analyse the broad-band emission of SNRs based on three-dimensional hydrodynamical simulations. From the hydrodynamical data, we estimate the distribution of non-thermal electrons accelerated at the shock wave and follow the subsequent evolution as they lose or gain energy by adiabatic expansion or compression and emit energy by radiation. As a first test case, we use a simulation of a bipolar supernova expanding into a cloudy medium. We find that our method qualitatively reproduces the main observational features of typical SNRs and produces fluxes that agree with observations to within a factor of a few. allowing for further use in more extended sets of models.

GRB 101225A is a prototype of the class of blackbody-dominated (BBD) gamma-ray bursts (GRBs). It has been suggested that BBD-GRBs result from the merger of a binary system formed by a neutron star and the helium core of an evolved star. We have modelled the propagation of ultrarelativistic jets through the environment left behind the merger by means of relativistic hydrodynamic simulations. In this paper, the output of our numerical models is post-processed to obtain the (thermal) radiative signature of the resulting outflow. We outline the most relevant dynamical details of the jet propagation and connect them to the generation of thermal radiation in GRB events akin to that of GRB 101225A. A comprehensive parameter study of the jet/environment interaction has been performed and synthetic light curves are confronted with the observational data. The thermal emission in our models originates from the interaction between the jet and the hydrogen envelope ejected during the neutron star/He core merger. We find that the lack of a classical afterglow and the accompanying thermal emission in BBD-GRBs can be explained by the interaction of an ultrarelativistic jet with a toroidally shaped ejecta whose axis coincides with the binary rotation axis. The spectral inversion and reddening happening at about 2 d in GRB 101225A can be related to the time at which the massive shell ejected in an early phase of the common envelope evolution of the progenitor system is completely ablated by the ultrarelativistic jet.

Blackbody-dominated (BBD) gamma-ray bursts (GRBs) are events characterized by long durations and the presence of a significant thermal component following the prompt emission, as well as by the absence of a typical afterglow. GRB 101225A is the most prominent member of this class. A plausible progenitor system for it and for BBD-GRBs is the merger a neutron star and a helium core of an evolved, massive star. Using relativistic hydrodynamic simulations we model the propagation of ultrarelativistic jets through the environments created by such mergers. In a previous paper we showed that the thermal emission in BBD-GRBs is linked to the interaction of an ultrarelativistic jet with the ejected envelope of the secondary star of the binary. Here we focus on explaining the emission properties of BBD-GRBs computing the whole radiative signature (both thermal and non-thermal) of the jet dynamical evolution. The non-thermal emission of the forward shock of the jet is dominant during the early phases of the evolution, when that shock is moderately relativistic. Our models do not produce a classical afterglow because the quick deceleration of the jet results primarily from the mass entrainment in the beam, and not from the process of plowing mass from the external medium in front of the GRB ejecta. The contribution of the reverse shock is of the same magnitude than that of the forward shock during the first 80 min after the GRB. Later, it quickly fades because the jet/environment interaction chocks the ultrarelativistic jet beam and effectively dumps the reverse shock. In agreement with observations, we obtain rather flat light curves during the first 2 d after the GRB, and a spectral evolution consistent with the observed reddening of the system.

We explore the signature imprinted by dynamically relevant magnetic fields on the spectral energy distribution (SED) of blazars. It is assumed that the emission from these sources originates from the collision of cold plasma shells, whose magnetohydrodynamic evolution we compute by numerically solving Riemann problems. We compute the SEDs including the most relevant radiative processes and scan a broad parameter space that encompasses a significant fraction of the commonly accepted values of not directly measurable physical properties. We reproduce the standard double hump SED found in blazar observations for unmagnetized shells, but show that the prototype double hump structure of blazars can also be reproduced if the dynamical source of the radiation field is very ultrarelativistic both, in a kinematically sense (namely, if it has Lorentz factors \gtrsim 50) and regarding its magnetization (e.g., with flow magnetizations \sigma \simeq 0.1). A fair fraction of the {\em blazar sequence} could be explained as a consequence of shell magnetization: negligible magnetization in FSRQs, and moderate or large (and uniform) magnetization in BL Lacs. The predicted photon spectral indices (\gph) in the \gamma-ray band are above the observed values (\Gamma_{\rm ph, obs} \lesssim 2.6 for sources with redshifts 0.4\le z \le 0.6) if the magnetization of the sources is moderate (\sigma \simeq 10^{-2}).

The internal-shocks scenario in relativistic jets has been used to explain the variability of blazars' outflow emission. Recent simulations have shown that the magnetic field alters the dynamics of these shocks producing a whole zoo of spectral energy density patterns. However, the role played by magnetization in such high-energy emission is still not entirely understood. With the aid of \emph{Fermi}'s second LAT AGN catalog, a comparison with observations in the $\gamma$-ray band was performed, in order to identify the effects of the magnetic field.

The temporal and spatial spectral evolution of the jets of AGN can be studied with multi-frequency, multi-epoch VLBI observations. The combination of both, morphological and spectral parameters can be used to derive source intrinsic physical properties such as the magnetic field and the non-thermal particle density. In the first two papers of this series, we analyzed the single-dish light curves and the VLBI kinematics of the blazar CTA 102 and suggested a shock-shock interaction between a traveling and a standing shock wave as a possible scenario to explain the observed evolution of the component associated to the 2006 flare. In this paper we investigate the core-shift and spectral evolution to test our hypothesis of a shock-shock interaction. We used 8 multi-frequency VLBA observations to analyze the temporal and spatial evolution of the spectral parameters during the flare. We observed CTA 102 between May 2005 and April 2007 using the VLBA at six different frequencies spanning from 2 - 86 GHz. After the calibrated VLBA images were corrected for opacity, we performed a detailed spectral analysis. From the derived values we estimated the magnetic field and the density of the relativistic particles. The detailed analysis of the opacity shift reveals that the position of the jet core is proportional to nu^-1 with some temporal variations. The value suggests possible equipartition between magnetic field energy and particle kinetic energy densities at the most compact regions. From the variation of the physical parameters we deduced that the 2006 flare in CTA 102 is connected to the ejection of a new traveling feature (t=2005.9) and the interaction between this shock wave and a stationary structure around 0.1 mas from the core. The source kinematics together with the spectral and structural variations can be described by helical motions in an over-pressured jet.

The 2006 radio flare in the blazar CTA102 is the largest ever reported for this source. The analysis of the single--dish light curves revealed a possible shock-shock interaction as a scenario for the 2006 outburst. In order to confirm this hypothesis we analyzed eight multi-frequency (2 GHz - 86 GHz) VLBI observations covering the 2006 flare. We detected a radio component ejected around 2005.9 and a stationary feature around 0.1 mas away from the core. We also performed a core-shift and spectral analysis on the fully calibrated VLBI maps. The frequency dependent position of the core slightly deviates from the typical $\nu^{-1}$-dependence during the flaring state. This behaviour may reflect a mismatch between magnetic energy density and energy density of the relativistic particles. From the core-shift corrected maps we obtained the spectral parameters and computed the temporal and spatial evolution of the jet intrinsic parameters such as the magnetic field and the particle density. The overall picture of the source could be best interpreted as an over-pressured jet perturbed by traveling shock waves.

Broadband emission from relativistic outflows (jets) of active galactic nuclei (AGN) and gamma-ray bursts (GRBs) contains valuable information about the nature of the jet itself, and about the central engine which launches it. Using special relativistic hydrodynamics and magnetohydronamics simulations we study the dynamics of the jet and its interaction with the surrounding medium. The observational signature of the simulated jets is computed using a radiative transfer code developed specifically for the purpose of computing multi-wavelength, time-dependent, non-thermal emission from astrophysical plasmas. We present results of a series of long-term projects devoted to understanding the dynamics and emission of jets in parsec-scale AGN jets, blazars and the afterglow phase of the GRBs.

We performed multifrequency multiepoch Very Long Baseline Array (VLBA) observations of the blazar CTA 102 during its 2006 radio flare, the strongest ever reported for this source. These observations provide an excellent opportunity to investigate the evolution of the physical properties of blazars, especially during these flaring events. We want to study the kinematic changes in the source during the strong radio outburst in April 2006 and test the assumption of a shock-shock interaction. This assumption is based on the analysis and modeling of the single-dish observations of CTA\,102 (Paper I). In this paper we study the kinematics of CTA 102 at several frequencies using VLBI observations. From the modeled jet features we derived estimates for the evolution of the physical parameters, such as the particle density and the magnetic field. Furthermore, we combined our observations during the 2006 flare with long-term VLBA monitoring of the source at 15 GHz and 43 GHz. We cross-identified seven features throughout our entire multifrequency observations and find evidence of two possible recollimation shocks around 0.1 mas (deprojected 18 pc at a viewing angle of 2.6 degrees) and 6.0 mas (deprojected 1 kpc) from the core. The 43 GHz observations reveal a feature ejected at epoch $t_\mathrm{ej}=2005.9\pm0.2$, which could be connected to the 2006 April radio flare. Furthermore, this feature might be associated with the traveling component involved in the possible shock-shock interaction, which gives rise to the observed double peak structure in the single-dish light curves reported in Paper I.

The recent transient event Swift J1644+57 has been interpreted as emission from a collimated relativistic jet, powered by the sudden onset of accretion onto a supermassive black hole following the tidal disruption of a star. Here we model the radio-microwave emission as synchrotron radiation produced by the shock interaction between the jet and the gaseous circumnuclear medium (CNM). At early times after the onset of the jet (t < 5-10 days) a reverse shock propagates through and decelerates the ejecta, while at later times the outflow approaches the Blandford-McKee self-similar evolution (possibly modified by additional late energy injection). The achromatic break in the radio light curve of J1644+57 is naturally explained as the transition between these phases. We show that the temporal indices of the pre- and post-break light curve are consistent with those predicted if the CNM has a wind-type radial density profile n ~ 1/r^2. The observed synchrotron frequencies and self-absorbed flux constrain the fraction of the post-shock thermal energy in relativistic electrons epsilon_e ~ 0.03-0.1; the CNM density at 1e18 cm ~ 1-10 1/cm^3; and the initial Lorentz factor Gamma_j ~ 10-20 and opening angle theta_j ~ (0.1-1)Gamma_j^(-1) ~ 0.01-0.1 of the jet. Radio modeling thus provides robust independent evidence for a narrowly collimated outflow. Extending our model to the future evolution of J1644+57, we predict that the radio flux at low frequencies (<~ few GHz) will begin to brighten more rapidly once the characteristic frequency crosses below the self-absorption frequency on a timescale of months (indeed, such a transition may already have begun). Our results demonstrate that relativistic outflows from tidal disruption events provide a unique probe of the conditions in distant, previously inactive galactic nuclei, complementing studies of normal AGN.

Common models of blazars and gamma-ray bursts assume that the plasma underlying the ob- served phenomenology is magnetized to some extent. Within this context, radiative signatures of dissipation of kinetic and conversion of magnetic energy in internal shocks of relativistic magnetized outflows are studied. We model internal shocks as being caused by collisions of homogeneous plasma shells. We compute the flow state after the shell interaction by solving Riemann problems at the contact surface between the colliding shells, and then compute the emission from the resulting shocks. Under the assumption of a constant flow luminosity we find that there is a clear difference between the models where both shells are weakly magne- tized ({\sigma}<\sim0.01) and those where, at least, one shell has a {\sigma}>\sim0.01. We obtain that the radiative efficiency is largest for models in which, regardless of the ordering, one shell is weakly and the other strongly magnetized. Substantial differences between weakly and strongly magne- tized shell collisions are observed in the inverse-Compton part of the spectrum, as well as in the optical, X-ray and 1GeV light curves. We propose a way to distinguish observationally between weakly magnetized from magnetized internal shocks by comparing the maximum frequency of the inverse-Compton and synchrotron part of the spectrum to the ratio of the inverse-Compton and synchrotron fluence. Finally, our results suggest that LBL blazars may correspond to barely magnetized flows, while HBL blazars could correspond to moderately magnetized ones. Indeed, by comparing with actual blazar observations we conclude that the magnetization of typical blazars is {\sigma} <\sim 0.01 for the internal shock model to be valid in these sources.

The analysis of the radio light curves of the blazar CTA102 during its 2006 flare revealed a possible interaction between a standing shock wave and a traveling one. In order to better understand this highly non-linear process, we used a relativistic hydrodynamic code to simulate the high energy interaction and its related emission. The calculated synchrotron emission from these simulations showed an increase in turnover flux density, $S_{m}$, and turnover frequency, $\nu_{m}$, during the interaction and decrease to its initial values after the passage of the traveling shock wave.

We study the dynamic and radiative efficiency of conversion of kinetic-to-thermal/magnetic energy by internal shocks in relativistic magnetized outflows. A parameter study of a large number of collisions of cylindrical shells is performed. We explore how, while keeping the total flow luminosity constant, the variable fluid magnetization influences the efficiency and find that the interaction of shells in a mildly magnetized jet yields higher dynamic, but lower radiative efficiency than in a non-magnetized flow. A multi-wavelength radiative signature of different shell magnetization is computed assuming that relativistic particles are accelerated at internal shocks.

The Compact Symmetric Objects (CSOs) are small (<1 kiloparsec) and powerful extragalactic radio sources showing emission on both sides of an active galactic nucleus and no signs of strong relativistic beaming. They may be young radio sources, progenitors of large FRII radio galaxies. We aim to study the statistical properties of CSOs by constructing and investigating a new large sample of CSO candidates on the basis of dual-frequency, parsec-scale morphology. For the candidate selection we utilized VLBI data for 4170 extragalactic objects obtained simultaneously at 2.3 and 8.6 GHz (S and X band) within the VLBA Calibrator Survey 1-6 and the Research and Development - VLBA projects. Properties of their broad-band radio spectra were characterized by using RATAN-600 observations. Numerical modeling was applied in an attempt to explain the observed effects. A sample of 64 candidate CSOs is identified. The median two-point S-X band spectral index of parsec-scale hotspots is found to be -0.52; with the median brightness temperature ~10^9 K at X band. Statistical analysis reveals a systematic difference between positions of brightest CSO components (associated with hotspots) measured in the S and X bands. The distance between these components is found to be on average 0.32+/-0.06 mas greater at 8.6 GHz than at 2.3 GHz. This difference in distances cannot be explained by different resolutions at the S and X bands. It is a manifestation of spectral index gradients across CSO components, which may potentially provide important physical information about them. Despite our detailed numerical modeling of a CSO hotspot, the model was not able to reproduce the magnitude of the observed positional difference. A more detailed modeling may shed light on the origin of the effect.

Long Gamma-Ray Bursts (GRBs) are the most dramatic examples of massive stellar deaths, usually as- sociated with supernovae (Woosley et al. 2006). They release ultra-relativistic jets producing non-thermal emission through synchrotron radiation as they interact with the surrounding medium (Zhang et al. 2004). Here we report observations of the peculiar GRB 101225A (the "Christmas burst"). Its gamma-ray emission was exceptionally long and followed by a bright X-ray transient with a hot thermal component and an unusual optical counterpart. During the first 10 days, the optical emission evolved as an expanding, cooling blackbody after which an additional component, consistent with a faint supernova, emerged. We determine its distance to 1.6 Gpc by fitting the spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a helium star-neutron star merger that underwent a common envelope phase expelling its hydrogen envelope. The resulting explosion created a GRB-like jet which gets thermalized by interacting with the dense, previously ejected material and thus creating the observed black-body, until finally the emission from the supernova dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy (Campana et al. 2011).

The afterglow emission that follows gamma-ray bursts (GRBs) contains valuable information about the circumburst medium and, therefore, about the GRB progenitor. Theoretical studies of GRB blast waves, however, are often limited to simple density profiles for the external medium (mostly constant density and power-law R^{-k} ones). We argue that a large fraction of long-duration GRBs should take place in massive stellar clusters where the circumburst medium is much more complicated. As a case study, we simulate the propagation of a GRB blast wave in a medium shaped by the collision of the winds of O and Wolf-Rayet stars, the typical distance of which is d /sim 0.1 - 1 pc. Assuming a spherical blast wave, the afterglow light curve shows a flattening followed by a shallow break on a timescale from hours up to a week after the burst, which is a result of the propagation of the blast wave through the shocked wind region. If the blast wave is collimated, the jet break may, in some cases, become very pronounced with the post break decline of the light curve as steep as t-5. Inverse Compton scattering of ultra-violet photons from the nearby star off energetic electrons in the blast wave leads to a bright \simGeV afterglow flare that may be detectable by Fermi.

The so-called internal shock model aims to explain the light-curves and spectra produced by non-thermal processes originated in the flow of blazars and gamma-ray bursts. A long standing question is whether the tenuous collisionless shocks, driven inside a relativistic flow, are efficient enough to explain the amount of energy observed as compared with the expected kinetic power of the outflow. In this work we study the dynamic efficiency of conversion of kinetic-to- thermal/magnetic energy of internal shocks in relativistic magnetized outflows. We find that the collision between shells with a non-zero relative velocity can yield either two oppositely moving shocks (in the frame where the contact surface is at rest), or a reverse shock and a forward rarefaction. For moderately magnetized shocks (magnetization {\sigma} ~ 0.1), the dynamic efficiency in a single two-shell interaction can be as large as 40%. Hence, the dynamic efficiency of moderately magnetized shocks is larger than in the corresponding unmagnetized two-shell interaction. We find that the efficiency is only weakly dependent on the Lorentz factor of the shells and, thus internal shocks in the magnetized flow of blazars and gamma-ray bursts are approximately equally efficient.

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

We use high-resolution relativistic MHD simulations coupled with a radiative transfer code to compute multiwavelength afterglow light curves of magnetized ejecta of gamma-ray bursts interacting with a uniform circumburst medium. The aim of our study is to determine how the magnetization of the ejecta at large distance from the central engine influences the afterglow emission, and to assess whether observations can be reliably used to infer the strength of the magnetic field. We find that, for typical parameters of the ejecta, the emission from the reverse shock peaks for magnetization $\sigma_0 \sim 0.01 - 0.1$ of the flow, and that it is greatly suppressed for higher $\sigma_0$. The emission from the forward shock shows an achromatic break shortly after the end of the burst marking the onset of the self-similar evolution of the blast wave. Fitting the early afterglow of GRB 990123 and 090102 with our numerical models we infer respective magnetizations of $\sigma_0 \sim 0.01$ and $\sigma_0 \sim 0.1$ for these bursts. We argue that the lack of observed reverse shock emission from the majority of the bursts can be understood if $\sigma_0 \simmore 0.1$, since we obtain that the luminosity of the reverse shock decreases significantly for $\sigma_0 \sim 1$. For ejecta with $\sigma_0 \simmore 0.1$ our models predict that there is sufficient energy left in the magnetic field, at least during an interval of ~10 times the burst duration, to produce a substantial emission if the magnetic energy can be dissipated (for instance, due to resistive effects) and radiated away.

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

We study the dynamic efficiency of conversion of kinetic-to-thermal/magnetic energy of internal shocks in relativistic magnetized outflows. We model internal shocks as being caused by collisions of shells of plasma with the same energy flux and a non-zero relative velocity. The contact surface, where the interaction between the shells takes place, can break up either into two oppositely moving shocks (in the frame where the contact surface is at rest), or into a reverse shock and a forward rarefaction. We find that for moderately magnetized shocks (magnetization $\sigma\simeq 0.1$), the dynamic efficiency in a single two-shell interaction can be as large as 40%. Thus, the dynamic efficiency of moderately magnetized shocks is larger than in the corresponding unmagnetized two-shell interaction. If the slower shell propagates with a sufficiently large velocity, the efficiency is only weakly dependent on its Lorentz factor. Consequently, the dynamic efficiency of shell interactions in the magnetized flow of blazars and gamma-ray bursts is effectively the same. These results are quantitatively rather independent on the equation of state of the plasma. The radiative efficiency of the process is expected to be a fraction $f_r<1$ of the estimated dynamic one, the exact value of $f_r$ depending on the particularities of the emission processes which radiate away the thermal or magnetic energy of the shocked states.

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

(Abridged) We aim to quantitatively understand the dynamical effect and observational signatures of magnetization of the GRB ejecta on the onset of the afterglow. We perform ultrahigh-resolution one-dimensional relativistic MHD simulations of the interaction of a radially expanding, magnetized ejecta with the interstellar medium. The need of ultrahigh numerical resolution derives from the extreme jump conditions in the region of interaction between the ejecta and the circumburst medium. We study the evolution of an ultrarelativistic shell all the way to a the self-similar asymptotic phase. Our simulations show that the complete evolution can be characterized in terms of two parameters, namely, the \xi parameter introduced by Sari & Piran (1995) and the magnetization \sigma_0. We exploit this property by producing numerical models where the shell Lorentz factor is \gamma_0 ~ tens and rescaling the results to arbitrarily large \gamma_0. We find that the reverse shock is typically very weak or absent for ejecta characterized by \sigma_0 >~ 1. The onset of the forward shock emission is strongly affected by the magnetization. On the other hand, the magnetic energy of the shell is transfered to the external medium on a short timescale (~several times the duration of the burst). The later forward shock emission does not contain information for the initial magnetization of the flow. The asymptotic evolution of strongly magnetized shells, after they have suffred a substantial deceleration, resembles that of hydrodynamic shells, i.e., they fully enter in the Blandford-McKee self-similar regime.

Recent theoretical work has pointed out that the transition layer between a jet an the medium surrounding it may be more complex than previously thought. Under physically realizable conditions, the transverse profile of the Lorentz factor in the boundary layer can be non-monotonic, displaying the absolute maximum where the flow is faster than at the jet spine, followed by an steep fall off. Likewise, the rest-mass density, reaches an absolute minimum (coincident with the maximum in Lorentz factor) and then grows until it reaches the external medium value. Such a behavior is in contrast to the standard monotonic decline of the Lorentz factor (from a maximum value at the jet central spine) and the corresponding increase of the rest-mass density (from the minimum reached at the jet core). We study the emission properties of the aforementioned anomalous shear layer structures in kiloparsec-scale jets aiming to show observable differences with respect to conventional monotonic and smooth boundary layers.

We study the afterglow phases of a GRB through relativistic magnetohydrodynamic simulations. The evolution of a relativistic shell propagating into a homogeneous external medium is followed. We focus on the effect of the magnetization of the ejecta on the initial phases of the ejecta-external medium interaction. In particular we are studying the condition for the existence of a reverse shock into the ejecta, the timescale for the transfer of the energy from the shell to the shocked medium and the resulting multiwavelength light curves. To this end, we have developed a novel scheme to include non-thermal processeses which is coupled to the relativistic magnetohydrodynamic code MRGENESIS in order to compute the non-thermal synchrotron radiation.

The role of magnetic fields in gamma-ray burst (GRB) flows remains controversial. The study of the early afterglow phases and, in particular, of the reverse shock dynamics and associated emission offers a promising probe of the magnetization of the ejecta. In this paper, we derive the conditions for the existence of a reverse shock in arbitrarily magnetized ejecta that decelerate and interact with the circumburst medium. Both constant and wind-like density profiles are considered. We show, in contrast to previous estimates, that ejecta with magnetization larger than unity are not crossed by a reverse shock for a large fraction of the parameter space relevant to GRB flows. Allowing for shell spreading, there is always a relativistic or mildly relativistic reverse shock forming in sigma_o < 0.3 ejecta. From this, we conclude that the paucity of optical flashes, believed to be a distinctive signature of a reverse shock, may be explained by the existence of dynamically important magnetic fields in the ejecta.

(Abridged): We study the collision of magnetized irregularities (shells) in relativistic outflows in order to explain the origin of the generic phenomenology observed in the non-thermal emission of both blazars and gamma-ray bursts. We focus on the influence of the magnetic field on the collision dynamics, and we further investigate how the properties of the observed radiation depend on the strength of the initial magnetic field and on the initial internal energy density of the flow. The collisions of magnetized shells and the radiation resulting from these collisions are calculated using the 1D relativistic magnetohydrodynamics code MRGENESIS. The interaction of the shells with the external medium prior to their collision is also determined using an exact solver for the corresponding 1D relativistic magnetohydrodynamic Riemann problem. Our simulations show that two magnetization parameters - the ratio of magnetic energy density and thermal energy density, \alpha_B, and the ratio of magnetic energy density and mass-energy density, \sigma - play an important role in the pre-collision phase, while the dynamics of the collision and the properties of the light curves depend mostly on the magnetization parameter \sigma. The interaction of the shells with the external medium changes the flow properties at their edges prior to the collision. For sufficiently dense shells moving at large Lorentz factors (\simgt 25) these properties depend only on the magnetization parameter \sigma. Internal shocks in GRBs may reach maximum efficiencies of conversion of kinetic into thermal energy between 6% and 10%, while in case of blazars, the maximum efficiencies are \sim 2%.

We show some of the most important reasons why the likely fate of the merger of a neutron star with another compact object may be to yield a short gamma-ray burst (sGRB). Emphasis is made on some robust results that general relativistic (magneto)hydrodynamic simulations have established regarding the aforementioned subject.

Light curves of the long-period Algols are known for their complex shape (asymmetry in the eclipse, light variations outside eclipse, changes from cycle-to-cycle), and their interpretation is not possible in the standard model of binary stars. Complex structures present in these active Algol systems could be studied with the eclipse-mapping method which was successfully applied to the new 7-color photometric observations in the Geneva system of W Cru, belonging to the isolated group of these active Algols. Several cycles of this long-period (198.5 days) eclipsing binary have been covered by observations. We have used a modified Rutten's approach to the eclipse-mapping. The optimization of the system's parameters and the recovery of the disk intensity distribution are performed using a genetic algorithm (GA). It is found that a hot component is hidden in the thick accretion disk which confirms previous findings. The mass of the component, M1 = 8.2 Ms indicates that it is a mid-B type star. The mass-losing component is filling its critical lobe which means it is a G-type supergiant with a mass M2 = 1.6 Ms. The disk is geometrically very extended and its outer radius is about 80% of the primary's critical lobe. A reconstructed image reveals a rather clumpy and nonuniform brightness distribution of an accretion disk rim in this almost edge-on seen system. This clumpyness accounts for light curve distortions and asymmetries, as well as for secular changes.

We present the results of a detailed temporal analysis of the bright BL Lac object Mrk 421 using the three available long timing mode observations by the EPIC PN camera. This detector mode is characterized by its long life time and is largely free of photon pile-up problems. The source was found in different intensity and variability states differing by up to more than a factor of three in count rate. A time resolved cross correlation analysis between the soft and hard energy bands revealed that the characteristics of the correlated emission, with lags of both signs, change on time scales of a few thousand seconds. Individual spectra, resolved on time scales of a few hundread seconds, can be quite well fitted by a broken power law. We find significant spectral variations on time scales as short as 500-1000 sec. Both the hard and the soft band spectral indices show a non-linear correlation with the source flux. A comparison of the observed light curves with numerical results from relativistic hydrodynamic computer simulations of the currently favored shock-in-jet models indicates that any determination of the jet's physical parameters from `simple' emission models must be regarded with caution: at any time we are seeing the emission from several emission regions distinct in space and time, which are connected by the complex hydrodynamic evolution of the non-uniform jet.

We present results of a detailed numerical study and theoretical analysis of the dynamics of internal shocks in relativistic jets and the non-thermal flares associated with these shocks. In our model internal shocks result from collisions of density inhomogeneities (shells) in relativistic jet flows. We find that the merged shell resulting from the inelastic collision of shells has a complicated internal structure due to the non-linear dynamics of the interaction. Furthermore, the instantaneous efficiency for converting kinetic energy into thermal energy is found to be almost twice as high as theoretically expected during the period of significant emission. The Lorentz factors of the internal shocks are correlated with the initial inertial masses of the shells. Because of the complexity of the non-linear evolution the merged shell becomes very inhomogeneous and simple one-zone models are inadequate to extract physical parameters of the emitting region from the resulting light curves. In order to improve on these one-zone approximations, we propose a novel way of analyzing the space-time properties of the emission. Based on these properties we construct an analytic model of non-thermal flares which can be used to constrain some (unobservable) physical parameters of the internal shocks. These are the ratio of the Lorentz factors between the forward and the reverse shock (caused by the shell collision), and the shell crossing times of these shocks. The analytic model is validated by applying it to the synthetic light curves computed from our models. It can equally well be applied to observations.

We present the results of relativistic hydrodynamic simulations of the collision of two dense shells in a uniform external medium, as envisaged in the internal shock model for BL Lac jets. The non-thermal radiation produced by highly energetic electrons injected at the relativistic shocks is computed following their temporal and spatial evolution. The acceleration of electrons at the relativistic shocks is parametrized using two different models and the corresponding X-ray light curves are computed. We find that the interaction time scale of the two shells is influenced by an interaction with the external medium. For the chosen parameter sets, the efficiency of the collision in converting dissipated kinetic energy into the observed X-ray radiation is of the order of one percent.