Measuring the spin of the primary black hole in OJ287

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.2Comment: 12 pages, 6 figure

Testing black hole no-hair theorem with OJ287

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.Comment: 14 pages, 14 figure

An Evaluation of Body-grip Trap Trigger Configurations for Reducing River Otter Take Incidental to Beaver Trapping

River otter (Lontra canadensis) populations in North America have been the focus of significant restoration efforts. Wildlife management agencies, concerned about the unintentional take of river otters incidental to beaver (Castor canadensis) trapping, may recommend techniques to avoid capturing river otters. River otter avoidance techniques that are ineffective or diminish trap performance for beavers are undesirable. We conducted a field evaluation in 2015 and 2016 in Wisconsin to assess how two trigger configurations (offset and center) on body-grip traps would affect the incidental capture rate of river otters during beaver trapping. We also evaluated effects of each configuration on beaver capture rates, body lengths, and anatomical locations of trap-jaw strikes. We used size 330 body-grip traps equipped with identical triggers and alternated between trigger configurations during beaver damage management activities. We captured 8 river otters with each trap trigger configuration. Trap-jaw strikes on beavers differed between trigger configurations, with offset triggers resulting in more abdomen strikes and center triggers causing more cervical vertebrae strikes. We found that an offset trigger configuration did not reduce incidental take of otters and was less effective for trapping beavers

Optical polarization angle and VLBI jet direction in the binary black hole model of OJ287

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.Comment: to appear in Monthly Notices of Royal Astronomical Societ

The Nature of Starburst Activity in M82

We present new evolutionary synthesis models of M82 based mainly on observations consisting of near-infrared integral field spectroscopy and mid-infrared spectroscopy. The models incorporate stellar evolution, spectral synthesis, and photoionization modeling, and are optimized for 1-45 micron observations of starburst galaxies. The data allow us to model the starburst regions on scales as small as 25 pc. We investigate the initial mass function (IMF) of the stars and constrain quantitatively the spatial and temporal evolution of starburst activity in M82. We find a typical decay timescale for individual burst sites of a few million years. The data are consistent with the formation of very massive stars (> 50-100 Msun) and require a flattening of the starburst IMF below a few solar masses assuming a Salpeter slope at higher masses. Our results are well matched by a scenario in which the global starburst activity in M82 occurred in two successive episodes each lasting a few million years, peaking about 10 and 5 Myr ago. The first episode took place throughout the central regions of M82 and was particularly intense at the nucleus while the second episode occurred predominantly in a circumnuclear ring and along the stellar bar. We interpret this sequence as resulting from the gravitational interaction M82 and its neighbour M81, and subsequent bar-driven evolution. The short burst duration on all spatial scales indicates strong negative feedback effects of starburst activity, both locally and globally. Simple energetics considerations suggest the collective mechanical energy released by massive stars was able to rapidly inhibit star formation after the onset of each episode.Comment: 48 pages, incl. 16 Postscript figures; accepted for publication in the Astrophysical Journa

Measuring Black Hole Spin in OJ287

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).Comment: 12 pages, 4 figures, IAU26

Order and Chaos in a Three-dimensional Binary System of Interacting Galaxies

In the present article, we present a galactic gravitational model of three degrees of freedom, in order to investigate and reveal the behavior of orbits in a binary quasar system. The two quasars are hosted in a pair of interacting disk galaxies. We study in detail the regular or chaotic character of motion, in two different cases: the time independent model in both 2D and 3D dynamical systems and the evolving 3D model. Our numerical calculations, indicate that a large fraction of orbits in the 2D system are chaotic in the time independent case. A careful analysis suggest that several Lindblad resonances are also responsible for the chaotic motion of stars in both hosts galaxies. In the time dependent system, we follow the evolution of 3D orbits in our dynamical model, as the two interacting host galaxies develop dense and massive quasars in their cores, by mass transportation from the disks to their nuclei. In this interesting case, there are orbits that change their orbital character from regular to chaotic and vise versa and also orbits that maintain their character during the galactic evolution. These results strongly indicate that the ordered or chaotic nature of 3D orbits depends, not only in the galactic interaction but also in the presence of the quasars in the galactic cores of the host galaxies. The outcomes derived from our dynamical model are compared with observational data. Some theoretical arguments to support the numerically derived outcomes are presented, both in 2D and 3D systems, while a comparison with similar earlier work is also made.Comment: Published in the Astrophysical Journal (ApJ

Simulations of the grand design galaxy M51: a case study for analysing tidally induced spiral structure

We present hydrodynamical models of the grand design spiral M51 (NGC 5194), and its interaction with its companion NGC 5195. Despite the simplicity of our models, our simulations capture the present day spiral structure of M51 remarkably well, and even reproduce details such as a kink along one spiral arm, and spiral arm bifurcations. We investigate the offset between the stellar and gaseous spiral arms, and find at most times (including the present day) there is no offset between the stars and gas to within our error bars. We also compare our simulations with recent observational analysis of M51. We compute the pattern speed versus radius, and like the observations, find no single global pattern speed. We also show that the spiral arms cannot be fitted well by logarithmic spirals. We interpret these findings as evidence that M51 does not exhibit a quasi-steady density wave, as would be predicted by density wave theory. The internal structure of M51 derives from the complicated and dynamical interaction with its companion, resulting in spiral arms showing considerable structure in the form of short-lived kinks and bifurcations. Rather than trying to model such galaxies in terms of global spiral modes with fixed pattern speeds, it is more realistic to start from a picture in which the spiral arms, while not being simple material arms, are the result of tidally induced kinematic density `waves' or density patterns, which wind up slowly over time.Comment: 23 pages, 20 figures, accepted for publication in MNRA

A search for periodicity in the light curves of selected blazars

We present an analysis of multifrequency light curves of the sources 2223-052 (3C 446), 2230+114 (CTA 102), and 2251+158 (3C 454.3), which had shown evidence of quasi-periodic activity. The analysis made use of data from the University of Michican Radio Astronomy Observatory (USA) at 4.8, 8, and 14.5 GHz, as well as the Metsahovi Radio Astronomy Observatory (Finland) at 22 and 37 GHz. Application of two different methods (the discrete autocorrelation function and the method of Jurkevich) both revealed evidence for periodicity in the flux variations of these sources at essentially all frequencies. The periods derived for at least two of the sources -- 2223-052 and 2251+158-- are in good agreement with the time interval between the appearance of successive VLBI components. The derived periods for 2251+158 (P = 12.4 yr and 2223-052 (P = 5.8 yr) coincide with the periods found earlier by other authors based on optical light curves.Comment: 27 pages, 11 figures, accepted for publication in Astronomy Report

HST detection of spiral structure in two Coma Cluster dwarf galaxies

We report the discovery of (stellar) spiral-like structure in Hubble Space Telescope images of two dwarf galaxies (GMP 3292 and GMP 3629) belonging to the Coma cluster. GMP 3629 is the faintest such galaxy detected in a cluster environment, and it is the first such galaxy observed in the dense Coma cluster. The large bulge and the faintness of the broad spiral-like pattern in GMP 3629 suggests that its disk may have been largely depleted. >We may therefore have found an example of the ``missing link'' in theories of galaxy evolution which have predicted that dwarf spiral galaxies, particularly in clusters, evolve into dwarf elliptical galaxies.Comment: Accepted for publication in AJ on July