Nuclear Outbursts in the Centers of Galaxies
A dissertation presented to the faculty of the College of Arts and Sciences of Ohio University
In partial fulfillment of the requirements for the degree Doctor of Philosophy
Reza Katebi December 2019
© 2019 Reza Katebi. All Rights Reserved. 2
This dissertation titled Nuclear Outbursts in the Centers of Galaxies
by REZA KATEBI
has been approved for the Department of Physics and Astronomy and the College of Arts and Sciences by
Ryan Chornock Assistant Professor of Physics and Astronomy
Florenz Plassmann Dean, College of Arts and Sciences 3 Abstract
KATEBI, REZA, Ph.D., December 2019, Physics Nuclear Outbursts in the Centers of Galaxies (182 pp.) Director of Dissertation: Ryan Chornock This dissertation consists of two parts. In the first part, we focus on studying the nuclear outbursts in the centers of galaxies and their nature in order to better understand the behavior of central Super Massive Black Holes (SMBHs) and their interaction with the surrounding environment, and to better understand the accretion disk structure. Nuclear outbursts can be better understood by studying the changes in the broad emission lines and the underlying continuum. We quantify the properties of these nuclear outbursts using multi-wavelength observations including optical, ultraviolet, and X-rays from MDM Observatory, the Sloan Digital Sky Survey, Swift, and Magellan. Some of these nuclear outbursts are linked to Tidal Disruption Events (TDEs) and nuclear supernovae (SNs), while a number of these events are proposed to be a rare phenomenon called “changing-look” Active Galactic Nuclei (AGN). These types of AGNs have been observed to optically transition from type 1 to type 2 and vice versa on timescales of months to years, where broad emission lines such as Hα and Hβ appeared or disappeared followed by an increase or decrease in the continuum light. We investigate two transient events called PS1-13cbe and PS1-10cdq that were observed during outburst by the PS1 survey in 2013 and 2010, respectively. We investigate TDE, SN, and AGN activity as the three possible scenarios for the nature of these events. In the case of PS1-13cbe, we conclude that “changing-look” behavior caused by thermal accretion disk instabilities is the most plausible explanation for the outburst. However, in the case of PS1-10cdq, we favor the tidal disruption scenario because of the structure of the lightcurve and spectral evolution. In the second part of this dissertation, we focus on galaxy morphology prediction using a newly designed neural network called “Capsule Networks”. We 4 automate the process of morphology prediction and eliminate the need for feature engineering and heavy data prepossessing prior to classification. We also reconstruct the galaxy images while preserving the brightness structure of the galaxies. This study provides one of the possible solutions for classifying objects for the era of large sky surveys where the amount of data for galaxies will dramatically increase and automated methods will play a very crucial role. 5 Dedication
To my father who taught to me to question everything To my mother who has always encouraged me to pursue my dreams To my sisters who have always supported me To my wife Volha who truly inspires me to be the best of myself and To everyone who is brave enough to wander in the realm of unknown and seek answers. 6 Acknowledgments
During my PhD, I had the privilege of meeting so many kind and smart people that I want to thank for their help and support. But at first I want to take a moment and express my deepest gratitude toward my advisor Professor Ryan Chornock who taught me the real meaning of science. Whenever I was struggling with a problem he was there to help me to move toward the right direction with his exceptionally vast knowledge. He taught me the correct path of scientific research and made me the physicist I am today. I would like to acknowledge him for his endless support throughout my PhD; no matter where I am or whoever I become, I will always be his lifelong student. I would like to thank physics and astronomy staff members for their effort to provide an efficient and proper atmosphere for research and learning. I would like to express my gratitude toward the members of my dissertation committee Dr. Razvan Bunescu, Dr. Madappa Prakash, and Dr. Joseph Shields for kindly accepting to be involved in my PhD. I would like to specially thank my wife Volha who has always supported me with love and inspired me to be better and stronger. At last but not least, I would like to thank my mother Sakineh, and my sisters Soheyla, Soraya, and Romina for always being there for me and their continuous support. 7 Table of Contents
Page
Abstract...... 3
Dedication...... 5
Acknowledgments...... 6
List of Tables...... 10
List of Figures...... 11
List of Acronyms...... 15
List of Sorted Acronyms...... 17
1 Introduction...... 19 1.1 Active Galactic Nuclei...... 19 1.1.1 AGN Classification...... 20 1.1.2 LINERs...... 24 1.1.3 AGN Unification...... 25 1.1.4 Accretion Process and Variability in AGNs...... 26 1.2 Changing Look AGNs...... 30 1.3 Tidal Disruption Events...... 33 1.4 Supernovae...... 38 1.4.1 Superluminous Supernovae...... 38 1.4.2 Sources that Power Superluminous Supernovae...... 39 1.4.2.1 Circumstellar Interaction...... 39 1.5 The Modular Open Source Fitter for Transients...... 41 1.6 Galaxy Morphology Prediction...... 42 1.6.1 Deep Learning...... 44 1.6.1.1 Feed-Forward Neural Networks...... 44 1.6.1.2 Convolutional Neural Networks...... 46
2 PS1-13cbe: The Rapid Transition of a Seyfert 2 to a Seyfert 1...... 49 2.1 Introduction...... 49 2.2 Observations of PS1-13cbe...... 51 2.2.1 Optical Photometry...... 52 2.2.2 Observations of the Host Galaxy...... 54 2.2.3 X-ray Photometry...... 55 2.2.4 Optical Spectroscopy...... 55 8
2.3 Observational Features of PS1-13cbe...... 57 2.3.1 Host galaxy of PS1-13cbe...... 57 2.3.2 Astrometry...... 60 2.3.3 Multi-band Light Curves of PS1-13cbe...... 61 2.3.4 Spectral Features of PS1-13cbe...... 65 2.4 Interpretation of the Features of PS1-13cbe...... 69 2.4.1 Type IIn Supernovae Interpretation of PS1-13cbe...... 70 2.4.2 PS1-13cbe as a TDE...... 72 2.4.3 PS1-13cbe as a “Changing Look” AGN...... 74 2.4.3.1 Obscuration of the AGN...... 75 2.4.3.2 Tidal Disruption Events...... 77 2.4.3.3 Accretion disk instabilities...... 77 2.4.4 Comparison to other Changing Look AGNs...... 79 2.5 Conclusions...... 80
3 PS1-10cdq...... 89 3.1 Observations of PS1-10cdq...... 89 3.1.1 Optical Photometry...... 89 3.1.2 Observations of the Host Galaxy...... 90 3.1.3 X-ray Photometry...... 90 3.1.4 Optical Spectroscopy...... 91 3.2 Observational Features of PS1-10cdq...... 92 3.2.1 Host Galaxy of PS1-10cdq...... 92 3.2.2 Multi-band Light Curves of PS1-10cdq...... 94 3.2.3 Spectral Features of PS1-10cdq...... 96 3.3 Interpretation of Features of PS1-10cdq...... 107 3.3.1 PS1-10cdq as a Variable AGN...... 108 3.3.2 PS1-10cdq as SLSN...... 111 3.3.3 PS1-10cdq as a TDE...... 113 3.4 Conclusions and Future Remarks...... 118
4 Galaxy morphology prediction using capsule networks...... 122 4.1 Introduction...... 122 4.2 Galaxy Zoo 2...... 125 4.3 Related Work...... 126 4.4 Approach...... 127 4.4.1 Experimental Setup...... 128 4.4.2 Data Preprocessing...... 128 4.4.3 Capsule Network...... 129 4.4.4 Network Architecture...... 131 4.4.4.1 Capsule Network...... 131 4.4.4.2 Baseline Network...... 132 4.4.5 Implementation and Resources...... 133 9
4.5 Results...... 134 4.5.1 Regression...... 134 4.5.2 Classification Based on Answers to Question 1 and Reconstruction of Galaxies...... 135 4.6 Conclusions...... 137
5 Conclusions and Outlook...... 142 5.1 Conclusions...... 142 5.2 Outlook...... 145
References...... 147
Appendix: Supplementary Material...... 178 10 List of Tables
Table Page
2.1 Photometry of PS1-13cbe...... 82 2.1 continued ...... 83 2.1 continued ...... 84 2.1 continued ...... 85 2.1 continued ...... 86 2.1 continued ...... 87 2.2 Luminosity of Broad Lines...... 88
3.1 Photometry of the host of PS1-10cdq...... 92 3.2 Late time spectra of the host of PS1-10cdq (Walter Baade (WB), and center (CN))...... 93 3.3 Luminosity of Broad Lines...... 105 3.4 Model parameters from fitting the CSM model from MOSFiT to the light-curve of PS10cdq...... 115 3.5 Model parameters from fitting the TDE model from MOSFiT to the light-curve of PS10cdq...... 120
4.1 Computed RMSE between the predictions and true crowd-sourced probabili- ties. A relative error reduction of 8.8% was achieved...... 134 4.2 Training and testing accuracy vs number of epochs for classification based on the answers to question one. A relative error reduction of 36.5% was achieved. 137 11 List of Figures
Figure Page
1.1 A representation of a feed forward neural network where HL stands for hidden layer. All of the layers are fully connected to each other...... 46 1.2 A schematic representation of a convolutional layer followed by a pooling layer. Each neuron in layer l + 1 is only connected to a local region of feature maps from layer l. The pooling layer applies a mean or max function on the local region of all feature maps from the convolutional layer l + 1...... 48
2.1 The observed transient luminosities of PS1-13cbe from the PS1 survey in grizyP1 filters after correction for Galactic extinction. S: marks the epoch of the LDSS spectrum (MJD 56570)...... 53 2.2 Light curves of PS1-13cbe in grizyP1 bands, corrected for Galactic extinction. Note the small-scale fluctuations in the light curve (such as near −20 d)..... 54 2.3 Optical spectra for PS1-13cbe, from bottom to top: host galaxy of PS1-13cbe from SDSS (grey), model of the host galaxy generated by FAST 1.0 (black), spectrum during outburst observed using LDSS3 (purple), spectra obtained with OSMOS in red and blue for the red and blue setups, respectively, and the most recent OSMOS spectrum (olive green)...... 57 2.4 The excitation diagrams using the [N II]λ6583/Hα, [O I]λ6300/Hα, [O III]λ5007/Hβ and [O III]λ5007/[O II]λ3726 line ratios (Baldwin et al., 1981; Kewley et al., 2006). The navy blue squares display the position of the host galaxy of PS1-13cbe. The shaded area represents the location of the SDSS galaxies calculated using MPA-JHU where darker regions represent higher number density of the galaxies. The extreme star formation line (solid red; Kewley et al., 2001, 2006), the revised star formation line (dotted black; Kauffmann et al., 2003) and the Seyfert-LINER classification line (dashed black; Kewley et al., 2006; Fernandes et al., 2010) are also plotted. Comp: AGN/star forming composites...... 59 2.5 SED of SDSS J2221+0030 in quiescence (purple circles), scaled SED templates of an Sb spiral galaxy (dotted orange line) and a Seyfert 2 galaxy (dashed blue lines) from the SWIRE template library (Polletta et al., 2007). The host exhibits a mid-infrared (MIR) excess relative to the star-forming template due to an AGN...... 61 2.6 The gP1 PS1 template image of SDSS J2221+0030 showing the position of the centroid of the galaxy (green “X”), the position of PS1-13cbe (red dot), and its position uncertainty is shown by a circle with 5σ radius (red circle). The position uncertainty is dominated by the outburst...... 62 12
2.7 Total optical luminosity light curve of PS1-13cbe integrated over the grizyP1 filters and relative to the baseline flux in the template images. The total optical luminosity was estimated using the spectral distribution at each epoch and the trapezoidal rule (blue triangles). The estimated luminosity at early and late times using yP1 (red circles) and zP1 (purple squares) bands assumed the same colour correction as measured from those epochs with all filters...... 63 2.8 Top: Rest-frame blackbody temperature from fitting the optical photometry. α Middle: Rest-frame spectral index from fitting a power law fν ∝ ν to the optical photometry, where α is the spectral index. Bottom: observed g-r colour diagram. All three panels show little or no evolution during the outburst..... 64 2.9 Continuum subtracted Hα line profiles. Top: Multiple-component Gaussian fit to the Hα+[NII] emission lines (blue), individual components (dashed green) and broad component of the Hα (dashed black). Bottom: The fit residuals. Top-left (on each panel): Observation date and numbers of days before/after the peak (navy blue)...... 67 2.10 Continuum subtracted Hβ line profiles. Top: Multiple-component Gaussian fit to the Hβ+[O III] emission lines (blue), individual components (dashed green) and broad component of the Hβ (dashed black). Bottom: The fit residuals. Top-left (on each panel): Observation date and numbers of days before/after the peak (navy blue). Strong residuals are present in the red wing of λ4959 due to poor subtraction of the 5577Å sky line...... 68 2.11 Optical spectrum of PS1-13cbe taken with LDSS3 during the outburst after subtracting the stellar continuum (blue). The spectra of the comparison objects include a QSO (SDSS QSO; Vanden Berk et al., 2001) template, a SN Type IIn (SN1994Y; Filippenko, 1997) and a TDE (ASASSN-14li; Holoien et al., 2015) from bottom to top (black), the emission lines are labeled (purple)...... 71
3.1 The observed Galactic extinction corrected luminosities of PS1-10cdq from the PS1 survey in grizyP1 filters. S: marks the epoch of the LDSS spectrum (MJD 55626)...... 90 3.2 Galactic extinction corrected light curves of PS1-10cdq in grizyP1 and Galex NUV bands...... 91 3.3 Optical spectra of PS1-10cdq. Spectrum ∼ 84 days after the peak of outburst observed using LDSS3 in dark red. Spectrum taken ∼ +352 days after the peak of outburst in purple. For detailed information about late time spectra shown in navy blue see Table 3.2. All the spectra are binned for the purpose of better visualization...... 94 3.4 Spectrum of the host galaxy generated by combining three late time spectra (spectra taken by IMACS on January, March, and May of 2017; see Table 3.2) (blue), the model of the host galaxy simulated with the FAST 1.0 code (Kriek et al., 2009) (orange)...... 95 13
3.5 SED of SDSS J100229+0130 in quiescence using observed and archival data (purple circles) (see Section 3.1.2), SED of PS1-13cbe at outburst from PS1 survey (green stars), scaled SED templates of an Sb spiral galaxy (dotted orange line) and a Seyfert 2 galaxy (dashed blue lines) from the SWIRE template library (Polletta et al., 2007). The observed MIR excess in the host relative to the star-forming template is due to an AGN...... 96 3.6 Total optical luminosity light curve of PS1-10cdq including the baseline that was calculated by integrating over grizp1 bands using trapezoidal method at each epoch...... 97 3.7 Top: Rest-frame blackbody temperature from fitting the optical photometry. α Middle: Rest-frame spectral index from fitting a power law fν ∝ ν to the optical photometry, where α is the spectral index. Red squares: represent fits where in addition to PS1 optical data, NUV data points from Galex are available. Bottom: observed g-r colour diagram. Strong evolution is evident in all panels...... 98 3.8 Continuum subtracted Hα+[N II] line profiles. Top: Multiple-component Gaussian fit to the Hα+[N II] emission lines (blue), individual components (dashed green) and broad component of the Hα (dashed black). Bottom: The fit residuals. Top-right (on each panel): Observation date and numbers of days after the peak (black)...... 101 3.9 Continuum subtracted Hβ line profile. Top: Multiple-component Gaussian fit to the Hβ emission line (blue), individual components (dashed green) and broad component of the Hβ (dashed black). Bottom: The fit residuals. Top-right (on each panel): Observation date and numbers of days after the peak (black).... 104 3.10 Continuum subtracted Mg II λλ2800. Top: Single Gaussian fit to the Mg II (blue). Bottom: The fit residuals. Top-right: Observation date and numbers of days after the peak (black)...... 106 3.11 Light-curve of PS1-10cdq (blue). Light curves of the comparison objects extremely luminous SLSN Type IIn/TDE candidate CSS100217 (Drake et al., 2011)(red), and TDE candidate PS16dtm (Blanchard et al., 2017)(green).... 109 3.12 Optical spectrum of PS1-10cdq taken with IMACS +84 days after the outburst after subtracting the stellar continuum (navy). Optical spectrum of TDE candidate PS16dtm +87 days after the outburst (Blanchard et al., 2017)(red). Optical spectrum of extremely luminous SLSN Type IIn/TDE candidate CSS100217 +84 days after the outburst (Drake et al., 2011)(purple). The emission lines are labeled (black)...... 110 3.13 light curves of PS1-10cdq in grizyPs and NUV GALAX bands (colored circles). Fits generated by MOSFiT’s CSM interaction model (colored lines)...... 113 3.14 Corner plot showing the posterior distributions of parameter realizations for MOSFiT’s CSM model. The 16th, 50th, and 84th quantiles of the posterior distribution of the model parameters (dash black lines). Median values with lower and upper error values of the model parameters are reported in the title of histograms (red)...... 114 14
3.15 Light curves of PS1-10cdq in grizyP1 and NUV Galex bands (colored circles). Fits generated by MOSFiT’s TDE model (colored lines)...... 118 3.16 Corner plot showing the posterior distributions of parameter realizations for MOSFiT’s CSM model. The 16th, 50th, and 84th quantiles the posterior distribution of the model parameters (dash black lines). Median values with lower and upper error values of the model parameters are reported in the title of histograms (red)...... 119
4.1 The architecture of the model used in this work. Top: represents architecture of the capsule layers where the GalaxCap layer has 2 or 37 capsules based on the different setups discussed in Section 4.4.4. Bottom: represents the structure of the decoder that acts as regularization during the training...... 132 4.2 Training and testing RMSE vs number of epochs for the regression scenario... 135 4.3 Training and testing accuracy vs number of epochs for classification based on the answers to question one...... 138 4.4 Original and reconstructed images of the galaxies...... 139 4.5 Estimated Sersic´ index of original images versus reconstructed images...... 140 4.6 Mean and 95% confidence interval of the difference between the estimated Sersic´ index for the original and reconstructed images (nOriginal − nReconstructed).. 141
A.1 A convolution on an 7×7 input image (green squares) with zero-padding P = 1 (white squares) with receptive field of size 3 × 3 (yellow squares) and a stride of S = 1 where the receptive field moved from left to right by one pixel..... 180 A.2 A max-pooling operation over a 3 × 3 convolutional matrix (blue squares) with receptive field of size 2 × 2 (red squares) and stride of 1...... 181 15 List of Acronyms In order of appearance.
AGN Active Galactic Nuclei
SMBH Super Massive Black Hole
BH Black Hole
UV ultraviolet
BLR Broad Line Region
NLR Narrow Line Region
GR General Theory of Relativity
Quasar Quasi Stellar Radio Source
QSO Quasar
RLQ Radio Loud Quasar
RQQ Radio Quiet Quasar
SSRQ Steep Spectrum Radio Quasars
FSRQ Flat Spectrum Radio Quasars
BLRG Broad-Line Radio Galaxies
NELG Narrow-Emission-Line X-ray Galaxy
NLRG Narrow-Line Radio Galaxies
FR Fanaroff-Riley
BL Lac BL Lacertae
LINER Low Ionization Nuclear Emission-Line Region
ISCO Innermost Stable Circular Orbit
CL Changing Look
TDE Tidal Disruption Event
SN Supernovae
CSI Circumstellar Interaction 16
CCSN Core Collapse Supernova
ZAMS zero-age Main Sequence
CSM Circumstellar Material
SLSN Super Luminous Supernova
FS Forward Shock
RS Backward Shock
CDS Cool Dense Shell
IC Inverse Compton
MCMC Markov Chain Monte Carlo
WAIC Watanabe-Akaike Information Criterion
PSRF Potential Scale Reduction Factor
FFNN Feed Forward Nerual Network
CNN Convolutional Neural Network
ReLU Rectified Linear Function
CapsNet Capsule Network 17 List of Sorted Acronyms Sorted in alphabetical order.
AGN Active Galactic Nuclei
BH Black Hole
BL Lac BL Lacertae
BLR Broad Line Region
BLRG Broad-Line Radio Galaxies
CapsNet Capsule Network
CCSN Core Collapse Supernova
CDS Cool Dense Shell
CL Changing Look
CNN Convolutional Neural Network
CSI Circumstellar Interaction
CSM Circumstellar Material
FFNN Feed Forward Nerual Network
FR Fanaroff-Riley
FS Forward Shock
FSRQ Flat Spectrum Radio Quasars
GR General Theory of Relativity
IC Inverse Compton
ISCO Innermost Stable Circular Orbit
LINER Low Ionization Nuclear Emission-Line Region
MCMC Markov Chain Monte Carlo
NELG Narrow-Emission-Line X-ray Galaxy
NLR Narrow Line Region
NLRG Narrow-Line Radio Galaxies 18
PSRF Potential Scale Reduction Factor
QSO Quasar
Quasar Quasi Stellar Radio Source
ReLU Rectified Linear Function
RLQ Radio Loud Quasar
RQQ Radio Quiet Quasar
RS Backward Shock
SLSN Super Luminous Supernova
SMBH Super Massive Black Hole
SN Supernovae
SSRQ Steep Spectrum Radio Quasars
TDE Tidal Disruption Event
UV ultraviolet
WAIC Watanabe-Akaike Information Criterion
ZAMS zero-age Main Sequence 19 1 Introduction
1.1 Active Galactic Nuclei
Active Galactic Nuclei are among the brightest objects in the Universe. Many galaxies have very bright central regions that can be more luminous than the light from the remainder of the galaxy. This bright central region leads to their naming as Active Galactic Nuclei (AGN). The origin of this emission is thought to be the accretion of the matter falling into a
6 9 Super Massive Black Hole (SMBH) with the mass of ≈ 10 − 10 M (Lynden-Bell, 1969). Historically, there were alternative suggestions for AGNs rather than SMBHs in the center of the galaxies. Some examples are multiple supernovae explosions caused by star formation bursts (Colgate and Cameron, 1963) or the presence of pulsars (Hoyle and Fowler, 1967). But shortly after that, it was shown that nuclear reactions would not be able to generate the observed energy and since all of the energy is coming from a very compact region (because of the observed rapid time variations), then accretion onto a SMBH at the center of the host galaxy would be a more suitable explanation (accretion is known to be one of the most efficient processes)(Lynden-Bell, 1969). Nowadays, it is an accepted idea that most of the galaxies contain a central Black Hole (BH) (Magorrian et al., 1998; Ho, 2004). One of the strongest pieces of evidence supporting this scenario is our own Milky Way which contains a BH at center; the presence of this BH was confirmed by tracing the orbital motions of the stars closest to Sgr. A* (a compact radio source which was a BH candidate) (e.g., Gillessen et al., 2009). As mentioned, accretion of the falling matter in SMBH is believed to be the origin of the emission from the center of the galaxies (Lynden-Bell, 1969). In the accretion process, matter emits ultraviolet (UV) and soft X-rays by losing angular momentum through viscous and turbulent processes. Hard X-rays are also produced very close to the BH because of the possible connection with a sea of hot electrons above the disk (Urry and Padovani, 20
1995). By realizing that the BH mass and the properties of the host galaxy have a strong correlation, researchers started linking BH mass growth with the formation of the host galaxy. Black hole growth that results from accretion and the energy feedback has affected most of the recent structure formation studies (Granato et al., 2004; Springel et al., 2005; Hopkins et al., 2006; Ho, 2008). All of the AGNs emit in a wide spectrum of wavelengths from γ-rays and X-rays to radio that are known to originate from different regions. Strong UV and optical emission lines are believed to be produced in the rapidly moving gas clouds in the gravitational potential of the BH known as Broad Line Region (BLR) clouds. UV and optical radiation experience obscuration by a torus, which is a warped disk of gas and dust particles outside of the accretion disk. Narrow emission lines, however, are produced in slowly moving gas clouds beyond the torus that are called the Narrow Line Region (NLR). Along the poles of the disk, the outflow of energetic articles generates jets that are the main source of radio emission. Outward high velocity streams of plasma in jets beams radiate relativistically in the forward direction (Urry and Padovani, 1995). Because of their high luminosity, AGNs are observable up to very high redshifts. For that reason they can provide very useful information about the structure of the Universe at early times and its evolution. Also, understanding more about characteristics of the AGNs such as BH mass, BH spin, the type of host galaxy and accretion rate will help us to find out more about the accretion process and interaction of the BH with its environment (Urry and Padovani, 1995). On the other hand, the SMBHs in their center can be a good probe to investigate the General Theory of Relativity (GR).
1.1.1 AGN Classification
The axisymmetric model of AGNs suggests that there can be different types of AGNs based on the angle toward the line of sight. Based on their optical spectra, luminosity, and 21 radio loudness, a brief summary of the classification of AGNs is provided in the following paragraphs (Urry and Padovani, 1995).
• Type 1: This type of AGN is known to have bright continua and broad emission lines that originate from fast moving clouds close to the BH in the BLR region. They are subdivisions of Type 1 galaxies which are significantly different only in luminosity and radio loudness.
– Seyfert 1 galaxies: The main classification of the Seyfert galaxies was done by Carl Seyfert (Seyfert, 1943) and they are divided to Seyfert 1 and Seyfert 2 sub- categories (Khachikian and Weedman, 1974). Seyfert 1 are radio-quiet Type 1 AGNs and they are seen nearby at low redshifts because of their low luminosity. Their spectra show both broad and narrow lines. As mentioned before, broad lines originate from high velocity, dense, and highly ionized gas clouds in the BLR which is close to the BH. In contrast, narrow lines are emitted from low- ionization, low velocity and less dense gas clouds in the NLR, which is beyond the torus and the accretion disk (Ho, 2008).
– Quasi Stellar Radio Source (Quasar): Quasars (QSOs) (Schmidt, 1963) are highly luminous and for that reason they can be observed up to high redshifts. Because of their high luminosity, only the star-like nucleus is observed, whereas a surrounding galaxy can be rarely observed. That is why they are called “Quasi Stellar Radio Sources”. QSOs can be both radio-loud (RLQ) and radio-quiet (RQQ). Based on the shape of their continuum, RLQs can be divided into a subgroup of Steep Spectrum Radio Quasars (SSRQ) or Flat Spectrum Radio Quasars (FSRQ) (Urry and Padovani, 1995) 22
– Broad-Line Radio Galaxies (BLRG): They are low luminosity Type 1 AGNs that have strong radio emission that probably comes from radio lobes which are powered by jets (Blandford and Rees, 1974; Urry and Padovani, 1995)
– Narrow Line Seyfert 1 (NLS1) galaxies: The most important characteristics
of NLS1s are the widths of the broad Balmer lines (specially, FWHMHβ < 2000
−1 km s ) and weakness of [O III]λ5007 emission relative to Hβ of [O III]/Hβtotl < 3 (Osterbrock and Pogge, 1985; Goodrich, 1989). Furthermore, they typically show strong Fe II features which anti-correlate with [O III] emission and broad Balmer emission lines. Also, the presence of Fe II features is suggested as an
extra measure for classification of NLS1s (Fe II/Hβtotl < 0.5;V eron-Cetty´ et al., 2001).
• Type 2: This type of the AGNs has weak continua and narrow emission lines only. It is believed that the broad line region is obscured by a torus of absorbing gas and dust around it. Like the Type 1 AGNs, the Type 2s can be also divided into subdivisions based on their luminosity and radio loudness (see Figure 1 in (Urry and Padovani, 1995)). The high luminosity Type 2 AGNs have been studied using samples selected from the SDSS survey in optical, radio, and X-rays (Zakamska et al., 2003, 2004). Using spectropolarimetry, Zakamska et al.(2005) revealed the presence of a hidden Type 1 nuclei in these high luminosity Type 2 AGNs.
– Seyfert 2 galaxies: Seyfert 2 galaxies are low luminosity and radio-quiet Type 2 AGNs. Only narrow emission lines can be seen in their spectra. Seyfert 2 galaxies have [OIII]λ5007/Hβ ratio of > 3, which can be useful for their identification (Shuder and Osterbrock, 1981). There are also other known low luminosity counterpart Type 2 AGNs that are called narrow-emission-line X- ray galaxies (NELGs) (Mushotzky, 1982). 23
– Narrow-Line Radio Galaxies (NLRGs): NLRGs are radio loud Type 2 AGNs and like BLRGs their radio emission comes from radio lobes which are powered by a jet of particles (Blandford and Rees, 1974; Urry and Padovani, 1995). The NLRGs are divided into different morphological categories based on the appearance of the radio jets (Fanaroff and Riley, 1974): The low luminosity Fanaroff-Riley Type 1 radio galaxies (FRI) are the first category, where the center is the brightest, the radio-jets are symmetric, and the intensity decreases by distance from the nucleus. The second category objects are high luminosity FRII radio galaxies, which have very well defined radio-lobes with hot spots and the brightness increasing by distance from the core. FRIIs are known to contain both NLRGs and BLRGs (Blandford and Rees, 1974; Urry and Padovani, 1995).
• Type 0: A small number of AGNs are called Type 0 AGNs because of their unusual characteristics (Urry and Padovani, 1995). These objects have close to zero degrees angles toward the line of sight.
– BL Lacertae (BL Lac): These objects are radio-loud AGNs with an unusual characteristic of not having strong emission or absorption features. Some other characteristics of them include surprisingly variable and high polarization, fast variability and high brightness temperatures (Urry and Padovani, 1995).
In summary, the AGNs are divided into the three main sub-classes of Type 1, Type 2 and Type 0, based on their optical spectra properties. Each of these types can have sub- classes based on radio-loudness and luminosity. As a last note, Low Ionization Nuclear Emission-Line Region galaxies (LINERs) will be discussed as a source of weak AGNs in the next section. 24
1.1.2 LINERs
Heckman(1980) studied nuclear activity of a class of objects that were di fferent from known AGNs and HII regions. These so called Low Ionization Nuclear Emission-Line Region galaxies, are known to have low-ionized narrow emission lines. The objects which had [OII]λ3727/[OIII]λ5007 ≥ 1 and [OI]λ6300/[OIII]λ5007 ≥ 1/3 were classified as LINERs (Filippenko, 1996). Numerous spectroscopic surveys concluded that LINERs exist in very large numbers and are more common than it was thought before (Heckman, 1980; Heckman et al., 1980; Filippenko, 1996). Moreover, there is a possibility that more than half of spiral galaxies contain LINERs (Ho et al., 1995). There is not too much activity observed in their nucleus and most of the emission comes from the stars in the galaxy. A number of other studies have shown that some of the LINERs contain broad Hα line similar to Seyfert 1 galaxies (Keel, 1983a,b; Filippenko and Sargent, 1985; Stauffer, 1982; Veron-Cetty´ and Veron´ , 1986; Filippenko, 1996). The LINERs are believed to have a continuum with non-stellar origin (Netzer, 1990). The main candidate for the excitation of the LINERs is photoionization by a central Low Luminosity Active Galactic Nuclei (LLAGN). Additionally, LINERs can be a radio (Wrobel and Heeschen, 1991) or an X-ray (Halpern and Steiner, 1983) emission source. ∼ 25% of the LINERs have UV emission from an unresolved core (Maoz et al., 1996) that can be the sign of accretion activity in their center. However, dust obscuration prevents UV detection in the other LINERs (Ho, 2008; Pogge et al., 2000). The LINERs can also fall into Type 1 and Type 2 categories, based on the observed obscuration. Understanding the nature of the LINERs is very important, since they make a large proportion of the nearby galaxies. Possible power sources for ionization of the LINERs are: shock heating created by the accretion process; winds generated by star-bursts or galaxy mergers and interactions; cooling gas flow and photoionization by evolved or main sequence very hot stars (Filippenko, 1996). 25
1.1.3 AGN Unification
The simple model of AGN suggests that there is a SMBH with an accretion disk around it. The BLR is located in a greater distance with the fast moving gas clumps. There is a warped torus of gas which surrounds the BLR and the NLR is located at larger distances. As mentioned in previous section, AGNs come with a great number of different sub-classes, but the classification of them between Type 1 and Type 2 is based on obscuration of the nucleus and the classification of them based on radio-loudness depends on the angle of the jets with the line of sight (Urry and Padovani, 1995). The main scenario for the classification of the AGNs is the orientation to the line of sight. By observing NGC 1068, Antonucci and Miller(1985) realized that this classic Seyfert 2 galaxy has very broad lines similar to the Seyfert 1 broad Balmer lines that are observed in polarized light. Later discovery of the other similar Seyfert 2s with Seyfert 1 like broad Balmer lines (Miller and Goodrich, 1990) inspired the model that states that Seyfert 2 and Seyfert 1 galaxies are the same objects which are observed from different angles. This model proposes that there is a warped torus of the absorbing gas and dust particles surrounding the BLR in a way that the BLR is hidden from some orientations which is classified as a Type 2 AGN. Based on this model, emission from the BLR that is scattered by high energy electrons outside the nucleus that permits the observation of the broad lines in the polarized light (Antonucci and Miller, 1985; Antonucci, 1993). As a part of this unification scheme, it is important to discuss broad and narrow lines and their origins in more details. Broad lines are permitted lines with Doppler-broadened
3 4 −1 10 widths of 10 −10 km s that originate from the BLR. The BLR contains dense (nH ≥ 10 atoms cm−3), fast moving and highly ionized gas clouds. The observed continuum radiation with wavelength λ < 912Å in most Seyfert galaxies (these photons did not pass through the BLR because if they did, they would be absorbed by hydrogen gas in the galaxy and along the line of sight to us) and Doppler-broadened emission lines and highly ionized 26 gas, suggest that the BLR is a compact region located deep in the gravitational well of the central BH. In contrast, narrow emission lines, such as [OIII]λ5007, are coming from the transitions which are mainly forbidden with widths of 102−103 km s−1 that are emitted from
8 the NLR. The forbidden lines are observed only where the density of the gas is nH ≤ 10 , which suggests that the NLR should contain low density gas clouds. These forbidden lines have not shown any variations in short timescales, demonstrating that the NLR is located far from the central BH (Sparke and Gallagher III, 2007; Bennert et al., 2002, chapter 9). However, in one of the recent studies Peterson et al. (2013) show that the fluxes of the narrow [OIII]λλ4959, 5007 emission lines in the spectra of NGC 5548 (a well-known Seyfert 1 galaxy) vary in the timescale of decades. Therefore, from this variability they show that the NLR is more compact and denser than previously known (Peterson et al., 2013). Although it is not complete, the simple unification model mentioned above that is only based on orientation effects has been the base model for the unification of Type 1 and Type 2 AGNs. However, multiple AGN types in some objects have been recently reported in different epochs, which challenges this model. These interesting objects and the possible explanations for this type of changes will be discussed in section 1.2.
1.1.4 Accretion Process and Variability in AGNs
The material that is far away from the SMBH in the central regions of the galaxies loses angular momentum while flowing inward and toward the central engine. The angular momentum is lost through the viscous processes and is converted to heat and the magnetic field. This rotationally supported flow forms what is known as an accretion disk (Koratkar and Blaes, 1999). Accretion disks are among the most efficient energy providers from a compact volume, which makes them one of the main components of the AGN models and a possible source of the power in them (Salpeter, 1964; Zeldovich and Novikov, 1964). 27
However, there is a limit to their luminosity which is known as Eddington luminosity, which happens when the radiative force is equal to gravitational force (Eddington, 1925)
4πGMmpc 38 −1 M Ledd = ≈ 1.26 × 10 ergs (1.1) σT M
where mp is the mass of the electron and σT is the Thompson scattering cross section. The accretion flow of the matter in the accretion disk is believed to be the source of X-ray, UV and optical emission which powers AGN by ionizing the surrounding gas in both the BLR and NLR and may also be the birth place of the jets and winds (Koratkar and Blaes, 1999). Observation of the central engine of AGN, which contains the accretion flow, is hard and it is difficult to resolve the small volume at the galaxy centers where it is located. However, one strong piece of observational evidence that shows that there is gas located deep in the potential well of the SMBH in AGNs is the detection of the very broad Fe Kα line (Tanaka et al., 1995; Mushotzky et al., 1995; Nandra et al., 1997). The simplest model of the accretion flow is that the plasma orbits the central BH in “Keplerian” orbits and the flow of gas is geometrically thin and optically thick. The vertical structure of the disk; however, results from balancing the gas pressure and the gravitational force provided by the BH (Koratkar and Blaes, 1999). The standard theoretical models proposed for this geometrically thin disk structure, (Shakura and Sunyaev, 1973; Novikov and Thorne, 1973; Koratkar and Blaes, 1999) make a couple of important assumptions: (1) the disk is stationary and axisymmetric and flows all the way to the Innermost Stable Circular Orbit (ISCO) where the shear stress is assumed to be zero; (2) the gravitational field of the hole is assumed to be Newtonian and matter orbits around the hole with Keplerian orbits; (3) the mass of the accretion disk is negligibly small compared to the mass of the BH; (4) the viscosity makes some material fall inward and some move outward as it orbits because of conservation of angular momentum. The lost gravitational energy during this process is converted to radiation due to conservation 28 of energy; (5) the disk is thin, so heat transfers vertically and at the same radius; therefore, the emission is a local blackbody (Pereyra et al., 2006; Koratkar and Blaes, 1999). Despite the fact that the above model is the base model used for accretion flow, there exist a number of limitations for it. The first one is that the model does not work for high accretion luminosities where the radiation pressure plays an important role. Also, the SED
1/3 of Fν ∝ ν resulting from this model is bluer than the observed UV and optical spectra of the AGNs. Furthermore, a local blackbody is not a good approximation even with an optically thick disk, which may not be always the case (Koratkar and Blaes, 1999). AGN luminosities have been observed to vary over time. Variability has been detected in a wide range of wavelengths with timescales from hours to years with the shortest in X-rays (Koratkar and Blaes, 1999). X-ray observations with variation time-scales of a few hours are observed in Seyfert galaxies (Green et al., 1993; Nandra et al., 1997) and the observed average variability time-scale in Seyfert galaxies is on the order of weeks to months (Kaspi et al., 1996; Giveon et al., 1999). A number of models have been proposed for the source of this variability, such as supernovae (Aretxaga et al., 1997) and microlensing (Hawkins, 2000), but accretion disk instability is the most promising one (Kawaguchi et al., 1998; Kelly et al., 2009). Reverberation mapping studies (Blandford and McKee, 1982) have shown that the broad emission lines respond to the variations with time lag which is strong evidence that the continuum variations are derived from accretion instabilities (Peterson et al., 2004; Kelly et al., 2009). Optical and UV emission are believed to be thermal emission originating from the accretion disk (Shields, 1978) and if this scenario is valid then thermal fluctuations would be the source of the variations (Kelly et al., 2009). One observational piece of evidence in favor of this scenario is that quasars have been observed to be bluer when brighter (Giveon et al., 1999; Trevese et al., 2001; Geha et al., 2003). 29
The model developed by Shakura and Sunyaev (1973) is also known as a standard α model. The key assumption in their model is that the viscosity is the source of thermal emission and angular momentum lost in accretion disks, which is proportional to the total pressure of the disk. The viscosity is
ν = αcsh (1.2)
where h is gas scale height and cs is speed of sound. A number of studies on quasar variability estimated the value of α to be α ≈ 0.001 (Siemiginowska and Czerny, 1989; Collier and Peterson, 2001; Starling et al., 2004; Kelly et al., 2009). Three important timescales in accretion disks are :
tdyn = 1/Ω (1.3)
−2 −2 −1 H H −1H tvisc ≈ α ≈ α tdyn (1.4) cs R R −1 tth = α tdyn (1.5) q GM where Ω = R3 is angular velocity, tdyn is dynamical timescale, tvisc is the viscous H timescale and tth is the thermal timescale and since in geometrically thin disks ( R ) 1 then the ordering of the timescales is tvisc tth > tdyn (Lasota, 2016). In radiation-dominated disks, which are also believed to be the characteristic of the inner disk environment, the α-disk is unstable both thermally and viscously (Shakura and Sunyaev, 1976; Lightman and Eardley, 1989). This thermal instability is believed to grow exponentially on thermal timescales which are on the order of months to years for AGNs and instability timescales in optical light curves of AGNs are not usually detected to be ≥ tth (Kelly et al., 2009), except some specific sources (Czerny et al., 2003; Lub and de Ruiter, 1992). In an optically-thick accretion disk, different wavelength ranges of the spectrum come from different regions of the flow that are spatially separated. Thus, time lags in variability can happen because of the finite speed of the propagating signal between different regions (Koratkar and Blaes, 1999). A possible explanation would be that optical 30 and UV variations are caused by reprocessing of X-ray or UV radiation from the inner regions of the disk (Krolik et al., 1991; Collin-Souffrin, 1991). This scenario has been proven to be correct by multiwavelength observations of the Seyfert galaxies NGC 4151 and NGC 5548 that showed time lags between X-ray and UV/optical which suggest that the UV/optical variability is driven by reprocessing of X-rays (Edelson et al., 1996; Clavel et al., 1992; McHardy et al., 2014).
1.2 Changing Look AGNs
As mentioned in section 1.1, there is a SMBH at the center of most galaxies that grows by accreting matter. AGNs are identified by the observed properties of their optical spectra, broad and narrow lines in Type 1s and just narrow lines in Type 2s. Also, intermediate AGNs such as Type 1.8 and 1.9, are classified by weak or absent broad Hβ while showing broad Hα in their spectra (Osterbrock, 1981). As discussed in section 1.1.3, the standard unification model suggests that the Type 1 and the Type 2 AGNs are the result of different viewing angles to the line of sight of an obscuring torus of gas and dust particles surrounding the BLR and the intermediate type AGNs are suggested to be the result of partial obscuration or reddening by optically thin dust (Stern and Laor, 2012). It is believed that the simple orientation-based model cannot explain the Type 1 and Type 2 classification scheme to the full extent. This model is challenged when different AGN Types are observed in the same object at different epochs of time. One alternative suggestion is that Type 2s were Type 1s and their engine is now not active anymore or basically turned off; therefore, the Type 2s are evolved version of the Type 1s (Penston and Perez´ , 1984; Runnoe et al., 2016). This alternative view was first suggested based on the observations of broad emission lines disappearance in NGC 4151 (Lyutyj et al., 1984; Penston and Perez´ , 1984) and 3C 390.3 (Penston and Perez´ , 1984). 31
The objects that change type at different epochs of time can help us to advance our understanding of the evolutionary side of unification model. Although rare, these objects have been observed in the past. The term “changing-look” was first used for the AGNs that showed X-ray absorption variations, which in simple terms means that they were observed to appear both “Compton-thin” or “reflection-dominated” at different epochs of time in timescale of years (Matt et al., 2003; Bianchi et al., 2009; Puccetti et al., 2007; Risaliti et al., 2009; Marchese et al., 2012). Additionally, the term “changing-look” (CL) AGNs has been also used to describe the objects that have been observed to optically transition from Type 1 to Type 1.8, 1.9 and 2 or vice versa, by sudden appearance or disappearance of broad emission lines such as the broad Hβ emission line (Runnoe et al., 2016). Such objects have been detected in the past four decades in a number of the AGNs. In some objects the broad emission lines disappeared and the continuum faded completely (Collin-Souffrin et al., 1973; Tohline and Osterbrock, 1976; Sanmartim et al., 2014; Denney et al., 2014; Barth et al., 2015; Runnoe et al., 2016), while in others broad emission lines appeared or, in other words, the AGN “turned-on” (Cohen et al., 1986; Storchi-Bergmann et al., 1993; Aretxaga et al., 1999; Eracleous and Halpern, 2001; Shappee et al., 2014). These objects are at low redshift with low luminosities. The “changing-look” behavior that is observed in some of the AGNs can possibly be explained by three suggested mechanisms. The first mechanism is by variation of obscuration where the obscuring material such as dust clouds moves in or out of the line of sight and obscures or gives a clear view to the BLR (Elitzur, 2012). Another possible explanation can be a change in the accretion rate. The change in the accretion rate transforms the structure of the BLR and the objects will transition from the high accretion rate and luminous Type 1 to a low accretion state Type 2 (Elitzur et al., 2014), because multiple models have been suggested that a radiatively efficient BLR cannot exist in low accretion rates. Some of the possible explanations are: the lack of sufficient ionizing 32 photons to power the BLR (Korista and Goad, 2004), the size of critical radius (at this radius disk changes from gas pressure dominated to radiation pressure dominated) becoming smaller than ISCO (Nicastro, 2000; Nicastro et al., 2003), a minimum required bolometric luminosity where the BLR cannot exist if it goes lower (Laor, 2003), and other reasons (Denney et al., 2014; Elitzur et al., 2014; Trump et al., 2011). It has also been argued that transient events, such as Tidal Disruption Events (TDEs), can cause type change (Eracleous et al., 1995). There are observational results that suggest that each of the three scenarios mentioned above is a possible physical process behind the CL behavior in different objects. For example, Alexander et al.(2013) observed a source at redshift z = 0.510 and argued that obscuring material moved into the line of the sight and obscured the nucleus of the object that changed the class of the source from a Type 1 to a Type 2. In another observation, Shappee et al.(2013) reported an interesting transition from Sy 1.8 to Sy 1 in NGC 2617 and they argued that X-ray irradiation of the disk, likely resulting from coronal activity, is causing this variability because the variability was first detected in X-rays and then with time lags in UV and NIR (Shappee et al., 2014). LaMassa et al.(2015) also favor changes in the accretion power as an explanation for the observed transition from a Type 1 quasar to a Type 1.9 AGN in J0159+0033 because they estimated the crossing time for an possible obscuring object outside the BLR with a Keplerian orbit to be tcross = 20 years that is too long to explain this transition which happened in shorter timescales. They are also skeptical about the existence of obscuring material with the required physical properties to obscure the BLR and whole continuum at such radii. However, Merloni et al.(2015) proposed an alternative scenario that the nature of flare in J0159+0033 is a tidal disruption event (TDE; will be discussed in more detail in section 1.3) by interpreting that this object was first in a low accretion state and suddenly showed an outburst and brightened dramatically. 33
Additionally, they showed that the very rapid rise and decay time of the flare is consistent with the t−5/3 behavior expected in TDEs. The number of objects that show CL behavior is increasing and there already have been a number of bulk searches in the SDSS and Pan-STARRS surveys that found new CL candidates (Ruan et al., 2016; MacLeod et al., 2016). As we can see, the reason behind CL behavior is not very well known yet and is still under debate. Thus, understanding the nature of these objects will help us to further study the structure of the accretion disk and the accretion process of BHs and consequently shed more light on the structure of the BLR. Furthermore, it will give us an opportunity to investigate the unification models of AGNs in more detail.
1.3 Tidal Disruption Events
When a star orbits a BH with a sufficiently close pericenter and the BH’s gravitational force exceeds the self-gravity of the star, the star is disrupted (Hills, 1975). This happens at a distance which is known as tidal radius:
M 1/3 M −1/3 r r ≈ × 12 BH ∗ ∗ cm t 7 10 6 (1.6) 10 M M r
In a Tidal Disruption Event (TDE), debris from the disrupted star is divided into almost two equal parts. One of the parts becomes unbound on hyperbolic orbits while the other part is bound to the hole with elliptical trajectories and a pericenter the same as original star’s pericenter. Eventually, bound debris will be accreted to the BH and produce a luminous transient event (Rees, 1988; Komossa, 2015). However, the tidal radius scales differently
2GMBH with the mass of the BH (as seen in eq.1.6) than the Schwarzschild radius (rs = c2 ), so solar type stars will pass the horizon without getting disrupted for SMBHs with masses
8 higher than ∼ 10 M . This limit is higher when spin is included (Beloborodov et al., 1992). 34
7 In TDEs involving a solar-type star and SMBHs with the mass of ≤ 10 M which accretes with Eddington limit, the spectrum is a blackbody with the characteristic blackbody temperatures of ∼ 2.5 × 105 K at the tidal radius (Ulmer, 1999; Strubbe and Quataert, 2009). The radiation from these events peaks in the UV and soft X-rays, declining on the timescale of months to years. One fundamental puzzle is that the UV/optical flux predicted by a standard BH accretion theory is orders of magnitude less than the observed values and has a blue constant color (Roth et al., 2016). A possible solution to this problem is the presence of gas which absorbs and re-emits the radiation with temperatures around ∼ 20000 − 30000 K. The source of this gas is thought to be the formation of an envelope around the BH (Loeb and Ulmer, 1997; Guillochon et al., 2014; Coughlin and Begelman, 2014); or supper-Eddington mass outflow (Strubbe and Quataert, 2009; Lodato and Rossi, 2010; Stone and Metzger, 2015; Vinko´ et al., 2014; Miller, 2015) or the circularization of the material at much larger distances than the tidal radius of the BH (Shiokawa et al., 2015; Piran et al., 2015; Hayasaki et al., 2016; Bonnerot et al., 2015; Guillochon and Ramirez- Ruiz, 2013, 2015; Dai et al., 2015). The light curves of TDEs are predicted to follow a rapid rise to a maximum peak and then decline by the fall back rate of t−5/3 (Rees, 1988; Evans and Kochanek, 1989). This fall-back rate was derived by assuming that the the spread of the specific energy for debris is constant with mass. However, recent detailed modeling has shown that the internal structure of the disrupted star can affect this assumption and predicts a faster decline rate than of t−5/3 for the light-curve of the TDEs after the peak (Ramirez-Ruiz and Rosswog, 2009; Lodato et al., 2008; Stone et al., 2013; Guillochon and Ramirez-Ruiz, 2013). Furthermore, the process which turns accreted mass into observable radiation is not simple to simulate. The simplest model to simulate the Spectral Energy Distribution (SED) of TDEs and their evolution is the thin disk model but the new and more complex models 35 have added a thick disk and a wind or outflow to represent super-Eddington accretion (Guillochon and Ramirez-Ruiz, 2013). The number of observed TDEs is rapidly increasing and the rate of these events are estimated to be ∼ 10−4 to 10−5 events per year per galaxy (e.g., Magorrian and Tremaine, 1999; Donley et al., 2002; Wang and Merritt, 2004; Kesden, 2012; Stone and Metzger, 2015; Van Velzen and Farrar, 2014). A large portion of TDE candidates are soft X-ray flares with large amplitudes from the galactic centers (Komossa and Bade, 1999; Donley et al., 2002; Komossa et al., 2004; Halpern et al., 2004; Esquej et al., 2007; Cappelluti et al., 2009a; Maksym et al., 2010; Saxton et al., 2012; Hryniewicz and Walter, 2016; Lin et al., 2015; Komossa, 2015). In these events, the luminosity peaks at 1044 erg s−1 and the SED peaks at soft X-rays with energies ≤ 0.1 kev. Additionally, the decay rate of the light curve is consistent with the predicted fall-back rate of the debris t−5/3 (Rees, 1988; Evans and Kochanek, 1989; Lodato et al., 2008; Guillochon and Ramirez-Ruiz, 2013; Roth et al., 2016). Van Velzen et al.(2011) reported the first optical TDEs (SDSS TDE1, SDSS TDE2) from multi-epoch imaging of Stripe 82 in SDSS. They selected these candidates because of their unusually blue colors with slow evolution which was very different from supernovae (SNe) and AGN variability. Gezari et al.(2012) reported the TDE PS1-10jh from PAN- STARRS1 (PS1) with optical light curves from the time of the rise and decline from the peak which became the subject of studies for the process that derives TDEs (Guillochon and Ramirez-Ruiz, 2013; Bogdanovic´ et al., 2014). There have been TDE candidates discovered in ultraviolet (UV) (Gezari et al., 2006, 2009) and other optical TDE candidates reported from PAN-STARRS1 (Chornock et al., 2013), ASASSN (Holoien et al., 2014, 2015), PTF (Cenko et al., 2012a; Arcavi et al., 2014) and ROTSE (Vinko´ et al., 2014). The light curves of these TDEs rise and peak at R-band luminosity 2 × 1043 erg s−1 and decline consistent with t−5/3(Arcavi et al., 2014). Additionally, TDEs like PTF10iya and 36
ASASSN-14li are observed both in optical and X-ray (Holoien et al., 2015; Miller et al., 2015). Swift J164449.3+573451 (Sw 1644+57; Bloom et al., 2011; Levan et al., 2011; Burrows et al., 2011; Zauderer et al., 2011), Swift J2058.4+0516 (Sw 2058+05; Cenko et al., 2012a,b) and Swift J1112.2-8238 (Brown et al., 2015) are three long lived (∼ 107 s) relativistic TDE candidates that first discovered by γ-ray triggers and then observed in X-rays. The decline rate of the X-ray light curves in Sw 1644+57 and Sw 2058+05 are approximately consistent with t−5/3. All of these three candidates are believed to have launched relativistic jets (Giannios and Metzger, 2011; Brown et al., 2015) which are very rare in TDEs and have not been observed in any other candidates yet. Recently, broad emission lines have been observed in some TDE candidates. SDSS TDE2 showed a broad Hα line (Van Velzen et al., 2011) and PS1-10jh exhibited broad He II lines with H-lines lacking (Gezari et al., 2012). These were the first objects to show spectral features among TDE candidates. This low ratio of H:He lines has been the matter of debate in recent studies. A number of explanations have been proposed, including the disruption of H-poor star with He-rich core (Gezari et al., 2012; Bogdanovic´ et al., 2014; Strubbe and Murray, 2015); the effects of the optical depth (Gaskell and Rojas Lobos, 2013); the presence of a H-rich BLR which is truncated (Guillochon et al., 2014); the presence of a H-rich envelope surrounding the accretion disk which absorbs and re-emits the radiation
from the central parts (Roth et al., 2016); H-burning in the stellar core of a M ≥ M star (Kochanek, 2016). Arcavi et al (2014) discovered new TDEs with wide range of H:He ratio whereas other candidates show no sign of emission lines (e.g., Chornock et al., 2013; Cenko et al., 2012a). Transients that are detected in the galactic nucleus can be related to the AGN activity. As mentioned in section 1.1.4 AGNs are known to show variability in their continuum and emission lines (Peterson, 2001) that is believed to be because of change in accretion rate 37 related to disk instability or thermal fluctuations (Kawaguchi et al., 1998; Pereyra et al., 2006; Kelly et al., 2009). Also, as discussed in section 1.2, there are newly discovered class of objects called CL AGNs that show appearing or disappearing broad emission lines that is usually followed by almost an order magnitude change (increase or decrease) in the continuum brightness (Storchi-Bergmann et al., 1993; Shappee et al., 2014; LaMassa et al., 2015; MacLeod et al., 2016; Gezari et al., 2017). The timescales of these changes are on the order of years. However, recently iPTF16bco exhibited a change from a LINER to a broad line quasar in a shorter timescale of . 1 yr and is thought to be due to intrinsic accretion changes (Gezari et al., 2017). These recently discovered class of objects can be interpreted as TDEs or vice versa. For example, Blanchard et al.(2017) favor a TDE as the source of the transient PS16dtm over AGN variability. They argue that the rise time for PS16dtm is 50 days, which is much shorter than viscous timescale (∼ 100 years). Also, rapid two order of magnitude change in UV/Optical flares in not typically seen in AGNs, since the order of magnitude variability in CL AGNs happens in timescale of years (Grupe et al., 2010; Drake et al., 2011). On the other hand, they argue that the dimming in the X-rays cannot be explained by CL AGN scenario, since as observed in NG 2617 and iPTF16bco, as a result of AGN accretion stage changes an increase in both X-ray and UV/optical is expected to be seen. The upcoming sky surveys will discover hundreds or thousands of TDEs. These surveys include SKA in radio (Donnarumma et al., 2015), Large Synoptic Survey Telescope (LSST) in optical (Van Velzen et al., 2011), LOFT in hard X-ray (Rossi et al., 2015) and possibly China’s Einstein Probe in soft X-ray (Yuan et al., 2014). These detections will help us to better understand the nature of TDEs and study SMBHs in the galaxies that are rather quiescent. 38
1.4 Supernovae
A supernova (SN) is the explosion of a star at the end if its life when the equilibrium between the gravitational pull from the mass of the star and radiation pressure caused by fusion is lost. SNe are classified to two major classes: Type I or hydrogen poor and Type II or Hydrogen rich based on the presence of the hydrogen features in the observed spectra of these events. Each of these general classes are again classified to different sub-classes based on the spectral features in the case of Type Is and light curve decay behaviour in the case of Type IIs. Type I is classified to Type Ia with strong silicon features, Type Ib with strong helium features, and Type Ic with no helium or silicon detected in their spectra. While, Type II is classified to Type IIL with linear decay, and Type IIP for a plateau in the light-curve after the peak of the outburst. Type II contains another subclass called Type IIn where strong narrow hydrogen emission and absorption lines are present in the spectra. Most of the Type II SNe are powered by circumstellar interaction (CSI) or hydrogen recombination while Type Is are powered by radioactivity of 56 Ni. Furthermore, majority of these events are caused by core collapse supernovae (CCSNe) as a part of the evolution of a star with zero-age main sequence (ZAMS) mass of > 8M where a compact object such as s neutron star or black hole is produced. On the other hand, Type Ia is the only exception that is caused by thermonuclear explosion of a white dwarf that gains mass and passes
Chandrasekhar limit (∼ 1.4M ) through binary interaction and accretion and explodes. Studying SNe can help us to better understand stellar evolution, properties of the circumstellar material (CSM), and the environment of the galaxy that hosts them.
1.4.1 Superluminous Supernovae
A number of outbursts have been observed that were 10-100 times brighter than classical SNe. Because of these intense brightness difference they are called superluminous supernovae (SLSNe). Their luminosity at the peak of the outburst is ∼ 1044 − 1045 erg s−1 39 that makes them detectable at redshifts z > 4 (Berger et al., 2012; Cooke et al., 2012; Quimby et al., 2011). These rare outbursts are observed in faint host galaxies (Angus et al., 2016; Lunnan et al., 2014; Schulze et al., 2018) with low metallicity and with the rate of 10−4 − 10−5 of the CCSN rate (McCrum et al., 2015). The classification scheme of SLSNe is similar to the SNe. Type I or hydrogen poor SLSNe are very similar to SN Ic with strong oxygen lines in their spectra (e.g., SN 2015bn; Jerkstrand et al., 2017; Nicholl et al., 2016b,a; Bhirombhakdi et al., 2018). Type II or hydrogen rich SLSNe typically show narrow emission and absorption with prominent Hα emission lines similar to SNe IIn that are similarly called SLSNe IIn (e.g., SN 2006gy; Quimby, 2006). However, in some rare cases the narrow lines are not detected that are called SLSNe II (e.g., SN 2008es; Gezari et al., 2008; Miller et al., 2008; Inserra et al., 2018; Bhirombhakdi et al., 2019).
1.4.2 Sources that Power Superluminous Supernovae
Based on their types SLSNe can be powered by different sources. The SLSNe IIn are thought to be powered by CSI very similar to SNe IIn while the power sources of the SLSNe I and II are not well known and under debate. Radioactive 56Ni, and magnetar spindown are other possible sources for both SLSNe I and II. In this section we will mainly discuss CSI as the power source of the SLSNe IIn, since we discuss an object detected by PS1 survey in chapter3 where SLSN IIn can be one of the possible interpretations for this outburst.
1.4.2.1 Circumstellar Interaction
The Type II SNe are caused by core-collapse of massive stars and come in different luminosity ranges. The more luminous ones are typically being interpreted as Type IIn SNe because of the presence of the relatively narrow Balmer emission lines that are believed to be originated from the strong interaction of the shocks with the CSM (Schlegel, 1990). 40
These shocks are generated from the interaction of the accelerated ejecta (with velocities of ∼ 104 km s−1) with CSM and consist of a forward (FS) and backward component (RS). The emission lines in the spectra of SNe IIn typically include multiple components: a narrow component with a width of ∼ 102 km s−1, an intermediate component with a width of ∼ 103 km s−1, and a broad component with a width of 104 km s−1. The narrow components of the Balmer emission lines are believed to be originated from the pre-shock CSM that was shed by a progenitor star while intermediate components and broad components of the emission lines are originated from the post-shock dense gas known as cool dense shell (CDS) and ejecta, respectively (Chugai and Danziger, 1994). The evolution of a limited number of SLSNe IIn candidates has been studied from the early time of the explosion to years after. They show blue spectra with few emission lines with simple profiles at early times. At the same time, later time spectra tend to have more complex continuum shape and emission line profiles and after several hundred days the spectra show very little continuum dominated by residual Hα emission lines (Smith et al., 2008a,b, 2010). When the CSM is sufficiently dense, it decelerates the shock waves thus converting their bulk kinetic energy to X-rays and then optical radiation that significantly increases the bolometric luminosity of the system (Fransson, 1984). The presence of such dense CSM suggests a more massive progenitor or a progenitor with a slow wind (e.g., Chugai et al., 2004; Smith et al., 2008a). Another source of radiation in the CSI is from synchrotron and inverse Compton (IC) caused by electrons that are accelerated by shocks to relativistic speeds. The radiation generated by IC is more significant at early times since the density of high energy photons is higher (Margutti et al., 2018). The radiation caused by upscattering of optical photons in the SN photosphere is in UV, X-rays and low energy γ rays. The peak of synchrotron radiation can range from X-rays to radio which depends on the time of generation after the explosion. It is notable that the synchrotron radiation is mostly self- 41 absorbed at early times and is observed to peak in radio only at late times (Pooley et al., 2002). The CSI has been studied through analytical models (Chatzopoulos et al., 2012; Chevalier, 1982; Chevalier and Fransson, 1994, 2017; Chevalier and Liang, 1989; Fransson, 1984; Nadyozhin, 1985) where the CSI is typically represented by self-similar solutions with a power-law density profile; however, these models’ assumption of symmetric CSM structure is questionable since there is evidence that in some cases CSM is in the shape of a disk or a torus with clumpy distribution (Andrews et al., 2010, 2011, 2016; Chugai et al., 1995; Fransson et al., 2002; Smith et al., 2009). Also, most of the parameters (e.g, energy, mass, velocity, and density profiles of the ejecta) in these models are uncertain and model sensitive (Chatzopoulos et al., 2012). The CSI is more accurately simulated through 3D hydrodynamic simulations (Balberg et al., 2000; Bersten et al., 2011; Blinnikov et al., 1998; Chieffi et al., 2003; Dessart and Hillier, 2010; Kasen and Woosley, 2009; Pumo and Zampieri, 2011; Utrobin, 2007; Woosley and Weaver, 1995; Zampieri et al., 1998).
1.5 The Modular Open Source Fitter for Transients
The Modular Open Source Fitter for Transients (MOSFiT) is a Python 2.7/3.x package that is installed in the Anaconda environment. This package is used for fitting, estimating, and sharing the parameters of transients such as but not limited to SLSNe, SN Ia and Ic, TDEs, and Kilonovas (Guillochon et al., 2018). It is designed to work with datasets with the data format of Astrocats: Open Astronomy Catalogs 1. The data can be downloaded directly by running MOSFiT on the name of a specific object. In order to use private data, the data should be transformed to a JSON 2 format that matches the data format of Astrocats. 1 https://github.com/astrocatalogs/astrocats 2 https://www.json.org/ 42
The MOSFiT package gives a lot of flexibility to the users. The package runs through the command line and users can specify the transient model, model parameters, and extra outputs such as luminosity and temperature. After running the package, each of the modules take the inputs and process them to return an output such as extinction correction, and SED that is going to be used by other modules. The MOSFiT contains an ensemble of Markov Chain Monte Carlo (MCMC) walkers that follow a Bayesian approach. The objective function is a Gaussian Process Likelihood where the Watanabe- Akaike information criterion (WAIC) is being applied to calculate this objective from the ensemble. The model converges when the potential scale reduction factor (PSRF) goes under a desired value. The standard and extra outputs of the model are then saved in JSON format and can be visualized through a provided Jupyter Notebook. The MOSFiT tries to fit each band pass separately. However, MOSFiT heavily relies on the SED to calculate fluxes through luminosity at each epoch in each band-pass because all of the models are represented as bolometric light curves. In all of the models a single blackbody SED is being used except the model for SLSNe where a modified blackbody SED is being applied. We use MOSFiT to fit for TDE and SLSN IIn scenarios in Chapter3.
1.6 Galaxy Morphology Prediction
Galaxies have a variety of shapes, colors, and sizes that are used to indicate their evolution and interactions during the course of cosmic history. Galaxy morphology studies can help us to better understand the evolution of galaxies with redshift and provide us information about the dynamical history of galaxies without the need of expensive spectroscopy. In order to carry out these studies, all sky surveys and an accurate method for morphological classification is crucial. Recently, large sky surveys such as SDSS made a large amount of data for objects in our universe available. Traditionally, the morphological 43 classification of the galaxies has been done by visual inspection by experts that is inefficient and impractical for the large data available from current surveys and even larger upcoming datasets from surveys such as LSST (Ivezic et al., 2008). Different automated methods have been proposed to solve this issue; however, historically these methods did not reach required reliability level for scientific studies (Clery, 2011). The Galaxy Zoo project was launched to help to accelerate finding a reliable classification method by introducing a crowd-sourcing method (Lintott et al., 2008). During the initial Galaxy Zoo project ∼ 900, 000 galaxies were annotated by online participants in the short time span of months. Since the successful initial project, there have been other iterations of the Galaxy Zoo projects with more fine grained classification scheme (e.g., Galaxy Zoo 2; Willett et al., 2013). For example, during the Galaxy Zoo 2 project participants were asked questions from a decision tree where the answer to each question would determine the next question to be asked from this decision tree (Figure 1; Willett et al., 2013). Then, the answers of the participants on the same image were transformed to vote fractions which were used to determine the confidence levels of the answers. The data from Galaxy Zoo projects have already been used for studies regarding galaxy structure, formation, and evolution (Skibba et al., 2009; Bamford et al., 2009; Schawinski et al., 2009; Lintott et al., 2009; Darg et al., 2010; Masters et al., 2010, 2011; Simmons et al., 2014; Melvin et al., 2014; Willett et al., 2015). Recent improvements in computer vision through deep neural networks (e.g., Krizhevsky et al., 2012), available large annotated datasets through Galaxy Zoo projects, and increased computing power through the GPUs have made the possibility of providing a more accurate automated approach more promising. Different machine learning methods such as Principal Component Analysis (PCA) (Naim et al., 1995; Lahav et al., 1995; De La Calleja and Fuentes, 2004), Support Vector Machines (SVM) (Tasca et al., 2009), feed forward neural network (Storrie-Lombardi et al., 1992), and convolutional neural 44 networks (Gravet et al., 2015; Dieleman et al., 2015; Dom´ınguez Sanchez´ et al., 2018) have been used for galaxy morphological predictions. In the the following section we will provide an overview of deep neural networks.
1.6.1 Deep Learning
Deep neural networks try to build models to extract hierarchical and abstract features and discover accurate symmetries directly from the data. These networks consist of multiple layers where each layer extracts more abstract features based on the information available from lower layers that represent the input data. These abstract features are typically calculated by a non-linear transformation function. Next, the parameters of this multi-layer structure are optimized by training the model over a large dataset. Although neural networks have been discovered decades ago (McCulloch and Pitts, 1943; Fukushima, 1980), only recent the advancements in GPU computing and advances in the underlying methods enabled training them on large-scale datasets.
1.6.1.1 Feed-Forward Neural Networks
One simple example of deep neural networks is a structure called Feed Forward Neural Network (FFFNN) that consists of three kinds of layers: input layer, hidden layer, and output layer (a representation of FFFNN can be found in Figure 1.1). Each layer is connected to the next layer via a set of initialized weight parameters. The input layer is just the layer that feeds the data to the network. Next, the network has multiple layers of what are called hidden layers where each of these layers is made of a number of neurons that compute a weighted summation of the incoming inputs followed by a non- linear transformation. Let al be the input of the layer l + 1, Wl to be matrix of weights that connect layer l to layer l + 1, and bl to be a vector of biases. Therefore, the output of layer l + 1 (al+1 ) can be 45 represented as, sl l+1 l+1 l+1 X l l l ai = f (Zi ), Zi = Wi ja j + bi (1.7) j=1
where i and j represent number of neurons in layer l and l + 1 respectively, sl represents number of layers being used in FFNN, and f is the activation non-linear function. The most common choices for activation function are Rectified Linear Unit (ReLU, f (x) = max(x, 0);
1 Nair and Hinton, 2010), sigmoidal function ( f (x) = 1+e−x ), and hyperbolic tangent ( f (x) = tanh(x)). During the training process, the objective is to optimize parameters of the network in a way that the output generated by network (anl ) closely resembles the true values of the desired target output (y). This prediction error is quantified by a cost function J(anl , y). Optimization of this error function is achieved by a process called gradient descent where the weights and biases are repeatedly updated by taking infinitesimal steps in the opposite direction of the gradient of the cost function with respect to updating parameters,
∂J(anl , y) Wl = Wl − η , (1.8) ∂Wl ∂J(anl , y) bl = bl − η (1.9) ∂bl where η is the learning rate which is a hyperparameter that determines the step size at each epoch of training. It is worth noting that models with many hidden layers were not commonly used, because the gradient information would vanish during the back propagation process which made lower level parameters a difficult task (Hochreiter et al., 2001). However, the introduction of the ReLU activation function (Nair and Hinton, 2010) significantly reduced the gradient vanishing problem and made training large networks computationally efficient. Also, the introduction of the dropout method (Hinton et al., 2012; Srivastava et al., 2014) made training even deeper and larger networks with many more parameters possible. The dropout is a regularizer that when applied on a layer l+1 randomly removes (or sets to zero) the output values of the previous layer (l) with probability p. Therefore, at each epoch of 46 training a different subset of output values from layer l is being removed. This will force each neuron in layer l+1 to learn useful features without heavily relying on other neurons in the same layer since they can be randomly removed (see Section A.1 for more information about FFFNNs).
HL2 (al+1) Output layer Input layer HL1 (al)
X1
X2
X3 anl
X4
Xn
Figure 1.1: A representation of a feed forward neural network where HL stands for hidden layer. All of the layers are fully connected to each other.
1.6.1.2 Convolutional Neural Networks
Convolutional neural networks (CNNs; Fukushima, 1980; LeCun et al., 1998) are another class of neural networks that extract hierarchical features from the input like FFNNs but with reduced connectivity and shared weights that decrease the number of weights substantially. They consist of three different types of layers: convolutional layers, pooling layers with shared weights, and fully connected layers. The convolutional layer extracts a local combination of features from a set of feature maps from the previous layer and convolves them using a set of filters to a new set of feature maps. Convolution operation 47 can be written as, K l+1 l+1 l+1 X l l l an = f (Zn ), Zn = Wmn ∗ am + bn (1.10) m=1 where K is the number of input feature maps in the form of matrices, ∗ represents the
l l linear convolution operation, Wmn represent the filters of layer l+1, and bn represent bias of feature map n. Each neuron in layer l + 1 is only connected to a local region (receptive field with size of F × F) in the input volume from layer l that is replicated to cover the entire input by shifting the receptive field (see Section A.1 for more information.). As we can see, the number of parameters reduces significantly when we replace the matrix multiplication in a fully connected layer shown in equation 1.7 with summation of convolutions of features map from previous layer in equation 1.10. This parameter reduction significantly improves the generalization performance of the network where the network performs better on unseen samples. The pooling layer is located after each convolutional layer that reduces the spatial dimensions of the input feature maps by applying a max or mean function (Boureau et al., 2010). This spatial reduction enables the convolutional layer to exploit larger regions of the input while reducing the number of parameters since pooling layers do not have trainable parameters. Therefore, they increase computation speed and reduce the chances of overfitting. State of the art CNNs are constructed of multiple convolutional layers followed by pooling layers which allow the deeper layers to extract more abstract representations of the input by allowing them to exploit a larger part of the input. We show a schematic representation of a convolutional layer followed by a pooling layer in Figure 1.2. The fully connected layers are added after multiple convolutional and pooling layers and the procedure is the same as mentioned in the Section 1.6.1.1. Because of the mentioned advantages and demonstrated strength for classification of large scale datasets (Krizhevsky et al., 2012), they have been used in many different computer vision problems. 48
Figure 1.2: A schematic representation of a convolutional layer followed by a pooling layer. Each neuron in layer l + 1 is only connected to a local region of feature maps from layer l. The pooling layer applies a mean or max function on the local region of all feature maps from the convolutional layer l + 1. 49 2 PS1-13cbe:The Rapid Transition of a Seyfert 2 to a
Seyfert 1
Originally published in (Katebi et al., 2019a):
Monthly Notices of the Royal Astronomical Society, Katebi, R., Chornock, R., Berger, E., Jones, D.O., Lunnan, R., Margutti, R., Rest, A., Scolnic, D.M., Burgett, W.S., Kaiser, N. and Kudritzki, R.P., 2019. PS1-13cbe: the rapid transition of a Seyfert 2 to a Seyfert 1. Monthly Notices of the Royal Astronomical Society, 487(3), 4057-4070. In this chapter, we present a nuclear transient event, PS1-13cbe, that was first discovered in the Pan-STARRS1 survey in 2013. The outburst occurred in the nucleus of the galaxy SDSS J222153.87+003054.2 at z = 0.12355, which was classified as a Seyfert 2 in a pre-outburst archival SDSS spectrum. PS1-13cbe showed the appearance of strong broad Hα and Hβ emission lines and a non-stellar continuum in a Magellan spectrum taken 57 days after the peak of the outburst that resembled the characteristics of a Seyfert 1. These broad lines were not present in the SDSS spectrum taken a decade earlier and faded away within two years, as observed in several late-time MDM spectra. We argue that the dramatic appearance and disappearance of the broad lines and factor of ∼ 8 increase in the optical continuum is most likely caused by variability in the pre-existing accretion disk than a tidal disruption event, supernova, or variable obscuration. The timescale for the turn-on of the optical emission of ∼ 70 days observed in this transient is among the shortest observed in a “changing look” active galactic nucleus.
2.1 Introduction
The axisymmetric unification model of active galactic nuclei (AGN) suggests that there can be different types of AGNs based on the angle toward the line of sight. However, 50 it is believed that the simple orientation-based model cannot fully explain the Type 1 (both narrow and broad emission lines are present in the spectra) and Type 2 (only narrow emission lines are present in the spectra) classification scheme because this model is challenged when different AGN types are observed in the same object at different epochs of time. One alternative suggestion is that some Type 2s were Type 1s and their engine is now not active anymore or basically turned off; therefore, the Type 2s are evolved version of the Type 1s (Penston and Perez´ , 1984; Runnoe et al., 2016). This alternative view was first suggested based on the observations of broad emission lines disappearance in the Type 2 NGC 4151 (Lyutyj et al., 1984; Penston and Perez´ , 1984) and 3C 390.3 (Penston and Perez´ , 1984). The term “changing-look” (CL) was first used for the AGNs that showed X-ray absorption variations (Matt et al., 2003; Bianchi et al., 2009; Puccetti et al., 2007; Risaliti et al., 2009; Marchese et al., 2012). Recently, this term has also been used to describe the type of AGNs that have been observed to optically transition from Type 1 to Type 1.8, 1.9 and 2 or vice versa, by the sudden appearance or disappearance of broad emission lines such as the broad Hβ emission line (MacLeod et al., 2016; Runnoe et al., 2016). In some objects, the broad emission lines disappeared and the continuum faded completely (Collin- Souffrin et al., 1973; Tohline and Osterbrock, 1976; Sanmartim et al., 2014; Denney et al., 2014; Barth et al., 2015; Runnoe et al., 2016), while in others broad emission lines appeared or, in other words, the AGN “turned-on” (Cohen et al., 1986; Storchi-Bergmann et al., 1993; Aretxaga et al., 1999; Eracleous and Halpern, 2001; Shappee et al., 2014). These objects are mostly at low redshift with low absolute luminosities. However, recently LaMassa et al.(2015) discovered a luminous changing look AGN with a redshift of z = 0.31 that transitioned from Type 1 quasar to a Type 1.9 AGN in ∼ 9 years. Subsequent studies (MacLeod et al., 2016; Runnoe et al., 2016; MacLeod et al., 2018) have started finding larger samples in the Sloan Digital Sky Survey (SDSS). Most recent, Gezari et al.(2017) 51 discovered a quasar with a rapid “turn-on” timescale of < 1 year that demonstrated one of the fastest changes of state to date. The luminosities of Seyfert galaxies have been observed to vary over time. Variability has been detected in a wide range of wavelengths with timescales from hours to years with the shortest in X-rays that show variation time-scales of a few hours in Seyfert galaxies (Green et al., 1993; Nandra et al., 1997). The observed average variability time-scale in the optical in Seyfert galaxies is on the order of weeks to months (Kaspi et al., 1996; Giveon et al., 1999). Accretion disk instabilities are the most promising model as the source of this variability (Kawaguchi et al., 1998; Kelly et al., 2009). Here, we report the rapid turn-on (∼ 70 days) of a nuclear transient, PS1-13cbe, from a galaxy at redshift z = 0.12355 classified as a Seyfert 2 in a pre-event SDSS spectrum which was accompanied by the appearance of broad Balmer lines and hence a transition to a Seyfert 1. This represents the most rapid “turn-on” in a changing look AGN to date. We discuss the observations in Sections 2.2 and 2.3. We discuss the possible scenarios for the origin of the variations observed in PS1-13cbe in Section 2.4. We summarize and present our conclusions in Section 2.5. Throughout this work we are assuming a standard
−1 −1 ΛCDM cosmology with H0 = 69.6 km s Mpc , Ωm = 0.286 and ΩΛ = 0.714 parameters
(Bennett et al., 2014) that yields a luminosity distance of dL = 582 Mpc. All magnitudes are in the AB system and all dates are UT. We assume a Galactic extinction value of E(B − V) = 0.06 (Schlafly and Finkbeiner, 2011).
2.2 Observations of PS1-13cbe
The Medium Deep Survey (MDS) of the Pan-STARRS1 (PS1) sky survey performed daily (in-season) deep monitoring of ten ∼ 7 sq. deg. fields over the years 2010–2014 to find transient and variable sources. The typical observation sequence was composed of
gP1 and rP1 bands on the first night, iP1 on the next night, and then zP1 on the third night. 52
This pattern was repeated during the ∼ 6 months of the observing season and was only interrupted by the weather and times near full moon, when observations in the yP1 filter were taken. A more complete description of the survey and photometric system were given by Tonry et al.(2012) and Chambers et al.(2016). On 2013 July 9, we detected a transient event, PS1-13cbe, at coordinates α = 22h21m53.86s, δ = +00◦30054.5600 (J2000) coincident with the nucleus of a galaxy using the photpipe transient discovery pipeline, described by Rest et al.(2014) and Scolnic et al.(2018). The host galaxy was observed as part of SDSS and given the name SDSS J222153.87+003054.2 (hereafter SDSS J2221+0030), with a spectroscopic redshift of z = 0.12355 (Ahn et al., 2014). We constructed a difference image light curve for PS1-13cbe using the PS1 transient pipeline (Scolnic et al., 2018). The template images were created from a stack of high-quality observations excluding the observing season containing the outburst (the year 2013) and were then subtracted from all observations of the transient. It is important to note that this photometry represents a flux difference relative to the host contribution present in the template (which is consistent with the SDSS photometry; Ahn et al. 2014). We include the PS1 photometry in Table 2.1.
2.2.1 Optical Photometry
As shown in Figure 2.1, the light curves of PS1-13cbe were constant and consistent with zero change in flux relative to the baseline in the template for three observational seasons and then showed a rise peaking at MJD 56512.6, decline, and a second rise
(gP1 = 19.5 mag at peak; see Figure 2.2). The gP1 band luminosity showed a rise of about ∼ 1.5 × 1043 erg s−1 from the base luminosity in the course of ∼ 70 days and then declined to ∼ 0.7 × 1043 erg s−1 in ∼ 50 days and rose back up again to ∼ 1 × 1043 erg s−1 in the course of next ∼ 50 days, at which point the MDS ended. 53
15 y
) z −1 10 i r g erg s
42 5 (10 ν L
ν 0 S −5 −1200 −1000 −800 −600 −400 −200 0 200 Time (MJD − 56512.6) Figure 2.1: The observed transient luminosities of PS1-13cbe from the PS1 survey in
grizyP1 filters after correction for Galactic extinction. S: marks the epoch of the LDSS spectrum (MJD 56570).
With the SDSS value for the quiescent host flux, g0 = 19.21 ± 0.01 mag (Ahn et al.,
2014), and gP1 = 19.5 mag for PS1-13cbe at the peak of the outburst, we can see that
the total luminosity of the galaxy increased by & 75%. However, the quiescent host flux value is dominated by star light. We estimate that the central AGN contributes .10% of the continuum flux in the quiescent spectrum taken by SDSS because the absorption lines from star light are not noticeably diluted by a non-stellar continuum, so the amplitude of the outburst from the AGN must be significantly larger, as discussed below. After the MDS ended, we obtained late-time photometry in g0 using the MDM4K and Templeton detectors on the 2.4 m Hiltner and 1.3 m McGraw-Hill telescopes at MDM Observatory on 2015 November 18 and 2017 June 18, respectively. The total magnitude of 54
15 16
17
18 y − 4 z − 3 i − 2 19 r − 1 g 20 Magnitude + constant 21 −100 −50 0 50 100 Time (MJD − 56512.6)
Figure 2.2: Light curves of PS1-13cbe in grizyP1 bands, corrected for Galactic extinction. Note the small-scale fluctuations in the light curve (such as near −20 d).
the host galaxy (including any possible transient contribution) was g0 = 19.23 ± 0.05 mag (g0 = 19.22 ± 0.13 mag) at 831 (1408) days after the peak, which are consistent with the SDSS pre-outburst photometry. Therefore, the system returned to the baseline flux value in
. 2 yr after the outburst.
2.2.2 Observations of the Host Galaxy
In addition to the SDSS observation of the host (∼ 10 years before the outburst), the source was detected by the AllWISE survey in the W1, W2, W3, and W4 bands with cataloged values of 17.82 ± 0.037, 18.20 ± 0.34, 16.65 ± 0.23 and 15.13 ± 0.37 magnitudes, respectively (∼ 3 years before the outburst; Chang et al., 2015). We also obtained photometry using the Neil Gehrels Swift Observatory (Gehrels et al., 2004) with 55 the UV Optical Telescope (UVOT; Roming et al., 2005) in the u, uvw1, and uvw2 filters on 2016 November 27–29 (∼ 3 years after the outburst) and we measure values of 21.10±0.13, 21.76 ± 0.21 and 22.41 ± 0.19 magnitudes, respectively.
2.2.3 X-ray Photometry
We obtained an X-ray observation of SDSS J2221+0030 using the X-ray telescope (XRT; Burrows et al., 2003) on Swift with a total exposure time of 7.1 ks between 2016 November 27 and November 29. The host was not detected with a 3σ upper limit of 2.13 × 10−3 cts s−1 (0.3 − 10 keV). Using the Galactic neutral hydrogen column density in
20 -2 the direction of PS1-13cbe of NH = 4.49 × 10 cm (Kalberla et al., 2005) and assuming no intrinsic absorption, and a typical photon index of Γ = 2, we calculate the unabsorbed X-
−14 −2 −1 42 ray flux to be fx(2 − 10 keV) ≤ 4.0 × 10 erg cm s , which translates to LX ≤ 1.6 × 10 erg s−1. The absorbing column densities of Seyfert 2 galaxies range from 1022 to 1025 cm-2 (e.g., Risaliti et al. 1999) and therefore, assuming a minimum intrinsic column density of 1022 cm-2 for the host galaxy in addition to the Galactic value, we recalculate the
−13 −2 −1 42 unabsorbed X-ray flux to be fx(2 − 10 keV) ≤ 1.0 × 10 erg cm s or LX ≤ 4.1 × 10
−1 erg s . Additionally, using LX and the empirical bolometric correction from Marconi et al. (2004), this corresponds to a minimal upper limit on the bolometric luminosity of the object
44 −1 Lbol ≤ 0.6 × 10 erg s 3 years after the outburst, although this value can increase if the
45 −1 intrinsic absorption is higher (e.g., Lbol ≤ 9.0 × 10 erg s for an intrinsic column density of 1024 cm-2). We did not find any archival X-ray observations of the object before or at the time of the outburst.
2.2.4 Optical Spectroscopy
A pre-outburst spectrum of SDSS J2221+0030 was obtained on 2003 May 26 by SDSS (Ahn et al., 2014). After the detection by PS1, we observed PS1-13cbe on 2013 56
October 5 (+57 days after peak) for 1200 s using the Low Dispersion Survey Spectrograph- 3 (LDSS3) on the 6.5 m Magellan Clay telescope. We used a 100-wide long slit with the VPH-all grism to cover the observed wavelength range 3700 − 10000 Å with a resolution of ∼ 9 Å. At late times, we obtained five epochs of spectroscopy using the Ohio State Multiple Object Spectrograph on the 2.4 m Hiltner telescope at MDM Observatory (OSMOS; Martini et al., 2011). The first two spectra were taken on 2015 October 3 and 2015 November 12 using the 1.200 center slit with a VPH-red grism and an OG530 filter that covered 5350 − 10200 Å (resolution = 5 Å) in the observed frame. No significant spectral differences were present, so we combined these two observations for all subsequent analysis. The next two spectra were obtained on 2015 December 1 and 2016 November 16 using the same 1.200 center slit and a VPH-blue grism that covered 3675−5945 Å (resolution = 2 Å) in the observed frame. The final spectrum was taken on 2017 June 18 using a 1.000 outer slit and VPH-red grism that covered 3930−9050 Å (resolution = 4 Å) in the observed frame. We preprocess our data using standard procedures such as flat-fielding, bias subtraction, and wavelength calibration using arc lamps in IRAF 3. Additionally, we removed cosmic rays using the L.A.Cosmic (Van Dokkum, 2001) task. Finally, we calibrate our data using our own IDL procedures and observations of the standard stars BD+174708 for red spectra and Feige110 for blue spectra. All of the spectra for PS1-13cbe are shown in Figure 2.3.
3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. 57
3 2017/06/18 ) −1 Å −1
s 2 2016/11/16 −2
2015/12/01 /
erg cm 2015/10/03
−16 1 (10 λ 2013/10/05 Logf 0 2003/06/26
−1 3000 4000 5000 6000 7000 8000 Rest frame wavelength (Å)
Figure 2.3: Optical spectra for PS1-13cbe, from bottom to top: host galaxy of PS1-13cbe from SDSS (grey), model of the host galaxy generated by FAST 1.0 (black), spectrum during outburst observed using LDSS3 (purple), spectra obtained with OSMOS in red and blue for the red and blue setups, respectively, and the most recent OSMOS spectrum (olive green).
2.3 Observational Features of PS1-13cbe
2.3.1 Host galaxy of PS1-13cbe
To isolate the emission lines of SDSS J2221+0030 and correct for stellar absorption lines, we simulated a galaxy model with the FAST 1.0 code (Kriek et al., 2009) using the archival SDSS spectrum and optical photometry in ugriz bands along with the Swift u, uvw1, and uvw2 photometry. We experimented with the initial parameters and generated the best fit with the stellar age of 3 × 109 years and e-folding timescale of τ ≈ 109 years, by 58 assuming the star-formation history to be exponentially declining, the stellar initial mass function from Chabrier(2003), the Bruzual and Charlot(2003) spectral library, a Milky Way dust law (Cardelli et al., 1989), and solar-like metallicity of Z = 0.02. The FAST model of the host galaxy is over-plotted on the original SDSS spectrum (grey) in the bottom of the Figure 2.3 (black). After subtracting this model from the pre-outburst optical spectrum of the host galaxy of PS1-13cbe (bottom of Figure 2.3 in grey), we fit the profiles of the narrow emission lines [S II]λλ6717, 6731, [O III]λλ4959, 5007, [O II]λ3727, [O I]λ6300, [N II]λλ6549, 6583, Hα, and Hβ. It has been shown that a model for narrow emission lines can be obtained from [S II] lines in most cases (Ho et al., 1997). By using this fact, we constrained parameters of the Gaussian profiles of the other narrow lines. The lines were not well modeled with single Gaussians, so we used double Gaussian profiles to fit the line profile of all the narrow lines except [O II]λ3726 and [O I]λ6300. In Figure 2.4, we plot the line ratios in excitation diagrams (Baldwin et al., 1981). We also show the extreme star formation line (Kewley et al., 2001, 2006), the pure star formation line (Kauffmann et al., 2003) and the Seyfert- LINER classification line (Kewley et al., 2006; Fernandes et al., 2010). Additionally, we show 30000 randomly selected galaxies (shaded area in Figure 2.4) with emission-line fluxes from the MPA-JHU DR7 catalog 4 (Aihara et al., 2011). The automatic Portsmouth pipeline from SDSS classified the host galaxy of PS1-13cbe as a LINER (Sarzi et al., 2006); however, our emission line ratios calculated from the SDSS spectrum after subtraction of the stellar continuum (navy squares) classify the host galaxy as a clear Seyfert (Figure 2.4). Next, we constructed the spectral energy distribution (SED) of the host galaxy using the observed and archival photometry data. We also scaled and plotted SEDs of an Sb spiral galaxy and a Seyfert 2 galaxy from the SWIRE template library (Polletta et al., 2007) alongside the SED of the host (see Figure 2.5). The SED of the Seyfert 2 galaxy
4 http://wwwmpa.mpa-garching.mpg.de/SDSS/DR7/ 59
1.0 AGN Seyfert 0.4 Seyfert
0.2 0.5
0.0 HII ) β
−0.2 0.0 HII LINER LOG([OIII/OII]) LOG([OIII]/H −0.4
HII −0.5 −0.6
−0.8 LINER Comp −1.0 −1.0 −1.0 −0.5 0.0 0.5 −2.0 −1.5 −1.0 −0.5 0.0 −2.0 −1.5 −1.0 −0.5 0.0 LOG([NII/Hα]) LOG([OI/Hα]) LOG([OI/Hα])
Figure 2.4: The excitation diagrams using the [N II]λ6583/Hα, [O I]λ6300/Hα, [O III]λ5007/Hβ and [O III]λ5007/[O II]λ3726 line ratios (Baldwin et al., 1981; Kewley et al., 2006). The navy blue squares display the position of the host galaxy of PS1-13cbe. The shaded area represents the location of the SDSS galaxies calculated using MPA-JHU where darker regions represent higher number density of the galaxies. The extreme star formation line (solid red; Kewley et al., 2001, 2006), the revised star formation line (dotted black; Kauffmann et al., 2003) and the Seyfert-LINER classification line (dashed black; Kewley et al., 2006; Fernandes et al., 2010) are also plotted. Comp: AGN/star forming composites.
fits well from the mid-IR to UV, including the contribution from star light. As we can see in Figure 2.5, the SED of the spiral galaxy fits the optical and UV part of the host SED as well; however, there is a broad excess in the mid-IR (W3 and W4). This excess is because of the contribution of dust heated by the AGN in the host galaxy. After subtracting the host contribution (using the Sb template) from the observed values, we estimate the luminosity 60
43 −1 of the AGN component to be νLν(W4) ≈ (1.63 ± 0.65) × 10 erg s . Then, using a bolometric correction factor of 10.1 ± 1.4 for the W4 band from Runnoe et al.(2012), we
44 estimate the total bolometric luminosity of the AGN to be Lbol(W4) ≈ (1.65 ± 0.65) × 10 erg s−1. Furthermore, narrow emission lines are possible indicators of the intrinsic bolometric luminosity (Netzer, 2009). We estimate the bolometric luminosity of the AGN from the luminosity of the narrow [O III] λ5007 emission line in the SDSS spectrum and used the conversion that Lbol = 3500L(λ5007), which has a variance of 0.38 dex by assuming a
44 standard AGN SED (Heckman et al., 2004). This results in an estimate of Lbol = 1.6 × 10 erg s−1, which is consistent with the bolometric luminosity calculated using the W4 band above. This luminosity can be an overestimate if some of the [OIII] emission originates from star formation. However, as shown in Figure 2.4, SDSS J2221+0030 is classified as a clear Seyfert galaxy and therefore narrow line excitation is dominated by the AGN. Moreover, as we discussed in Section 2.2.3, we estimate the upper limit on the bolometric luminosity in quiescence from the X-ray non-detection to be ≤ 0.6 × 1044 erg s−1, which is consistent with the bolometric luminosity estimated before the outburst in quiescence. This is only an estimate because the X-ray flux is a non-detection, whereas having an X- ray detection would help us to have a better estimate of the intrinsic absorption and column density.
2.3.2 Astrometry
To find the location of PS1-13cbe in its host galaxy, we performed relative astrometry between the PS1/MDS template images of the host galaxy and the position of the transient reported by photpipe. First, we fit a 2-dimensional Gaussian function to the templates in all filters to find the centroid of SDSS J2221+0030 and then, using the weighted average centroid coordinates of the transient, we calculate the offset of PS1-13cbe from the center 61
44 10 Sy2 Sb
) data − 43 10 (ergs ν ν
42 Log 10
12 13 14 15 16 10 10 10 10 10 Log ν (Hz)
Figure 2.5: SED of SDSS J2221+0030 in quiescence (purple circles), scaled SED templates of an Sb spiral galaxy (dotted orange line) and a Seyfert 2 galaxy (dashed blue lines) from the SWIRE template library (Polletta et al., 2007). The host exhibits a mid-infrared (MIR) excess relative to the star-forming template due to an AGN.
of its host galaxy to be 0.03600 ± 0.03500 (101 ± 100 pc), consistent with the nucleus. We used a systematic astrometric error floor of 0.1 pixel to calculate the uncertainty in the offset (Scolnic et al., 2018). We show the position of PS1-13cbe relative to the center of its host in gP1 band in Figure 2.6.
2.3.3 Multi-band Light Curves of PS1-13cbe
We calculate the total optical luminosity light curve of PS1-13cbe using multi-band observations shown in Figure 2.1. This has been done by first finding the epochs where 62
Nucleus Outburst
1′′ = 2.8 kpc
Figure 2.6: The gP1 PS1 template image of SDSS J2221+0030 showing the position of the centroid of the galaxy (green “X”), the position of PS1-13cbe (red dot), and its position uncertainty is shown by a circle with 5σ radius (red circle). The position uncertainty is dominated by the outburst.
gP1 and rP1 bands were observed simultaneously and then interpolating the iP1, zP1, and yP1 bands at these epochs using Legendre polynomials. Next, by integrating the spectral 63 distribution at each epoch and using the trapezoidal rule, we calculate the total flux and thereafter optical luminosity over 3685 − 8910 Å in the rest frame. At late and early times, we had data only from the yP1 and zP1 bands, so we calculate the luminosity assuming the same colour correction as measured from those epochs with all filters. As seen in Figure 2.7, the total optical luminosity of PS1-13cbe rises to a peak value of (1.06 ± 0.01) × 1043 erg s−1 in the course of ∼ 70 days, declines in next ∼ 50 days, and then rises back up again.
12 )
−1 10 z y 8 g−r−i−z−y Trapezoidal erg s 42 6 4 2
Luminosity (10 0 −300 −200 −100 0 100 Time (Rest frame days after peak)
Figure 2.7: Total optical luminosity light curve of PS1-13cbe integrated over the grizyP1 filters and relative to the baseline flux in the template images. The total optical luminosity was estimated using the spectral distribution at each epoch and the trapezoidal rule (blue triangles). The estimated luminosity at early and late times using yP1 (red circles) and zP1 (purple squares) bands assumed the same colour correction as measured from those epochs with all filters. 64
Furthermore, to study the evolution in continuum colour and temperature while remaining agnostic about the overall SED of the transient, we fit both power law and blackbody models. We also calculate the g − r colour and show it alongside the spectral index and blackbody temperatures for all epochs in Figure 2.8. The reason we chose gP1 and rP1 bands as a proxy for the colour is that they were taken on the same night in the PS1 survey and do not need to be interpolated. No strong colour evolution during the outburst is evident in Figure 2.8.
9500 9000 8500 8000 7500 7000 Temperature (K) 6500 0.0 −0.2 −0.4
index −0.6 Spectral −0.8 −1.0 0.3 0.2 0.1 0.0 g−r (mag) −0.1 −0.2 −50 0 50 100 Time (Rest frame days after peak)
Figure 2.8: Top: Rest-frame blackbody temperature from fitting the optical photometry.
α Middle: Rest-frame spectral index from fitting a power law fν ∝ ν to the optical photometry, where α is the spectral index. Bottom: observed g-r colour diagram. All three panels show little or no evolution during the outburst. 65
2.3.4 Spectral Features of PS1-13cbe
In Figure 2.3, we show the optical spectra for PS1-13cbe. The spectra of PS1-13cbe show remarkable evolution over the course of fourteen years. The spectra contain narrow and broad emission line profiles, including hydrogen Balmer lines, [S II]λλ6717, 6731, [O III]λλ4959, 5007, [O II]λ3726, [O I]λ6300 and [N II]λλ6549, 6583. We note that the spectra were taken with different effective apertures (slit and seeing) and different spectral resolutions. We fit the Hα+[N II] and Hβ+[O III] complexes in order to investigate the presence of broad Hα and Hβ lines. First, we scaled our OSMOS (blue setup) and LDSS3 spectra using the flux of [O III] λ5007 from the SDSS spectrum and scaled the OSMOS red setup using the [N II] flux from SDSS. We did these scalings with the assumptions that these narrow lines are centrally concentrated in the host and that the fluxes do not change in such a short timescale because it has been shown that narrow emission lines only slowly vary over decades (Peterson et al., 2013). In addition, the LDSS3 spectrum was taken during the outburst and contains
α transient flux. In that case, first we fit a power law fν ∝ ν to the optical photometry from PS1 at the epoch when the spectrum was taken, where α = −0.58. After that, we constructed a model from a linear combination of the power-law continuum and the host galaxy model shown in Figure 2.3. Additionally, we smoothed this model by a Gaussian with full width at half maximum (FWHM) = 5 Å to better match the resolution of our data and subtracted it from the LDSS3 spectrum, isolating the emission-line spectrum. Next, we fit the [S II] lines using double Gaussian profiles and use this model to constrain the multi-component Gaussian profiles that were used to fit the Hα+[N II] complex to reduce the number of free parameters (Ho et al., 1997). Specifically, we model emission lines using two Gaussian profiles for each narrow line and one broad component for Hα and simultaneously fit for parameters of the broad and and narrow lines. We fix the narrow components to the wavelengths of Hα and [N II] λλ6549, 6583. Also, we fix the 66 widths and relative amplitudes of the two Gaussian components of the narrow lines using the values from the [S II] model, leaving only the overall amplitudes of the narrow Hα and [NII] lines as free parameters. The parameters of the broad component of Hα were allowed to vary freely. We fit the [O III] λλ4959, 5007 lines using two Gaussian profiles for each narrow line without constraining the parameters. For Hβ, however, we use two Gaussian profiles, a single one for the narrow component with just the centroid fixed and allowing the other parameters to vary, and one for the broad component with no constraints. The SDSS and OSMOS spectra lack transient flux, so we only subtract a scaled galaxy model and perform the same procedure as for the LDSS3 spectrum to fit the narrow and broad emission lines, with the exception of fixing the width and the the centroid of the broad component of both Hα and Hβ to the values derived from LDSS3 spectra, allowing only the normalization to vary. We show the resultant Gaussian fits to the Hα+[NII] and Hβ+[O III] emission lines in Figure 2.9 and 2.10 respectively. As we show in Figure 2.9 and Figure 2.10, the existence of the broad Hα and Hβ components is clear in the LDSS3 spectrum (a weak broad Hγ emission line was also detected). Ho et al.(1997) found that they could reliably extract a weak broad H α component that comprised ≥ 20% of the Hα + [NII] blend from spectra with spectral resolution close to ours. As shown in Figure 2.9, there is marginal evidence for a broad component of Hα in the SDSS and OSMOS data, which is 21% and 17% of the flux in the Hα + [NII] complex, respectively. However, it is notable that there are visible wiggles in the residuals of the fit from the spectra that suggest the broad emission line profiles are not purely Gaussian. For this reason, we are not confident that broad Hβ line in the SDSS data and broad lines in the OSMOS are real and are not result of deviations from a Gaussian profile. We report the measured luminosities of the broad components of the Hα and Hβ emission lines in Table 2.2. 67
SDSS Data LDSS3 Data OSMOS Data OSMOS Data Best Fit Best Fit Best Fit Best Fit 4 Individual Individual Individual Individual
) Gaussians Gaussians Gaussians Gaussians
−1 Broad Hα Broad Hα Broad Hα Broad Hα
Å 2003/06/26 2013/10/05 2015/10/03 2017/06/18 −1 3 − 3697 d + 57 days + 785 d + 1409 d s −2
2 erg cm −16 1 (10 λ f
0
0.4 0.2 0.0
Residual −0.2 −0.4 6500 6550 6600 6500 6550 6600 6500 6550 6600 6500 6550 6600 Rest wavelength (Å)
Figure 2.9: Continuum subtracted Hα line profiles. Top: Multiple-component Gaussian fit to the Hα+[NII] emission lines (blue), individual components (dashed green) and broad component of the Hα (dashed black). Bottom: The fit residuals. Top-left (on each panel): Observation date and numbers of days before/after the peak (navy blue).
As mentioned in Section 2.3.1, we classified the SDSS spectrum to be that of a Seyfert. Moreover, the galaxy can be classified as a Seyfert 1.9/2 galaxy since it just shows narrow emission lines and possible weak broad Hα without any sign of a broad Hβ line, as shown in the left panels of Figures 2.9 and 2.10. At the time of the outburst, the spectrum of PS1- 13cbe taken with LDSS3 showed broad Hα and Hβ components with FWHM velocities of 3385 ± 32 km s−1 and 3277 ± 110 km s−1, respectively. These high values of FWHM are the sign of the high velocity, dense, and highly ionised gas clouds in the Broad Line Region (BLR) close to the central black hole (BH) where the broad emission lines originate from. The presence of these broad lines suggests that the galaxy transformed from a Seyfert 1.9/2 68
SDSS Data LDSS3 Data OSMOS Data OSMOS Data Best Fit Best Fit Best Fit Best Fit Individual Individual Individual Individual 1.5 Gaussians ) Gaussians Gaussians Gaussians
−1 Broad Hβ Broad Hβ Broad Hβ Broad Hβ
Å 2003/06/26 2013/10/05 2015/10/03 2017/06/18 −1 − 3697 days + 57 days + 785 days + 1409 d s
−2 1.0 erg cm
−16 0.5 (10 λ f
0.0
0.3 0.2 0.1 0.0
Residual −0.1 −0.2 4800 4900 5000 4800 4900 5000 4800 4900 5000 4800 4900 5000 Rest wavelength (Å)
Figure 2.10: Continuum subtracted Hβ line profiles. Top: Multiple-component Gaussian fit to the Hβ+[O III] emission lines (blue), individual components (dashed green) and broad component of the Hβ (dashed black). Bottom: The fit residuals. Top-left (on each panel): Observation date and numbers of days before/after the peak (navy blue). Strong residuals are present in the red wing of λ4959 due to poor subtraction of the 5577Å sky line.
to a Seyfert 1 galaxy. On the other hand, we did not detect broad emission lines in the spectra taken with OSMOS and that means the host galaxy of PS1-13cbe transformed back from a Seyfert 1 to a Seyfert 2 galaxy in less than 2 years and continues in that state as there are no signs of broad lines in spectra taken with OSMOS ∼ 4 years after the peak of the outburst. We estimate the mass of the central BH using two methods. The first method uses the stellar kinematics and the revised scaling relation between the SMBH mass and stellar velocity dispersion (McConnell and Ma, 2013), where σ = 93.52 km s−1 was provided 69
+3.5 6 from the SDSS spectrum of the host galaxy, and results in an estimate of MBH = 2.9−1.7×10
M . However, this estimation is subject to uncertainties because of the large scatter and lack
of constraints at the low σ and MBH region of the MBH − σ scaling relation, particularly in a case where we do not have a decomposition of the galaxy that separates the bulge component (McConnell and Ma, 2013; Kormendy and Ho, 2013). The mass of the SMBH can also be estimated using photoionisation equilibrium, by applying a mass-scaling relationship based on the FWHM of the broad Hβ emission line and continuum luminosity. Therefore, using the measured FWHM ≈ 3277 km s−1, intrinsic
42 −1 luminosity at 5100Å (λL5100) of (6.2 ± 0.2) × 10 erg s at the time of the spectrum, and the mass-scaling relationship of Vestergaard and Peterson(2006), we calculate the
7 mass of the central BH to be (2.2 ± 0.1) × 10 M . We note that the values used in this calculation are measured during the outburst and can be affected by the changes in the accretion structure and also rely on the assumption that the (unknown) shape of the ionising continuum is similar to those objects used to calibrate the scaling relationship. We prefer this mass estimate in our calculations below because we believe the assumptions behind this photoionisation calculation to be more robust.
45 −1 The Eddington luminosity for this BH mass is LEdd = (2.7±0.1)×10 erg s . We also
43 calculate the intrinsic luminosity at 5100Å (λL5100) for PS1-13cbe that is (1.16±0.01)×10
−1 erg s at the time of the peak and then we estimate the bolometric luminosity using λLλ and a conversion factor of 8.1 (Runnoe et al., 2012) to convert from monochromatic to bolometric luminosity to be (9.4±0.1)×1043 erg s−1. This results in an Eddington parameter of λEdd = Lbol/LEdd ≈ 0.03 at the peak of the outburst.
2.4 Interpretation of the Features of PS1-13cbe
In Figure 2.11, we show the spectrum of PS1-13cbe at the time of the outburst (after subtraction of the star light component from the host galaxy model) alongside the 70 comparison objects, including a QSO (SDSS QSO template), a Type IIn SN (SN1994Y; Filippenko, 1997), and a tidal disruption event (TDE: ASASSN-14li; Holoien et al. 2015). All of these objects show the presence of broad Hα and Hβ lines with emission line profiles similar to the PS1-13cbe spectrum. It is notable that the PS1-13cbe spectrum closely resembles that of the QSO template, which is consistent with AGN activity. However, we estimate the spectral index at the epoch of the LDSS3 spectrum to be α = −0.58, which is redder than the highly variable QSOs studied in Wilhite et al.(2005) that had spectral index of α = −2 for the average difference spectrum (bright phase minus faint phase). In this section, we discuss SNe, TDE, and AGN variability as three possible interpretations for PS1-13cbe. However, before discussing the details about these scenarios we have summarized the key features of PS1-13cbe:
• PS1-13cbe occurred in the nucleus of a Seyfert 2 galaxy with a central SMBH with
7 mass of ∼ 2 × 10 M .
• PS1-13cbe brightened in the course of ∼ 70 days and reached a peak total optical luminosity of (1.06 ± 0.01) × 1043 erg s−1.
• The temperature of PS1-13cbe roughly stayed constant and g − r (transient component) did not show any colour evolution.
• The spectra of PS1-13cbe show significant evolution over the course of ∼ 12 years where broad Hα,Hβ, and Hγ lines appear and disappear from the spectra.
2.4.1 Type IIn Supernovae Interpretation of PS1-13cbe
Type IIn SNe show broad H Balmer lines and blue continua in their optical spectra that can resemble the spectra of AGNs in specific phases of their evolution (Filippenko, 1989). As we show in Figure 2.11, the lack of P Cygni profiles in the spectrum of PS1-13cbe is similar to the spectral features of Type IIn SNe such as SN1994Y (Filippenko, 1997) 71
5 γ Hβ H He II Hα He I 4 ASASSN−14li He I )
−1 β Hα Hγ H Å
−1 3 He I s SN1994Y He I −2
2 Hα β Hγ H erg cm Fe II Fe II He I He I
−16 QSO 1 SDSS template
(10 Hα λ Hβ
Logf γ 0 H Fe II PS1−13cbe Fe II He I 2013/10/05 He I
−1
3000 4000 5000 6000 7000 8000 Rest Frame Wavelength (Å)
Figure 2.11: Optical spectrum of PS1-13cbe taken with LDSS3 during the outburst after subtracting the stellar continuum (blue). The spectra of the comparison objects include a QSO (SDSS QSO; Vanden Berk et al., 2001) template, a SN Type IIn (SN1994Y; Filippenko, 1997) and a TDE (ASASSN-14li; Holoien et al., 2015) from bottom to top (black), the emission lines are labeled (purple).
that is used as an example here for the purpose of the comparison. Additionally, the peak
luminosity of PS1-13cbe, Mr ≈ −19.4, is not substantially brighter than the typical Type
IIn SNe, which have a broad peak luminosity range of −18.5 ≤ Mr ≤ −17 (e.g., Kiewe et al. 2011). The Type-IIn SN hypothesis is however disfavored for the following reasons. First, the narrow emission line ratios in Figure 2.4 clearly identify the galaxy as a Seyfert, which increases the chances that this variability was caused by the existing AGN rather than a SN. Furthermore, the temperature of PS1-13cbe is roughly constant and it does not show 72 any colour evolution, which is not consistent with the cooling of the ejecta typically seen in SNe. Additionally, the double hump behavior seen in the light curves (Figure 2.2 and Figure 2.7) of PS1-13cbe is not commonly observed in Type IIn SNe (e.g., Kiewe et al., 2011; Taddia et al., 2013). While Type II SNe have been very rarely observed to exhibit multiple rebrightening bumps after the first peak (e.g., iPTF14hls and iPTF13z; Arcavi et al., 2017; Nyholm et al., 2017), the small fluctuations observed in the light curve of PS1- 13cbe are not present in their light curves and the rise and rebrightening timescales in the case of PS1-13cbe are much shorter in comparison. Based on the reasons provided in this section, we disfavor the Type IIn SN origin of PS1-13cbe, but this possibility cannot be completely ruled out.
2.4.2 PS1-13cbe as a TDE
Most of the optically detected TDEs show a lack of colour evolution and constant blackbody temperature, consistent with our observations of PS1-13cbe (Gezari et al., 2012; Chornock et al., 2013). Furthermore, broad emission lines such as Hα,Hβ and He II have been detected in the spectra of the TDEs and disappeared at later times (Arcavi et al., 2014; Van Velzen et al., 2011). Additionally, broad Hα and Hβ lines are detected in the spectrum close to the time of the peak and either were not present (Hβ) or only weakly present (Hα) in earlier spectra (SDSS) and disappear in the later spectra (OSMOS). Therefore, it is plausible that PS1-13cbe has a TDE origin. However, we do not have optical photometry data on the decline to check whether the light-curve decays at the predicted rate of t−5/3 for TDEs (Rees, 1988; Evans and Kochanek, 1989). The lack of any X-ray observations at the time of the outburst is another limitation we face to investigate the TDE and SN scenarios. As seen in Figure 2.4, the host of PS1-13cbe is classified as a Seyfert galaxy, which points to the presence of a pre-existing AGN. However, TDEs can happen in AGN galaxies, and in fact it has been suggested that TDEs may prefer galaxies with pre-existing steady 73 accretion to their central BH (Hills, 1975; Blanchard et al., 2017) where the dense star- formation clouds and the pre-existing accretion disk can increase the chance of the tidal encounter with the stars and reduce the relaxation time (Perets et al., 2006; Blanchard et al., 2017). For example, the spectrum of the hosts of the optical/UV TDE SDSS J0748 (Wang et al., 2011) and optical TDE candidate PS16dtm (Blanchard et al., 2017) show weak and narrow emission lines that might be the result of AGN presence in the core. Similar to them, the host of the TDE ASASSN-14li (Holoien et al., 2015) also shows traces of ongoing weak AGN activity at the center (van Velzen et al., 2015; Alexander et al., 2016). Another possibility for brightening and rebrightening of nuclear transients could be repeat tidal stripping of stars (Campana et al., 2015; Ivanov and Chernyakova, 2006; Komossa et al., 2016), although the timescale for rebrightening in PS1-13cbe is very short. However, we disfavor the TDE interpretation of PS1-13cbe because of the following reasons. The blackbody temperatures in optically-selected TDEs typically range from ∼ 20000 K (e.g., PS1-11af; Chornock et al., 2013) to ∼ 35000 K and higher (e.g., ASASSN-14li; Holoien et al., 2015), while the continuum of PS1-13cbe at the time of the outburst is not as blue as a typical TDE and, as we show in Figure 2.8, the inferred blackbody temperature of PS1-13cbe is no more than ∼ 9000 K. Moreover, the He II λ4686 line that is frequently seen in the spectra of TDEs (e.g., ASASSN-14li in Figure 2.11), is not visible in the spectrum of PS1-13cbe that was taken close to the time of outburst. However, it is worth mentioning that there are TDE candidates that lack the presence of the He II λ4686. Additionally, the rebrightening of PS1-13cbe by almost 75% is not typically seen in TDEs. However, it is notable that the TDE candidate ASASSN-15lh showed re- brightening in the UV +60 days after the peak (Dong et al., 2016; Leloudas et al., 2016; Margutti et al., 2017). The optical TDE candidate PS16dtm also showed dimming and rebrightening ∼ 150 days after the beginning of the rise (Blanchard et al., 2017). 74
Furthermore, the light curves of the PS1-13cbe have noticeable fluctuations on few days timescales (such as the one visible in griP1 near −20 days in Figure 2.2) that are inconsistent with the observed smooth optical light curves of most of the other known optically-selected TDEs (e.g., Van Velzen et al., 2011; Gezari et al., 2012; Chornock et al., 2013; Arcavi et al., 2014; Holoien et al., 2014). However, noticeable fluctuations have been observed in the light curves of optical TDE candidate PS16dtm (Blanchard et al., 2017). The optical and spectral features of PS1-13cbe and the reasons provided in this section show that the TDE origin of PS1-13cbe is a possible but not a likely scenario.
2.4.3 PS1-13cbe as a “Changing Look” AGN
Recently, a new type of AGN variability was discovered in objects called “changing look” AGNs that show the appearance or disappearance of broad emission lines followed by an order of magnitude increase or decrease in the continuum and change type from Type 1 to Type 1.8, 1.9 or vice versa (e.g., Runnoe et al., 2016; Shappee et al., 2014; LaMassa et al., 2015; Gezari et al., 2017; MacLeod et al., 2016). This CL behavior observed in some AGNs can be caused by at least three mechanisms. In the first scenario, variation of the obscuration when material such as dust clouds outside of the BLR move in or out of the line of sight that can obscure or clear the view to the BLR (Elitzur, 2012). Another mechanism can be variations of accretion rate that transforms the structure of the BLR. The AGN will transfer from Type 1 when the accretion rate is high and the broad lines are visible to a Type 2 when the accretion rate is low and the broad lines cannot exist or vice versa (Elitzur et al., 2014). Furthermore, it has been suggested that transient events such as TDEs can cause this CL behavior (Eracleous et al., 1995). As shown in Figure 2.9 and 2.10, we only detected the presence of weak broad Hα and no broad Hβ emission line in the spectrum taken by SDSS. However, the strong broad Hα line and Hβ lines appeared at the time of the outburst that followed an observed increase 75 in the flux and disappeared again in the spectra taken at later times. This shows that the AGN changed type from a Type 1.9/2 to a Type 1 and then to a Type 2 because no broad emission lines were detected in later observations.
We estimate the luminosity at 5100 Å (λL5100) based on our estimate that the non- stellar continuum contributes ≤ 10% in the SDSS spectrum to be ≤ 0.15 × 1043 erg s−1 in quiescence. Additionally, we measured the intrinsic 5100 Å (λL5100) luminosity at the
43 −1 peak of the outburst to be 1.16 × 10 erg s , which shows a factor & 8 increase in optical luminosity. Also, during this time the broad Hα emission varied by a factor of ∼ 4 in flux. The observed luminosity changes and appearance and disappearance of the Hα and Hβ lines can be caused by one of the three mentioned mechanisms that are discussed in more detail in the following sections.
2.4.3.1 Obscuration of the AGN
The CL behavior seen in PS1-13cbe can be caused by intervening material which is located outside of the BLR and orbits on a Keplerian orbit that can obscure or give a clear view to the BLR by moving in or out of the line of sight. The bolometric luminosity of ∼ 9.4 × 1043 erg s−1 estimated at the time of outburst using bolometric corrections for unobscured AGNs is less than the bolometric luminosity of ∼ 1.6×1044 erg s−1 at the time of quiescence estimated from the narrow λ5007 and mid-IR excess, which possibly indicates that the AGN is not fully unobscured at the peak of the outburst. In addition, the redder spectrum compared to the QSO template suggests such a scenario to be possible. However, we disfavor the changing obscuration scenario in the case of PS1-13cbe for the following reasons. First, we estimate the characteristic radius of the BLR using the R-L relation calibrated presented by (Bentz et al., 2013) to be RBLR ∼ 11 light days. Then, using the relation for crossing time presented by LaMassa et al.(2015), we estimate the crossing time for an obscuring object orbiting on a Keplerian orbit outside of BLR. In the most ideal case 76 that minimizes the crossing time, the obscuring object should be at rorb ≥ RBLR = 11 light days so that it can intercept a substantial amount of the broad Balmer flux from the BLR.
Even in the extreme case of assuming rorb = RBLR, the crossing time for the obscuring object is ∼ 23 years. In a more realistic scenario where rorb ≥ 3RBLR the crossing time is even higher. Not only are these timescales too long to explain the CL behavior of PS1-13cbe, but also the existence of intervening material with physical properties that can obscure the whole region of the continuum and the BLR at such radii is not obvious. One of the possible candidates for obscuration is the torus that lies just outside of the BLR beyond the dust sublimation radius (e.g., Suganuma et al., 2006; Koshida et al., 2014; LaMassa et al., 2015). However, the sublimation radius is itself 4-5 times larger than the RBLR where the crossing timescale is ∼ 29 years, which is again too long to explain the observed event. Another possible obscuring scenario is that the dimming of the continuum itself promotes the formation of the dust that is able to cover the BLR. However, the timescale for such dust formation with the gas density of ∼ 105 cm−3 in the torus (e.g., Nenkova et al., 2002) is ∼ 103 years (Draine, 2009; LaMassa et al., 2015) that is far too long for this scenario to be true in the case of PS1-13cbe. In addition, if the dust obscuration scenario is true, then the change in luminosity should follow the colour evolution and thus the colour evolution and luminosity change should be correlated. By contrast, there is no colour evolution while the luminosity is changing in the case of PS1-13cbe. We fit a linear regression model to luminosity versus g − r colour index and found a slope of −0.01 ± 0.01 with a coefficient of determination R2 ∼ 0.056. These results confirm the lack of correlation between luminosity and g − r colour index. 77
2.4.3.2 Tidal Disruption Events
As we mentioned in the the Section 2.4.2, the lack of the colour evolution, constant blackbody temperature, and appearance/disappearance of the broad Hα and Hβ lines suggest that the optical outburst and the appearance of the apparent blue continuum may have been powered by a TDE in the pre-existing AGN. Based on the reasons that we provided in the Section 2.4.2, we conclude that TDE origin of the outburst in PS1-13cbe is not favored, but it cannot be ruled out completely.
2.4.3.3 Accretion disk instabilities
In the light curves of PS1-13cbe shown in Figure 2.2, there are small undulations that can be seen in all of the filters (e.g., near −20 days). The amplitudes of these changes are consistent with the variability observed in AGNs (e.g., MacLeod et al., 2012) and point to fluctuations in ongoing accretion activity. It has been shown that in AGNs the optical/UV emission is generated in the accretion disk, with two possible classes of models for the propagation of fluctuations. One is that the locally generated viscous perturbations can produce local blackbody emission (e.g., Krolik, 1999; Liu et al., 2008). In the case of these so-called “outside-in” variations that are produced by changes in accretion rate, the accretion flow fluctuations propagating inward and across the accretion disk first affect the optical region located at outer radii and then affect the UV and X-ray emitting regions. Another origin can be the reprocessing of the UV or X-rays (e.g., Krolik et al., 1991; Cackett et al., 2007; Liu et al., 2008). In this case, the X-rays from the central source irradiate the disk and produce “inside-out” variations from short to long wavelengths (LaMassa et al., 2015; Shappee et al., 2013). The timescale over which the accretion changes happen that might produce “outside- in” variations is known as the inflow timescale. More accurately, changes in the accretion responsible for the “outside-in” variations happen on the inflow time scale, which is the 78 time it takes a parcel of gas in a given radius in the accretion disk to radially move to the center. Assuming the optical continuum emission radius of R ≈ 200 rS (e.g., Morgan et al., 2010; Fausnaugh et al., 2016) and using the relationship between radius and inflow timescale presented by LaMassa et al.(2015), we calculate the inflow timescale to be
6 tinfl ∼ 10 years, which is much longer than the observed change in the continuum flux of PS1-13cbe. However, it is notable that the optical continuum originates in a part of the disk where the radius is an order of magnitude larger than the UV-emitting region, which results in a several thousand times larger inflow timescale than for the UV-emitting region (LaMassa et al., 2015). Thus, we conclude that the rapid continuum flux change in PS1-13cbe is too short to be generated by outside-in variations (perturbations in a given radius of the disk that propagate radially inward), but might be more characteristic of a disturbance in the inner accretion flow that propagates outward. Assuming a standard thin disk model (Shakura and Sunyaev, 1973) with an optically thick and geometrically thin accretion disk, we calculate the dynamical (orbital) time-scale q GM tdyn = 1/Ω of the SMBH where Ω = R3 . We can rewrite the dynamical time-scale to be 3/2 3/2 GM R 2GM tdyn = 2 3 where c is the speed of light and rS = 2 is the Schwarzschild radius c rS c 7 of the central BH. The dynamical timescale around a SMBH with mass of 2.2 × 10 M 3/2 R can be written as tdyn ≈ 310 s. Assuming an optical emission distance of R ≈ 200 rS rS (e.g., Morgan et al., 2010; Fausnaugh et al., 2016), we calculate the dynamical timescales
−1 of ≈ 10 days. Then, using tth = α tdyn we calculate the thermal timescale to be ≈ 99 days. We note that these values are not strongly dependent on the uncertainties in the BH mass for this object, so the ordering of timescales is robust. Among all of the calculated timescales, only the thermal timescale of ∼ 99 days is reasonably similar to the observed timescale of ∼ 70 days in the case of PS1-13cbe and suggests another possible scenario where reprocessing of the UV or X-rays can produce 79 the optical variations. In this case, an increase in the X-rays that originate in the smaller hot corona that is closer to the central BH can heat the inner part of the accretion disk first, move outward, and generate inside-out variations by irradiating the disk and driving an increase in the blue and then red emission. This scenario has been observed and well studied in NGC 2617 (Shappee et al., 2013; Oknyansky et al., 2017). Shappee et al. (2013) detected the variability first in X-rays and then with time lags in UV and NIR and concluded that the observed continuum flux variability resulted from inside-out variations. In another example, NGC 4151, it was also observed that X-ray variability drove UV- optical variability (Edelson et al., 2017).
2.4.4 Comparison to other Changing Look AGNs
PS1-13cbe is one of the few CL AGNs that has been observed during the turn-on phase (Cohen et al., 1986; Storchi-Bergmann et al., 1993; Aretxaga et al., 1999; Eracleous and Halpern, 2001; Shappee et al., 2014; Gezari et al., 2017) by suddenly demonstrating the appearance of broad Hα and Hβ emission lines. The CL behavior has been observed in other candidates where the broad lines appeared or disappeared in spectra taken a couple of years to decades apart (e.g., LaMassa et al., 2015; Runnoe et al., 2016; MacLeod et al., 2018). By contrast, in PS1-13cbe the outburst timescale is very short. PS1-13cbe “turned- on” in only 70 days and the broad lines were observed in a spectrum taken 57 days after the peak of the outburst. Furthermore, the broad lines disappeared before two years after the time of the peak, which is again very short, and they never re-appeared in later spectra (see Figure 2.3). However, we should note that by “turn-on” we refer to brightening of a pre-existing AGN, but the presence of the strong narrow emission lines with a Seyfert 2 spectrum implies that narrow-line region was exposed to AGN accretion activity in the recent past. 80
Gezari et al.(2017) also presented the CL quasar iPTF 16bco that had a “turn-on”
8 timescale of ≤ 1 year which is very short for a BH with MBH ∼ 10 M compared to others. PS1-13cbe and iPTF 16bco are the only CLs that demonstrate extremely short turn- on timescales which push the limit of the accretion disk theory. One more interesting fact about PS1-13cbe is that the light curves showed a double peak behavior observed in all of the optical bands (see Figure 2.2). This behavior was also observed in the X-rays/UV in the case of NGC 2617 (Shappee et al., 2014), in the UV/optical in the case of the ASASSN- 15lh TDE candidate (Dong et al., 2016; Leloudas et al., 2016; Margutti et al., 2017), and in the optical band in the case of the PS16dtm TDE candidate (Blanchard et al., 2017).
2.5 Conclusions
We present a transient event that was discovered in the PS1/MDS survey, PS1-13cbe, at redshift z = 0.12355. The outburst happened in the nucleus of a galaxy that is classified as a Seyfert 2 (see Figure 2.4) using the SDSS archival data that was taken a decade before the outburst. At the time of the outburst, the galaxy changed type to a Seyfert 1 as broad Hα and Hβ appeared and the continuum brightened in the spectrum taken with LDSS3 +57 days after the peak and then changed its type back to a Seyfert 2 as the broad Hα and Hβ disappeared in the spectrum taken with OSMOS 2 years later and did not reappear in spectra taken 3 and 4 years after the outburst. The optical photometry shows that the continuum flux increased by a factor of ∼8 on a timescale of ∼ 70 days and declined for next ∼ 50 days and then rose again over the course of the next ∼ 50 days. Observational evidence presented in this work argues against the Type IIn SN and TDE interpretations. The constant colour evolution and blackbody temperature during the outburst and also the presence of a pre-existing AGN disfavour the SN Type IIn scenario. As mentioned, TDEs have been observed in the galaxies with pre-existing AGNs (e.g., SDSS J0748; Wang et al., 2011); however, the lack of a broad He IIλ4686 emission line, 81 low blackbody temperature at the time of peak, a light curve that has small fluctuations, and unusual re-brightening by 75% are inconsistent with properties of known TDEs. We conclude that PS1-13cbe is a changing-look AGN that has been powered by instabilities in the accretion disk. We argued against the obscuration scenario and TDE origin of these accretion disk instabilities by showing that the crossing and viscous timescales are longer than the timescale observed in the case of PS1-13cbe. Furthermore, we also argued against outside-in variations by calculating the inflow timescale which is too long to explain the observed timescale here. We also conclude that the thermal instabilities in the accretion disk are most likely the source of the outburst and CL behavior. These thermal instabilities may have caused inside-out variations that have generated the observed optical variability in the light curves of PS1-13cbe. One very interesting point about PS1-13cbe is that the observed turn-on timescale pushes the limits of viscous accretion disk theory which predicts much longer timescales and it might be one of the CLs that have shown the most rapid change of the state (iPTF 16bco; Gezari et al., 2017) compared to the other CLs. Other CLs have been observed over a timespan of years and sometimes decades apart; however, the short timescale observed here suggests that other candidates may have been through these short timescale outbursts. Therefore, more frequent observations with a higher cadence and multiwavelength coverage to overcome limitations such as the one we faced here with the lack of X-ray observations at the time of outburst, will provide us more insight to better understand the changing look behavior of AGNs. 82
Table 2.1: Photometry of PS1-13cbe.
Data of Observation Epoch Filter Mag Mag (MJD) (days) (observed) Uncertainty
56217.3 -295.3 gp1 (23.5)
56475.5 -37.1 gp1 19.89 0.03
56480.5 -32.1 gp1 19.79 0.03
56486.6 -26.0 gp1 19.76 0.03
56489.5 -23.1 gp1 19.91 0.03
56501.4 -11.2 gp1 19.79 0.03
56507.6 -5.0 gp1 19.76 0.03
56512.5 -0.1 gp1 19.66 0.03
56518.5 5.9 gp1 19.85 0.03
56531.5 18.9 gp1 20.16 0.04
56534.5 21.9 gp1 20.20 0.04
56539.4 26.8 gp1 20.30 0.04
56545.3 32.7 gp1 20.38 0.04
56548.3 35.7 gp1 20.42 0.06
56563.3 50.7 gp1 20.49 0.05
56571.2 58.6 gp1 20.50 0.04
56574.2 61.6 gp1 20.43 0.06
56590.3 77.7 gp1 20.35 0.04
56593.3 80.7 gp1 20.24 0.04
56596.2 83.6 gp1 20.22 0.04
56601.2 88.6 gp1 20.34 0.07
56628.3 115.7 gp1 20.11 0.17
56285.2 -227.4 rp1 (23.5) 83
Table 2.1: continued
Data of Observation Epoch Filter Mag Mag (MJD) (days) (observed) Uncertainty
56475.6 -37.0 rp1 19.85 0.04
56479.5 -33.1 rp1 19.75 0.04
56486.6 -26.0 rp1 19.66 0.03
56489.5 -23.1 rp1 19.80 0.03
56501.4 -11.2 rp1 19.70 0.03
56507.6 -5.0 rp1 19.64 0.03
56508.5 -4.1 rp1 19.63 0.03
56512.6 0.0 rp1 19.60 0.03
56518.5 5.9 rp1 19.72 0.03
56531.5 18.9 rp1 20.04 0.04
56534.5 21.9 rp1 20.08 0.04
56538.3 25.7 rp1 20.03 0.05
56545.3 32.7 rp1 20.11 0.04
56571.2 58.6 rp1 20.44 0.06
56574.3 61.7 rp1 19.81 0.37
56590.3 77.7 rp1 20.29 0.05
56593.3 80.7 rp1 20.11 0.05
56596.3 83.7 rp1 20.05 0.04
56597.3 84.7 rp1 20.08 0.04
56626.2 113.6 rp1 20.11 0.07
56284.2 -228.4 ip1 (23.5)
56473.5 -39.1 ip1 19.65 0.04
56476.6 -36.0 ip1 19.60 0.04 84
Table 2.1: continued
Data of Observation Epoch Filter Mag Mag (MJD) (days) (observed) Uncertainty
56481.6 -31.0 ip1 19.51 0.05
56484.5 -28.1 ip1 19.42 0.03
56491.5 -21.1 ip1 19.54 0.03
56505.4 -7.2 ip1 19.46 0.03
56508.6 -4.0 ip1 19.44 0.03
56513.6 1.0 ip1 19.44 0.03
56516.5 3.9 ip1 19.46 0.03
56520.5 7.9 ip1 19.54 0.03
56532.4 19.8 ip1 19.79 0.04
56538.4 25.8 ip1 19.72 0.04
56540.4 27.8 ip1 19.70 0.04
56550.5 37.9 ip1 20.18 0.06
56558.4 45.8 ip1 20.07 0.05
56564.4 51.8 ip1 19.96 0.05
56567.2 54.6 ip1 20.10 0.06
56569.3 56.7 ip1 20.02 0.05
56575.3 62.7 ip1 20.01 0.05
56588.4 75.8 ip1 19.79 0.04
56597.2 84.6 ip1 19.88 0.04
56599.3 86.7 ip1 19.89 0.04
56626.2 113.6 ip1 19.85 0.06
56216.2 -296.4 zp1 (23.5)
56453.5 -59.1 zp1 19.86 0.05 85
Table 2.1: continued
Data of Observation Epoch Filter Mag Mag (MJD) (days) (observed) Uncertainty
56456.6 -56.0 zp1 19.71 0.04
56474.6 -38.0 zp1 19.58 0.05
56477.5 -35.1 zp1 19.67 0.04
56480.5 -32.1 zp1 19.44 0.04
56482.5 -30.1 zp1 19.41 0.03
56485.6 -27.0 zp1 19.54 0.04
56488.4 -24.2 zp1 19.55 0.03
56490.5 -22.1 zp1 19.48 0.05
56492.5 -20.1 zp1 19.54 0.04
56506.6 -6.0 zp1 19.53 0.04
56511.6 -1.0 zp1 19.47 0.04
56517.4 4.8 zp1 19.49 0.04
56519.6 7.0 zp1 19.60 0.04
56521.6 9.0 zp1 19.65 0.04
56530.4 17.8 zp1 19.81 0.04
56536.5 23.9 zp1 19.84 0.05
56539.4 26.8 zp1 19.91 0.05
56541.3 28.7 zp1 19.79 0.05
56544.3 31.7 zp1 19.87 0.06
56549.3 36.7 zp1 19.90 0.06
56557.5 44.9 zp1 20.07 0.06
56568.3 55.7 zp1 20.35 0.14
56570.3 57.7 zp1 20.03 0.06 86
Table 2.1: continued
Data of Observation Epoch Filter Mag Mag (MJD) (days) (observed) Uncertainty
56584.3 71.7 zp1 19.92 0.05
56589.3 76.7 zp1 19.94 0.05
56592.3 79.7 zp1 19.87 0.05
56595.2 82.6 zp1 19.84 0.05
56600.3 87.7 zp1 19.81 0.05
56613.3 100.7 zp1 19.87 0.05
56628.2 115.6 zp1 20.03 0.06
56266.3 -246.0 yp1 (23.5)
56434.6 -78.0 yp1 20.63 0.17
56436.6 -76.0 yp1 20.42 0.10
56437.6 -75.0 yp1 20.73 0.29
56454.6 -58.0 yp1 19.70 0.11
56467.5 -45.0 yp1 19.82 0.04
56484.5 -28.0 yp1 19.34 0.05
56493.6 -19.0 yp1 19.50 0.03
56494.6 -18.0 yp1 19.58 0.03
56511.5 -1.0 yp1 19.37 0.03
56517.4 5.0 yp1 19.32 0.06
56522.6 10.0 yp1 19.66 0.04
56532.4 20.0 yp1 19.66 0.04
56556.4 44.0 yp1 19.73 0.03
56559.3 47.0 yp1 19.94 0.09
56562.3 50.0 yp1 19.88 0.11 87
Table 2.1: continued
Data of Observation Epoch Filter Mag Mag (MJD) (days) (observed) Uncertainty
56563.2 51.0 yp1 19.89 0.04
56566.5 54.0 yp1 20.04 0.17
56567.3 55.0 yp1 19.77 0.14
56569.2 57.0 yp1 20.03 0.19
56573.3 61.0 yp1 20.00 0.08
56574.2 62.0 yp1 19.89 0.19
56584.4 72.0 yp1 19.72 0.04
56585.3 73.0 yp1 19.96 0.07
56586.3 74.0 yp1 20.15 0.05
56613.3 101.0 yp1 19.79 0.04
56616.2 104.0 yp1 19.88 0.08
56620.2 108.0 yp1 19.80 0.04
56626.2 114.0 yp1 19.70 0.17
56638.2 126.0 yp1 19.77 0.04
Note: Magnitudes provided here are not reddening corrected and the epochs are relative to the peak in the rest frame. 3σ upper limit values are represented in parentheses. 88
Table 2.2: Luminosity of Broad Lines
Date (UT) Instrument Broad Hα Broad Hβ
2003 June 26 SDSS 4.87 ± 0.37 1.25 ± 0.39 2013 Oct 05 LDSS3 16 ± 0.28 4.39 ± 0.22 2015 Oct 03 OSMOS 3.92 ± 0.45 < 1.2 2016 Nov 16 OSMOS – < 0.85 2017 June 18 OSMOS < 0.86 < 1.14
Note: Luminosity of broad Hα and Hβ lines in SDSS, LDSS3, and OSMOS spectra reported in units of 1040 erg s−1. Non-detections are reported as 3σ upper limits. 89 3 PS1-10cdq
In this chapter, we discuss a transient that was discovered in PS1 survey in 2010. We discuss TDE, SLSN, and AGN activity as the possible sources of this nuclear outburst.
3.1 Observations of PS1-10cdq
On 2010 December 10, we discovered a transient event consistent with the nucleus of a galaxy, PS1-10cdq, at coordinates α = 10h02m29.92s, δ = +01◦30006.100 (J2000) using the Photpipe transient discovery pipeline, described by Rest et al.(2014) and Scolnic et al.(2018). The host galaxy was observed by SDSS and was given the name SDSS J100229.92+013006.1 (hereafter SDSS J1002+013) with a redshift of z = 0.373. The difference image light curve of PS1-10cdq is constructed by the PS1 transient pipeline. A stack of high quality images (excluding the outburst observing season; year 2010) was used to create template images that were subtracted from all of the observations of the transient PS1-10cdq. This flux difference is relative to the host galaxy contribution in the templates.
3.1.1 Optical Photometry
As we show in Figure 3.1, the light curves of PS1-10cdq were constant with zero change relative to baseline in one season before the outburst. Then, the light curves showed a peak at MJD 55545.5 and declined. The luminosity in gp1 band (νLν) increased from a base luminosity of 0.25 × 1044 erg s−1 to a peak luminosity of 2.01 × 1044 erg s−1 in the
44 course of ∼ 200 days (gp1 = 19.50 at peak; see Figure 3.2) and then declined to 0.35 × 10 erg s−1 in the course of ∼ 135 days. The magnitude of these variations are ∼ 100 times larger than the AGN variability of RMS ∼ 0.02 × 1044 erg s−1 observed in other seasons excluding the one that contains the outburst. Also, the transient was observed by Galex in the Near Ultra Violet (NUV) band shown by olive bars in Figure 3.1. 90
2.0 y z i r
) g
−1 1.5 NUV erg s 44 1.0 (10 ν L ν 0.5
0.0 S −1000 −500 0 500 1000 1500 Time (MJD − 55545.5)
Figure 3.1: The observed Galactic extinction corrected luminosities of PS1-10cdq from the
PS1 survey in grizyP1 filters. S: marks the epoch of the LDSS spectrum (MJD 55626).
3.1.2 Observations of the Host Galaxy
In addition to the PS1 observations of the host in gp1, rp1, ip1, zp1, and yp1 with values of 21.61 ± 0.03, 20.85 ± 0.01, 20.67 ± 0.03, 20.22 ± 0.02, and 19.21 ± 0.22 magnitudes, we found archival observations from XMM, GALEX, SDSS, IRAC, and MIPS that are summarized in Table 3.1.
3.1.3 X-ray Photometry
As we show in Table 3.1, the host galaxy of the PS1-10cdq has been observed with XMM Newton in Cosmic Evolution Survey (COSMOS) survey with total exposure time of 14.8 ks in November 24 2008 (∼ 2 years before the outburst; Cappelluti et al., 2009b).
−14 −2 −1 The host galaxy was detected with fx(0.5 − 2 keV) = (0.46 ± 0.07) × 10 erg cm s . 91
16
18
20 z − 4 i − 3
Magnitude + constant r − 2 22 g − 1 NUV
0 20 40 60 80 100 120 140 Time (MJD − 55545.5)
Figure 3.2: Galactic extinction corrected light curves of PS1-10cdq in grizyP1 and Galex NUV bands.
20 −2 Using galactic neutral hydrogen column density of nH ≈ 2.68 × 10 cm (Kalberla et al., 2005) in the direction of PS1-10cdq, and a photon index of Γ = 2, we calculated
−14 −2 −1 the X-ray flux to be fx(2 − 10 keV) = (0.58 ± 0.09) × 10 erg cm s which gives
42 −1 Lx(2 − 10keV) = (2.83 ± 0.44) × 10 erg s . Additionally, using the empirical bolometric
correction from Marconi et al.(2004) and LX, we calculated a bolometric luminosity of
43 −1 Lbol ≈ (3.71 ± 0.41) × 10 erg s .
3.1.4 Optical Spectroscopy
We observed PS1-10cdq using the Inamori Magellan Areal Camera and Spectrograph (IMACS) on the 6.5 m Magellan Baade (MC) telescope for 5400 s on 2011 March 9 (∼ +84 days after the peak of the outburst). We used a 0.900 long slit with the 200-15.0 grism to 92
Table 3.1: Photometry of the host of PS1-10cdq.
Observed Passband Photometry Uncertainty Units No. Measurement
1 2-10 keV (XMM) – < 1.75E-14 erg cm−2 s−1 2 0.5-2 keV (XMM) 0.46E-14 ±0.07E-14 erg cm−2 s−1 3 FUV (GALEX) AB 23.03 ±0.22 mag 4 NUV (GALEX) AB 22.32 ±0.13 mag 5 u (SDSS Model) AB 21.80 ±0.25 asinh mag 6 g (SDSS Model) AB 21.10 ±0.05 asinh mag 7 r (SDSS Model) AB 20.18 ±0.037 asinh mag 8 i (SDSS Model) AB 20.01 ±0.05 asinh mag 9 z (SDSS Model) AB 19.55 ±0.12 asinh mag 10 3.6 microns (IRAC) 69.79 ±0.30 microJy 11 4.5 microns (IRAC) 73.41 ±0.41 microJy 12 5.8 microns (IRAC) 74.67 ±1.44 microJy 13 8.0 microns (IRAC) 136.10 ±2.55 microJy 14 24 microns (MIPS) 522.88 ±11.78 microJy
cover wavelength range of 3900 − 10150 Å with a resolution of ∼ 8Å. We have many late-time observations that we summarize in Table 3.2 and show in Figure 3.3.
3.2 Observational Features of PS1-10cdq
3.2.1 Host Galaxy of PS1-10cdq
We generated a quiescent spectrum by combining three post-outburst spectra of SDSS J100229.92+013006.1 (spectra taken with IMACS on January, March, and May of 2017; 93
Table 3.2: Late time spectra of the host of PS1-10cdq (Walter Baade (WB), and center (CN)).
Date Observed INST EXP Telescope Grism Slit Wavelength Resolution Seconds Range (Å) (Å)
2013 Dec 30 LDSS3 1800 6.5 m MC VPH-all 100 CN 3580 − 10700 ∼ 9 2014 Dec 19 IMACS 1800 6.5 m WB 300-17.5 0.900 top 3700 − 9450 ∼ 5 2015 Feb 18 IMACS 3000 6.5 m WB 300-17.5 0.900 3700 − 9450 ∼ 5 2015 Apr 23 IMACS 3000 6.5 m WB 300-17.5 0.900 3700 − 9450 ∼ 5 2017 Jan 21 LDDS3-C 900 6.5 m MC VPH-all 100 CN 3580 − 10700 ∼ 9 2017 May 1 LDSS3-C 2400 6.5 m MC VPH-all 100 CN 3580 − 10700 ∼ 9 2018 Mar 14 LDSS3-C 3600 6.5 m MC VPH-all 100 CN 3580 − 10700 ∼ 9
see Table 3.2). In order to isolate emission lines and subtract stellar absorption lines, we used this combined spectrum, optical photometry in SDSS griz bands to simulate a model of the host galaxy using the FAST 1.0 code (Kriek et al., 2009). We experimented with different initial parameters and were able to generate the best fit using a stellar age of ∼ 109 years and e-folding timescale of τ ≈ 109 years, by assuming the star-formation history to be exponentially declining, the stellar initial mass function from Chabrier(2003), the Bruzual and Charlot(2003) spectral library, a Milky Way dust law (Cardelli et al., 1989), and solar- like metallicity of Z = 0.02. The model of the host galaxy simulated by FAST is plotted alongside the combined spectrum of the host in the Figure 3.4. As we see in Table 3.1, the host galaxy of PS1-10cdq was detected in X-rays that points to presence of an AGN. Additionally, as we show in Figure 3.5, we constructed the SED of the host galaxy using observed values from PS1 and archival photometry data (see Table 3.1 for more details). 94
8 Mar 2018
6 May 2017
Jan 2017
Apr 2015 ) + constant 4 −1 Feb 2015 Å −1
s 2 Dec 2014 −2 Dec 2014
0 Dec 2013 (erg cm λ Dec 2011
Logf −2 Mar 2011 Hβ α Mg II H 3000 4000 5000 6000 7000 8000 Rest frame wavelength (Å)
Figure 3.3: Optical spectra of PS1-10cdq. Spectrum ∼ 84 days after the peak of outburst observed using LDSS3 in dark red. Spectrum taken ∼ +352 days after the peak of outburst in purple. For detailed information about late time spectra shown in navy blue see Table 3.2. All the spectra are binned for the purpose of better visualization.
The MIR, UV, and NUV excess of the host relative to the star-forming template observed in the SED is another evidence of AGN activity.
3.2.2 Multi-band Light Curves of PS1-10cdq
We calculate the total optical luminosity of the PS1-10cdq using multi-band observations (see Figure 3.1). In order to do this, first, we find the epochs where gp1 and rp1 bands were observed simultaneously. Next, using Legendre polynomials, we interpolated for ip1 and zp1 bands at these epochs. Finally, we calculate the total flux which translates to total optical luminosity over 3685 − 8910 Å in the rest frame, by integrating 95
Combined Host Spectrum