We homogeneously reanalyse $124$ transit light curves for the WASP-4 b hot Jupiter. This set involved new observations secured in 2019 and nearly all observations mentioned in the literature, including high-accuracy GEMINI/GMOS transmission spectroscopy of 2011-2014 and TESS observations of 2018. The analysis confirmed a nonlinear TTV trend with $P/|\dot P|\sim (17-30)$ Myr (1-sigma range), implying only half of the initial decay rate estimation. The trend significance is at least $3.4$-sigma in the agressively conservative treatment. Possible radial acceleration due to unseen companions is not revealed in Doppler data covering seven years 2007-2014, and radial acceleration of $-15$ m s$^{-1}$yr$^{-1}$ reported in a recent preprint by another team is not confirmed. If present, it is a very nonlinear RV variation. Assuming that the entire TTV is tidal in nature, the tidal quality factor $Q_\star'\sim (4.5-8.5)\cdot 10^4$ does not reveal a convincing disagreement with available theory predictions.

We homogeneously analyse $\sim 3.2\times 10^5$ photometric measurements for $\sim 1100$ transit lightcurves belonging to $17$ exoplanet hosts. The photometric data cover $16$ years 2004--2019 and include amateur and professional observations. Old archival lightcurves were reprocessed using up-to-date exoplanetary parameters and empirically debiased limb-darkening models. We also derive self-consistent transit and radial-velocity fits for $13$ targets. We confirm the nonlinear TTV trend in the WASP-12 data at a high significance, and with a consistent magnitude. However, Doppler data reveal hints of a radial acceleration about $(-7.5\pm 2.2)$~m/s/yr, indicating the presence of unseen distant companions, and suggesting that roughly $10$ per cent of the observed TTV was induced via the light-travel (or Roemer) effect. For WASP-4, a similar TTV trend suspected after the recent TESS observations appears controversial and model-dependent. It is not supported by our homogeneus TTV sample, including $10$ ground-based EXPANSION lightcurves obtained in 2018 simultaneously with TESS. Even if the TTV trend itself does exist in WASP-4, its magnitude and tidal nature are uncertain. Doppler data cannot entirely rule out the Roemer effect induced by possible distant companions.

We present a self-consistent framework to perform the wavelet analysis of two-dimensional statistical distributions. The analysis targets the 2D probability density function (p.d.f.) of an input sample, in which each object is characterized by two peer parameters. The method performs a probabilistic detection of various `patterns', or `structures' related to the behaviour of the p.d.f. Laplacian. These patterns may include regions of local convexity or local concavity of the p.d.f., in particular peaks (groups of objects) or gaps. In the end, the p.d.f. itself is reconstructed based on the least noisy (most economic) superposition of such patterns. Among other things, our method also involves optimal minimum-noise wavelets and minimum-noise reconstruction of the distribution density function. The new 2D algorithm is now implemented and released along with an improved and optimized 1D version. The code relies on the C++11 language standard, and is fully parallelized. The algorithm has a rich range of applications in astronomy: Milky Way stellar population analysis, investigations of the exoplanets diversity, Solar System minor bodies statistics, etc.

We investigate 1D exoplanetary distributions using a novel analysis algorithm based on the continuous wavelet transform. The analysis pipeline includes an estimation of the wavelet transform of the probability density function (p.d.f.) without pre-binning, use of optimized wavelets, a rigorous significance testing of the patterns revealed in the p.d.f., and an optimized minimum-noise reconstruction of the p.d.f. via matching pursuit iterations. In the distribution of orbital periods, $P$, our analysis revealed a narrow subfamily of exoplanets within the broad family of "warm jupiters", or massive giants with $P\gtrsim 300$~d, which are often deemed to be related with the iceline accumulation in a protoplanetary disk. We detected a p.d.f. pattern that represents an upturn followed by an overshooting peak spanning $P\sim 300-600$~d, right beyond the "period valley". It is separated from the other planets by p.d.f. concavities from both sides. It has at least two-sigma significance. In the distribution of planet radii, $R$, and using the California Kepler Survey sample properly cleaned, we confirm the hints of a bimodality with two peaks about $R=1.3 R_\oplus$ and $R=2.4 R_\oplus$, and the "evaporation valley" between them. However, we obtain just a modest significance for this pattern, two-sigma only at the best. Besides, our follow-up application of the Hartigan & Hartigan dip test for unimodality returns $3$ per cent false alarm probability (merely $2.2$-sigma significance), contrary to $0.14$ per cent (or $3.2$-sigma), as claimed by Fulton et al. (2017).

We present an attempt to improve models of the Rossiter-McLaughlin effect by relaxing several restrictive assumptions. We consider the entire multiline stellar spectrum rather than just a single line, use no assumptions about the shape of the lines profiles, and allow arbitrary size ratio for the star and its eclipser. However, we neglect the effect of macro-turbulence and differential rotation. We construct our model as a power series in the stellar rotation velocity, $V\sin i$, giving a closed set of analytic formulae for up to three terms, and assuming quadratic limb-darkening law. We consider three major approaches of determining the Doppler shift: cross-correlation with a predefined template, cross-correlation with an out-of-transit stellar spectrum, and parametric modelling of the spectrum. A numerical testcase revels that our model preserves good accuracy for the rotation velocity of up to the limit of $2-3$ times the average linewidth in the spectrum. We also apply our approach to the Doppler data of HD 189733, for which we obtain an improved model of the Rossiter-McLaughlin effect with two correction terms, and derive a reduced value for $V\sin i$.

We perform an analysis of ~80000 photometric measurements for the following 10 stars hosting transiting planets: WASP-2, -4, -5, -52, Kelt-1, CoRoT-2, XO-2, TrES-1, HD 189733, GJ 436. Our analysis includes mainly transit lightcurves from the Exoplanet Transit Database, public photometry from the literature, and some proprietary photometry privately supplied by other authors. Half of these lightcurves were obtained by amateurs. From this photometry we derive 306 transit timing measurements, as well as improved planetary transit parameters. Additionally, for 6 of these 10 stars we present a set of radial velocity measurements obtained from the spectra stored in the HARPS, HARPS-N, and SOPHIE archives using the HARPS-TERRA pipeline. Our analysis of these TTV and RV data did not reveal significant hints of additional orbiting bodies in almost all of the cases. In the WASP-4 case, we found hints of marginally significant TTV signals having amplitude 10-20 sec, although their parameters are model-dependent and uncertain, while radial velocities did not reveal statistically significant Doppler signals.