BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Aneta Siemiginowska

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.COastro-ph.IMastro-phastro-ph.GAApplicationscs.CYDigital LibrariesMethodologyphysics.data-an

Catalog footprint

What is connected

38works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

38 published item(s)

preprint2026arXiv

The LIRA-Ising Model: Estimating the boundaries of irregularly shaped X-ray sources

Mapping the boundary of an extended source is a key step in the study of its morphology. The background contamination and statistical fluctuations of typical astronomical images make this a challenging statistical task, particularly for X-ray images with low surface brightness. We develop a three-step Bayesian procedure to identify the boundaries of irregularly shaped sources. We first apply a Bayesian multiscale reconstruction algorithm known as LIRA to obtain posterior pixelwise probability distributions of the source intensity that properly account for known structures, astrophysical background, and the effect of the telescope point spread function. Next, we adopt an Ising model to group pixels with similar intensities into cohesive regions corresponding to background and source. Finally, the boundary is derived on the basis of the most likely aggregation of pixels into the source region. Because the overall model combines LIRA and the Ising model, we call it LIRA-Ising. We verify the proposed method using a set of simulation studies. We then apply it to the Chandra X-ray Observatory images of two high redshift quasars, PKS J1421-0643 and 0730+257, to determine the extent and morphology of X-ray jets. Our method shows a uniform X-ray surface brightness of PKS J1421-0643 jet, and identifies knotty structure in the X-ray jet of 0730+257.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

The Origin of High Energy Emission in the Young Radio Source PKS 1718-649

We present a model for the broadband radio-to-$γ$-ray spectral energy distribution of the compact radio source, PKS 1718-649. Because of its young age (100 years) and proximity ($z=0.014$), PKS 1718-649 offers a unique opportunity to study nuclear conditions and the jet/host galaxy feedback process at the time of an initial radio jet expansion. PKS 1718-649 is one of a handful of young radio jets with $γ$-ray emission confirmed with the Fermi/LAT detector. We show that this $γ$-ray emission can be successfully explained by Inverse Compton scattering of the ultraviolet photons, presumably from an accretion flow, off non-thermal electrons in the expanding radio lobes. The origin of the X-ray emission in PKS 1718-649 is more elusive. While Inverse Compton scattering of the infrared photons emitted by a cold gas in the vicinity of the expanding radio lobes contributes significantly to the X-ray band, the data require that an additional X-ray emission mechanism is at work, e.g. a weak X-ray corona or a radiatively inefficient accretion flow, expected from a LINER type nucleus such as that of PKS 1718-649. We find that the jet in PKS 1718-649 has low power, $L_j \simeq 2.2 \times 10^{42}$ erg s$^{-1}$, and expands in an environment with density $n_0 \simeq 20$ cm$^{-3}$. The inferred mass accretion rate and gas mass reservoir within 50-100 pc are consistent with estimates from the literature obtained by tracing molecular gas in the innermost region of the host galaxy with SINFONI and ALMA.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Late-Time X-ray Observations of the Transient Source Cygnus A-2

We examine Chandra observations of the powerful Fanaroff-Riley class II (FR II) radio galaxy Cygnus A for an X-ray counterpart to the radio transient Cygnus A-2 that was first detected in 2011. Observations are performed using the High-Resolution Camera (HRC) instrument in order to spatially resolve Cygnus A-2 and the central Active Galactic Nucleus (AGN) at a separation of 0.42 arcseconds. Simulated images are generated of the emission region, and radial profiles for the region of interest are extracted. A comparison between the simulations and observations reveals no X-ray detection of Cygnus A-2 to a 0.5-7.0 keV flux upper limit of $1.04 \times 10^{-12}\rm\,erg\,cm^{-2}\,s^{-1}$, or a rest-frame 2-10 keV luminosity of $8.6\times 10^{42}\rm\,erg\,s^{-1}$. We estimate the black hole mass of Cygnus A-2 based on our X-ray flux limit and find it to be consistent with a flaring black hole rather than a steadily accreting source. The HRC observations are additionally compared with archival ACIS data from 2016-2017, and both the overall morphology and the flux limits of the AGN complex agree between the two datasets. This consistency is despite the pile-up effect in ACIS which was previously considered to bias the observed morphology of the AGN. The agreement between the datasets demonstrates the viability of utilizing the archival Chandra data of Cygnus A to analyze its AGN at an unprecedented level of precision.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Studies on slow radio transients

We present a brief overview of a very extensive studies of the group of active galactic nuclei (AGNs) that transitioned to radio-loud state over the past few decades. The sample consists of twelve sources, both quasars and galaxies, showing the characteristics of gigahertz-peaked spectrum (GPS) objects undergoing relatively rapid changes, due to the evolution of their newly-born radio jets. Discussed objects also show a wide range of physical parameters such as bolometric luminosity, black hole mass and jet power, suggesting a great diversity among young active galactic nuclei and their hosts. Furthermore, we introduce a new observational project, the aim of which will be to investigate and gain a more in-depth understanding of the phenomenon of slow radio transients.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

1321+045: a Compact Steep Spectrum radio source in a cool-core galaxy cluster

Cluster-central gigahertz peak and compact steep spectrum (CSS) sources offer an opportunity to study the earliest phases of AGN feedback, but few have yet been examined in detail. We present results from radio and X-ray observations of 1321+045, a CSS source in a 4.4 keV cluster at z=0.263. The cluster has a strongly cooling core, and disturbances from a minor cluster merger may have triggered a period of jet activity which formed the 16 kpc radio lobes 2.0 [+0.3,-0.2] Myr ago. However, new VLBA imaging shows a ~20 pc jet on a different projected axis, which is probably only a few hundred years old. We consider possible histories for the system, with either one or two periods of jet activity. While this single system is informative, a broader study of the youngest cluster-central radio sources is desirable.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Role of Machine Learning in the Next Decade of Cosmology

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Powerful ionized gas outflows in the interacting radio galaxy 4C +29.30

We investigate the ionized gas excitation and kinematics in the inner $4.3 \times 6.2$ kpc$^{2}$ of the merger radio galaxy 4C +29.30. Using optical integral field spectroscopy with the Gemini North Telescope, we present flux distributions, line-ratio maps, peak velocities and velocity dispersion maps as well as channel maps with a spatial resolution of $\approx 955$ pc. We observe high blueshifts of up to $\sim -650$ km s$^{-1}$, in a region $\sim 1''$ south of the nucleus (the southern knot, SK), which also presents high velocity dispersions ($\sim 250$ km s$^{-1}$), which we attribute to an outflow. A possible redshifted counterpart is observed north from the nucleus (the northern knot, NK). We propose that these regions correspond to a bipolar outflow possibly due to the interaction of the radio jet with the ambient gas. We estimate a total ionized gas mass outflow rate of $\dot{M}_{out} = 25.4 \substack{+11.5 \\ -7.5}$ M$_\odot$ yr$^{-1}$ with a kinetic power of $\dot{E} = 8.1 \substack{+10.7 \\ -4.0} \times 10^{42}$ erg s$^{-1}$, which represents $5.8 \substack{+7.6 \\ -2.9} \%$ of the AGN bolometric luminosity. These values are higher than usually observed in nearby active galaxies with the same bolometric luminosities and could imply a significant impact of the outflows in the evolution of the host galaxy. The excitation is higher in the NK (that correlates with extended X-ray emission, indicating the presence of hotter gas) than in the SK, supporting a scenario in which an obscuring dust lane is blocking part of the AGN radiation to reach the southern region of the galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

X-ray Properties of Young Radio Quasars at z > 4.5

We present a comprehensive analysis of Chandra X-ray observations of 15 young radio quasars at redshifts $4.5 < z < 5.0$. All sources are detected in the $0.5-7.0$ keV energy band. Emission spectra are extracted, and the average photon index for the sample is measured to be $1.5\pm0.1$. Unabsorbed rest-frame $2-10$ keV luminosities are found to range between $(0.5-23.2) \times 10^{45}$ erg s$^{-1}$. The optical-X-ray power-law spectral index $α_{ox}$ is calculated for each source using optical/UV data available in the literature. The $α_{ox}$-UV relationship is compared with other quasar surveys, and an anticorrelation is observed that agrees with independent estimates. Rest-frame radio and X-ray luminosities are established for the sample, and a correlation between the luminosities is detected. These multiwavelength results reinforce a lack of spectral evolution for quasars over a broad redshift range. We additionally identify three quasars from our multiwavelength analysis that are statistically significant outliers, with one source being a Compton-thick candidate in the early universe, and discuss each in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detecting Relativistic X-ray Jets in High-Redshift Quasars

We analyze Chandra X-ray images of a sample of 11 quasars that are known to contain kiloparsec scale radio jets. The sample consists of five high-redshift (z >= 3.6) flat-spectrum radio quasars, and six intermediate redshift (2.1 < z < 2.9) quasars. The dataset includes four sources with integrated steep radio spectra and seven with flat radio spectra. A total of 25 radio jet features are present in this sample. We apply a Bayesian multi-scale image reconstruction method to detect and measure the X-ray emission from the jets. We compute deviations from a baseline model that does not include the jet, and compare observed X-ray images with those computed with simulated images where no jet features exist. This allows us to compute p-value upper bounds on the significance that an X- ray jet is detected in a pre-determined region of interest. We detected 12 of the features unambiguously, and an additional 6 marginally. We also find residual emission in the cores of 3 quasars and in the background of 1 quasar that suggest the existence of unresolved X-ray jets. The dependence of the X-ray to radio luminosity ratio on redshift is a potential diagnostic of the emission mechanism, since the inverse Compton scattering of cosmic microwave background photons (IC/CMB) is thought to be redshift dependent, whereas in synchrotron models no clear redshift dependence is expected. We find that the high-redshift jets have X-ray to radio flux ratios that are marginally inconsistent with those from lower redshifts, suggesting that either the X-ray emissions is due to the IC/CMB rather than the synchrotron process, or that high redshift jets are qualitatively different.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

First detection in gamma-rays of a young radio galaxy: Fermi-LAT observations of the Compact Symmetric Object PKS 1718-649

We report the $γ$-ray detection of a young radio galaxy, PKS 1718$-$649, belonging to the class of Compact Symmetric Objects (CSOs), with the Large Area Telescope (LAT) on board the {\it Fermi} satellite. The third {\it Fermi} Gamma-ray LAT catalog (3FGL) includes an unassociated $γ$-ray source, 3FGL J1728.0$-$6446, located close to PKS 1718$-$649. Using the latest Pass 8 calibration, we confirm that the best fit $1 σ$ position of the $γ$-ray source is compatible with the radio location of PKS 1718$-$649. Cross-matching of the $γ$-ray source position with the positions of blazar sources from several catalogs yields negative results. Thus, we conclude that PKS 1718$-$649 is the most likely counterpart to the unassociated LAT source. We obtain a detection test statistics TS$\sim 36$ ($>$5$σ$) with a best fit photon spectral index $Γ=$2.9$\pm$0.3 and a 0.1-100 GeV photon flux density $F_{\rm 0.1-100GeV}=$(11.5$\pm$0.3)$\times{\rm 10^{-9}}$ ph cm$^{-2}$ s$^{-1}$. We argue that the linear size ($\sim$2 pc), the kinematic age ($\sim$100 years), and the source distance ($z=0.014$) make PKS 1718$-$649 an ideal candidate for $γ$-ray detection in the framework of the model proposing that the most compact and the youngest CSOs can efficiently produce GeV radiation via inverse-Compton scattering of the ambient photon fields by the radio lobe non-thermal electrons. Thus, our detection of the source in $γ$-rays establishes young radio galaxies as a distinct class of extragalactic high-energy emitters, and yields an unique insight on the physical conditions in compact radio lobes interacting with the interstellar medium of the host galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

X-ray properties of the Youngest Radio Sources and their Environments

We present the results of the first X-ray study of a sample of 16 young radio sources classified as Compact Symmetric Objects (CSOs). We observed six of them for the first time in X-rays using {\it Chandra}, re-observed four with the previous {\it XMM-Newton} or {\it Beppo-SAX} data, and included six other with the archival data. All the sources are nearby, $z<1$ with the age of their radio structures ($<3000$~years) derived from the hotspots advance velocity. Our results show heterogeneous nature of the CSOs indicating a complex environment associated with young radio sources. The sample covers a range in X-ray luminosity, $L_{2-10\,\rm keV} \sim 10^{41}$-$10^{45}$\,erg\,s$^{-1}$, and intrinsic absorbing column density of $N_H \simeq 10^{21}$--10$^{22}$\,cm$^{-2}$. In particular, we detected extended X-ray emission in 1718$-$649; a hard photon index of $Γ\simeq 1$ in 2021$+$614 and 1511$+$0518 consistent with either a Compton thick absorber or non-thermal emission from compact radio lobes, and in 0710$+$439 an ionized iron emission line at $E_{rest}=(6.62\pm0.04)$\,keV and EW $\sim 0.15-$1.4\,keV, and a decrease by an order of magnitude in the 2-10 keV flux since the 2008 {\it XMM-Newton} observation in 1607$+$26. We conclude that our pilot study of CSOs provides a variety of exceptional diagnostics and highlights the importance of deep X-ray observations of large samples of young sources. This is necessary in order to constrain theoretical models for the earliest stage of radio source evolution and study the interactions of young radio sources with the interstellar environment of their host galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detecting Unspecified Structure in Low-Count Images

Unexpected structure in images of astronomical sources often presents itself upon visual inspection of the image, but such apparent structure may either correspond to true features in the source or be due to noise in the data. This paper presents a method for testing whether inferred structure in an image with Poisson noise represents a significant departure from a baseline (null) model of the image. To infer image structure, we conduct a Bayesian analysis of a full model that uses a multiscale component to allow flexible departures from the posited null model. As a test statistic, we use a tail probability of the posterior distribution under the full model. This choice of test statistic allows us to estimate a computationally efficient upper bound on a p-value that enables us to draw strong conclusions even when there are limited computational resources that can be devoted to simulations under the null model. We demonstrate the statistical performance of our method on simulated images. Applying our method to an X-ray image of the quasar 0730+257, we find significant evidence against the null model of a single point source and uniform background, lending support to the claim of an X-ray jet.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optical Detection of the Pictor A Jet and Tidal Tail: Evidence against an IC/CMB jet

New images from the Hubble Space Telescope of the FRII radio galaxy Pictor A reveal a previously undiscovered tidal tail, as well as a number of jet knots coinciding with a known X-ray and radio jet. The tidal tail is approximately 5" wide (3 kpc projected), starting 18" (12 kpc) from the center of Pictor A, and extends more than 90" (60 kpc). The knots are part of a jet observed to be about 4' (160 kpc) long, extending to a bright hotspot. These images are the first optical detections of this jet, and by extracting knot flux densities through three filters we set constraints on emission models. While the radio and optical flux densities are usually explained by synchrotron emission, there are several emission mechanisms which might be used to explain the X-ray flux densities. Our data rule out Doppler boosted inverse Compton scattering as a source of the high energy emission. Instead, we find that the observed emission can be well described by synchrotron emission from electrons with a low energy index ($p\sim2$) that dominates the radio band, while a high energy index ($p\sim3$) is needed for the X-ray band and the transition occurs in the optical/infrared band. This model is consistent with a continuous electron injection scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

10 Simple Rules for the Care and Feeding of Scientific Data

This article offers a short guide to the steps scientists can take to ensure that their data and associated analyses continue to be of value and to be recognized. In just the past few years, hundreds of scholarly papers and reports have been written on questions of data sharing, data provenance, research reproducibility, licensing, attribution, privacy, and more, but our goal here is not to review that literature. Instead, we present a short guide intended for researchers who want to know why it is important to "care for and feed" data, with some practical advice on how to do that.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets

We present the use of continuous-time autoregressive moving average (CARMA) models as a method for estimating the variability features of a light curve, and in particular its power spectral density (PSD). CARMA models fully account for irregular sampling and measurement errors, making them valuable for quantifying variability, forecasting and interpolating light curves, and for variability-based classification. We show that the PSD of a CARMA model can be expressed as a sum of Lorentzian functions, which makes them extremely flexible and able to model a broad range of PSDs. We present the likelihood function for light curves sampled from CARMA processes, placing them on a statistically rigorous foundation, and we present a Bayesian method to infer the probability distribution of the PSD given the measured lightcurve. Because calculation of the likelihood function scales linearly with the number of data points, CARMA modeling scales to current and future massive time-domain data sets. We conclude by applying our CARMA modeling approach to light curves for an X-ray binary, two AGN, a long-period variable star, and an RR-Lyrae star, in order to illustrate their use, applicability, and interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stochastic Modeling of the Fermi/LAT Gamma-ray Blazar Variability

We study the gamma-ray variability of 13 blazars observed with the Fermi Large Area Telescope (LAT). These blazars have the most complete light curves collected during the first 4 years of the Fermi sky survey. We model them with the Ornstein-Uhlenbeck (OU) process or a mixture of the OU processes. The OU process has power spectral density (PSD) proportional to 1/f^alpha with alpha changing at a characteristic time scale, tau_0, from 0 (tau>>tau_0) to 2 (tau<<tau_0). The PSD of the mixed OU process has two characteristic time scales and an additional intermediate region with 0<alpha<2. We show that the OU model provides a good description of the Fermi/LAT light curves of three blazars in our sample. For the first time we constrain a characteristic gamma-ray time scale of variability in two BL Lac sources, 3C 66A and PKS 2155-304 (tau_0=25 day and tau_0=43 day, respectively, in the observer's frame), which are longer than the soft X-ray time scales detected in blazars and Seyfert galaxies. We find that the mixed OU process approximates the light curves of the remaining 10 blazars better than the OU process. We derive limits on their long and short characteristic time scales, and infer that their Fermi/LAT PSDs resemble a power-law function. We constrain the PSD slopes for all but one source in the sample. We find hints for sub-hour Fermi/LAT variability in four flat spectrum radio quasars. We discuss the implications of our results for theoretical models of blazar variability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strong Lens Time Delay Challenge: II. Results of TDC1

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $χ^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $χ^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Jet Emission in Young Radio Sources: a Fermi-LAT Gamma-ray View

We investigate the contribution of the beamed jet component to the high energy emission in young and compact extragalactic radio sources, focusing for the first time on the gamma-ray band. We derive predictions on the gamma-ray luminosities associated with the relativistic jet assuming a leptonic radiative model. The high energy emission is produced via Compton scattering by the relativistic electrons in a spherical region at the considered scales ($\lesssim$10 kpc). Simulations show a wide range of gamma-ray luminosities, with intensities up to $\sim10^{46}-10^{48}$ erg s$^{-1}$ depending on the assumed jet parameters. We find a highly linear relation between the simulated X-ray and gamma-ray luminosities that can be used to select candidates for a gamma-ray detection. We compare the simulated luminosity distributions in the radio, X-ray and gamma-ray regimes with observations for the largest sample of X-ray detected young radio quasars. Our analysis of $\sim$4-year Fermi Large Area Telescope (LAT) data does not give any statistically significant detection. However, the majority of the model-predicted gamma-ray fluxes for the sample are near or below the current Fermi-LAT flux threshold and compatible with the derived upper limits. Our study gives constraints on the minimum jet power ($L_{jet,kin}/L_{disk}>0.01$), on a potential jet contribution to the X-ray emission in the most compact sources ($\lesssim1$ kpc) and on the particles to magnetic field energy density ratio in broad agreement with equipartition assumption.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Deep Chandra X-ray Imaging of a Nearby Radio Galaxy 4C+29.30: X-ray/Radio Connection

We report results from our deep Chandra X-ray observations of a nearby radio galaxy, 4C+29.30 (z=0.0647). The Chandra image resolves structures on sub-arcsec to arcsec scales, revealing complex X-ray morphology and detecting the main radio features: the nucleus, a jet, hotspots, and lobes. The nucleus is absorbed (N(H)=3.95 (+0.27/-0.33)x10^23 atoms/cm^2) with an unabsorbed luminosity of L(2-10 keV) ~ (5.08 +/-0.52) 10^43 erg/s characteristic of Type 2 AGN. Regions of soft (<2 keV) X-ray emission that trace the hot interstellar medium (ISM) are correlated with radio structures along the main radio axis indicating a strong relation between the two. The X-ray emission beyond the radio source correlates with the morphology of optical line-emitting regions. We measured the ISM temperature in several regions across the galaxy to be kT ~ 0.5 with slightly higher temperatures (of a few keV) in the center and in the vicinity of the radio hotspots. Assuming these regions were heated by weak shocks driven by the expanding radio source, we estimated the corresponding Mach number of 1.6 in the southern regions. The thermal pressure of the X-ray emitting gas in the outermost regions suggest the hot ISM is slightly under-pressured with respect to the cold optical-line emitting gas and radio-emitting plasma, which both seem to be in a rough pressure equilibrium. We conclude that 4C+29.30 displays a complex view of interactions between the jet-driven radio outflow and host galaxy environment, signaling feedback processes closely associated with the central active nucleus.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Models of Accretion Disks

An accretion flow onto a supermassive black hole is the primary process powering quasars. However, a geometry of this flow is not well constrained. Both global MHD simulations and observations suggest that there are several emission components present in the nucleus: an accretion disk, hot plasma (corona or sphere) with electrons scattering the optical and UV photons, and an outflow (wind/jet). The relative location and size of these emission components, as well as their "interplay" affect the emerging quasar spectrum. I review briefly standard accretion disk models and the recent progress, point out discrepancies between the predicted and observed spectra and discuss some issues in fitting these models to the broad-band spectral energy distribution of quasars. I present examples of models fitted simultaneously to the optical-UV-X-ray data and possible constraints on the parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nuclear X-ray properties of the peculiar radio-loud hidden AGN 4C+29.30

We present results from a study of a nuclear emission of a nearby radio galaxy, 4C+29.30, over a broad 0.5-200 keV X-ray band. This study used new XMM-Newton (~17 ksec) and Chandra (~300 ksec) data, and archival Swift/BAT data from the 58-month catalog. The hard (>2 keV) X-ray spectrum of 4C+29.30 can be decomposed into an intrinsic hard power-law (Gamma ~ 1.56) modified by a cold absorber with an intrinsic column density N_{H,z} ~ 5x10^{23} cm^{-2}, and its reflection (|Omega/2pi| ~ 0.3) from a neutral matter including a narrow iron Kalpha emission line at the rest frame energy ~6.4 keV. The reflected component is less absorbed than the intrinsic one with an upper limit on the absorbing column of N^{refl}_{H,z} < 2.5x10^{22} cm^{-2}. The X-ray spectrum varied between the XMM-Newton and Chandra observations. We show that a scenario invoking variations of the normalization of the power-law is favored over a model with variable intrinsic column density. X-rays in the 0.5-2 keV band are dominated by diffuse emission modeled with a thermal bremsstrahlung component with temperature ~0.7 keV, and contain only a marginal contribution from the scattered power-law component. We hypothesize that 4C+29.30 belongs to a class of `hidden' AGN containing a geometrically thick torus. However, unlike the majority of them, 4C+29.30 is radio-loud. Correlations between the scattering fraction and Eddington luminosity ratio, and the one between black hole mass and stellar velocity dispersion, imply that 4C+29.30 hosts a black hole with ~10^8 M_{Sun} mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Stochastic Model for the Luminosity Fluctuations of Accreting Black Holes

In this work we have developed a new stochastic model for the fluctuations in lightcurves of accreting black holes. The model is based on a linear combination of stochastic processes and is also the solution to the linear diffusion equation perturbed by a spatially correlated noise field. This allows flexible modeling of the power spectral density (PSD), and we derive the likelihood function for the process, enabling one to estimate the parameters of the process, including break frequencies in the PSD. Our statistical technique is computationally efficient, unbiased by aliasing and red noise leak, and fully accounts for irregular sampling and measurement errors. We show that our stochastic model provides a good approximation to the X-ray lightcurves of galactic black holes, and the optical and X-ray lightcurves of AGN. We use the estimated time scales of our stochastic model to recover the correlation between characteristic time scale of the high frequency X-ray fluctuations and black hole mass for AGN, including two new `detections' of the time scale for Fairall 9 and NGC 5548. We find a tight anti-correlation between the black hole mass and the amplitude of the driving noise field, which is proportional to the amplitude of the high frequency X-ray PSD, and we estimate that this parameter gives black hole mass estimates to within ~ 0.2 dex precision, potentially the most accurate method for AGN yet. We also find evidence that ~ 13% of AGN optical PSDs fall off flatter than 1 / f^2, and, similar to previous work, find that the optical fluctuations are more suppressed on short time scales compared to the X-rays, but are larger on long time scales, suggesting the optical fluctuations are not solely due to reprocessing of X-rays.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accounting for Calibration Uncertainties in X-ray Analysis: Effective Areas in Spectral Fitting

While considerable advance has been made to account for statistical uncertainties in astronomical analyses, systematic instrumental uncertainties have been generally ignored. This can be crucial to a proper interpretation of analysis results because instrumental calibration uncertainty is a form of systematic uncertainty. Ignoring it can underestimate error bars and introduce bias into the fitted values of model parameters. Accounting for such uncertainties currently requires extensive case-specific simulations if using existing analysis packages. Here we present general statistical methods that incorporate calibration uncertainties into spectral analysis of high-energy data. We first present a method based on multiple imputation that can be applied with any fitting method, but is necessarily approximate. We then describe a more exact Bayesian approach that works in conjunction with a Markov chain Monte Carlo based fitting. We explore methods for improving computational efficiency, and in particular detail a method of summarizing calibration uncertainties with a principal component analysis of samples of plausible calibration files. This method is implemented using recently codified Chandra effective area uncertainties for low-resolution spectral analysis and is verified using both simulated and actual Chandra data. Our procedure for incorporating effective area uncertainty is easily generalized to other types of calibration uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statistical Characterization of the Chandra Source Catalog

The first release of the Chandra Source Catalog (CSC) contains ~95,000 X-ray sources in a total area of ~0.75% of the entire sky, using data from ~3,900 separate ACIS observations of a multitude of different types of X-ray sources. In order to maximize the scientific benefit of such a large, heterogeneous data-set, careful characterization of the statistical properties of the catalog, i.e., completeness, sensitivity, false source rate, and accuracy of source properties, is required. Characterization efforts of other, large Chandra catalogs, such as the ChaMP Point Source Catalog (Kim et al. 2007) or the 2 Mega-second Deep Field Surveys (Alexander et al. 2003), while informative, cannot serve this purpose, since the CSC analysis procedures are significantly different and the range of allowable data is much less restrictive. We describe here the characterization process for the CSC. This process includes both a comparison of real CSC results with those of other, deeper Chandra catalogs of the same targets and extensive simulations of blank-sky and point source populations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraints on Black Hole Growth, Quasar Lifetimes, and Eddington Ratio Distributions from the SDSS Broad Line Quasar Black Hole Mass Function

We present an estimate of the black hole mass function (BHMF) of broad line quasars (BLQSOs) that self-consistently corrects for incompleteness and the statistical uncertainty in the mass estimates, based on a sample of 9886 quasars at 1 < z < 4.5 drawn from the Sloan Digital Sky Survey. We find evidence for `cosmic downsizing' of black holes in BLQSOs, where the peak in their number density shifts to higher redshift with increasing black hole mass. The cosmic mass density for black holes seen as BLQSOs peaks at z ~ 2. We estimate the completeness of the SDSS as a function of black hole mass and Eddington ratio, and find that at z > 1 it is highly incomplete at M_BH < 10^9 M_Sun and L / L_Edd < 0.5. We also estimate a lower limit on the lifetime of a single BLQSO phase and we place constraints on the maximum mass of a black hole in a BLQSO. Our estimated distribution of BLQSO Eddington ratios peaks at L / L_Edd ~ 0.05 and has a dispersion of ~ 0.4 dex, implying that most BLQSOs are not radiating at or near the Eddington limit; however the location of the peak is subject to considerable uncertainty. The steep increase in number density of BLQSOs toward lower Eddington ratios is expected if the BLQSO accretion rate monotonically decays with time. Furthermore, our estimated lifetime and Eddington ratio distributions imply that the majority of the most massive black holes spend a significant amount of time growing in an earlier obscured phase, a conclusion which is independent of the unknown obscured fraction. These results are consistent with models for self-regulated black hole growth, at least for massive systems at z > 1, where the BLQSO phase occurs at the end of a fueling event when black hole feedback unbinds the accreting gas, halting the accretion flow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High redshift X-ray cooling-core cluster associated with the luminous radio loud quasar 3C186

We present the first results from a new, deep (200ks) Chandra observation of the X-ray luminous galaxy cluster surrounding the powerful (L ~10^47 erg/s), high-redshift (z=1.067), compact-steep-spectrum radio-loud quasar 3C186. The diffuse X-ray emission from the cluster has a roughly ellipsoidal shape and extends out to radii of at least ~60 arcsec (~500 kpc). The centroid of the diffuse X-ray emission is offset by 0.68(+/-0.11) arcsec (5.5+/-0.9 kpc) from the position of the quasar. We measure a cluster mass within the radius at which the mean enclosed density is 2500 times the critical density, r_2500=283(+18/-13)kpc, of 1.02 (+0.21/-0.14)x10^14 M_sun. The gas mass fraction within this radius is f_gas=0.129(+0.015/-0.016). This value is consistent with measurements at lower redshifts and implies minimal evolution in the f_gas(z) relation for hot, massive clusters at 0<z<1.1. The measured metal abundance of 0.42(+0.08/-0.07) Solar is consistent with the abundance observed in other massive, high redshift clusters. The spatially-resolved temperature profile for the cluster shows a drop in temperature, from kT~8 keV to kT~3 keV, in its central regions that is characteristic of cooling core clusters. This is the first spectroscopic identification of a cooling core cluster at z>1. We measure cooling times for the X-ray emitting gas at radii of 50 kpc and 25 kpc of 1.7(+/-0.2)x10^9 years and 7.5(+/-2.6)x 10^8 years, as well as a nominal cooling rate (in the absence of heating) of 400(+/-190)M_sun/year within the central 100 kpc. In principle, the cooling gas can supply enough fuel to support the growth of the supermassive black hole and to power the luminous quasar. The radiative power of the quasar exceeds by a factor of 10 the kinematic power of the central radio source, suggesting that radiative heating may be important at intermittent intervals in cluster cores.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Intermittent Activity of Jets in AGN

Large scale X-ray jets that extend to >100 kpc distances from the host galaxy indicate the importance of jets interactions with the environment on many different physical scales. Morphology of X-ray clusters indicate that the radio-jet activity of a cD galaxy is intermittent. This intermittency might be a result of a feedback and/or interactions between galaxies within the cluster. Here we consider the radiation pressure instability operating on short timescales (<10^5 years) as the origin of the intermittent behaviour. We test whether this instability can be responsible for short ages (< 10^4 years) of Compact Symmetric Objects measured by hot spots propagation velocities in VLBI observations. We model the accretion disk evolution and constrain model parameters that may explain the observed compact radio structures and over-abundance of GPS sources. We also describe effects of consequent outbursts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On Computing Upper Limits to Source Intensities

A common problem in astrophysics is determining how bright a source could be and still not be detected. Despite the simplicity with which the problem can be stated, the solution involves complex statistical issues that require careful analysis. In contrast to the confidence bound, this concept has never been formally analyzed, leading to a great variety of often ad hoc solutions. Here we formulate and describe the problem in a self-consistent manner. Detection significance is usually defined by the acceptable proportion of false positives (the TypeI error), and we invoke the complementary concept of false negatives (the TypeII error), based on the statistical power of a test, to compute an upper limit to the detectable source intensity. To determine the minimum intensity that a source must have for it to be detected, we first define a detection threshold, and then compute the probabilities of detecting sources of various intensities at the given threshold. The intensity that corresponds to the specified TypeII error probability defines that minimum intensity, and is identified as the upper limit. Thus, an upper limit is a characteristic of the detection procedure rather than the strength of any particular source and should not be confused with confidence intervals or other estimates of source intensity. This is particularly important given the large number of catalogs that are being generated from increasingly sensitive surveys. We discuss the differences between these upper limits and confidence bounds. Both measures are useful quantities that should be reported in order to extract the most science from catalogs, though they answer different statistical questions: an upper bound describes an inference range on the source intensity, while an upper limit calibrates the detection process. We provide a recipe for computing upper limits that applies to all detection algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Probing the accretion disk - jet connection via instabilities in the inner accretion flow. From microquasars to quasars

We present various instability mechanisms in the accreting black hole systems which might indicate at the connection between the accretion disk and jet. The jets observed in microquasars can have a persistent or blobby morphology. Correlated with the accretion luminosity, this might provide a link to the cyclic outbursts of the disk. Such duty-cycle type of behavior on short timescales results from the thermal instability caused by the radiation pressure domination. The same type of instability may explain the cyclic radioactivity of the supermassive black hole systems. The somewhat longer timescales are characteristic for the instability caused by the partial hydrogen ionization. The distortions of the jet direction and complex morphology of the sources can be caused by precession of the disk-jet axis.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Submillimeter Variability and the Gamma-ray Connection in Fermi Blazars

We present multi-epoch observations from the Submillimeter Array (SMA) for a sample of 171 bright blazars, 43 of which were detected by Fermi during the first three months of observations. We explore the correlation between their gamma-ray properties and submillimeter observations of their parsec-scale jets, with a special emphasis on spectral index in both bands and the variability of the synchrotron component. Subclass is determined using a combination of Fermi designation and the Candidate Gamma-Ray Blazar Survey (CGRaBS), resulting in 35 BL Lac objects and 136 flat-spectrum radio quasars (FSRQs) in our total sample. We calculate submillimeter energy spectral indices using contemporaneous observations in the 1 mm and 850 micron bands during the months August-October 2008. The submillimeter light curves are modeled as first-order continuous autoregressive processes, from which we derive characteristic timescales. Our blazar sample exhibits no differences in submillimeter variability amplitude or characteristic timescale as a function of subclass or luminosity. All of the the light curves are consistent with being produced by a single process that accounts for both low and high states, and there is additional evidence that objects may be transitioning between blazar class during flaring epochs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Chandra Source Catalog

The Chandra Source Catalog (CSC) is a general purpose virtual X-ray astrophysics facility that provides access to a carefully selected set of generally useful quantities for individual X-ray sources, and is designed to satisfy the needs of a broad-based group of scientists, including those who may be less familiar with astronomical data analysis in the X-ray regime. The first release of the CSC includes information about 94,676 distinct X-ray sources detected in a subset of public ACIS imaging observations from roughly the first eight years of the Chandra mission. This release of the catalog includes point and compact sources with observed spatial extents <~ 30''. The catalog (1) provides access to the best estimates of the X-ray source properties for detected sources, with good scientific fidelity, and directly supports scientific analysis using the individual source data; (2) facilitates analysis of a wide range of statistical properties for classes of X-ray sources; and (3) provides efficient access to calibrated observational data and ancillary data products for individual X-ray sources, so that users can perform detailed further analysis using existing tools. The catalog includes real X-ray sources detected with flux estimates that are at least 3 times their estimated 1 sigma uncertainties in at least one energy band, while maintaining the number of spurious sources at a level of <~ 1 false source per field for a 100 ks observation. For each detected source, the CSC provides commonly tabulated quantities, including source position, extent, multi-band fluxes, hardness ratios, and variability statistics, derived from the observations in which the source is detected. In addition to these traditional catalog elements, for each X-ray source the CSC includes an extensive set of file-based data products that can be manipulated interactively.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Accretion Disk Model of Short-Timescale Intermittent Activity in Young Radio Sources

We associate the existence of short-lived compact radio sources with the intermittent activity of the central engine caused by a radiation pressure instability within an accretion disk. Such objects may constitute a numerous sub-class of Giga-Hertz Peaked Spectrum sources, in accordance with the population studies of radio-loud active galaxies, as well as detailed investigations of their radio morphologies. We perform the model computations assuming the viscosity parametrization as proportional to a geometrical mean of the total and gas pressure. The implied timescales are consistent with the observed ages of the sources. The duration of an active phase for a moderate accretion rate is short enough (< 10^3-10^4 years) that the ejecta are confined within the host galaxy and thus these sources cannot evolve into large size radio galaxies unless they are close to the Eddington limit.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Are the Variations in Quasar Optical Flux Driven by Thermal Fluctuations?

We analyze a sample of optical light curves for 100 quasars, 70 of which have black hole mass estimates. Our sample is the largest and broadest used yet for modeling quasar variability. The sources in our sample have z < 2.8 and 10^6 < M_BH < 10^10. We model the light curves as a continuous time stochastic process, providing a natural means of estimating the characteristic time scale and amplitude of quasar variations. We employ a Bayesian approach to estimate the characteristic time scale and amplitude of flux variations; our approach is not affected by biases introduced from discrete sampling effects. We find that the characteristic time scales strongly correlate with black hole mass and luminosity, and are consistent with disk orbital or thermal time scales. In addition, the amplitude of short time scale variations is significantly anti-correlated with black hole mass and luminosity. We interpret the optical flux fluctuations as resulting from thermal fluctuations that are driven by an underlying stochastic process, such as a turbulent magnetic field. In addition, the intranight variations in optical flux implied by our empirical model are < 0.02 mag, consistent with current microvariability observations of radio-quiet quasars. Our stochastic model is therefore able to unify both long and short time scale optical variations in radio-quiet quasars as resulting from the same underlying process, while radio-loud quasars have an additional variability component that operates on time scales < 1 day.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Bayesian Estimation of Hardness Ratios: Modeling and Computations

A commonly used measure to summarize the nature of a photon spectrum is the so-called Hardness Ratio, which compares the number of counts observed in different passbands. The hardness ratio is especially useful to distinguish between and categorize weak sources as a proxy for detailed spectral fitting. However, in this regime classical methods of error propagation fail, and the estimates of spectral hardness become unreliable. Here we develop a rigorous statistical treatment of hardness ratios that properly deals with detected photons as independent Poisson random variables and correctly deals with the non-Gaussian nature of the error propagation. The method is Bayesian in nature, and thus can be generalized to carry out a multitude of source-population--based analyses. We verify our method with simulation studies, and compare it with the classical method. We apply this method to real world examples, such as the identification of candidate quiescent Low-mass X-ray binaries in globular clusters, and tracking the time evolution of a flare on a low-mass star.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Chandra Discovery of a 300 kpc X-ray Jet in the GPS Quasar PKS1127-145

We have discovered an X-ray jet with Chandra imaging of the z=1.187 radio-loud quasar PKS1127-145. In this paper we present the Chandra X-ray data, follow-up VLA observations, and optical imaging using the HST WFPC2. The X-ray jet contains 273+/-5 net counts in 27ksec and extends ~30 arcsec, from the quasar core, corresponding to a minimum projected linear size of ~330/h_50 kpc. The evaluation of the X-ray emission processes is complicated by the observed offsets between X-ray and radio brightness peaks. We discuss the problems posed by these observations to jet models. In addition, PKS1127-145 is a Giga-Hertz Peaked Spectrum radio source, a member of the class of radio sources suspected to be young or ``frustrated'' versions of FRI radio galaxies. However the discovery of an X-ray and radio jet extending well outside the host galaxy of PKS1127-145 suggests that activity in this and other GPS sources may be long-lived and complex.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Statistics: Handle with Care, Detecting Multiple Model Components with the Likelihood Ratio Test

The likelihood ratio test (LRT) and the related $F$ test, do not (even asymptotically) adhere to their nominal $χ^2$ and $F$ distributions in many statistical tests common in astrophysics, thereby casting many marginal line or source detections and non-detections into doubt. Although there are many legitimate uses of these statistics, in some important cases it can be impossible to compute the correct false positive rate. For example, it has become common practice to use the LRT or the $F$ test for detecting a line in a spectral model or a source above background despite the lack of certain required regularity conditions. In these and other settings that involve testing a hypothesis that is on the boundary of the parameter space, {\it contrary to common practice, the nominal $χ^2$ distribution for the LRT or the $F$ distribution for the $F$ test should not be used}. In this paper, we characterize an important class of problems where the LRT and the $F$ test fail and illustrate this non-standard behavior. We briefly sketch several possible acceptable alternatives, focusing on Bayesian posterior predictive probability-values. We present this method in some detail, as it is a simple, robust, and intuitive approach. This alternative method is illustrated using the gamma-ray burst of May 8, 1997 (GRB 970508) to investigate the presence of an Fe K emission line during the initial phase of the observation.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Chandra Detection of X-ray Absorption Associated with a Damped Lyman Alpha System

We have observed three quasars, PKS 1127-145, Q 1331+171 and Q0054+144, with the ACIS-S aboard the Chandra X-ray Observatory, in order to measure soft X-ray absorption associated with intervening 21-cm and damped Ly$α$ absorbers. For PKS 1127-145, we detect absorption which, if associated with an intervening z_{abs}=0.312 absorber, implies a metallicity of 23% solar. If the absorption is not at z_{abs}=0.312, then the metallicity is still constrained to be less than 23% solar. The advantage of the X-ray measurement is that the derived metallicity is insensitive to ionization, inclusion of an atom in a molecule, or depletion onto grains. The X-ray absorption is mostly due to oxygen, and is consistent with the oxygen abundance of 30% solar derived from optical nebular emission lines in a foreground galaxy at the redshift of the absorber. For Q1331+171 and Q 0054+144, only upper limits were obtained, although the exposure times were intentionally short, since for these two objects we were interested primarily in measuring flux levels to plan for future observations. The imaging results are presented in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Analysis of Energy Spectra with Low Photon Counts via Bayesian Posterior Simulation

Over the past 10 years Bayesian methods have rapidly grown more popular as several computationally intensive statistical algorithms have become feasible with increased computer power. In this paper, we begin with a general description of the Bayesian paradigm for statistical inference and the various state-of-the-art model fitting techniques that we employ (e.g., Gibbs sampler and Metropolis- Hastings). These algorithms are very flexible and can be used to fit models that account for the highly hierarchical structure inherent in the collection of high-quality spectra and thus can keep pace with the accelerating progress of new space telescope designs. The methods we develop, which will soon be available in the CIAO software package, explicitly model photon arrivals as a Poisson process and, thus, have no difficulty with high resolution low count X-ray and gamma-ray data. We expect these methods to be useful not only for the recently launched Chandra X-ray observatory and XMM but also new generation telescopes such as Constellation X, GLAST, etc. In the context of two examples (Quasar S5 0014+813 and Hybrid-Chromosphere Supergiant Star alpha TrA) we illustrate a new highly structured model and how Bayesian posterior sampling can be used to compute estimates, error bars, and credible intervals for the various model parameters.

0 citations0 reviewsNo code linkedOpen access