Binary black hole (BH) central engine description for the unique blazar OJ 287 predicted that the next secondary BH impact-induced bremsstrahlung flare should peak on 2019 July 31. This prediction was based on detailed general relativistic modeling of the secondary BH trajectory around the primary BH and its accretion disk. The expected flare was termed the Eddington flare to commemorate the centennial celebrations of now-famous solar eclipse observations to test general relativity by Sir Arthur Eddington. We analyze the multi-epoch Spitzer observations of the expected flare between 2019 July 31 and 2019 September 6, as well as baseline observations during 2019 February-March. Observed Spitzer flux density variations during the predicted outburst time display a strong similarity with the observed optical pericenter flare from OJ 287 during 2007 September. The predicted flare appears comparable to the 2007 flare after subtracting the expected higher base-level Spitzer flux densities at 3.55 and 4.49 $\mu$m compared to the optical R-band. Comparing the 2019 and 2007 outburst lightcurves and the previously calculated predictions, we find that the Eddington flare arrived within 4 hours of the predicted time. Our Spitzer observations are well consistent with the presence of a nano-Hertz gravitational wave emitting spinning massive binary BH that inspirals along a general relativistic eccentric orbit in OJ 287. These multi-epoch Spitzer observations provide a parametric constraint on the celebrated BH no-hair theorem.

Chaos is present in most stellar dynamical systems and manifests itself through the exponential growth of small perturbations. Exponential divergence drives time irreversibility and increases the entropy in the system. A numerical consequence is that integrations of the N-body problem unavoidably magnify truncation and rounding errors to macroscopic scales. Hitherto, a quantitative relation between chaos in stellar dynamical systems and the level of irreversibility remained undetermined. In this work we study chaotic three-body systems in free fall initially using the accurate and precise N-body code Brutus, which goes beyond standard double-precision arithmetic. We demonstrate that the fraction of irreversible solutions decreases as a power law with numerical accuracy. This can be derived from the distribution of amplification factors of small initial perturbations. Applying this result to systems consisting of three massive black holes with zero total angular momentum, we conclude that up to five percent of such triples would require an accuracy of smaller than the Planck length in order to produce a time-reversible solution, thus rendering them fundamentally unpredictable.

In the binary black hole model of OJ 287 the secondary black hole orbits a much more massive primary, and impacts on the primary accretion disk at predictable times. We update the parameters of the disk, the viscosity $\alpha$ and the mass accretion rate $\dot m$. We find $\alpha=0.26 \pm 0.1$ and $\dot m = 0.08 \pm 0.04$ in Eddington units. The former value is consistent with Coroniti (1981) and the latter with Marscher and Jorstad (2011). Predictions are made for the 2019 July 30 superflare in OJ 287. We expect that it will take place simultaneously at the Spitzer infrared channels as well as in the optical and that therefore the timing of the flare in optical can be accurately determined from Spitzer observations. We also discuss in detail the light curve of the 2015 flare and find that the radiating volume has regions where bremsstrahlung dominates as well as regions that radiate primarily in synchrotron radiation. The former region produces the unpolarised first flare while the latter region gives rise to a highly polarized second flare.

The bright blazar OJ 287 is the best-known candidate for hosting a nanohertz gravitational wave (GW) emitting supermassive binary black hole (SMBBH) in the present observable universe. The binary black hole (BBH) central engine model, proposed by Lehto and Valtonen in 1996, was influenced by the two distinct periodicities inferred from the optical light curve of OJ 287. The current improved model employs an accurate general relativistic description to track the trajectory of the secondary black hole (BH) which is crucial to predict the inherent impact flares of OJ 287. The successful observations of three predicted impact flares open up the possibility of using this BBH system to test general relativity in a hitherto unexplored strong field regime. Additionally, we briefly describe an on-going effort to interpret observations of OJ 287 in a Bayesian framework.

Blazar OJ287 exhibits large thermal flares at least twice every 12 years. The times of these flares have been predicted successfully using the model of a quasi-Keplerian eccentric black hole binary where the secondary impacts the accretion disk of the primary, creating the thermal flares. New measurements of the historical light curve have been combined with the observations of the 2015 November/December flare to identify the impact record since year 1886, and to constrain the orbit of the binary. The orbit solution shows that the binary period, now 12.062 year, is decreasing at the rate of 36 days per century. This corresponds to an energy loss to gravitational waves that is 6.5 +- 4 % less than the rate predicted by the standard quadrupolar gravitational wave (GW) emission. We show that the difference is due to higher order gravitational radiation reaction terms that include the dominant order tail contributions.

Results from regular monitoring of relativistic compact binaries like PSR 1913+16 are consistent with the dominant (quadrupole) order emission of gravitational waves (GWs). We show that observations associated with the binary black hole central engine of blazar OJ 287 demand the inclusion of gravitational radiation reaction effects beyond the quadrupolar order. It turns out that even the effects of certain hereditary contributions to GW emission are required to predict impact flare timings of OJ 287. We develop an approach that incorporates this effect into the binary black hole model for OJ~287. This allows us to demonstrate an excellent agreement between the observed impact flare timings and those predicted from ten orbital cycles of the binary black hole central engine model. The deduced rate of orbital period decay is nine orders of magnitude higher than the observed rate in PSR 1913+16, demonstrating again the relativistic nature of OJ 287's central engine. Finally, we argue that precise timing of the predicted 2019 impact flare should allow a test of the celebrated black hole "no-hair theorem" at the 10% level.

OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.

We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

The current cosmological paradigm sees the formation and evolution of the cosmic large-scale structure as governed by the gravitational attraction of the Dark Matter (DM) and the repulsion of the Dark Energy (DE). We characterize the relative importance of uniform and constant dark energy, as given by the Lambda term in the standard LCDM cosmology, in galaxy systems of different scales, from groups to superclusters. An instructive "Lambda significance graph" is introduced where the matter-DE density ratio <rho_M>/rho_Lambda for different galaxy systems is plotted against the radius R. This presents gravitation and DE dominated regions and shows directly the zero velocity radius, the zero-gravity radius, and the Einstein-Straus radius for any fixed value of mass. Example galaxy groups and clusters from the local universe illustrate the use of the Lambda significance graph. These are generally located deep in the gravity-dominated region <rho_M}>/rho_Lambda > 2, being virialized. Extended clusters and main bodies of superclusters can reach down near the border line between gravity-dominated and DE dominated regions <rho_M>/rho_Lambda = 2. The scale--mass relation from the standard 2-point correlation function intersects this balance line near the correlation lenght. The log <rho_M>/rho_Lambda vs. log R diagram is a useful and versatile way to characterize the dynamical state of systems of galaxies within the Lambda dominated expanding universe.

We have shown previously (Bobylev et al 2011) that some of the stars in the Solar neighborhood today may have originated in the same star cluster as the Sun, and could thus be called Solar Siblings. In this work we investigate the sensitivity of this result to Galactic models and to parameters of these models, and also extend the sample of orbits. There are a number of good candidates for the Sibling category, but due to the long period of orbit evolution since the break-up of the birth cluster of the Sun, one can only attach probabilities of membership. We find that up to 10% (but more likely around 1 %) of the members of the Sun's birth cluster could be still found within 100 pc from the Sun today.

Einstein's General Relativity theory has been tested in many ways during the last hundred years as reviewed in this chapter. Two tests are discussed in detail in this article: the concept of a zero gravity surface, the roots of which go back to J\"arnefelt, Einstein and Straus, and the no-hair theorem of black holes, first proposed by Israel, Carter and Hawking. The former tests the necessity of the cosmological constant Lambda, the latter the concept of a spinning black hole. The zero gravity surface is manifested most prominently in the motions of dwarf galaxies around the Local Group of galaxies. The no-hair theorem is testable for the first time in the binary black hole system OJ287. These represent stringent tests at the limit of weak and strong gravitational fields, respectively. In this article we discuss the current observational situation and future possibilities.

We used VIMOS on VLT to perform $V$ and $I$ band imaging of the outermost halo of NGC 5128 / Centaurus A ($(m-M)_0=27.91\pm0.08$), 65 kpc from the galaxy's center and along the major axis. The stellar population has been resolved to $I_0 \approx 27$ with a $50\%$ completeness limit of $I_0 = 24.7$, well below the tip of the red-giant branch (TRGB), which is seen at $I_0 \approx 23.9$. The surface density of NGC 5128 halo stars in our fields was sufficiently low that dim, unresolved background galaxies were a major contaminant in the source counts. We isolated a clean sample of red-giant-branch (RGB) stars extending to $\approx 0.8$ mag below the TRGB through conservative magnitude and color cuts, to remove the (predominantly blue) unresolved background galaxies. We derived stellar metallicities from colors of the stars via isochrones and measured the density falloff of the halo as a function of metallicity by combining our observations with HST imaging taken of NGC 5128 halo fields closer to the galaxy center. We found both metal-rich and metal-poor stellar populations and found that the falloff of the two follows the same de Vaucouleurs' law profiles from $\approx 8$ kpc out to $\approx$ 70 kpc. The metallicity distribution function (MDF) and the density falloff agree with the results of two recent studies of similar outermost halo fields in NGC 5128. We found no evidence of a "transition" in the radial profile of the halo, in which the metal-rich halo density would drop rapidly, leaving the underlying metal-poor halo to dominate by default out to greater radial extent, as has been seen in the outer halo of two other large galaxies. If NGC 5128 has such a transition, it must lie at larger galactocentric distances.

Context. OJ287 is a quasar with a quasi-periodic optical light curve, with the periodicity observed for over 120 years. This has lead to a binary black hole model as a common explanation of the quasar. The radio jet of OJ287 has been observed for a shorter time of about 30 years. It has a complicated structure that varies dramatically in a few years time scale. Aims. Here we propose that this structure arises from a helical jet being observed from a small and varying viewing angle. The viewing angle variation is taken to be in tune with the binary orbital motion. Methods. We calculate the effect of the secondary black hole on the inner edge of the accretion disk of the primary using particle simulations. We presume that the axis of the helix is perpendicular to the disk. We then follow the jet motion on its helical path and project the jet to the sky plane. This projection is compared with observations both at mm waves and cm waves. Results. We find that this model reproduces the observations well if the changes in the axis of the conical helix propagate outwards with a relativistic speed of about 0.85c. In particular, this model explains at the same time the long-term optical brightness variations as varying Doppler beaming in a component close to the core, i.e. at parsec scale in real linear distance, while the mm and cm radio jet observations are explained as being due to jet wobble at much larger (100 parsec scale) distances from the core.

We have studied short-term variations of the blazar OJ 287, suspected to host a supermassive black hole binary. In this study, we use a two-season optical R-band dataset from 2004--2006 which consists of 3991 data points from the OJ 287 observation campaign. It has sections of dense time coverage, and is largely independent from previously published data. We find that this data confirms the existence of a ~50 day periodic component, presumably related to the half-period of the innermost stable circular orbit (ISCO) of the primary black hole. In addition we find several pseudo-periodic components in the 1 to 7 day range, most prominently at 3.5 days, which are likely Lorentz contracted jet re-emission of the 50 day component. The typical 50 day cycle exhibits a slow rise of brightness and a rapid dimming before the start of the new cycle. We explain this as being due to a spiral wave in the accretion disc which feeds the central black hole in this manner.

{We consider the Coma cluster of galaxies as a gravitationally bound physical system embedded in the perfectly uniform static dark energy background as implied by the $\Lambda$CDM cosmology.} {We ask if the density of dark energy is high enough to affect the structure of a large rich cluster of galaxies?} {We use recent observational data on the cluster together with our theory of local dynamical effects of dark energy, including the zero-gravity radius $R_{\rm ZG}$ of the local force field as the key parameter.} {1) Three masses are defined which characterize the structure of a regular cluster: the matter mass $M_{\rm M}$, the dark-energy effective mass $M_{\rm DE}$ ($<0$) and the gravitating mass $M_{\rm G}$ ($= M_{\rm M} + M_{\rm DE}$). 2) A new matter density profile is suggested which reproduces well the observational data for the Coma cluster in the radius range from 1.4 Mpc to 14 Mpc and takes into account the dark energy background. 3) Using this profile, we calculate upper limits for the total size of the Coma cluster, $R \le R_{\rm ZG} \approx 20$ Mpc, and its total matter mass, $M_{\rm M} \la M_{\rm M}(R_{\rm ZG}) = 6.2 \times 10^{15} M_{\odot}$.} {The dark energy antigravity affects strongly the structure of the Coma cluster at large radii $R \ga 14$ Mpc and should be taken into account when its total mass is derived.}

We have studied three most recent precursor flares in the light curve of the blazar OJ 287 while invoking the presence of a precessing binary black hole in the system to explain the nature of these flares. Precursor flare timings from the historical light curves are compared with theoretical predictions from our model that incorporate effects of an accretion disk and post-Newtonian description for the binary black hole orbit. We find that the precursor flares coincide with the secondary black hole descending towards the accretion disk of the primary black hole from the observed side, with a mean z-component of approximately z_c = 4000 AU. We use this model of precursor flares to predict that precursor flare of similar nature should happen around 2020.96 before the next major outburst in 2022.

The blazar OJ287 is the most promising (and the only) case for an extragalactic binary black hole system inspiralling under the action of gravitational radiation reaction. At present, though it is not possible to directly observe the binary components, it is possible to observe the jet emanating form the primary black hole. We argue that the orbital motion of the secondary black hole is reflected in the wobble of the jet and demonstrate that the wobble is orbital position dependent. The erratic wobble of the jet, reported in Agudo et al. (2012), is analyzed by taking into account the binary nature of the system and we find that the erratic component of jet wobble is very small.

Two multifrequency campaigns were carried out on OJ287 in 2005: in April when it was in its pre-outburst state, and in November, during the main 12 yr cycle outburst. The wavelength coverage was from radio to X-rays. In the optical-to-UV range the differential spectrum between the observations has a bremsstrahlung spectral shape, consistent with gas at $3 \times 10^{5}K$ temperature. Our result supports the hydrogen column density of the OJ287 host galaxy of $\sim9.3\times 10^{20} cm^{-2}$, the average value found by Gosh & Soundararajaperumal. The $3 \times 10^{5}K$ bremsstrahlung radiation was predicted in the binary black hole model of OJ287, and it arises from a hot bubble of gas which is torn off the accretion disc by the impact of the secondary. As this radiation is not Doppler boosted, the brightness of the outburst provides an estimate for the mass of the secondary black hole, $\sim1.4\times10^{8}$ solar mass. In order to estimate the mass of the primary black hole, we ask what is the minimum mass ratio in a binary system which allows the stability of the accretion disc. By using particle simulations, we find that the ratio is $\sim1.3\times10^{2}$. This makes the minimum mass of the primary $\sim1.8\times10^{10}$ solar mass, in agreement with the mass determined from the orbit solution, $1.84 \times 10^{10}$ solar mass. With this mass value and the measured K-magnitude of the bulge of the host galaxy of OJ287, the system lies almost exactly on the previously established correlation in the black hole mass vs. K-magnitude diagramme. It supports the extension of this correlation to brighter magnitudes and to more massive black holes than has been done previously.

The light curve of the quasar OJ 287 extends from 1891 up today without major gaps. Here we summarize the results of the 2005 - 2010 observing campaign. The main results are the following: (1) The 2005 October optical outburst came at the expected time, thus confirming the general relativistic precession in the binary black hole system. This result disproved the model of a single black hole system with accretion disk oscillations, as well as several toy models of binaries without relativistic precession. In the latter models the main outburst would have been a year later. (2) The nature of the radiation of the 2005 October outburst was expected to be bremsstrahlung from hot gas at the temperature of $3\times 10^{5}$ $^{\circ}$K. This was confirmed by combined ground based and ultraviolet observations using the XMM-Newton X-ray telescope. (3) A secondary outburst of the same nature was expected at 2007 September 13. Within the accuracy of observations (about 6 hours), it started at the correct time. Thus the prediction was accurate at the same level as the prediction of the return of Halley's comet in 1986. (4) Further synchrotron outbursts were expected following the two bremsstrahlung outbursts. They came as scheduled between 2007 October and 2009 December. (5) Due to the effect of the secondary on the overall direction of the jet, the parsec scale jet was expected to rotate in the sky by a large angle around and after 2009. This rotation may have been seen at high frequency radio observations. OJ 287 binary black hole system is currently our best laboratory for testing theories of gravitation. Using OJ 287, the correctness of General Relativity has now been demonstrated at higher order than has been possible using the binary pulsars.

The light curve of quasar OJ287 extends from 1891 up today without major gaps. This is partly due to extensive studies of historical plate archives by Rene Hudec and associates, and partly due to several observing campaigns in recent times. Here we summarize the results of the 2005 - 2010 observing campaign, in which several hundred scientists and amateur astronomers took part. The main results are the following: (1) The 2005 October optical outburst came at the expected time, thus confirming the General Relativistic precession in the binary black hole system. At the same time, this result disproved the model of a single black hole system with accretion disk oscillations, as well as several toy models of binaries without relativistic precession. (2) The nature of the radiation of the 2005 October outburst was expected to be bremsstrahlung from hot gas at a temperature of 3 10^5 degrees K. This was confirmed by combined ground based and ultraviolet observations using the XMM-Newton X-ray telescope. (3) A secondary outburst of the same nature was expected at 2007 September 13. Within the accuracy of the observations (about 6 hours), it started at the correct time. Due to the bremsstrahlung nature of the outburst, the radiation was unpolarized, as expected. (4) Further synchrotron outbursts were expected following the two bremsstrahlung outbursts. They came as scheduled between 2007 October and 2009 December. (5) Due to the effect of the secondary on the overall direction of the jet, the parsec scale jet is expected to rotate in the sky by a large angle, which has been confirmed. The OJ287 binary black hole system is currently our best laboratory for testing theories of gravitation. Using OJ287, the correctness of General Relativity has now been demonstrated up to the second Post-Newtonian order, higher than has been possible using binary pulsars.

We study the variation of the optical polarization angle in the blazar OJ287 and compare it with the precessing binary black hole model with a 'live' accretion disk. First, a model of the variation of the jet direction is calculated, and the main parameters of the model are fixed by the long term optical brightness evolution. Then this model is compared with the variation of the parsec scale radio jet position angle in the sky. Finally, the variation of the polarization angle is calculated using the same model, but using a magnetic field configuration which is at a constant angle relative to the optical jet. It is found that the model fits the data reasonably well if the field is almost parallel to the jet axis. This may imply a steady magnetic field geometry, such as a large-scale helical field.

A new estimator of the local density of dark energy is suggested which comes from the virial theorem for non-relativistic gravitating systems embedded in the uniform dark energy background.

We examine the ability to test the black hole no-hair theorem at the 10% level in this decade using the binary black hole in OJ287. In the test we constrain the value of the dimensionless parameter q that relates the scaled quadrupole moment and spin of the primary black hole: q2 = -q 2 . At the present we can say that q = 1 \pm 0.3 (one), in agreement with General Relativity and the no-hair theorems. We demonstrate that this result can be improved if more observational data is found in historical plate archives for the 1959 and 1971 outbursts. We also show that the predicted 2015 and 2019 outbursts will be crucial in improving the accuracy of the test. Space-based photometry is required in 2019 July due the proximity of OJ287 to the Sun at the time of the outburst. The best situation would be to carry out the photometry far from the Earth, from quite a different vantage point, in order to avoid the influence of the nearby Sun. We have considered in particular the STEREO space mission which would be ideal if it has a continuation in 2019 or LORRI on board the New Horizons mission to Pluto.

We propose a kinematic approach to searching for the stars that could be formed with the Sun in a common "parent" open cluster. The approach consists in preselecting suitable candidates by the closeness of their space velocities to the solar velocity and analyzing the parameters of their encounters with the solar orbit in the past in a time interval comparable to the lifetime of stars. We consider stars from the Hipparcos catalog with available radial velocities. The Galactic orbits of stars have been constructed in the Allen--Santillan potential by taking into account the perturbations from the spiral density wave. We show that two stars, HIP 87382 and HIP 47399, are of considerable interest in our problem. Their orbits oscillate near the solar orbit with an amplitude of about 250 pc; there are short-term close encounters to distances <10 pc. Both stars have an evolutionary status and metallicity similar to the solar ones.

We explore the origin of a population of distant companions (~1000 - 5000 AU) to Class I protostellar sources recently found by Connelley and co-workers, who noted that the companion fraction diminished as the sources evolved. Here we present N-body simulations of unstable triple systems embedded in dense cloud cores. Many companions are ejected into unbound orbits and quickly escape, but others are ejected with insufficient momentum to climb out of the potential well of the cloud core and associated binary. These loosely bound companions reach distances of many thousands of AU before falling back and eventually being ejected into escapes as the cloud cores gradually disappear. We use the term orphans to denote protostellar objects that are dynamically ejected from their placental cloud cores, either escaping or for a time being tenuously bound at large separations. Half of all triple systems are found to disintegrate during the protostellar stage, so if multiple systems are a frequent outcome of the collapse of a cloud core, then orphans should be common. Bound orphans are associated with embedded close protostellar binaries, but escaping orphans can travel as far as ~0.2 pc during the protostellar phase. The steep climb out of a potential well ensures that orphans are not kinematically distinct from young stars born with a less violent pre-history. The identification of orphans outside their heavily extincted cloud cores will allow the detailed study of protostars high up on their Hayashi tracks at near-infrared and in some cases even at optical wavelengths.

The standard \LambdaCDM cosmological model implies that all celestial bodies are embedded in a perfectly uniform dark energy background, represented by Einstein's cosmological constant, and experience its repulsive antigravity action. Can dark energy have strong dynamical effects on small cosmic scales as well as globally? Continuing our efforts to clarify this question, we focus now on the Virgo Cluster and the flow of expansion around it. We interpret the Hubble diagram, from a new database of velocities and distances of galaxies in the cluster and its environment, using a nonlinear analytical model which incorporates the antigravity force in terms of Newtonian mechanics. The key parameter is the zero-gravity radius, the distance at which gravity and antigravity are in balance. Our conclusions are: 1. The interplay between the gravity of the cluster and the antigravity of the dark energy background determines the kinematical structure of the system and controls its evolution. 2. The gravity dominates the quasi-stationary bound cluster, while the antigravity controls the Virgocentric flow, bringing order and regularity to the flow, which reaches linearity and the global Hubble rate at distances \ga 15 Mpc. 3. The cluster and the flow form a system similar to the Local Group and its outflow. In the velocity-distance diagram, the cluster-flow structure reproduces the group-flow structure with a scaling factor of about 10; the zero-gravity radius for the cluster system is also 10 times larger. The phase and dynamical similarity of the systems on the scales of 1-30 Mpc suggests that a two-component pattern may be universal for groups and clusters: a quasi-stationary bound central component and an expanding outflow around it, due to the nonlinear gravity-antigravity interplay with the dark energy dominating in the flow component.

We model the binary black hole system OJ287 as a spinning primary and a non-spinning secondary. It is assumed that the primary has an accretion disk which is impacted by the secondary at specific times. These times are identified as major outbursts in the light curve of OJ287. This identification allows an exact solution of the orbit, with very tight error limits. Nine outbursts from both the historical photographic records as well as from recent photometric measurements have been used as fixed points of the solution: 1913, 1947, 1957, 1973, 1983, 1984, 1995, 2005 and 2007 outbursts. This allows the determination of eight parameters of the orbit. Most interesting of these are the primary mass of $1.84\cdot 10^{10} M_\odot$, the secondary mass $1.46\cdot 10^{8} M_\odot$, major axis precession rate $39^\circ.1$ per period, and the eccentricity of the orbit 0.70. The dimensionless spin parameter is $0.28\:\pm\:0.01$ (1 sigma). The last parameter will be more tightly constrained in 2015 when the next outburst is due. The outburst should begin on 15 December 2015 if the spin value is in the middle of this range, on 3 January 2016 if the spin is 0.25, and on 26 November 2015 if the spin is 0.31. We have also tested the possibility that the quadrupole term in the Post Newtonian equations of motion does not exactly follow Einstein's theory: a parameter $q$ is introduced as one of the 8 parameters. Its value is within 30% (1 sigma) of the Einstein's value $q = 1$. This supports the $no-hair theorem$ of black holes within the achievable precision. We have also measured the loss of orbital energy due to gravitational waves. The loss rate is found to agree with Einstein's value with the accuracy of 2% (1 sigma).

The compact binary system in OJ287 is modelled to contain a spinning primary black hole with an accretion disk and a non-spinning secondary black hole. Using Post Newtonian (PN) accurate equations that include 2.5PN accurate non-spinning contributions, the leading order general relativistic and classical spin-orbit terms, the orbit of the binary black hole in OJ287 is calculated and as expected it depends on the spin of the primary black hole. Using the orbital solution, the specific times when the orbit of the secondary crosses the accretion disk of the primary are evaluated such that the record of observed outbursts from 1913 up to 2007 is reproduced. The timings of the outbursts are quite sensitive to the spin value. In order to reproduce all the known outbursts, including a newly discovered one in 1957, the Kerr parameter of the primary has to be $0.28 \pm 0.08$. The quadrupole-moment contributions to the equations of motion allow us to constrain the `no-hair' parameter to be $1.0\:\pm\:0.3$ where 0.3 is the one sigma error. This supports the `black hole no-hair theorem' within the achievable precision. It should be possible to test the present estimate in 2015 when the next outburst is due. The timing of the 2015 outburst is a strong function of the spin: if the spin is 0.36 of the maximal value allowed in general relativity, the outburst begins in early November 2015, while the same event starts in the end of January 2016 if the spin is 0.2

Supermassive black holes are common in centers of galaxies. Among the active galaxies, quasars are the most extreme, and their black hole masses range as high as to $6\cdot 10^{10} M_\odot$. Binary black holes are of special interest but so far OJ287 is the only confirmed case with known orbital elements. In OJ287, the binary nature is confirmed by periodic radiation pulses. The period is twelve years with two pulses per period. The last four pulses have been correctly predicted with the accuracy of few weeks, the latest in 2007 with the accuracy of one day. This accuracy is high enough that one may test the higher order terms in the Post Newtonian approximation to General Relativity. The precession rate per period is $39^\circ.1 \pm 0^\circ.1$, by far the largest rate in any known binary, and the $(1.83\pm 0.01)\cdot 10^{10} M_\odot$ primary is among the dozen biggest black holes known. We will discuss the various Post Newtonian terms and their effect on the orbit solution. The over 100 year data base of optical variations in OJ287 puts limits on these terms and thus tests the ability of Einstein's General Relativity to describe, for the first time, dynamic binary black hole spacetime in the strong field regime. The quadrupole-moment contributions to the equations of motion allows us to constrain the `no-hair' parameter to be $1.0\:\pm\:0.3$ which supports the black hole no-hair theorem within the achievable precision.

Dark energy must be taken into account to estimate more reliably the amount of dark matter and how it is distributed in the local universe. For systems several Mpc across like the Local Group, we introduce three self-consistent independent mass estimators. These account for the antigravity effect of dark energy treated as Einstein's cosmological constant Lambda. The first is a modified Kahn-Woltjer model which gives a value of the Local Group mass via the particular motions of the two largest members, the Milky Way and M31. Inclusion of dark energy in this model increases the minimum mass estimate by a factor of three compared to the "classical estimate". The increase is less but still significant for different ways of using the timing argument. The second estimator is a modified virial theorem which also demonstrates how dark energy can "hide" from detection a part of the gravitating mass of the system. The third is a new zero-gravity method which gives an upper limit to the group mass which we calculate with high precision HST observations. In combination, the estimators lead to a robust and rather narrow range for a group's mass, M. For the Local Group, 3.2 < M < 3.7 x 10^{12} M_sun. Our result agrees well with the Millennium Simulation based on the LambdaCDM cosmology.

Several authors have found a statistically significant excess of galaxies with higher redshifts relative to the group centre, so-called discordant redshifts, in particular in groups where the brightest galaxy, identified in apparent magnitudes, is a spiral. Our aim is to explain the observed redshift excess. We use a semi-analytical galaxy catalogue constructed from the Millennium Simulation to study redshift asymmetries in spiral-dominated groups in the Lambda cold dark matter (LambdaCDM) cosmology. We show that discordant redshifts in small galaxy groups arise when these groups are gravitationally unbound and the dominant galaxy of the group is misidentified. The redshift excess is especially significant when the apparently brightest galaxy can be identified as a spiral, in full agreement with observations. On the other hand, the groups that are gravitationally bound do not show a significant redshift asymmetry. When the dominant members of groups in mock catalogues are identified by using the absolute B-band magnitudes, our results show a small blueshift excess. This result is due to the magnitude limited observations that miss the faint background galaxies in groups. When the group centre is not correctly identified it may cause the major part of the observed redshift excess. If the group is also gravitationally unbound, the level of the redshift excess becomes as high as in observations. There is no need to introduce any "anomalous" redshift mechanism to explain the observed redshift excess. Further, as the Friends-of-Friends percolation algorithm picks out the expanding parts of groups, in addition to the gravitationally bound group cores, group catalogues constructed in this way cannot be used as if the groups are purely bound systems.

Tests of Einstein's general theory of relativity have mostly been carried out in weak gravitational fields where the space-time curvature effects are first-order deviations from Newton's theory. Binary pulsars provide a means of probing the strong gravitational field around a neutron star, but strong-field effects may be best tested in systems containing black holes. Here we report such a test in a close binary system of two candidate black holes in the quasar OJ287. This quasar shows quasi-periodic optical outbursts at 12 yr intervals, with two outburst peaks per interval. The latest outburst occurred in September 2007, within a day of the time predicted by the binary black-hole model and general relativity. The observations confirm the binary nature of the system and also provide evidence for the loss of orbital energy in agreement (within 10 per cent) with the emission of gravitational waves from the system. In the absence of gravitational wave emission the outburst would have happened twenty days later.

We investigate whether it is possible that viable microbes could have been transported to Earth from the planets in extra-solar systems by means of natural vehicles such as ejecta expelled by comet or asteroid impacts on such planets. The probabilities of close encounters with other solar systems are taken into account as well as the limitations of bacterial survival times inside ejecta in space, caused by radiation and DNA decay. The conclusion is that no potentially DNA/RNA life-carrying ejecta from another solar system in the general Galactic star field landed on Earth before life already existed on Earth, not even if microbial survival time in space is as long as tens of millions of years. However, if the Sun formed initially as a part of a star cluster, as is commonly assumed, we cannot rule out the possibility of transfer of life from one of the sister systems to us. Likewise, there is a possibility that some extra-solar planets carry life that originated in our solar system. It will be of great interest to identify the members of the Sun's birth cluster of stars and study them for evidence for planets and life on the planets. The former step may be accomplished by the GAIA mission, the latter step by the SIM and DARWIN missions. Therefore it may not be too long until we have experimental knowledge on the question whether the natural transfer of life from one solar system to another has actually taken place.

In this paper we carry out a quantitative analysis of the three-body systems and map them as a function of decaying time and intial conguration, look at this problem as an example of a simple deterministic system, and ask to what extent the orbits are really predictable. We have investigated the behavior of about 200 000 general Newtonian three body systems using the simplest initial conditions. Within our resolution these cover all the possible states where the objects are initially at rest and have no angular momentum. We have determined the decay time-scales of the triple systems and show that the distribution of this parameter is fractal in appearance. Some areas that appear stable on large scales exhibit very narrow strips of instability and the overall pattern, dominated by resonances, reminds us of a traditional Maasai warrior shield. Also an attempt is made to recover the original starting conguration of the three bodies by backward integration. We find there are instances where the evolution to the future and to the past lead to different orbits, in spite of time symmetric initial conditions. This implies that even in simple deterministic systems there exists an Arrow of Time.

We have compared numerical simulations to observations for the nearby (< 40 Mpc) groups of galaxies (Huchra & Geller 1982 and Ramella et al. 2002). The group identification is carried out using a group-finding algorithm developed by Huchra and Geller (1982). Using cosmological N-body simulation code with the LambdaCDM cosmology, we show that the dynamical properties of groups of galaxies identified from the simulation data are, in general, in a moderate, within 2sigma, agreement with the observational catalogues of groups of galaxies. As simulations offer more dynamical information than observations, we used the N-body simulation data to calculate whether the nearby groups of galaxies are gravitationally bound objects by using their virial ratio. We show that in a LambdaCDM cosmology about 20 per cent of nearby groups of galaxies, identified by the same algorithm as in the case of observations, are not bound, but merely groups in a visual sense. This is quite significant, specifically because estimations of group masses in observations are often based on an assumption that groups of galaxies found by the friends-of-friends algorithm are gravitationally bound objects. Simulations with different resolutions show the same results. We also show how the fraction of gravitationally unbound groups varies when the apparent magnitude limit of the sample and the value of the cosmological constant is changed. In general, a larger value of the Omega_Lambda generates slightly more unbound groups.

The Hubble Space Telescope observations of the nearby galaxy group M 81/M 82 and its vicinity indicate that the expansion outflow around the group is dominated by the antigravity of the dark energy background. The local density of dark energy in the area is estimated to be near the global dark energy density or perhaps exactly equal to it. This conclusion agrees with our previous results for the Local group vicinity and the vicinity of the Cen A/M 83 group.

We report the detection of dark energy near the Milky Way made with precision observations of the local Hubble flow of expansion. We estimate the local density of dark energy and find that it is near, if not exactly equal to, the global dark energy density. The result is independent of, compatible with, and complementary to the horizon-scale observations in which dark energy was first discovered. Together with the cosmological concordance data, our result forms direct observational evidence for the Einstein antigravity as a universal phenomenon -- in the same sense as the Newtonian universal gravity.

A structure with a massive group in its center and a cool expansion outflow outside is studied around the Cen A galaxy with the use of the Hubble Space Telescope observations. It is demonstrated that the dynamics of the flow is dominated by the antigravity of the dark energy background. The density of dark energy in the cell is estimated to be near the global cosmological density. This agrees with our previous result from the neighborhood of the Local group. A notion of the ``Hubble cell'' is introduced as a building block of the local structure of the universe.

The physical nature of the very local (<3 Mpc) Hubble flow is studied on the basis of the recent high precision observations in the Local Volume. A model including both analytical treatment and computer simulations describes the flow dynamical evolution from a chaotic Little Bang initial state to the present-day state of a quasi-regular expansion. The two major observational parameters of the flow, which are the rate of expansion and the velocity dispersion, find a clear qualitative and quantitative explanation in the model.

The phenomenon of the very local ($\le3$ Mpc) Hubble flow is studied on the basis of the data of recent precision observations. A set of computer simulations is performed to trace the trajectories of the flow galaxies back in time to the epoch of the formation of the Local Group. It is found that the `initial conditions' of the flow are drastically different from the linear velocity-distance relation. The simulations enable also to recognize the major trends of the flow evolution and identify the dynamical role of universal antigravity produced by cosmic vacuum.

We obtained B and R-band CCD images for the dwarf elliptical (dE) galaxies DDO44, UGC4998, KK98_77, DDO71, DDO113, and UGC7356 at the NOT. Using Fourier analysis technique we measure stellar R-band surface brightness fluctuations (SBFs) and magnitudes in 29 different fields of the galaxies. Independent tip of the red giant branch distances for DDO44, KK98_77, DDO71 are used to convert their set of apparent into absolute SBF magnitudes. The results are combined with the local (B-R) colours and compared with the (B-R)-\bar{M}_R relation for mainly old, metal-poor stellar populations as predicted by Worthey's population synthesis models using Padova isochrones. While the colour dependency of the theoretical relation is confirmed by the empirical data, we find a systematic zero point offset between observations and theory in the sense that models are too faint by 0.13+-0.02 mag. Based on these findings we establish a new semiempirical calibration of the SBF method as distance indicator for dE galaxies with an estimated uncertainty of ~10 percent. Taking first advantage of the improved calibration, we determine SBF distances for the other three early-type dwarfs. Although found in the M81 group region, previous velocity measurements suggested UGC4998 is in the background. This picture is confirmed by our SBF distance of 10.5+-0.9 Mpc. We can further identify DDO113 as a faint stellar system at the near side of the Canes Venatici I (CVn I) cloud at a distance of 3.1+-0.3 Mpc. The second CVn I member in our sample, UGC7356, lies at 6.7+-0.6 Mpc and spatially close to M106. We derive BR surface brightness profiles and colour gradients for all dwarfs and determine photometric and Sersic parameters. Finally, we discuss two objects in DDO71 and UGC7356 which may resemble GCs.