Engineering Papers⌕ Search

Engineering topics

Siemiginowska, Aneta

Publications and source records attributed to Siemiginowska, Aneta.

At least 19 records

Algorithms and Statistical Models for Scientific Discovery in the Petabyte Era

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Correlation between X-Ray and Radio Absorption in Compact Radio Galaxies

Compact radio galaxies with a GHz-peaked spectrum (GPS) and/or compact-symmetric-object (CSO) morphology (GPS/CSOs) are increasingly detected in the X-ray domain. Their radio and X-ray emissions are affected by significant absorption. However, the locations of the X-ray and radio absorbers are still debated. We investigated the relationship between the column densities of the total (N{sub H}) and neutral (N{sub HI}) hydrogen to statistically constrain the picture. We compiled a sample of GPS/CSOs including both literature data and new radio data that we acquired with the Westerbork Synthesis Radio Telescope for sources whose X-ray emission was either established or under investigation. In this sample, we compared the X-ray and radio hydrogen column densities, and found that N{sub H} and N{sub HI} display a significant positive correlation with N{sub HI} ∝ N{sub H} {sup b}, where b = 0.47 and b = 0.35, depending on the subsample. The N{sub H}–N{sub HI} correlation suggests that the X-ray and radio absorbers are either co-spatial or different components of a continuous structure. The correlation displays a large intrinsic spread that we suggest to originate from fluctuations, around a mean value, of the ratio between the spin temperature and the covering factor of the radio absorber, T{sub s}/C{sub f}.

79 ASTRONOMY AND ASTROPHYSICS↗

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 data set 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 six marginally. We also find residual emission in the cores of three quasars and in the background of one 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 are due to the IC/CMB rather than the synchrotron process, or that high-redshift jets are qualitatively different.

79 ASTRONOMY AND ASTROPHYSICS↗

X-RAY PROPERTIES OF THE YOUNGEST RADIO SOURCES AND THEIR ENVIRONMENTS

We present the first results from our X-ray study of young radio sources classified as compact symmetric objects (CSOs). Using the Chandra X-ray Observatory we observed six CSOs for the first time in X-rays, and re-observed four CSOs already observed with XMM-Newton or BeppoSAX. We also included six other CSOs with archival data to built a pilot study of a sample of the 16 CSO sources observed in X-rays to date. All the sources are nearby, z<1, and the age of their radio structures (<3000 yr) has been derived from the expansion velocity of their hot spots. Our results show the heterogeneous nature of the CSOs’ X-ray emission, indicating a complex environment associated with young radio sources. The sample covers a range in X-ray luminosity, L{sub 2--10keV}∼10{sup 41}–10{sup 45} erg s{sup −1}, and intrinsic absorbing column density of N{sub H}≃10{sup 21}–10{sup 22} cm{sup −2}. In particular, we detected extended X-ray emission in 1718−649; a hard photon index of Γ≃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{sub rest}=(6.62±0.04) keV and EW ∼0.15–1.4 keV, and a decrease by an order of magnitude in the 2–10 keV flux since the 2008 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 to study the interactions of young radio sources with the interstellar environment of their host galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

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.

79 ASTRONOMY AND ASTROPHYSICS↗

Optical detection of the Pictor A jet and tidal tail: Evidence against an IC/CMB JET

New images of the FR II radio galaxy Pictor A from the Hubble Space Telescope 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 that 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∼2) that dominates the radio band, while a high-energy index (p∼3) 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.

79 ASTRONOMY AND ASTROPHYSICS↗

A fully Bayesian method for jointly fitting instrumental calibration and X-ray spectral models

Owing to a lack of robust principled methods, systematic instrumental uncertainties have generally been ignored in astrophysical data analysis despite wide recognition of the importance of including them. Ignoring calibration uncertainty can cause bias in the estimation of source model parameters and can lead to underestimation of the variance of these estimates. We previously introduced a pragmatic Bayesian method to address this problem. The method is 'pragmatic' in that it introduced an ad hoc technique that simplified computation by neglecting the potential information in the data for narrowing the uncertainty for the calibration product. Following that work, we use a principal component analysis to efficiently represent the uncertainty of the effective area of an X-ray (or γ-ray) telescope. Here, however, we leverage this representation to enable a principled, fully Bayesian method that coherently accounts for the calibration uncertainty in high-energy spectral analysis. In this setting, the method is compared with standard analysis techniques and the pragmatic Bayesian method. The advantage of the fully Bayesian method is that it allows the data to provide information not only for estimation of the source parameters but also for the calibration product—here the effective area, conditional on the adopted spectral model. In this way, it can yield more accurate and efficient estimates of the source parameters along with valid estimates of their uncertainty. Provided that the source spectrum can be accurately described by a parameterized model, this method allows rigorous inference about the effective area by quantifying which possible curves are most consistent with the data.

79 ASTRONOMY AND ASTROPHYSICS↗

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 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 light curve. 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 active galactic nuclei, a long-period variable star, and an RR Lyrae star in order to illustrate their use, applicability, and interpretation.

79 ASTRONOMY AND ASTROPHYSICS↗

NUCLEAR X-RAY PROPERTIES OF THE PECULIAR RADIO-LOUD HIDDEN AGN 4C+29.30

We present results from a study of nuclear emission from a nearby radio galaxy, 4C+29.30, over a broad 0.5-200 keV X-ray band. This study used new XMM-Newton ({approx}17 ks) and Chandra ({approx}300 ks) 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} {approx} 1.56) modified by a cold absorber with an intrinsic column density N {sub H,z} {approx} 5 Multiplication-Sign 10{sup 23} cm{sup -2}, and its reflection (|{Omega}/2{pi}| {approx} 0.3) from a neutral matter including a narrow iron K{alpha} emission line at a rest-frame energy {approx}6.4 keV. The reflected component is less absorbed than the intrinsic one with an upper limit on the absorbing column of N {sup refl} {sub H,z} < 2.5 Multiplication-Sign 10{sup 22} cm{sup -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 {approx}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' active galactic nuclei containing a geometrically thick torus. However, unlike the majority of hidden AGNs, 4C+29.30 is radio-loud. Correlations between the scattering fraction and Eddington luminosity ratio, and between black hole mass and stellar velocity dispersion, imply that 4C+29.30 hosts a black hole with {approx}10{sup 8} M {sub Sun} mass.

79 ASTRONOMY AND ASTROPHYSICS↗

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{sub H} {approx_equal} 3.95{sup +0.27}{sub -0.33} Multiplication-Sign 10{sup 23} cm{sup -2}) with an unabsorbed luminosity of L{sub 2-10keV} {approx_equal} (5.08 {+-} 0.52) Multiplication-Sign 10{sup 43} erg s{sup -1} characteristic of Type 2 active galactic nuclei. 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 extends beyond the radio source and correlates with the morphology of optical-line-emitting regions. We measured the ISM temperature in several regions across the galaxy to be kT {approx_equal} 0.5 keV, with slightly higher temperatures (of a few keV) in the center and in the vicinity of the radio hotspots. Assuming that 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 suggests that the hot ISM is slightly underpressured 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.

79 ASTRONOMY AND ASTROPHYSICS↗

BROADBAND JET EMISSION IN YOUNG AND POWERFUL RADIO SOURCES: THE CASE OF THE COMPACT STEEP SPECTRUM QUASAR 3C 186

We present the X-ray analysis of a deep ({approx}200 ks) Chandra observation of the compact steep spectrum radio-loud quasar 3C 186 (z = 1.06) and investigate the contribution of the unresolved radio jet to the total X-ray emission. The spectral analysis is not conclusive on the origin of the bulk of the X-ray emission. In order to examine the jet contribution to the X-ray flux, we model the quasar spectral energy distribution, adopting several scenarios for the jet emission. For the values of the main physical parameters favored by the observables, a dominant role of the jet emission in the X-ray band is ruled out when a single-zone (leptonic) scenario is adopted, even including the contribution of the external photon fields as seed photons for inverse Compton emission. We then consider a structured jet, with the blazar component that-although not directly visible in the X-ray band-provides an intense field of seed synchrotron photons Compton-scattered by electrons in a mildly relativistic knot. In this case, the whole X-ray emission can be accounted for if we assume a blazar luminosity within the range observed from flat spectrum radio quasars. The X-ray radiative efficiency of such a (structured) jet is intimately related to the presence of a complex velocity structure. The jet emission can provide a significant contribution in X-rays if it decelerates within the host galaxy on kiloparsec scales. We discuss the implications of this model in terms of jet dynamics and interaction with the ambient medium.

79 ASTRONOMY AND ASTROPHYSICS↗

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.

79 ASTRONOMY AND ASTROPHYSICS↗

Handbook of X-Ray Astronomy

X-ray astronomy was born in the aftermath of World War II as military rockets were repurposed to lift radiation detectors above the atmosphere for a few minutes at a time. These early flights detected and studied X-ray emission from the Solar corona. The first sources beyond the Solar System were detected during a rocket flight in 1962 by a team headed by Riccardo Giaccom at American Science and Engineering, a company founded by physicists from MIT. The rocket used Geiger counters with a system designed to reduce non-X-ray backgrounds and collimators limiting the region of sky seen by the counters. As the rocket spun, the field of view (FOV) happened to pass over what was later found to be the brightest non-Solar X-ray source; later designated See X-1. It also detected a uniform background glow which could not be resolved into individual sources. A follow-up campaign using X-ray detectors with better spatial resolution and optical telescopes identified See X-1 as an interacting binary with a compact (neutron star) primary. This success led to further suborbital rocket flights by a number of groups. More X-ray binaries were discovered, as well as X-ray emission from supernova remnants, the radio galaxies M87 and Cygnus-A, and the Coma cluster. Detectors were improved and Geiger counters were replaced by proportional counters, which provided information about energy spectra of the sources. A constant challenge was determining precise positions of sources as only collimators were available.

Arnaud, Keith A.↗

HIGH-REDSHIFT X-RAY COOLING-CORE CLUSTER ASSOCIATED WITH THE LUMINOUS RADIO-LOUD QUASAR 3C 186

We present the first results from a new, deep (200 ks) Chandra observation of the X-ray luminous galaxy cluster surrounding the powerful (L {approx} 10{sup 47} erg s{sup -1}), high-redshift (z = 1.067), compact-steep-spectrum radio-loud quasar 3C 186. The diffuse X-ray emission from the cluster has a roughly ellipsoidal shape and extends out to radii of at least {approx}60 arcsec ({approx}500 kpc). The centroid of the diffuse X-ray emission is offset by 0.68 {+-} 0.''11 ({approx}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{sub 2500} = 283{sup +18}{sub -13} kpc, of 1.02{sup +0.21}{sub -0.14} x 10{sup 14} M{sub sun}. The gas-mass fraction within this radius is f{sub gas} = 0.129{sup +0.015}{sub -0.016}. This value is consistent with measurements at lower redshifts and implies minimal evolution in the f{sub gas}(z) relation for hot, massive clusters at 0 < z < 1.1. The measured metal abundance of 0.42{sup +0.08}{sub -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 {approx} 8 keV to kT {approx} 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 x 10{sup 9} years and 7.5 {+-} 2.6 x 10{sup 8} years, as well as a nominal cooling rate (in the absence of heating) of 400 {+-} 190 M{sub sun} year{sup -1} 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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 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 (SDSS). 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 {approx} 2. We estimate the completeness of the SDSS as a function of the black hole mass and Eddington ratio, and find that at z > 1 it is highly incomplete at M {sub BH} {approx}< 10{sup 9} M {sub sun} and L/L{sub Edd} {approx}< 0.5. We estimate a lower limit on the lifetime of a single BLQSO phase to be t {sub BL} > 150 {+-} 15 Myr for black holes at z = 1 with a mass of M {sub BH} = 10{sup 9} M{sub sun}, and we constrain the maximum mass of a black hole in a BLQSO to be {approx}3 x 10{sup 10} M{sub sun}. Our estimated distribution of BLQSO Eddington ratios peaks at L/L {sub Edd} {approx} 0.05 and has a dispersion of {approx}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.

79 ASTRONOMY AND ASTROPHYSICS↗

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 in an observation. Despite the simplicity with which the problem can be stated, the solution involves complicated statistical issues that require careful analysis. In contrast to the more familiar 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 (background fluctuations that are claimed as detections, or Type I error), and we invoke the complementary concept of false negatives (real sources that go undetected, or Type II 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 Type II 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. It 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, with examples, 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.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

ORIGIN OF THE COMPLEX RADIO STRUCTURE IN BAL QSO 1045+352

We present new, more sensitive, high-resolution radio observations of a compact broad absorption line (BAL) quasar, 1045+352, made with the EVN+MERLIN at 5 GHz. These observations allowed us to trace the connection between the arcsecond structure and the radio core of the quasar. The radio morphology of 1045+352 is dominated by a knotty jet showing several bends. We discuss possible scenarios that could explain such a complex morphology: galaxy merger, accretion disk instability, precession of the jet, and jet-cloud interactions. It is possible that we are witnessing an ongoing jet precession in this source due to internal instabilities within the jet flow; however, a dense environment detected in the submillimeter band and an outflowing material suggested by the X-ray absorption could strongly interact with the jet. It is difficult to establish the orientation between the jet axis and the observer in 1045+352 because of the complex structure. Nevertheless, taking into account the most recent inner radio structure, we conclude that the radio jet is oriented close to the line of sight, which can mean that the opening angle of the accretion disk wind can be large in this source. We also suggest that there is no direct correlation between the jet-observer orientation and the possibility of observing BALs.

79 ASTRONOMY AND ASTROPHYSICS↗

The X-Ray Weakness of GPS Radio Galaxies: A Volume-Limited Complete Sample

The XMM-Newton observations of the GPS sample was completed last summer. We are in process of finalizing the paper describing the data and the results. The main goal of the project was to determine the X-ray spectra of the GPS galaxies in comparison to regular radio galaxies. Our XMM observations show evidence that the GPS galaxies are heavily obscured with the large absorbing columns exceeding N(H)greater than le22 cm^-2. Taking into account the obscuration we determined that the intrinsic X-ray luminosities of GPS galaxies are of order le43-le44 erg/s, comparable to low luminosity radio loud quasars. The large GPS samples can confirm the result, as at this moment our evidence is based only on 7 GPS galaxies observed with good S/N/ in X-rays. The first paper summarizing the results of the XMM observation of Mkn 668 has been published Astronomy & Astrophysics. We found soft X-ray signatures of a hot plasma (kT approximately 10^7~K) and a hard X-ray emission from the nucleus. The X-ray spectrum above 2.5-keV is characterized by a very flat (observed photon index, Gamma-0.5) power-law continuum, alongside with a strong Fe-K-alpha neutral iron fluorescent line (EW~600-eV). The best explanation for the origin of this high energy X-ray emission is in terms of the Compton-reflection of the nuclear emission. The primary X-ray emission is obscured by a Compton-thick (N_H ~ to 10^24 ~cm-2) matter which becomes transparent at higher energies. The observed above 2.5-keV X-rays are mostly due to reflection which is indicated by a strong Fe-K-alpha line. This source represented the second hard X-ray detection of the GPS galaxy ever (the first one being 1345+125; O'Dea et al. 2000). The observations of the other GPS galaxies in our sample confirmed the trend of the large obscuration present in the spectra. However, we do not have a compelling evidence for a hot gas in the nucleus. The two other GPS galaxies observed with Chandra were added to the total of 7 GPS galaxies. This GPS sample was compared against the sample of radio loud galaxies and quasar. The paper is about to be submitted. Our preliminarily results based on the XMM-Newton data for the entire sample were presented at the meeting in July 2004.

Mushotsky, Richard↗