Summer Research Fellowship Programme of India's Science Academies

Spectral analysis of the Narrow Line Seyfert 1 galaxy PG 1244+026

Prerna Rao

Msc Physics, Manipal Institute of Technology, Eshwar Nagar, Manipal, Karnataka 576104

Guided by:

Dr. C S Stalin

Indian Institute of Astrophysics, Sarjapur Main Road, 2nd Block, Koramangala, Bengaluru, Karnataka, 560034


Active galactic nuclei (AGN) are one of the brightest objects in the sky and they are believed to be powered by the accretion of matter onto super massive black holes situated at the centres of galaxies. About three decades ago Osterbrock and his group identified a new group of active galactic nuclei called Narrow Line Seyfert 1 (NLSy1) galaxies. They have narrow emission lines in their optical spectra compared to broad line Seyfert 1 galaxies. In addition to that, they have excess emission in their soft X-ray spectra and they are also known to show strong X-ray variability. In this project, we aim to understand the peculiar X-ray properties of Narrow Line Seyfert 1 galaxies, using data from the XMM-Newton satellite which was launched by the European Space Agency (ESA) in 1999. For the purpose of this project, I have selected a source named PG 1244+026. We have reduced the data and extracted the lightcurves and spectra corresponding to both source and background regions. We have also done spectral fitting using the clean spectra and tried to characterize the X-ray variability and spectral properties of PG 1244+026.

Keywords: Active galaxies, Accretion, X-rays.


AGNActive Galactic Nuclei
NLSy1Narrow Line Seyfert1
XMM X-ray Multi Mirror 
EPIC European Photon Imaging Camera 
 MOSMetal Oxide Semiconductor 
CCD Charge Coupled Device 
HEASOFT  High Energy Astrophysics Software
 HEASARCHigh Energy Astrophyscs Archive Research Center 
 SASScience Analysis System 
 ODF Observation Data Files
 RMF Redistribution Matrix File
 ARF Ancillary Response File
PPS Processing Pipeline Files 
 CIFCalibration Index File 


1.1 Background/Rationale

Active galaxies are galaxies with an extremely luminous central region emitting radiations over the entire range of the electromagnetic spectrum. The central region is often more luminous than the total luminosity of the rest of the galaxy. Active galactic nuclei are powered by central supermassive black holes (BH) whose mass lies in the range 1016 -1010 MM_\odot. The superamassive black hole is surrounded by a gaseous accretion disc which is a few light days across. The broad line region of AGN corresponds to fast moving gas clouds at a distance of ~100 light days from the centre. The narrow line region corresponds to small, low density gas clouds moving at lower velocities at a distance of ~1000 lightyears from the centre. In many AGNs emission is collimated into jets whose direction of propogation is perpendicular to the accretion disc plane. AGNs are classified into several categories, one of which is based on radio emissions. The radio loud AGNs produce radio jets and lobes and a significantly large part of the total emission is due to the jets. Radio loud AGNs are typically found in elliptical host galaxies whereas radio quiets are associated with spiral hosts.

Radio quiet galaxies consist of Seyfert galaxies, Quasars and LINERs. Carl Seyfert discovered the first class of AGN that are now named after him. Seyfert galaxies are further classified into type 1 and type 2 based on the presence or absence of broad, permitted emission lines. Seyfert 1 spectra consists of both broad and narrow emission lines, while Seyfert 2 spectra consists of only narrow emission lines.

Quasars are the most luminous among all the subclasses of AGNs. Their spectra is quite similar to those of Seyferts, except that quasar spectra have weak or absent stellar absorption features and their narrow lines are weaker than the broad lines of Seyferts. Low-ionization Nuclear Emission-line Region galaxies (LINERs) are spectroscopically similar to Seyfert 2 galaxies except that they show relatively strong low ionization lines.

    : Unified model of AGN

    About 15-20% of all AGNs have been observed to be radio loud. Radio loud AGNs are further classified into Blazars and Radio galaxies. The appearance of active galactic nuclei depends so strongly on orientation that our current classification schemes are dominated by random pointing directions instead of the more interesting physical properties (Urry and Padovany 1995). ​​

    For this project we have studied the X-ray spectral and variability properties of the narrow line Seyfert 1 galaxy PG 1244+026 using the data collected by the X-ray space observatory XMM Newton. XMM Newton short for X-ray Multi mirror mission was built by the European Space Agency and was launched on December 10, 1999. The European Photo Imaging camera (EPIC) carries two MOS-CCD cameras used to detect low energy X-rays and one PN-CCD camera used to detect high energy X-rays.

    1.2 Statement of the Problems

    We selected the source PG 1244+026 which is a narrow line Seyfert 1 galaxy to understand the origin of the steeper soft X-ray continuum slope in the spectrum.

    This source has been previously selected as one of the brightest unobscured type 1 AGN with both XMM–Newton and SDSS spectra (Jin et al 2012).

    From an analysis of optical spectra, PG 1244+026 can also be seen to have the narrowest Hβ line (830 km s −1 ) and the fourth smallest Hβ equivalent width (41 Å) of all of the PG quasars (Boroson & Green 1992).

    An observation corresponding to 120 ks from data obtained by XMM-Newton has been studied and the spectrum and variability of the source corresponding to the observation have been discussed ( ​Jin et al 2013​).

    In this project we have tried to understand the spectral and variability properties corresponding to all the eight observations as recorded by XMM-Newton.

    1.3 Objectives of the Research

    1. To characterize the X-ray variability properties of the narrow line Seyfert 1 galaxy PG 1244+026.

    2. To characterize the X-ray spectral properties of the NLSy1 PG 1244+026.


    Narrow Line Seyfert 1 (NLSy1) galaxies were classified as a special class of active galactic nuclei by Osterbrock and collaborators about 30 years ago. NLSy1 spectra show narrow emission lines and generally steeper soft X-ray continuum slopes. The full width half maximum of Hβ line is less than 2200. For NLSy1 galaxies the balmer lines of hydrogen are only slightly broader than the forbidden lines, emission lines of Fe II or higher ionisation iron emission lines are often seen in the spectra and the ratio of O3 5007A to Hβ is less than 3.

    NLSy1 are found in spiral host galaxies powered by low mass black holes. They were originally believed to be radio-quiet but it is now known that about 7% of all observed NLSy1 are radio loud.

    Interest in the study of NLSy1 galaxies has been renewed in the past few years owing to the detection of gamma rays in about 6 NLSy1 galaxies as studied from observations collected by the Fermi gamma-ray space telescope. Detailed studies of these observations have produced enough evidence to support the claim that the gamma ray emitting NLSy1 galaxies are similar to blazars in their general properties.

    This has led to the development of new arenas in the study of NLSy1 gaaxies.

    Through the study of data collected by ROSAT, it has been found that NLSy1 galaxies have remarkable X-ray properties. (Boller, Brandt and Fink 1996).

    The X-ray spectra of NLSy1 galaxies are often much softer than those of Seyfert 1 galaxies with broader lines. NLSy1 galaxies can exhibit rapid and giant amplitude variability. (​Boller 1997​)



    We selected the narrow line Seyfert 1 galaxy PG1244+026 as our sample and the data was obtained from NASA’s High Energy Astrophysics Archive Research Center (HEASARC). The data consisted of eight observation ids corresponding to different exposure durations and we have done the data reduction for seven of the eight observation ids. The two softwares I used for the purpose are HEASOFT and SAS. HEASOFT 6.24 provided by NASA’S HEASARC is an unified release of the FTOOLS and XANDAU software packages which are used for X-ray astronomical spectral, timing, imaging analysis and manipulaton of FITS files. The Science Analysis System software version 17.0.0 was used for basic XMM-Newton data reduction. Data from XMM consists of ODF and PPS. The Observation Data Files (ODF) contain raw data and need to be reprocessed. The Pipeline Processing System(PPS) contains pre-processed data. Before creating event files correct environment is set for analysis, by creating cif and by defining SAS_CCFPATH, SAS_CCF and SAS_ODF. The data from all three cameras (PN, MOS1 and MOS2) have been used reduced, although only the data from PN camera was used for analysis. The source and background regions were extracted from the image. The source extraction region was chosen to be circular with a radius of 20 arcsec for each EPIC camera. The background region was chosen to be circular with a radius of 25 arcsec in the case of PN and 20 arsec for MOS as far as possible from the source. The corresponding lightcurves were obtained and in the case of high flares in the background lightcurve, rate filters were applied. We selected data with PATTERN ≤ 12 for MOS1 and MOS2, and PATTERN ≤ 4 for PN. Spectra were extracted from source and background regions separately. The areas of source and background were calcualated using backscale. The detector response was reformatted by making a redistribution matrix file (RMF) and ancillary response file (ARF) using RMFGEN and ARFGEN. The source and background spectra were then grouped together using SPECGROUP. The clean spectra generated was used for spectral fitting.

    In the case of 0675320201, data reduction could not be performed as we were only interested in the analysis of PN data and the PN data files were missing.

    : The different observations corresponding to the source PG1244+026 as obtained by XMM-Newton
    0744440401 26-12-2014 128500
    0675320101 25-12-2011 123869
    0744440301 24-12-2014 122002
    0744440501 15-01-2015 120100
    0744440201 22-12-2014 120000
    0744440101 20-12-2014 119499
    0051760101 17-06-2001 13061
    0675320201 25-12-2011 7829

    3.2 ANALYSIS


    Lightcurves display the variation of the brightness of a source with time. In order to calculate Fvariance, lightcurves were plotted in three energy bands: 0.3-1.0 kev, 1.0-2.0 kev and 2.0-10.0 kev. The source extraction region was again selected to be circular with a radius of 20 arcsec and the background extraction region was a circle of radius 25 arcsec . LCMATH in FTOOLS was used to subtract the background lightcurve from that of the source. The corrected lightcurve was used to calculate the variability amplitude. We have only used the data from PN camera for analysis.


    To characterize the flux variations we used the fractional root mean square (rms) variability amplitude Fvar . This gives an estimate of the intrinsic variability amplitude relative to the mean count rate exceeding the measurement noise in the lightcurves. Fvar is defined as:

    Fvar  =  S2σerr2x2\displaystyle F_{var\;=\;}\sqrt{\frac{S^2-\overline{\sigma_{err}^2}}{\overline x^2}}

    where S2 is the sample variance, x\overline x the arithmetic mean of xi and σ2err\overline{{\sigma^2}_{err}}the mean square error, given by:

    S2=1N1i1N(xix)2\displaystyle S^2=\frac1{N-1}\sum_{i-1}^N\left(x_i-\overline x\right)^2
    σerr2=  1Ni=1Nσerr,i2\displaystyle \overline{\sigma_{err}^2}=\;\frac1N\sum_{i=1}^N\sigma_{err,i}^2

    The uncertainty in Fvar is given by:

    err(Fvar)=(12Nσerr2x2Fvar)+(σerr2N1x)\displaystyle err\left(F_{var}\right)=\sqrt{\left(\sqrt{\frac1{2N}}\frac{\overline{\sigma_{err}^2}}{\overline x^2F_{var}}\right)+\left(\sqrt{\frac{\overline{\sigma_{err}^2}}N}\frac1{\overline x}\right)}

    (​Rani P, Stalin and Rakshit 2017​)


    Spectral fitting for the clean spectra was done on xspec v.12.10.0. We have tried the following combinations of models for fitting:

    1) wabs✴zwabs✴pow





    Wabs uses Wisconsin cross-sections and accounts for photoelectric absorption. Zwabs permits us to specify the redshift of the source. For our source, the redshift was found to be 4.81780x10-2 and was frozen at that value. Pow is a simple power law model and bbody is a blackbody spectrum model. Nthcomp gives the shape of the continuum from thermal comptonisation. It accounts for both high and low energy rollovers and is also parametrized by the temperature of electron(kTe) and seed photons(kTbb). Comptt is an analytical model developed by Titarchuk. It describes the comptonization of soft photons in a hot plasma and also includes relativistic effects. (Titarchuk 1994)


    4.1 Lightcurves


                  4.2 Fvariance

                  Fvariance calculated using corrected lightcurves
                  OBSERVATION ID Fvar ±err(Fvar)
                  0.3-1.0 kev 1.0-2.0 kev 2.0-10.0 kev
                  0744440401 0.237±0.002 0.752±0.005 0.965±0.010
                  0675320101 0.212±0.001 0.199±0.003 0.190±0.006
                  0744440301 0.231±0.002 0.227±0.004 0.280±0.009
                  0744440501 0.172±0.002 0.182±0.004 0.701±0.009
                  0744440201 0.173±0.002 0.192±0.004 0.468±0.008
                  0744440101 0.192±0.002 0.206±0.004 0.695±0.009
                  0051760101 0.136±0.006 0.171±0.013 0.194±0.025

                  4.3 Spectra













                                                                              4.4 Details of spectral fitting

                                                                              MODEL COMPONENT PARAMETER VALUE±ERROR REDUCED χ2
                                                                              1)wabs*zwabs*pow wabs nH*(1022) 1.72000x10-2 14.7030
                                                                              zwabs nH(1022) 2.28682x10-7±1.72939x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.01895±1.03787x10-2
                                                                              norm 2.58182x10-3±1.49937x10-5
                                                                              2)wabs*zwabs*(pow+bbody) wabs nH*(1022) 1.72000x10-2 3.7588
                                                                              zwabs nH(1022) 6.78123x10-2±3.86284x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.63088±2.95110x10-2
                                                                              norm 2.95940x10-3±2.90458x10-5
                                                                              bbody kT(keV) 1.31393±2.72855x10-2
                                                                              norm 1.77245x10-5±3.98979x10-7
                                                                              3)wabs*zwabs*(bbody+nthcomp) wabs nH*(1022) 1.72000x10-2 3.2144
                                                                              zwabs nH(1022) 4.57195x10-10±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 0.132134±8.84000x10-4
                                                                              norm 6.65534x10-5±1.09255x10-6
                                                                              nthcomp Gamma 2.45877±4.68492x10-2
                                                                              kTe(keV) 1000.00±2.15271x10+5
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 1.66985x10-3±3.28833x10-5
                                                                              4)wabs*zwabs*(bbody+comptt) wabs nH*(1022) 1.72000x10-2 3.2609
                                                                              zwabs nH(1022) 6.85925x10-8±4.16595x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 1.42029±3.70398x10-2
                                                                              norm 1.56498x10-5±2.82521x10-6
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 9.70939x10-2±2.27574x10-3
                                                                              kT(keV) 90.8407±2278.92
                                                                              taup 1.73134x10-2±1.34060
                                                                              approx* 0.0
                                                                              norm 2.31987x10-4±5.83202x10-3
                                                                              5)wabs*zwabs*(bbody+nthcomp+comptt) wabs nH*(1022) 1.72000x10-2 2.6866
                                                                              zwabs nH(1022) 4.97783x10-3±2.14162x10-2
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 2.58079±23.3003
                                                                              norm 5.14265x10-6±1.91838x10-4
                                                                              nthcomp Gamma 4.67517±16.2787
                                                                              kTe(keV) 1.62813±21.6613
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 1.34427x10-3±1.04180x10-2
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 0.148355±0.130947
                                                                              kT(keV) 52.7752±3.60331x10+4
                                                                              taup 0.308383±289.939
                                                                              approx* 0.0
                                                                              norm 7.22215x10-5±4.90055x10-2
                                                                              MODEL COMPONENT PARAMETER VALUE±ERROR REDUCED χ2
                                                                              1)wabs*zwabs*pow wabs nH*(1022) 1.72000x10-2 24.8781
                                                                              zwabs nH(1022) 3.33424x10-9±2.97897
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.05625±3.88920x10-3
                                                                              norm 3.74562x10-3±6.13079x10-6
                                                                              2)wabs*zwabs*(pow+bbody) wabs nH*(1022) 1.72000x10-2 4.6252
                                                                              zwabs nH(1022) 8.35991x10-2±3.00305x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.78896±2.28523x10-2
                                                                              norm 4.48674x10-3±3.40643x10-5
                                                                              bbody kT(keV) 1.24392±1.72413x10-2
                                                                              norm 2.50135x10-5±4.00541x10-7
                                                                              3)wabs*zwabs*(bbody+nthcomp) wabs nH*(1022) 1.72000x10-2 3.2080
                                                                              zwabs nH(1022) 5.84361x10-10±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 0.134098±6.86813x10-4
                                                                              norm 1.00233x10-4±1.08691x10-6
                                                                              nthcomp Gamma 2.48241±1.61054x10-2
                                                                              kTe(keV) 1000.00±1.05529x10+5
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 2.35632x10-3±2.41520x10-5
                                                                              4)wabs*zwabs*(bbody+comptt) wabs nH*(1022) 1.72000x10-2 3.5244
                                                                              zwabs nH(1022) 6.98194x10-10±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 1.32167±5.72284x10-2
                                                                              norm 2.19974x10-5±3.64577x10-6
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 0.101827±4.72888x10-3
                                                                              kT(keV) 98.6492±5896.86
                                                                              taup 1.00000x10-2±4.07929
                                                                              approx* 0.0
                                                                              norm 2.95481x10-4±1.76900x10-2
                                                                              5)wabs*zwabs*(bbody+nthcomp+comptt) wabs nH*(1022) 1.72000x10-2 2.9049
                                                                              zwabs nH(1022) 8.59785x10-10±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 0.132927±1.41675x10-3
                                                                              norm 9.56510x10-5±2.58082x10-6
                                                                              nthcomp Gamma 2.54913±4.46431x10-2
                                                                              kTe(keV) 997.892±2.69894x10+5
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 2.45990x10-3±5.22694x10-5
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 19.8176±5745.83
                                                                              kT(keV) 500.000±3.63377x10+7
                                                                              taup 41.0030±1.69412x10+6
                                                                              approx* 0.0
                                                                              norm 2.15289x10-2±1542.34
                                                                              MODEL COMPONENT PARAMETER VALUE±ERROR REDUCED χ2
                                                                              1)wabs*zwabs*pow wabs nH*(1022) 1.72000x10-2 16.9736
                                                                              zwabs nH(1022) 4.04581x10-8±1.59923x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.03536±8.41937x10-3
                                                                              norm 3.01226x10-3±1.50876x10-5
                                                                              2)wabs*zwabs*(pow+bbody) wabs nH*(1022) 1.72000x10-2 4.6959
                                                                              zwabs nH(1022) 7.84638x10-2±3.50139x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.71992±2.75307x10-2
                                                                              norm 3.53915x10-3±3.05403x10-5
                                                                              bbody kT(keV) 1.23142±2.29903x10-2
                                                                              norm 2.15061x10-5±4.35623x10-7
                                                                              3)wabs*zwabs*(bbody+nthcomp) wabs nH*(1022) 1.72000x10-2 4.4590
                                                                              zwabs nH(1022) 1.06154x10-9±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 1.24314±2.32916x10-2
                                                                              norm 2.03566x10-5±4.39303x10-7
                                                                              nthcomp Gamma 3.51369±1.15980x10-2
                                                                              kTe(keV) 995.223±x1.2064310+5
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 2.99579x10-3±9.85313x10-6
                                                                              4)wabs*zwabs*(bbody+comptt) wabs nH*(1022) 1.72000x10-2 4.5557
                                                                              zwabs nH(1022) 6.52668x10-9±1.67823x10-4
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 1.32723±3.09316x10-2
                                                                              norm 1.90674x10-5±3.24671x10-6
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 0.100974±1.18370x10-3
                                                                              kT(keV) 75.3737±1866.07
                                                                              taup 2.65705x10-2±1.90094
                                                                              approx* 0.0
                                                                              norm 3.07655x10-4±7.63538x10-3
                                                                              5)wabs*zwabs*(bbody+nthcomp+comptt) wabs nH*(1022) 1.72000x10-2 3.4627
                                                                              zwabs nH(1022) 1.16749x10-6±1.21427x10-2
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 0.165557±0.557379
                                                                              norm 3.80793x10-5±1.44709x10-3
                                                                              nthcomp Gamma 8.14933±170.336
                                                                              kTe(keV) 1000.000±1.87847x10+6
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 5.84499x10-4±1.85886x10-2
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 0.191681±2.88912
                                                                              kT(keV) 104.599±3781.30
                                                                              taup 0.106258±7.68170
                                                                              approx* 0.0
                                                                              norm 3.37542x10-5±1.83796x10-3
                                                                              MODEL COMPONENT PARAMETER VALUE±ERROR REDUCED χ2
                                                                              1)wabs*zwabs*pow wabs nH*(1022) 1.72000x10-2 15.6867
                                                                              zwabs nH(1022) 2.61464x10-9±3.24397
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.00622±5.08263x10-3
                                                                              norm 2.36875x10-3±5.06479x10-6
                                                                              2)wabs*zwabs*(pow+bbody) wabs nH*(1022) 1.72000x10-2 166.5531
                                                                              zwabs nH(1022) 2.30013x10-8±9.80580x10-4
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex -0.713758±6.66585x10-2
                                                                              norm 6.12615x10-6±7.54522x10-7
                                                                              bbody kT(keV) 0.194105±3.71728x10-4
                                                                              norm 8.51399x10-5±2.50620x10-7
                                                                              3)wabs*zwabs*(bbody+nthcomp) wabs nH*(1022) 1.72000x10-2 2.9554
                                                                              zwabs nH(1022) 1.21663x10-9±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 0.129244±9.71287x10-4
                                                                              norm 6.07928x10-5±1.76214x10-6
                                                                              nthcomp Gamma 2.46998±4.14084x10-2
                                                                              kTe(keV) 1000.00±1.28953x10+5
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 1.57162x10-3±3.46909x10-5
                                                                              4)wabs*zwabs*(bbody+comptt) wabs nH*(1022) 1.72000x10-2 78.4056
                                                                              zwabs nH(1022) 7.22590x10-11±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 0.172975±1.27386x10-3
                                                                              norm 9.17977x10-5±8.64895x10-7
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 0.369034±0.118386
                                                                              kT(keV) 500.000±7.75400x10+4
                                                                              taup 31.9789±3243.60
                                                                              approx* 0.0
                                                                              norm 1.16726x10-3±0.206033
                                                                              5)wabs*zwabs*(bbody+nthcomp+comptt) wabs nH*(1022) 1.72000x10-2 3.0141
                                                                              zwabs nH(1022) 1.58463x10-10±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 0.128568±1.06868x10-3
                                                                              norm 5.95123x10-5±1.83031x10-6
                                                                              nthcomp Gamma 2.50071±5.80128x10-2
                                                                              kTe(keV) 1000.000±2.08079x10+5
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 1.60167x10-3±4.11917x10-5
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 98.9588±2.34835x10+5
                                                                              kT(keV) 500.000±4.98044x10+7
                                                                              taup 52.3542±2.92371x10+6
                                                                              approx* 0.0
                                                                              norm 1.16871±1.19581x10+5
                                                                              MODEL COMPONENT PARAMETER VALUE±ERROR REDUCED χ2
                                                                              1)wabs*zwabs*pow wabs nH*(1022) 1.72000x10-2 13.6436
                                                                              zwabs nH(1022) 8.37707x10-8±1.56797x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 2.99942±8.93618x10-3
                                                                              norm 2.84004x10-3±1.51325x10-5
                                                                              2)wabs*zwabs*(pow+bbody) wabs nH*(1022) 1.72000x10-2 2.9527
                                                                              zwabs nH(1022) 6.03218x10-2±3.89099x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.57065±2.96906x10-2
                                                                              norm 3.17723x10-3±3.10220x10-5
                                                                              bbody kT(keV) 1.27039±2.55482x10-2
                                                                              norm 1.93139x10-5±4.55528x10-7
                                                                              3)wabs*zwabs*(bbody+nthcomp) wabs nH*(1022) 1.72000x10-2 2.6213
                                                                              zwabs nH(1022) 2.00722x10-2±7.52258x10-4
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 1.24515±2.21716x10-2
                                                                              norm 1.90447x10-5±1.11517x10-6
                                                                              nthcomp Gamma 3.46345±5.63836x10-2
                                                                              kTe(keV) 463.251±4.76986x10+4
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 2.80471x10-3±1.01276x10-5
                                                                              4)wabs*zwabs*(bbody+comptt) wabs nH*(1022) 1.72000x10-2 2.3534
                                                                              zwabs nH(1022) 5.45947x10-10±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 1.33310±3.34597x10-2
                                                                              norm 1.72507x10-5±3.04018x10-6
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 9.74030x10-2±1.18549x10-3
                                                                              kT(keV) 100.242±2559.10
                                                                              taup 1.30389x10-2±1.07023
                                                                              approx* 0.0
                                                                              norm 2.27032x10-4±5.81346x10-3
                                                                              5)wabs*zwabs*(bbody+nthcomp+comptt) wabs nH*(1022) 1.72000x10-2 5.3155
                                                                              zwabs nH(1022) 2.22404x10-3±x5.0034810-2
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 2.15805x10-2±4.09431
                                                                              norm 2.27886x10-2±78.2571
                                                                              nthcomp Gamma 4.67813±2.94354
                                                                              kTe(keV) 12.9446±1704.37
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 1.38756x10-3±1.82598x10-3
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 0.155120±7.59645x10-2
                                                                              kT(keV) 52.7752±3.60331x10+4
                                                                              taup 0.341165±11.7323
                                                                              approx* 0.0
                                                                              norm 7.88281x10-5±1.93502x10-3
                                                                              MODEL COMPONENT PARAMETER VALUE±ERROR REDUCED χ2
                                                                              1)wabs*zwabs*pow wabs nH*(1022) 1.72000x10-2 14.6654
                                                                              zwabs nH(1022) 1.64745x10-9±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.03270±4.99826x10-3
                                                                              norm 2.26829x10-3±4.93389x10-6
                                                                              2)wabs*zwabs*(pow+bbody) wabs nH*(1022) 1.72000x10-2 14.9191
                                                                              zwabs nH(1022) 4.09184x10-9±7.98099x10-3
                                                                              Redshift* 4.81780x10-2
                                                                              powerlaw PhoIndex 3.02973±6.16764x10-3
                                                                              norm 2.26026x10-3±5.61326x10-6
                                                                              bbody kT(keV) 9.93761x10-3±1.57247x10+16
                                                                              norm 9.68636x10-10±6.66407x10+10
                                                                              3)wabs*zwabs*(bbody+nthcomp) wabs nH*(1022) 1.72000x10-2 2.0912
                                                                              zwabs nH(1022) 8.97276x10-10±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 0.132267±1.30321x10-3
                                                                              norm 6.07940x10-5±1.01894x10-6
                                                                              nthcomp Gamma 2.43516±2.27917x10-2
                                                                              kTe(keV) 1000.00±1.35185x10+5
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 1.42632x10-3±3.57143x10-5
                                                                              4)wabs*zwabs*(bbody+comptt) wabs nH*(1022) 1.72000x10-2 2.2135
                                                                              zwabs nH(1022) 1.44329x10-9±1.00000
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 1.34161±3.12829x10-2
                                                                              norm 1.44837x10-5±2.46322x10-6
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 9.93462x10-2±1.27580x10-3
                                                                              kT(keV) 94.4659±2703.27
                                                                              taup 1.30878x10-2±1.21951
                                                                              approx* 0.0
                                                                              norm 1.90986x10-4±5.48190x10-3
                                                                              5)wabs*zwabs*(bbody+nthcomp+comptt) wabs nH*(1022) 1.72000x10-2 1.6401
                                                                              zwabs nH(1022) 3.30653x10-2±7.74116x10-2
                                                                              Redshift* 4.81780x10-2
                                                                              bbody kT(keV) 3.26570x10-2±6.18345x10-2
                                                                              norm 6.74603x10-3±0.128110
                                                                              nthcomp Gamma 4.84156±4.44563
                                                                              kTe(keV) 1.62383±38.8928
                                                                              kT*bb(keV) 0.100000
                                                                              inp*type 0.0
                                                                              Redshift* 4.81780x10-2
                                                                              norm 1.41492x10-3±1.54098x10-3
                                                                              comptt Redshift* 4.81780x10-2
                                                                              T0(keV) 0.139449±0.172681
                                                                              kT(keV) 53.6914±1424.76
                                                                              taup 0.370376±13.0472
                                                                              approx* 0.0
                                                                              norm 6.11173x10-5±1.56782x10-3

                                                                              * indicates frozen parameter value

                                                                              5. CONCLUSIONS

                                                                              ∙ Through lightcurve analysis we have found that the source is variable for all the observattion ids analysed.

                                                                              ∙ Through spectral fitting using the different models specified, we have calculaed the coronal temperaturre.


                                                                              I would like to thank Indian Academy of Sciences for giving me this wonderful opportunity of working at Indian Institute of Astrophysics. This research internship has been a wonderful learning experience.

                                                                              I am extremely grateful to my guide Prof. C. S. Stalin, whose insight, encouragement and supervision have been invaluable.

                                                                              I would also like to thank Ms. Priyanka Rani, Indian Institute of Astrophysics for her constant support and for patiently answering my questions even though they were innumerable.


                                                                              • Urry C.M and Padovani P., (1995), Unified Schemes for Radio-Loud Active Galactic Nuclei, Publications of the Astronomical Society of the Pacific, 107, 803-845.

                                                                              • C. Jin, M.J. Ward, C. Done, and J.M. Gelbord, (2012), A combined optical and X-ray study of Seyfert 1 active galactic nuclei. I. Optical spectra and spectral energy distribution modeling, Monthly Notices of the Royal Astronomical Society (MNRAS), 420, 1825.

                                                                              • Boroson T. A., Green R. F., (1992), The emission-line properties of low-redshift quasi-stellar objects, Astrophysical Journal Supplement Series, 80, 109

                                                                              • C. Jin, C. Done, M. Middleton and M. Ward, (2013), A long XMM–Newton observation of an extreme narrow-line Seyfert 1: PG 1244+026, Monthly Notices of the Royal Astronomical Society (MNRAS), 436, 3173–3185.

                                                                              • Th. Boller, W. N. Brandt, H. Fink, (1996), Soft X-ray properties of "narrow-line" Seyfert 1 galaxies, Journal of Astronomy and Astrophysics, 305, 53.

                                                                              • Th Boller, (1997), Narrow‐line Seyfert 1 galaxies, Astronomische Nachrichten, 318, 209.​

                                                                              • P. Rani, C. S. Stalin and S. Rakshit, (2017), X-ray flux variability of active galactic nuclei observed using NuSTAR, Monthly Notices of the Royal Astronomical Society, 466, 3309–3322.

                                                                              • L. Titarchuk, (1994), Generalized Comptonization models and application to the recent high energy observations, American Journal of Physics conference proceedings, 304, 380-384.

                                                                              Written, reviewed, revised, proofed and published with