emcee is a Python library implementing a class of affine-invariant ensemble samplers for Markov chain Monte Carlo (MCMC). This package has been widely applied to probabilistic modeling problems in astrophysics where it was originally published, with some applications in other fields. When it was first released in 2012, the interface implemented in emcee was fundamentally different from the MCMC libraries that were popular at the time, such as PyMC, because it was specifically designed to work with "black box" models instead of structured graphical models. This has been a popular interface for applications in astrophysics because it is often non-trivial to implement realistic physics within the modeling frameworks required by other libraries. Since emcee's release, other libraries have been developed with similar interfaces, such as dynesty (Speagle 2019). The version 3.0 release of emcee is the first major release of the library in about 6 years and it includes a full re-write of the computational backend, several commonly requested features, and a set of new "move" implementations.

The largest Key Science Program of the RadioAstron space VLBI mission is a survey of active galactic nuclei (AGN). The main goal of the survey is to measure and study the brightness of AGN cores in order to better understand the physics of their emission while taking interstellar scattering into consideration. In this paper we present detection statistics for observations on ground-space baselines of a complete sample of radio-strong AGN at the wavelengths of 18, 6, and 1.3 cm. Two-thirds of them are indeed detected by RadioAstron and are found to contain extremely compact, tens to hundreds of $\mu$as structures within their cores.

Physical parameters of AGN jets observed with Very Long Baseline Interferometry (VLBI) are usually inferred from the core shift measurements or flux and size measured at a peak frequency of the synchrotron spectrum. Both are preceded by modelling of the observed VLBI jet structure with a simple Gaussian templates. We propose to infer the jets parameters using the inhomogeneous jet model directly - bypassing the modelling of the source structure with a Gaussian templates or image deconvolution. We applied Bayesian analysis to multi-frequency VLBA observations of radio galaxy NGC 315 and found that its parsec-scale jet is well described by the inhomogeneous conical model. Our results favour electron-positron jet. We also detected a component in a counter jet. Its position implies the presence of an external absorber with a steep density gradient at close ($r=0.1$ pc) distance from the central engine.

Synchrotron self-absorption in active galactic nuclei (AGN) jets manifests itself as a time delay between flares observed at high and low radio frequencies. It is also responsible for the observing frequency dependent change in size and position of the apparent base of the jet, aka the core shift effect, detected with very long baseline interferometry (VLBI). We measure the time delays and the core shifts in 11 radio-loud AGN to estimate the speed of their jets without relying on multi-epoch VLBI kinematics analysis. The 15$-$8 GHz total flux density time lags are obtained using Gaussian process regression, the core shift values are measured using VLBI observations and adopted from the literature. A strong correlation is found between the apparent core shift and the observed time delay. Our estimate of the jet speed is higher than the apparent speed of the fastest VLBI components by the median coefficient of 1.4. The coefficient ranges for individual sources from 0.5 to 20. We derive Doppler factors, Lorentz factors and viewing angles of the jets, as well as the corresponding de-projected distance from the jet base to the core. The results support evidence for acceleration of the jets with bulk motion Lorentz factor $\Gamma\propto R^{0.52\pm0.03}$ on de-projected scales $R$ of 0.5$-$500 parsecs.

Very Long Baseline Interferometric (VLBI) observations of quasar jets enable one to measure many theoretically expected effects. Estimating the significance of observational findings is complicated by the correlated noise in the image plane. A reliable and well justified approach to estimate the uncertainties of VLBI results is needed as well as significance testing criteria. We propose to use bootstrap for both tasks. Using simulations we find that bootstrap-based errors for the full intensity, rotation measure, and spectral index maps have coverage closer to the nominal values than conventionally obtained errors. The proposed method naturally takes into account heterogeneous interferometric arrays (such as Space VLBI) and can be easily extended to account for instrumental calibration factors.

Photometric variability detection is often considered as a hypothesis testing problem: an object is variable if the null-hypothesis that its brightness is constant can be ruled out given the measurements and their uncertainties. Uncorrected systematic errors limit the practical applicability of this approach to high-amplitude variability and well-behaving data sets. Searching for a new variability detection technique that would be applicable to a wide range of variability types while being robust to outliers and underestimated measurement uncertainties, we propose to consider variability detection as a classification problem that can be approached with machine learning. We compare several classification algorithms: Logistic Regression (LR), Support Vector Machines (SVM), k-Nearest Neighbors (kNN) Neural Nets (NN), Random Forests (RF) and Stochastic Gradient Boosting classifier (SGB) applied to 18 features (variability indices) quantifying scatter and/or correlation between points in a light curve. We use a subset of OGLE-II Large Magellanic Cloud (LMC) photometry (30265 light curves) that was searched for variability using traditional methods (168 known variable objects identified) as the training set and then apply the NN to a new test set of 31798 OGLE-II LMC light curves. Among 205 candidates selected in the test set, 178 are real variables, 13 low-amplitude variables are new discoveries. We find that the considered machine learning classifiers are more efficient (they find more variables and less false candidates) compared to traditional techniques that consider individual variability indices or their linear combination. The NN, SGB, SVM and RF show a higher efficiency compared to LR and kNN.

Clues to the physical conditions in radio cores of blazars come from measurements of brightness temperatures as well as effects produced by intrinsic opacity. We study the properties of the ultra compact blazar AO 0235+164 with RadioAstron ground-space radio interferometer, multi-frequency VLBA, EVN and single-dish radio observations. We employ visibility modeling and image stacking for deriving structure and kinematics of the source, and use Gaussian process regression to find the relative multi-band time delays of the flares. The multi-frequency core size and time lags support prevailing synchrotron self absorption. The intrinsic brightness temperature of the core derived from ground-based VLBI is close to the equipartition regime value. In the same time, there is evidence for ultra-compact features of the size of less than 10 $\mu$as in the source, which might be responsible for the extreme apparent brightness temperatures of up to $10^{14}$ K as measured by RadioAstron. In 2007--2016 the VLBI components in the source at 43 GHz are found predominantly in two directions, suggesting a bend of the outflow from southern to northern direction. The apparent opening angle of the jet seen in the stacked image at 43 GHz is two times wider than that at 15 GHz, indicating a collimation of the flow within the central 1.5 mas. We estimate the Lorentz factor $\Gamma = 14$, the Doppler factor $\delta=21$, and the viewing angle $\theta = 1.7^\circ$ of the apparent jet base, derive the gradients of magnetic field strength and electron density in the outflow, and the distance between jet apex and the core at each frequency.

Using standard cosmological model we show that the cores of ultra-compact radio sources observed with ground-based Very Long Baseline Interferometry (VLBI) on the angular scales of milliarcseconds cannot be used as a reasonable standard unit of linear size. "Luminosity - linear size" correlation obtained by many authors for ultra compact radio sources has different origin than that for the radio galaxies and quasars on the angular scales of arcminutes. It is just the manifestation of the fact that ground-based VLBI networks are unable to resolve VLBI-cores and the Malmquist bias presents. Thus, the cores of compact radio sources can't be used as "standard rods" at least with resolution offered by ground-based VLBI. This conclusion is illustrated on 15 GHz VLBA sample of radio sources.

Based on the miscellaneous published radio and optical data, SDSS and APM catalogue we consider the various properties of the giant radio sources (gRS) with the aim of refining the conditions leading to the formation of these objects. We compare gRSs with the regular-sized radio sources in radio and optical bands, yielding the following results: 1. The fraction of broad line objects among gRSs with high excitation spectrum is the same as for the RSs from isotropic samples. According to Unified Scheme this leads to the isotropic angle distribution of gRSs jets, thus gRSs cannot be characterized as objects with jets lying in the plane of sky. 2. gRSs do not differ from normal sized RSs in apparent asymmetry distribution of their extended radio components (ERC). However the fact that asymmetry distributions for gRSs and giant radio quasars (gQSS) are essentially the same leads within the Unification Scheme to the conclusion that the origin of this asymmetry is in the non-uniform environment. 3. The richness of the environment for gRSs is the same as for normal sized RSs. This contradicts the opinion that the low density of the environment is the single reason for gRSs formation. 4. About 10% of FRII RSs have by order of magnitude longer lifetimes and eventually evolve to gRSs. 5. The observed relative quantity of radio quasars in gRSs population (~0.1) can be interpreted as the presence of long-living population of radio loud QSSs as ~0.1 of all radio quasars.