The supermassive binary black hole system OJ287

Abstract This doctoral thesis concerns the active galactic nucleus (AGN) most often referred to with the catalogue number OJ287. The publications in the thesis present new discoveries of the system in the context of a supermassive binary black hole model. In addition, the introduction discusses general characteristics of the OJ287 system and the physical fundamentals behind these characteristics. The place of OJ287 in the hierarchy of known types of AGN is also discussed. The introduction presents a large selection of fundamental physics required to have a basic understanding of active galactic nuclei, binary black holes, relativistic jets and accretion disks. Particularly the general relativistic nature of the orbits of close binaries of supermassive black holes is explored with some detail. Analytic estimates of some of the general relativistic effects in such a binary are presented, as well as numerical methods to calculate the effects more precisely. It is also shown how these results can be applied to the OJ287 system. The binary orbit model forms the basis for models of the recurring optical outbursts in the OJ287 system. In the introduction, two physical outburst models are presented in some detail and compared. The radiation hydrodynamics of the outbursts are discussed and optical light curve predictions are derived. The precursor outbursts studied in Paper III are also presented, and tied into the model of OJ287. To complete the discussion of the observable features of OJ287, the nature of the relativistic jets in the system, and in active galactic nuclei in general, is discussed. Basic physics of relativistic jets are presented, with additional detail added in the form of helical jet models. The results of Papers II, IV and V concerning the jet of OJ287 are presented, and their relation to other facets of the binary black hole model is discussed. As a whole, the introduction serves as a guide, though terse, for the physics and numerical methods required to successfully understand and simulate a close binary of supermassive black holes. For this purpose, the introduction necessarily combines a large number of both fundamental and specific results from broad disciplines like general relativity and radiation hydrodynamics. With the material included in the introduction, the publications of the thesis, which present new results with a much narrower focus, can be readily understood. Of the publications, Paper I presents newly discovered optical data points for OJ287, detected on archival astronomical plates from the Harvard College Observatory. These data points show the 1900 outburst of OJ287 for the first time. In addition, new data points covering the 1913 outburst allowed the determination of the start of the outburst with more precision than was possible before. These outbursts were then successfully numerically modelled with an N-body simulation of the OJ287 binary and accretion disc. In Paper II, mechanisms for the spin-up of the secondary black hole in OJ287 via interaction with the primary accretion disc and the magnetic fields in the system are discussed. Timescales for spin-up and alignment via both processes are estimated. It is found that the secondary black hole likely has a high spin. Paper III reports a new outburst of OJ287 in March 2013. The outburst was found to be rather similar to the ones reported in 1993 and 2004. All these outbursts happened just before the main outburst season, and are called precursor outbursts. In this paper, a mechanism was proposed for the precursor outbursts, where the secondary black hole collides with a gas cloud in the primary accretion disc corona. From this, estimates of brightness and timescales for the precursor were derived, as well as a prediction of the timing of the next precursor outburst. In Paper IV, observations from the 2004–2006 OJ287 observing program are used to investigate the existence of short periodicities in OJ287. The existence of a _50 day quasiperiodic component is confirmed. In addition, statistically significant 250 day and 3.5 day periods are found. Primary black hole accretion of a spiral density wave in the accretion disc is proposed as the source of the 50 day period, with numerical simulations supporting these results. Lorentz contracted jet re-emission is then proposed as the reason for the 3.5 day timescale. Paper V fits optical observations and mm and cm radio observations of OJ287 with a helical jet model. The jet is found to have a spine–sheath structure, with the sheath having a much lower Lorentz gamma factor than the spine. The sheath opening angle and Lorentz factor, as well as the helical wavelength of the jet are reported for the first time. Tiivistelmä Tässä väitöskirjatutkimuksessa on keskitytty tutkimaan aktiivista galaksiydintä OJ287. Väitöskirjan osana olevat tieteelliset julkaisut esittelevät OJ287-systeemistä saatuja uusia tuloksia kaksoismusta-aukkomallin kontekstissa. Väitöskirjan johdannossa käsitellään OJ287:n yleisiä ominaisuuksia ja niitä fysikaalisia perusilmiöitä, jotka näiden ominaisuuksien taustalla vaikuttavat. Johdanto selvittää myös OJ287-järjestelmän sijoittumisen aktiivisten galaksiytimien hierarkiassa. Johdannossa käydään läpi joitakin perusfysiikan tuloksia, jotka ovat tarpeen aktiivisten galaksiydinten, mustien aukkojen binäärien, relativististen suihkujen ja kertymäkiekkojen ymmärtämiseksi. Kahden toisiaan kiertävän mustan aukon keskinäisen radan suhteellisuusteoreettiset perusteet käydään läpi yksityiskohtaisemmin. Johdannossa esitetään joitakin analyyttisiä tuloksia tällaisessa binäärissä havaittavista suhteellisuusteoreettisista ilmiöistä. Myös numeerisia menetelmiä näiden ilmiöiden tarkempaan laskemiseen esitellään. Tuloksia sovelletaan OJ287-systeemiin, ja verrataan havaintoihin. OJ287:n mustien aukkojen ratamalli muodostaa pohjan systeemin toistuvien optisten purkausten malleille. Johdannossa esitellään yksityiskohtaisemmin kaksi fysikaalista purkausmallia, ja vertaillaan niitä. Purkausten säteilyhydrodynamiikka käydään läpi, ja myös ennusteet purkausten valokäyrille johdetaan. Johdannossa esitellään myös Julkaisussa III johdettu prekursoripurkausten malli, ja osoitetaan sen sopivan yhteen OJ287:n binäärimallin kanssa. Johdanto esittelee myös relativististen suihkujen fysiikkaa sekä OJ287- systeemiin liittyen että aktiivisten galaksiydinten kontekstissa yleisesti. Relativististen suihkujen perusfysiikka esitellään, kuten myös malleja kierteisistä suihkuista. Julkaisujen II, IV ja V OJ287-systeemin suihkuja koskevat tulokset esitellään binäärimallin kontekstissa. Kokonaisuutena johdanto palvelee suppeana oppaana, joka esittelee tarvittavan fysiikan ja tarpeelliset numeeriset menetelmät mustien aukkojen binäärijärjestelmän ymmärtämiseen ja simulointiin. Tätä tarkoitusta varten johdanto yhdistää sekä perustuloksia että joitakin syvällisempiä tuloksia laajoilta fysiikan osa-alueilta kuten suhteellisuusteoriasta ja säteilyhydrodynamiikasta. Johdannon sisältämän materiaalin avulla väitöskirjan julkaisut, ja niiden esittämät tulokset, ovat hyvin ymmärrettävissä. Väitöskirjan julkaisuista ensimmäinen esittelee uusia OJ287-systeemistä saatuja havaintopisteitä, jotka on paikallistettu Harvardin yliopiston observatorion arkiston valokuvauslevyiltä. OJ287:n vuonna 1900 tapahtunut purkaus nähdään ensimmäistä kertaa näissä havaintopisteissä. Uudet havaintopisteet mahdollistivat myös vuoden 1913 purkauksen alun ajoittamisen tarkemmin kuin aiemmin oli mahdollista. Havaitut purkaukset mallinnettiin onnistuneesti simuloimalla OJ287-järjestelmän mustien aukkojen paria ja kertymäkiekkoa. Julkaisussa II käsitellään mekanismeja OJ287:n sekundäärisen mustan aukon spinin kasvamiseen vuorovaikutuksessa primäärin kertymäkiekon ja systeemin magneettikenttien kanssa. Julkaisussa arvioidaan maksimispinin saavuttamisen ja spinin suunnan vakiintumisen aikaskaalat kummallakin mekanismilla. Tutkimuksessa havaitaan sekundäärin spinin olevan todennäköisesti suuri. Julkaisu III esittelee OJ287-systeemissä maaliskuussa 2013 tapahtuneen purkauksen. Purkauksen havaittiin muistuttavan vuosina 1993 ja 2004 tapahtuneita purkauksia, joita kutsutaan yhteisnimityksellä prekursoripurkaus (precursor outburst). Julkaisussa esitellään purkauksen synnylle mekanismi, jossa OJ287-systeemin sekundäärinen musta aukko osuu primäärisen mustan aukon kertymäkiekon koronassa olevaan kaasupilveen. Mekanismin avulla johdetaan arviot prekursoripurkausten kirkkaudelle ja aikaskaalalle. Julkaisussa johdetaan myös ennuste seuraavan prekursoripurkauksen ajankohdalle. Julkaisussa IV käytetään vuosina 2004–2006 kerättyjä havaintoja OJ287- systeemistä lyhyiden jaksollisuuksien etsintään. Julkaisussa varmennetaan systeemissä esiintyvä n. 50 päivän kvasiperiodisuus. Lisäksi tilastollisesti merkittävät 250 päivän ja 3,5 päivän jaksollisuudet havaitaan. Julkaisussa esitetään malli, jossa primäärisen mustan aukon kertymäkiekossa oleva spiraalitiheysaalto aiheuttaa 50 päivän jaksollisuuden. Mallista tehty numeerinen simulaatio tukee tulosta. Systeemin relativistisen suihkun emittoima aikadilatoitunut säteily esitetään aiheuttajaksi 3,5 päivän jaksollisuusaikaskaalalle. Julkaisussa V sovitetaan kierresuihkumalli OJ287-systeemistä tehtyihin optisiin havaintoihin ja millimetri- sekä senttimetriaallonpituuden radiohavaintoihin. Suihkun rakenteen havaitaan olevan kaksijakoinen ja koostuvan ytimestä ja kuoresta. Suihkun kuorella on merkittävästi pienempi Lorentzin gamma-tekijä kuin suihkun ytimellä. Kuoren avautumiskulma ja Lorentztekijä sekä suihkun kierteen aallonpituus raportoidaan julkaisussa ensimmäistä kertaa.

Behavior of black liquor nitrogen in combustion : formation of cyanate

The Kraft pulping process is the dominant chemical pulping process in the world. Roughly 195 million metric tons of black liquor are produced annually as a by-product from the Kraft pulping process. Black liquor consists of spent cooking chemicals and dissolved organics from the wood and can contain up to 0.15 wt% nitrogen on dry solids basis. The cooking chemicals from black liquor are recovered in a chemical recovery cycle. Water is evaporated in the first stage of the chemical recovery cycle, so the black liquor has a dry solids content of 65-85% prior to combustion. During combustion of black liquor, a portion of the black liquor nitrogen is volatilized, finally forming N2 or NO. The rest of the nitrogen remains in the char as char nitrogen. During char conversion, fixed carbon is burned off leaving the pulping chemicals as smelt, and the char nitrogen forms mostly smelt nitrogen (cyanate, OCN-). Smelt exits the recovery boiler and is dissolved in water. The cyanate from smelt decomposes in the presence of water, forming NH3, which causes nitrogen emissions from the rest of the chemical recovery cycle. This thesis had two focuses: firstly, to determine how the nitrogen chemistry in the recovery boiler is affected by modification of black liquor; and secondly, to find out what causes cyanate formation during thermal conversion, and which parameters affect cyanate formation and decomposition during thermal conversion of black liquor. The fate of added biosludge nitrogen in chemical recovery was determined in Paper I. The added biosludge increased the nitrogen content of black liquor. At the pulp mill, the added biosludge did not increase the NO formation in the recovery boiler, but instead increased the amount of cyanate in green liquor. The increased cyanate caused more NH3 formation, which increased the NCG boiler’s NO emissions. Laboratory-scale experiments showed an increase in both NO and cyanate formation after biosludge addition. Black liquor can be modified, for example by addition of a solid biomass to increase the energy density of black liquor, or by separation of lignin from black liquor by precipitation. The precipitated lignin can be utilized in the production of green chemicals or as a fuel. In Papers II and III, laboratory-scale experiments were conducted to determine the impact of black liquor modification on NO and cyanate formation. Removal of lignin from black liquor reduced the nitrogen content of the black liquor. In most cases NO and cyanate formation decreased with increasing lignin removal; the exception was NO formation from lignin lean soda liquors. The addition of biomass to black liquor resulted in a higher nitrogen content fuel mixture, due to the higher nitrogen content of biomass compared to black liquor. More NO and cyanate were formed from the fuel mixtures than from pure black liquor. The increased amount of formed cyanate led to the hypothesis that black liquor is catalytically active and converts a portion of the nitrogen in the mixed fuel to cyanate. The mechanism behind cyanate formation during thermal conversion of black liquor was not clear before this thesis. Paper IV studies the cyanate formation of alkali metal loaded fuels during gasification in a CO2 atmosphere. The salts K2CO3, Na2CO3, and K2SO4 all promoted char nitrogen to cyanate conversion during gasification, while KCl and CaCO3 did not. It is now assumed that cyanate is formed when alkali metal carbonate or an active intermediate of alkali metal carbonate (e.g. -CO2K) reacts with the char nitrogen forming cyanate. By testing different fuels (bark, peat, and coal), each of which had a different form of organic nitrogen, it was concluded that the form of organic nitrogen in char also has an impact on cyanate formation. Cyanate can be formed during pyrolysis of black liquor, but at temperatures 900°C or above, the formed cyanate will decompose. Cyanate formation in gasifying conditions with different levels of CO2 in the atmosphere was also studied. Most of the char nitrogen was converted to cyanate during gasification at 800-900°C in 13-50% CO2 in N2, and only 5% of the initial fuel nitrogen was converted to NO during char conversion. The formed smelt cyanate was stable at 800°C 13% CO2, while it decomposed at 900°C 13% CO2. The cyanate decomposition was faster at higher temperatures and in oxygen-containing atmospheres than in an inert atmosphere. The presence of CO2 in oxygencontaining atmospheres slowed down the decomposition of cyanate. This work will provide new information on how modification of black liquor affects the nitrogen chemistry during thermal conversion of black liquor and what causes cyanate formation during thermal conversion of black liquor. The formation and decomposition of cyanate was studied in order to provide new data, which would be useful in modeling of nitrogen chemistry in the recovery boiler.

Regime-Switching in the Impact of Oil Price Shocks on Stock Market Volatility: Evidence from Oil-Importing and Oil-Exporting Countries

Research has highlighted the adequacy of Markov regime-switching model to address dynamic behavior in long term stock market movements. Employing a purposed Extended regime-switching GARCH(1,1) model, this thesis further investigates the regime dependent nonlinear relationship between changes in oil price and stock market volatility in Saudi Arabia, Norway and Singapore for the period of 2001-2014. Market selection is prioritized to national dependency on oil export or import, which also rationalizes the fitness of implied bivariate volatility model. Among two regimes identified by the mean model, high stock market return-low volatility regime reflects the stable economic growth periods. The other regime characterized by low stock market return-high volatility coincides with episodes of recession and downturn. Moreover, results of volatility model provide the evidence that shocks in stock markets are less persistent during the high volatility regime. While accelerated oil price rises the stock market volatility during recessions, it reduces the stock market risk during normal growth periods in Singapore. In contrast, oil price showed no significant notable impact on stock market volatility of target oil-exporting countries in either of the volatility regime. In light to these results, international investors and policy makers could benefit the risk management in relation to oil price fluctuation.

Interdependence of Macroeconomics Variables, Stock Market Returns and Stock Market Volatility in Latin America

This thesis examines the interdependence of macroeconomic variables, stock market returns and stock market volatility in Latin America between 2000 and 2015. Argentina, Brazil, Chile, Colombia, Mexico and Peru were chosen as the sample markets, while inflation, interest rate, exchange rate, money supply, oil and gold were chosen as the sample macroeconomic variables. Bivariate VAR (1) model was applied to examine the mean return spillovers between the variables, whereas GARCH (1, 1) – BEKK model was applied to capture the volatility spillovers. The sample was divided into two smaller sub-periods, where the first sub-period covers from 2000 to 2007, and the second sub-period covers from 2007 to 2015. The empirical results report significant shock transmissions and volatility spillovers between inflation, interest rate, exchange rate, money supply, gold, oil and the selected markets, which suggests interdependence between the variables.

Return and Volatility spillovers among stock and Foreign Exchange Markets: Empirical Evidence from selected African Markets

Return and volatility dynamics in financial markets across the world have recently become important for the purpose of asset pricing, portfolio allocation and risk management. However, volatility, which come about as a result of the actions of market participants can help adapt to different situations and perform when it really matters. With recent development and liberalization among financial markets in emerging and frontier markets, the need for how the equity and foreign exchange markets interact and the extent to which return and volatility spillover are spread across countries is of importance to investors and policy makers at large. Financial markets in Africa have received attention leading to investors diversifying into them in times of crisis and contagion effects in developed countries. Regardless of the benefits these markets may offer, investors must be wary of issues such as thin trading, volatility that exists in the equity and currency markets and its related fluctuations. The study employs a VAR-GARCH BEKK model to study the return and volatility dynamics between the stock and foreign exchange sectors and among the equity markets of Egypt, Kenya, Nigeria, South Africa and Tunisia. The main findings suggest a higher dependence of own return in the stock markets and a one way return spillover from the currencies to the equity markets except for South Africa which has a weaker interrelation among the two markets. There is a relatively limited integration among the equity markets. Return and volatility spillover is mostly uni-directional except for a bi-directional relationship between the equity markets of Egypt and Tunisia. The study implication still proves a benefit for portfolio managers diversifying in these African equity markets, since they are independent of each other and may not be highly affected by the influx of negative news from elsewhere. However, there is the need to be wary of return and volatility spillover between the equity and currency markets, hence devising better hedging strategies to curb them.

Return and volatility linkages between the S&P 500 and emerging stock markets in Eastern Europe

The aim of this thesis is to research mean return spillovers as well as volatility spillovers from the S&P 500 stock index in the USA to selected stock markets in the emerging economies in Eastern Europe between 2002 and 2014. The sample period has been divided into smaller subsamples, which enables taking different market conditions as well as the unification of the World’s capital markets during the financial crisis into account. Bivariate VAR(1) models are used to analyze the mean return spillovers while the volatility linkages are analyzed through the use of bivariate BEKK-GARCH(1,1) models. The results show both constant volatility pooling within the S&P 500 as well as some statistically significant spillovers of both return and volatility from the S&P 500 to the Eastern European emerging stock markets. Moreover, some of the results indicate that the volatility spillovers have increased as time has passed, indicating unification of global stock markets.