Source author record

Pauli Pihajoki

Pauli Pihajoki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2020arXiv

MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates

We present the novel algorithmically regularised integration method MSTAR for high accuracy ($|ΔE/E| \gtrsim 10^{-14}$) integrations of N-body systems using minimum spanning tree coordinates. The two-fold parallelisation of the $\mathcal{O}(N_\mathrm{part}^2)$ force loops and the substep divisions of the extrapolation method allows for a parallel scaling up to $N_\mathrm{CPU} = 0.2 \times N_\mathrm{part}$. The efficient parallel scaling of MSTAR makes the accurate integration of much larger particle numbers possible compared to the traditional algorithmic regularisation chain (AR-CHAIN) methods, e.g. $N_\mathrm{part} = 5000$ particles on $400$ CPUs for $1$ Gyr in a few weeks of wall-clock time. We present applications of MSTAR on few particle systems, studying the Kozai mechanism and N-body systems like star clusters with up to $N_\mathrm{part} =10^4$ particles. Combined with a tree or a fast multipole based integrator the high performance of MSTAR removes a major computational bottleneck in simulations with regularised subsystems. It will enable the next generation galactic-scale simulations with up to $10^9$ stellar particles (e.g. $m_\star = 100 M_\odot$ for a $M_\star = 10^{11} M_\odot$ galaxy) including accurate collisional dynamics in the vicinity of nuclear supermassive black holes.

preprint2020arXiv

Spitzer Observations of the Predicted Eddington Flare from Blazar OJ 287

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 $μ$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.

preprint2016arXiv

Black hole accretion disc impacts

We present an analytic model for computing the luminosity and spectral evolution of flares caused by a supermassive black hole impacting the accretion disc of another supermassive black hole. Our model includes photon diffusion, emission from optically thin regions and relativistic corrections to the observed spectrum and time-scales. We test the observability of the impact scenario with a simulated population of quasars hosting supermassive black hole binaries. The results indicate that for a moderate binary mass ratio of 0.3, and impact distances of 100 primary Schwarzschild radii, the accretion disc impacts can be expected to equal or exceed the host quasar in brightness at observed wavelength λ = 510 nm up to z = 0.6. We conclude that accretion disc impacts may function as an independent probe for supermassive black hole binaries. We release the code used for computing the model light curves to the community.

preprint2014arXiv

Explicit methods in extended phase space for inseparable Hamiltonian problems

We present a method for explicit leapfrog integration of inseparable Hamiltonian systems by means of an extended phase space. A suitably defined new Hamiltonian on the extended phase space leads to equations of motion that can be numerically integrated by standard symplectic leapfrog (splitting) methods. When the leapfrog is combined with coordinate mixing transformations, the resulting algorithm shows good long term stability and error behaviour. We extend the method to non-Hamiltonian problems as well, and investigate optimal methods of projecting the extended phase space back to original dimension. Finally, we apply the methods to a Hamiltonian problem of geodesics in a curved space, and a non-Hamiltonian problem of a forced non-linear oscillator. We compare the performance of the methods to a general purpose differential equation solver LSODE, and the implicit midpoint method, a symplectic one-step method. We find the extended phase space methods to compare favorably to both for the Hamiltonian problem, and to the implicit midpoint method in the case of the non-linear oscillator.

preprint2013arXiv

A helical jet model for OJ287

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.

preprint2013arXiv

Short time-scale periodicity in OJ 287

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.