We explore the unique supernova (SN) "Refsdal" - the first discovered gravitationally lensed SN with multiple images. A large magnification provided by the galactic-scale lens, augmented by the cluster lens, gave us the opportunity to perform a detailed modelling of a distant SN at $z \simeq 1.5$. % We present results of radiation hydrodynamics modelling of \snref. According to our calculations, the \snref progenitor is likely to be a more massive and energetic version of \sna, i.e. a blue supergiant star with a radius of ${\rm R} = 50 R_{\odot}$ before an explosion and a total mass of ${\rm M}= 26 M_{\odot}$. The total energy release is ${\rm E}=5\times10^{51}$ erg. Reconstruction of SN light curves allowed us to obtain time delays and magnifications for the images S2-S4 relative to S1 with higher accuracy than previous estimates of Rodney et al. We estimate the Hubble constant $H_0 = 68.6^{+13.6}_{-9.7}$ km s$^{-1}$ Mpc$^{-1}$ via re-scaling the time delays predicted by different lens models to match the values obtained in this work. With more photometric data on the fifth image SX, we will be able to further refine the time delay and magnification estimates for SX and obtain competitive constraints on $H_0$.

We consider a possibility of detecting the jitter effect of apparent celestial positions of distant sources due to local fluctuations of the Galaxy gravitational field. It is proposed to observe two samples of extragalactic sources (target and control) in different sky directions using the high-precision radio interferometry. It is shown that on a scale of ~2 years, it is possible to detect a systematic increase in the standard deviation of measured arc lengths of pairs of target sources compared to the control ones at the $3\sigma$-level if the accuracy of differential astrometric observations is around 10 $\mu$as. For the current state-of-the-art accuracy of 30 $\mu$as achieved at the KVN or VERA interferometers, which have shorter baselines in comparison with VLBI, the target and control samples will differ only at the 2$\sigma$-level on the scale of 10 years. To achieve the $3\sigma$-level on this time interval, it is necessary to improve the accuracy up to ~20 $\mu$as. Other possible effects that can also affect the arc length measurements between two sources are discussed, and an observational strategy to minimize them is suggested.

In this work, we present the first detailed analysis of the supernova remnant RX J1713.7-3946 in the hard X-ray energy range with the IBIS coded-mask telescope on board the INTEGRAL observatory. The shell-type morphology of the entire remnant is mapped in hard X-rays for the first time and significantly detected up to 50 keV. The IBIS sky image of RX J1713.7-3946, accumulated over 14 years of operations, demonstrates two extended hard X-ray sources. These sources are spatially consistent with northwest and southwest rims of RX J1713.7-3946 and are also clearly visible at energies below 10 keV with XMM-Newton. This points to a single emission mechanism operating in soft and hard X-rays. The INTEGRAL 17-120 keV spectrum of RX J1713.7-3946 is characterized by a power-law continuum with the photon index of $\Gamma\approx3$, that is significantly softer than $\Gamma\approx2$ determined by XMM-Newton in the 1-10 keV energy range, suggesting a progressive steepening of the spectrum with the energy.

Moderately strong shocks arise naturally when two subclusters merge. For instance, when a smaller subcluster falls into the gravitational potential of a more massive cluster, a bow shock is formed and moves together with the subcluster. After pericenter passage, however, the subcluster is decelerated by the gravity of the main cluster, while the shock continues moving away from the cluster center. These shocks are considered as promising candidates for powering radio relics found in many clusters. The aim of this paper is to explore the fate of such shocks when they travel to the cluster outskirts, far from the place where the shocks were initiated. In a uniform medium, such a "runaway" shock should weaken with distance. However, as shocks move to large radii in galaxy clusters, the shock is moving down a steep density gradient that helps the shock to maintain its strength over a large distance. Observations and numerical simulations show that, beyond $R_{500}$, gas density profiles are as steep as, or steeper than, $\sim r^{-3}$, suggesting that there exists a "Habitable zone" for moderately strong shocks in cluster outskirts where the shock strength can be maintained or even amplified. A characteristic feature of runaway shocks is that the strong compression, relative to the initial state, is confined to a narrow region just behind the shock. Therefore, if such a shock runs over a region with a pre-existing population of relativistic particles, then the boost in radio emissivity, due to pure adiabatic compression, will also be confined to a narrow radial shell.

We study a merger of the NGC 4839 group with the Coma cluster using X-ray observations from the XMM-Newton and Chandra telescopes. X-ray data show two prominent features: (i) a long (~600 kpc in projection) and bent tail of cool gas trailing (towards south-west) the optical center of NGC 4839, and ii) a "sheath" region of enhanced X-ray surface brightness enveloping the group, which is due to hotter gas. While at first glance the X-ray images suggest that we are witnessing the first infall of NGC 4839 into the Coma cluster core, we argue that a post-merger scenario provides a better explanation of the observed features and illustrate this with a series of numerical simulations. In this scenario, the tail is formed when the group, initially moving to the south-west, reverses its radial velocity after crossing the apocenter, the ram pressure ceases and the ram-pressure-displaced gas falls back toward the center of the group and overshoots it. Shortly after the apocenter passage, the optical galaxy, dark matter and gaseous core move in a north-east direction, while the displaced gas continues moving to the south-west. The "sheath" is explained as being due to interaction of the re-infalling group with its own tail of stripped gas mixed with the Coma gas. In this scenario, the shock, driven by the group before reaching the apocenter, has already detached from the group and would be located close to the famous relic to the south-west of the Coma cluster.

We show that there is a new class of gas tails - slingshot tails - which form as a subhalo (i.e. a subcluster or early-type cluster galaxy) moves away from the cluster center towards the apocenter of its orbit. These tails can point perpendicular or even opposite to the subhalo direction of motion, not tracing the recent orbital path. Thus, the observed tail direction can be misleading, and we caution against naive conclusions regarding the subhalo's direction of motion based on the tail direction. A head-tail morphology of a galaxy's or subcluster's gaseous atmosphere is usually attributed to ram pressure stripping and the widely applied conclusion is that gas stripped tail traces the most recent orbit. However, during the slingshot tail stage, the subhalo is not being ram pressure stripped (RPS) and the tail is shaped by tidal forces more than just the ram pressure. Thus, applying a classic RPS scenario to a slingshot tail leads not only to an incorrect conclusion regarding the direction of motion, but also to incorrect conclusions in regard to the subhalo velocity, expected locations of shear flows, instabilities and mixing. We describe the genesis and morphology of slingshot tails using data from binary cluster merger simulations, discuss their observable features and how to distinguish them from classic RPS tails. We identify three examples from the literature that are not RPS tails but slingshot tails and discuss other potential candidates.

We discuss constraints on the mass density distribution (parameterized as $\rho\propto r^{-\gamma}$) in early-type galaxies provided by strong lensing and stellar kinematics data. The constraints come from mass measurements at two `pinch' radii. One `pinch' radius $r_1=2.2 R_{Einst}$ is defined such that the Einstein (i.e. aperture) mass can be converted to the spherical mass almost independently of the mass-model. Another `pinch' radius $r_2=R_{opt}$ is chosen so that the dynamical mass, derived from the line-of-sight velocity dispersion, is least sensitive to the anisotropy of stellar orbits. We verified the performance of this approach on a sample of simulated elliptical galaxies and on a sample of 15 SLACS lens galaxies at $0.01 \leq z \leq 0.35$, which have already been analysed in Barnabe et al. (2011) by the self-consistent joint lensing and kinematic code. For massive simulated galaxies the density slope $\gamma$ is recovered with an accuracy of $\sim 13\%$, unless $r_1$ and $r_2$ happen to be close to each other. For SLACS galaxies, we found good overall agreement with the results of Barnabe et al. (2011) with a sample-averaged slope $\gamma=2.1\pm0.05$. While the two-pinch-radii approach has larger statistical uncertainties, it is much simpler and uses only few arithmetic operations with directly observable quantities.

We investigate the influence of random variations of the Galactic gravitational field on the apparent celestial positions of extragalactic sources. The basic statistical characteristics of a stochastic process (first-order moments, an autocorrelation function and a power spectral density) are used to describe a light ray deflection in a gravitational field of randomly moving point masses as a function of the source coordinates. We map a 2D distribution of the standard deviation of the angular shifts in positions of distant sources (including reference sources of the International Celestial Reference Frame) with respect to their true positions. For different Galactic matter distributions the standard deviation of the offset angle can reach several tens of $\mu as$ (microarcsecond) toward the Galactic center, decreasing down to 4-6 $\mu as$ at high galactic latitudes. The conditional standard deviation (`jitter') of 2.5 $\mu as$ is reached within 10 years at high galactic latitudes and within a few months toward the inner part of the Galaxy. The photometric microlensing events are not expected to be disturbed by astrometric random variations anywhere except the inner part of the Galaxy as the Einstein--Chvolson times are typically much shorter than the jittering timescale. While a jitter of a single reference source can be up to dozens of $\mu as$ over some reasonable observational time, using a sample of reference sources would reduce the error in relative astrometry. The obtained results can be used for estimating the physical upper limits on the time-dependent accuracy of astrometric measurements.

We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $\sigma_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{\sigma_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.

We consider an influence of a non-stationary gravitational field of the Galaxy on the visible positions of extragalactic sources. A contribution of the baryonic component of the galactic matter as well as of the hidden matter (including a population of brown dwarfs) were took into account. The observed variations of the deflection angle of light rays in a gravitational field of randomly moving point-like masses can be considered as a stochastic process. Using such an approach we constructed an autocorrelation function of studied stochastic process and found that its relative changes are about 15% for one year and about 35% for ten years.

We discuss a simple and fast method for estimating masses of early-type galaxies from optical data and compare the results with X-ray derived masses. The optical method relies only on the most basic observables such as the surface brightness $I(R)$ and the line-of-sight velocity dispersion $\sigma_p(R)$ profiles and provides an anisotropy-independent estimate of the galaxy circular speed $V_c$. The mass-anisotropy degeneracy is effectively overcome by evaluating $V_c$ at a characteristic radius $R_{\rm sweet}$ defined from {\it local} properties of observed profiles. The sweet radius $R_{\rm sweet}$ is expected to lie close to $R_2$, where $I(R) \propto R^{-2}$, and not far from the effective radius $R_{\rm eff}$. We apply the method to a sample of five X-ray bright elliptical galaxies observed with the 6-m telescope BTA-6 in Russia. We then compare the optical $V_c$-estimate with the X-ray derived value, and discuss possible constraints on the non-thermal pressure in the hot gas and configuration of stellar orbits. We find that the average ratio of the optical $V_c$-estimate to the X-ray one is equal to $\approx 0.98$ with $11 \%$ scatter, i.e. there is no evidence for the large non-thermal pressure contribution in the gas at $\sim R_{\rm sweet}$. From analysis of the Lick indices H$\beta$, Mgb, Fe5270 and Fe5335, we calculate the mass of the stellar component within the sweet radius. We conclude that a typical dark matter fraction inside $R_{\rm sweet}$ in the sample galaxies is $\sim 60\%$ for the Salpeter IMF and $\sim 75 \%$ for the Kroupa IMF.

We study the dark and luminous mass distributions, circular velocity curves (CVC), line-of-sight kinematics, and angular momenta for a sample of 42 cosmological zoom simulations of massive galaxies. Using a temporal smoothing technique, we are able to reach large radii. We find that: (i)The dark matter halo density profiles outside a few kpc follow simple power-law models, with flat dark matter CVCs for lower-mass systems, and rising CVCs for high-mass haloes. The projected stellar density distributions at large radii can be fitted by Sersic functions with n>10, larger than for typical ETGs. (ii)The massive systems have nearly flat total CVCs at large radii, while the less massive systems have mildly decreasing CVCs. The slope of the CVC at large radii correlates with v_circ itself. (iii)The dark matter fractions within Re are in the range 15-30% and increase to 40-65% at 5Re. Larger and more massive galaxies have higher dark matter fractions. (iv)The short axes of simulated galaxies and their host dark matter haloes are well aligned and their short-to-long axis ratios are correlated. (v)The stellar vrms(R) profiles are slowly declining, in agreement with planetary nebulae observations in the outer haloes of most ETGs. (vi)The line-of-sight velocity fields v show that rotation properties at small and large radii are correlated. Most radial profiles for the cumulative specific angular momentum parameter lambda(R) are nearly flat or slightly rising, with values in [0.06,0.75] from 2Re to 5Re. (vii)Stellar mass, ellipticity at 5Re, and lambda(5Re) are correlated: the more massive systems have less angular momentum and are rounder, as for observed ETGs. (viii)More massive galaxies with a large fraction of accreted stars have radially anisotropic velocity distributions outside Re. Tangential anisotropy is seen only for galaxies with high fraction of in-situ stars. (Full abstract in PDF)

The possibility to use observations of lensed relativistic jets on very small angular scales to construct proper models of spiral lens galaxies and to independently determine the Hubble constant is considered. On the example of the system B0218+357 it is shown that there exists a great choice of model parameters adequately reproducing its observed large-scale properties but leading to a significant spread in the Hubble constant. The jet image position angle is suggested as an additional parameter that allows the range of models under consideration to be limited. It is shown that the models for which the jet image position angles differ by at least 40 deg can be distinguished between themselves during observations on very small angular scales. Observations with the space interferometer RadioAstron-Earth is proposed to determine the jet orientation.

The accuracy and robustness of a simple method to estimate the total mass profile of a galaxy is tested using a sample of 65 cosmological zoom-simulations of individual galaxies. The method only requires information on the optical surface brightness and the projected velocity dispersion profiles and therefore can be applied even in case of poor observational data. In the simulated sample massive galaxies ($\sigma \simeq 200-400$ $\kms$) at redshift $z=0$ have almost isothermal rotation curves for broad range of radii (RMS $\simeq 5%$ for the circular speed deviations from a constant value over $0.5R_{\rm eff} < r < 3R_{\rm eff}$). For such galaxies the method recovers the unbiased value of the circular speed. The sample averaged deviation from the true circular speed is less than $\sim 1%$ with the scatter of $\simeq 5-8%$ (RMS) up to $R \simeq 5R_{\rm eff}$. Circular speed estimates of massive non-rotating simulated galaxies at higher redshifts ($z=1$ and $z=2$) are also almost unbiased and with the same scatter. For the least massive galaxies in the sample ($\sigma < 150$ $\kms$) at $z=0$ the RMS deviation is $\simeq 7-9%$ and the mean deviation is biased low by about $1-2%$. We also derive the circular velocity profile from the hydrostatic equilibrium (HE) equation for hot gas in the simulated galaxies. The accuracy of this estimate is about RMS $\simeq 4-5%$ for massive objects ($M > 6.5\times 10^{12} M_\odot$) and the HE estimate is biased low by $\simeq 3-4%$, which can be traced to the presence of gas motions. This implies that the simple mass estimate can be used to determine the mass of observed massive elliptical galaxies to an accuracy of $5-8 %$ and can be very useful for galaxy surveys.

The possibility of using lensed relativistic jets on very small angular scales to construct proper models of spiral lens galaxies and to independently determine the Hubble constant is considered. The system B0218+357 is used as an example to illustrate that there exists a great choice of model parameters adequately reproducing its observed large-scale properties but leading to a significant spread in the Hubble constant. The jet image position angle is suggested as an additional parameter that allows the range of models under consideration to be limited. It is shown that the models for which the jet image position angles differ by at least $40^o$ can be distinguished between themselves during observations on very small angular scales. The possibility of observing the geometric properties of lensed relativistic jets and measuring the superluminal velocities of knot images on time scales of several months with very long baseline space interferometers is discussed.

The images of relativistic jets from extragalactic sources produced by gravitational lensing by galaxies with different mass surface density distributions are modeled. In particular, the following models of the gravitational lens mass distribution are considered: a singular isothermal ellipsoid, an isothermal ellipsoid with a core, two- and three-component models with a galactic disk, halo, and bulge. The modeled images are compared both between themselves and with available observations. Different sets of parameters are shown to exist for the gravitationally lensed system B0218+357 in multicomponent models. These sets allow the observed geometry of the system and the intensity ratio of the compact core images to be obtained, but they lead to a significant variety in the Hubble constant determined from the modeling results.

The jets image modelling of gravitationally lensed sources have been performed. Several basic models of the lens mass distribution were considered, in particular, a singular isothermal ellipsoid, an isothermal ellipsoid with the core, different multi-components models with the galactic disk, halo and bulge. The obtained jet images were compared as with each other as with results of observations. A significant dependence of the Hubble constant on the model parameters was revealed for B0218+357, when the circular structure was took into account.