Solar flares are known to generate seismic waves in the Sun. We present a detailed analysis of seismic emission in sunspots accompanying M- and X-class solar flares. For this purpose, we have used high-resolution Dopplergrams and line-of-sight magnetograms at a cadence of 45 s, along with vector magnetograms at a cadence of 135 s obtained from Helioseismic and Magnetic Imager (HMI) instrument aboard the Solar Dynamic Observatory (SDO) space mission. In order to identify the location of flare ribbons and hard X-ray foot-points, we have also used H-alpha chromospheric intensity observations obtained from Global Oscillation Network Group (GONG) instruments and hard X-ray images in 12-25 KeV band from the Reuvan Ramaty High Energy Solar Spectroscopic Imager (RHESSI) spacecraft. The Fast Fourier Transform (FFT) technique is applied to construct the acoustic velocity power map in 2.5-4 mHz band for pre-flare, spanning flare, and post flare epochs for the identification of seismic emission locations in the sunspots. In the power maps, we have selected only those locations which are away from the flare ribbons and hard X-rays foot-points. These regions are believed to be free from any flare related artefacts in the observational data. We have identified concentrated locations of acoustic power enhancements in sunspots accompanying major flares. Our investigation provides evidence that abrupt changes in the magnetic fields and associated impulsive changes in the Lorentz force could be the driving source for these seismic emissions in the sunspots during solar flares.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

We present radio and optical afterglow observations of the TeV-bright long Gamma Ray Burst (GRB) 190114C, which was detected by the MAGIC telescope at a redshift of z=0.425. Our observations with ALMA, ATCA, and uGMRT were obtained by our low frequency observing campaign and range from ~1 to ~140 days after the burst and the optical observations were done with the 0.7-m GROWTH-India telescope up to ~25 days after the burst. Long term radio/mm observations reveal the complex nature of the afterglow, which does not conform to the predictions of the standard afterglow model. We find that the microphysical parameters of the external forward shock, representing the share of shock-created energy in the non-thermal electron population and magnetic field, are evolving with time. The kinetic energy in the blast-wave is almost an order of magnitude higher than that measured in the prompt emission.

Solar X-ray Monitor (XSM) is one of the scientific instruments on-board Chandrayaan-2 orbiter. The XSM along with instrument CLASS (Chandras Large Area Soft x-ray Spectrometer) comprise the remote X-ray fluorescence spectroscopy experiment of Chandrayaan-2 mission with an objective to determine the elemental composition of the lunar surface on a global scale. XSM instrument will measure the solar X-rays in the energy range of 1-15 keV using state-of-the-art Silicon Drift Detector (SDD). The Flight Model (FM) of the XSM payload has been designed, realized and characterized for various operating parameters. XSM provides energy resolution of 180 eV at 5.9 keV with high time cadence of one second. The X-ray spectra of the Sun observed with XSM will also contribute to the study of solar corona. The detailed description and the performance characteristics of the XSM instrument are presented in this paper.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 1987A. A $\rm^{56}Ni$ mass of ~0.087$\pm$0.004 $\rm M_{\odot}$ is inferred for SN 2018hna from its bolometric light curve. Hydrodynamical modelling of the cooling phase suggests a progenitor with a radius ~50 $\rm R_{\odot}$, a mass of ~14-20 $\rm M_{\odot}$ and explosion energy of ~1.7-2.9$\rm \times$ $\rm 10^{51}\ erg$. The smaller inferred radius of the progenitor than a standard red supergiant is indicative of a blue supergiant progenitor of SN 2018hna. A sub-solar metallicity (~0.3 $\rm Z_{\odot}$) is inferred for the host galaxy UGC 07534, concurrent with the low-metallicity environments of 1987A-like events.

Galaxy clusters are known to be reservoirs of Cosmic Rays (CRs), as inferred from theoretical calculations or detection of CR-derived observables. CR acceleration in clusters is mostly attributed to the dynamical activity that produces shocks. Shocks in clusters emerge out of merger or accretion, but which one is more effective in producing CRs? at which dynamical phase? and why? To this aim, we study the production or injection of CRs through shocks and its evolution in the galaxy clusters using cosmological simulations with the {\sc enzo} code. Particle acceleration model considered here is primarily the Diffusive Shock Acceleration (DSA) of thermal particles, but we also report a tentative study with pre-existing CRs. Defining appropriate dynamical states using the concept of virialization, we studied a sample of merging and non-merging clusters. We report that the merger shocks (with Mach number $\mathcal{M}\sim2-5$) are the most effective CR producers, while high-Mach peripheral shocks (i.e. $\mathcal{M}>5$) are mainly responsible for the brightest phase of CR injection in clusters. Clusters once merged, permanently deviate from CR and X-ray mass scaling of non-merging systems, enabling us to use it as a tool to determine the state of merger. Through a temporal and spatial evolution study, we found a strong correlation between cluster merger dynamics and CR injection. We observed that the brightest phase of X-ray and CR injection from clusters occur respectively at about 1.0 and 1.5 Gyr after every mergers, and CR injection peaks near to the cluster virial radius (i.e $r_{200}$). Delayed CR injection peaks found in this study deserve further investigation for possible impact on the evolution of CR-derived observables from galaxy clusters.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.

`Galaxy groups' have hardly been realised as a separate class of objects with specific characteristics in the structural hierarchy. The presumption that the self-similarity of dark matter structures is a valid prescription for the baryonic universe at all scales has rendered smaller structures undetectable by current observational facilities, leading to lesser dedicated studies on them. Some recent reports that indicate a deviation from $\rm{L_x}$-T scaling in groups compared to clusters have motivated us to study their physical properties in depth. In this article, we report the extensive study on physical properties of groups in comparison to the clusters through cosmological hydrodynamic plus N-body simulations using ENZO 2.2 code. As additional physics, radiative cooling, heating due to supernova and star motions, star formation and stellar feedback has been implemented. We have produced a mock sample of 362 objects with mass ranging from $5\times10^{12}\; \rm{M_{\odot}}$ to 2.5$\times 10^{15}\; \rm{M_{\odot}}$. Strikingly, we have found that objects with mass below $\sim$ $8\times 10^{13}\;\rm{M_{\odot}}$ do not follow any of the cluster self-similar laws in hydrostatics, not even in thermal and non-thermal energies. Two distinct scaling laws are observed to be followed with breaks at $\sim$ $8\times 10^{13}\;\rm{M_{\odot}}$ for mass, $\sim$1 keV for temperature and $\sim$1 Mpc for radius. This places groups as a distinct entity in the hierarchical structures, well demarcated from clusters. This study reveals that groups are mostly far away from virialization, suggesting the need for formulating new models for deciphering their physical parameters. They are also shown to have high turbulence and more non-thermal energy stored, indicating better visibility in the non-thermal regime.

We present a Chandra and XMM-Newton study of the supernova remnant (SNR) Kes 73 hosting the anomalous X-ray pulsar 1E 1841-045. The Chandra image reveals clumpy structures across the remnant with enhanced emission along the western rim. The X-ray emission fills the radio shell and spatially correlates with the infrared image. The global X-ray spectrum is described by a two-component thermal model with a column density N_H ~ 2.6e22 cm^{-2} and a total luminosity of L_X ~ 3.3e37 ergs/s (0.5-10 keV, at an assumed distance of 8.5 kpc). The soft component is characterized by a temperature kT_s ~ 0.5 keV, a high ionization timescale, and enhanced Si and S abundances suggesting emission that is dominated by shocked ejecta. The hard component has a temperature kT_h ~ 1.6 keV, a relatively low ionization timescale, and mostly solar abundances suggesting emission that is dominated by interstellar/circumstellar shocked material. A spatially resolved spectroscopy study reveals no significant variations in the spectral properties. We infer an SNR age ranging between 750 yr and 2100 yr, an explosion energy of ~0.3e51 ergs and a shock velocity of 1200 km/s (under the Sedov phase assumption). We also discuss the possible scenario for Kes 73 expanding into the late red supergiant wind phase of its massive progenitor. Comparing the inferred metal abundances to core-collapse nucleosynthesis model yields, we estimate a progenitor mass >20 solar masses, adding a candidate to the growing list of highly magnetized neutron stars proposed to be associated with very massive progenitors.

We explore the low energy cosmological dynamics of feebly self-interacting cold dark matter and propose a new simple explanation for the rotation curves of the core-halo model in massive LSB (Low Surface brightness)galaxies. We argue in favor of the truly collisionless nature of cold dark matter,which is feebly,self-interacting at small scales between epochs of equality and recombination.For this, we assume a model, wherein strongly coupled baryon-radiation plasma ejects out of small regions of concentrated cold dark matter without losing its equilibrium. We use the Merscerskii equation i.e. the variable mass formalism of classical dynamics.We obtain new results relating the oscillations in the CMB anisotropy to the ejection velocity of the baryon-radiation plasma,which can be useful tool for numerical work for exploring the second peak of CMB. Based on this model, we discuss the growth of perturbations in such a feebly self-interacting,cold dark matter both in the Jeans theory and in the expanding universe using Newton's theory.We obtain an expression for the growth of fractional perturbations in cold dark matter,which reduce to the standard result of perturbation theory for late recombination epochs. We see the effect of the average of the perturbations in the cold dark matter potential on the cosmic microwave background temperature anisotropy that originated at redshifts between equality and recombination i.e. 1100 < z < z_{eq}. Also we obtain an expression for the Sachs-Wolfe effect,i.e. the CMB temperature anisotropy at decoupling in terms of the average of the perturbations in cold dark matter potential.

With an aim to argue for the truly collisionless nature of cold dark matter between epochs of equality and recombination, we assume a model, wherein strongly coupled baryon-radiation plasma ejects out of small regions of concentrated cold dark matter without losing its equilibrium. We use the Meschersky equation to describe the dynamics of cold dark matter in the presence of varying mass of strongly coupled baryon-radiation plasma. Based on this model, we discuss the growth of perturbations in cold dark matter both in the Jeans theory and in the expanding universe using Newton's theory. We see the effect of the perturbations in the cold dark matter potential on the cosmic microwave background anisotropy that originated at redshifts between equality and recombination i.e. $1100 < z < z_{eq}$. Also we obtain an expression for the Sachs-Wolfe effect, i.e. the CMB temperature anisotropy at decoupling in terms of the perturbations in cold dark matter potential. We obtain similar solutions both in the static and in the expanding universe, for epochs of recombination. From this, we infer about the time scale when the dark energy starts to dominate.

The distinction between the high-magnetic field pulsars (HBPs, thought to be mainly rotation-powered) and magnetars (commonly believed to be powered by their super-strong magnetic fields) has been recently blurred with the discovery of magnetar-like activity from the HBP J1846-0258 in the Supernova Remnant (SNR) Kes 75. What determines the spin properties of a neutron star at birth and its manifestation as a magnetar-like or more classical pulsar is still not clear. Furthermore, although a few studies have suggested very massive progenitors for magnetars, there is currently no consensus on the progenitors of these objects. To address these questions, we examine their environments by studying or revisiting their securely associated SNRs. Our approach is by: 1) inferring the mass of their progenitor stars through X-ray spectroscopic studies of the thermally emitting supernova ejecta, and 2) investigating the physical properties of their hosting SNRs and ambient conditions. We here highlight our detailed studies of two SNRs: G292.2-0.5 associated with the HBP J1119-6127 and Kes 73 associated with the AXP 1E 1841-045, and summarize the current view of the other (handful) HBP/magnetar-SNR associations.

We present the first detailed imaging and spatially resolved spectroscopic study of the Galactic supernova remnant G292.2-0.5, associated with the high-magnetic field radio pulsar J1119-6127, using Chandra and XMM-Newton. The high-resolution X-ray images reveal a partially limb-brightened morphology in the west, with diffuse emission concentrated towards the interior of the remnant unlike the complete shell-like morphology observed at radio wavelengths. The spectra of most of the diffuse emission regions within the remnant are best described by a thermal+non-thermal model. The thermal component is described by a plane-parallel, non-equilibrium ionization plasma model with a temperature ranging from 1.3 keV in the western side of the remnant to 2.3 keV in the east, a column density increasing from 1.0e22 cm^-2 in the west to 1.8e22 cm^-2 in the east, and a low ionization timescale ranging from 5.7e9 cm^-3 s in the SNR interior to 3.6e10 cm^-3 s in the western side - suggestive of expansion of a young remnant in a low-density medium. The spatial and spectral differences across the SNR are consistent with the presence of a dark cloud in the eastern part of the SNR, absorbing the soft X-ray emission. The spectra from some of the regions also show slightly enhanced metal abundances from Ne, Mg and Si, hinting at the first evidence for ejecta heated by the reverse shock. Comparing our inferred metal abundances to core-collapse nucleosynthesis models yields, we estimate a high progenitor mass of ~30 solar mass suggesting a type Ib/c supernova. We confirm the presence of non-thermal hard X-ray emission from regions close to the pulsar. We estimate an SNR age range between 4.2 kyr (free expansion phase) and 7.1 kyr (Sedov phase), a factor of a few higher than the measured pulsar's age upper limit of 1.9 kyr.

We report the study of the short (32 ms) and first SGR-like burst observed from the anomalous X-ray pulsar (AXP) 1E 1841-045 associated with the supernova remnant (SNR) Kes 73, discovered on 2010 May 6 by the Burst Alert Telescope (BAT) onboard the Swift gamma-ray observatory. The 15-100 keV time-averaged burst spectrum is modeled by a single power-law (PL) with a photon index Gamma=3.2$^{+1.8}_{-1.0}$, and has a fluence of 1.1$^{+0.4}_{-0.6}$$\times$10$^{-8}$ ergs cm$^{-2}$, luminosity of 2.9$^{+1.1}_{-1.6}$$\times$10$^{39}$ ergs s$^{-1}$, and energy of 7.2$^{+0.4}_{-0.6}$$\times$10$^{36}$ ergs. The prompt after-burst 0.5-10 keV quiescent spectrum obtained with the Swift X-ray Telescope (XRT) is best-fit by an absorbed PL model with Gamma=2.6$\pm$0.2 and an unabsorbed flux of 9.1$^{+1.2}_{-1.4}$$\times$10$^{-11}$ ergs cm$^{-2}$ s$^{-1}$. To investigate the pre-burst 0.5-10 keV persistent emission, we analyzed the archival XMM-Newton observations and the spectra are well fitted by a two-component blackbody (BB) plus PL model with a temperature kT=0.45$\pm$0.03 keV, Gamma=1.9$\pm$0.2, and an unabsorbed flux of 4.3$^{+0.9}_{-1.2}$$\times$10$^{-11}$ ergs cm$^{-2}$ s$^{-1}$. Comparing the Swift-XRT spectrum with the XMM-Newton spectrum, spectral softening post-burst is evident with a 2.1 times increase in the unabsorbed flux. We discuss the burst activity and the persistent emission properties of AXP 1E 1841-045 in comparison with other magnetars and in the context of the magnetar model.

We present results on the Soft Gamma Repeater (SGR) 0501+4516, discovered by the SWIFT Burst Alert Telescope (BAT) on 2008 August 22. More than 50 bursts were identified from this source, out of which 18 bursts had enough counts to carry out spectral analysis. We performed time-averaged spectral analysis on these 18 bursts using 8 models, among which the cut-off powerlaw and the two-blackbody models provided the best fit in the 15-150 keV energy range. The cut-off powerlaw model fit yields a mean photon index Gamma_{CPL} = 0.54+/-0.11 and a cut-off energy E_C = 19.1+/-1.8 keV for the bursts. The mean hard and soft blackbody temperatures are found to be kT_{BB_h} = 12.8+/-0.7 keV and kT_{BB_s} = 4.6+/-0.5 keV, respectively, and are anti-correlated with the square of the radii of the hard and soft emitting regions (R_{BB_h} and R_{BB_s}) as R_{BB_h}^2 $\propto$ kT^{-5.8} and R_{BB_s}^2 $\propto$ kT^{-2.7}, respectively. The soft and hard component temperatures with different indices support the idea of two distinct emitting regions with the hard component corresponding to a smaller radius and the soft component corresponding to a larger radius, which further corroborate the idea of the propagation of extraordinary (E) and ordinary (O) mode photons across the photosphere, as predicted in the magnetar model. We notice strong burst fluence-duration correlation as well as hardness ratio-duration and hardness ratio-fluence anti-correlations for the SGR 0501+4516 bursts. The burst fluences range from ~ 4.4\times10^{-9} ergs/cm^-2 to ~ 2.7\times10^{-6} ergs/cm^{-2}, consistent with those observed for typical short SGR bursts.

(shortened) PSR J1119-6127 is a high magnetic field (B=4.1E13 Gauss), young (<=1,700 year-old), and slow (P=408 ms) radio pulsar associated with the supernova remnant (SNR) G292.2-0.5. In 2003, Chandra allowed the detection of the X-ray counterpart of the radio pulsar, and provided the first evidence for a compact pulsar wind nebula (PWN). We here present new Chandra observations which allowed for the first time an imaging and spectroscopic study of the pulsar and PWN independently of each other. The PWN is only evident in the hard band and consists of jet-like structures extending to at least 7" from the pulsar, with the southern `jet' being longer than the northern `jet'. The spectrum of the PWN is described by a power law with a photon index~1.1 for the compact PWN and ~1.4 for the southern long jet (at a fixed column density of 1.8E22/cm2), and a total luminosity of 4E32 ergs/s (0.5-7 keV), at a distance of 8.4 kpc. The pulsar's spectrum is clearly softer than the PWN's spectrum. We rule out a single blackbody model for the pulsar, and present the first evidence of non-thermal (presumably magnetospheric) emission that dominates above ~3keV. A two-component model consisting of a power law component (with photon index ~1.5--2.0) plus a thermal component provides the best fit. The thermal component can be fit by either a blackbody model with a temperature kT~0.21 keV, or a neutron star atmospheric model with a temperature kT~0.14 keV. The efficiency of the pulsar in converting its rotational power, Edot, into non-thermal X-ray emission from the pulsar and PWN is ~5E-4, comparable to other rotation-powered pulsars with a similar Edot. We discuss our results in the context of the X-ray manifestation of high-magnetic field radio pulsars in comparison with rotation-powered pulsars and magnetars.

We present results from the archival Chandra observations of the 0.3 s X-ray pulsar PSR J1846-0258 associated with the supernova remnant (SNR) Kes 75. The pulsar has the highest spin-down luminosity (Edot = 8.3e36 erg/s) among all the high magnetic field pulsars (HBPs) and has been classified as a Crab-like pulsar despite its magnetic field (5e13 G) being above the quantum critical field. It is the only HBP described by a non-thermal Crab-like spectrum, powering a bright pulsar wind nebula (PWN). Our spectroscopic study shows evidence of spectral softening (photon index = 1.32 to 1.97) and temporal brightening (unabsorbed flux = 4.3e-12 to 2.7e-11 erg/cm^2/s) of the pulsar by ~6 times from 2000 to 2006. The 0.5-10 keV luminosity of the pulsar at the revised distance of 6 kpc has also increased from 1.85e34 to 1.16e35 erg/s, and the X-ray efficiency increased from 0.2% to 1.4%. The observed X-ray brightening and softening of the pulsar suggests for the first time that this HBP is revealing itself as a magnetar.