Source author record

William G. Mathews

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

28 published item(s)

preprint2015arXiv

Hot gaseous atmospheres in galaxy groups and clusters are both heated and cooled by X-ray cavities

Expanding X-ray cavities observed in hot gas atmospheres of many galaxy groups and clusters generate shock waves and turbulence that are primary heating mechanisms required to avoid uninhibited radiatively cooling flows which are not observed. However, we show here that the evolution of buoyant cavities also stimulates radiative cooling of observable masses of low-temperature gas. During their early evolution, radiative cooling occurs in the wakes of buoyant cavities in two locations: in thin radial filaments parallel to the buoyant velocity and more broadly in gas compressed beneath rising cavities. Radiation from these sustained compressions removes entropy from the hot gas. Gas experiencing the largest entropy loss cools first, followed by gas with progressively less entropy loss. Most cooling occurs at late times, $\sim 10^8-10^9$ yrs, long after the X-ray cavities have disrupted and are impossible to detect. During these late times, slightly denser low entropy gas sinks slowly toward the centers of the hot atmospheres where it cools intermittently, forming clouds near the cluster center. Single cavities of energy $10^{57}-10^{58}$ ergs in the atmosphere of the NGC 5044 group create $10^8 - 10^9$ $M_{\odot}$ of cooled gas, exceeding the mass of extended molecular gas currently observed in that group. The cooled gas clouds we compute share many attributes with molecular clouds recently observed in NGC 5044 with ALMA: self-gravitationally unbound, dust-free, quasi-randomly distributed within a few kpc around the group center.

preprint2014arXiv

First Dynamic Computation of Synchrotron Emission from the Cygnus A Radio Cavity; Evidence for Electron Pair Plasma in Cavity

Cosmic rays, thermal gas and magnetic fields in FRII radio cavities are assumed to come entirely from winds flowing from just behind the jet shocks. Combining analytic and computational methods, it is shown that the computed radio-electron energy distribution and synchrotron emissivity spectra everywhere in the Cygnus A radio cavity agrees with radio observations of the Cygnus A lobes. The magnetic field energy density is small everywhere and evolves passively in the post-shock wind. Most synchrotron emission arises in recent post-shock material as it flows back along the radio cavity wall. Because it experienced less adiabatic expansion, the magnetic field in this young backflow is larger than elsewhere in the radio lobe, explaining the observed radio synchrotron limb-brightening. The boundary backflow decelerates due to small cavity pressure gradients, causing large-scale fields perpendicular to the backflow (and synchrotron emission) to grow exponentially unlike observations. However, if the field is random on subgrid (sub-kpc) scales, the computed field reproduces both the magnitude and slowly decreasing radio synchrotron emissivity observed along the backflow. The radio synchrotron spectrum and image computed with a small-scale random field agree with VLA observations. The total relativistic energy density in the post-jet shock region required in computations to inflate the radio cavity matches the energy density of relativistic electrons observed in the post-shock region of Cygnus A. This indicates that the component in the jet and cavity that dominates the dynamical evolution is a relativistic pair plasma.

preprint2013arXiv

Hot versus Cold: the Dichotomy in Spherical Accretion of Cooling Flows onto Supermassive Black Holes in Elliptical Galaxies, Galaxy Groups and Clusters

Feedback heating from active galactic nuclei (AGNs) has been commonly invoked to suppress cooling flows predicted in hot gas in elliptical galaxies, galaxy groups and clusters. Previous studies have focused on if and how AGN feedback heats the gas, but little paid attention to its triggering mechanism. Using spherically symmetric simulations, we investigate how large-scale cooling flows are accreted by central supermassive black holes (SMBHs) in eight well-observed systems and find an interesting dichotomy. In massive clusters, the gas develops a central cooling catastrophe within about the cooling time (typically ~ 100 - 300 Myr), resulting in a cold-mode accretion onto SMBHs. However, in our four simulated systems on group and galaxy scales at a low metallicity Z=0.3Z_{sun}, the gas quickly settles into a long-term state which has a cuspy central temperature profile extending to several tens to about 100 pc. At the more realistic solar metallicity, two groups (with R_{e} ~ 4 kpc) still host the long-term hot-mode accretion. Both accretion modes naturally appear in our idealized calculations where only cooling, gas inflow, and compressional heating are considered. The long-term hot-mode accretion is maintained by the quickly-established closeness between the timescales of these processes, preferably in systems with low gas densities, low gas metallicities, and importantly, compact central galaxies, which result in strong gravitational acceleration and compressional heating at the intermediate radii. Our calculations predict that central cuspy temperature profiles appear more often in smaller systems than galaxy clusters, which instead often host significant cold gas and star formation.

preprint2013arXiv

Variation of Mid and Far-IR Luminosities among Early-Type Galaxies: Relation to Stellar Metallicity and Cold Dust

The Hubble morphological sequence from early to late galaxies corresponds to an increasing rate of specific star formation. The Hubble sequence also follows a banana-shaped correlation between 24 and 70 micron luminosities, both normalized with the K-band luminosity. We show that this correlation is significantly tightened if galaxies with central AGN emission are removed, but the cosmic scatter of elliptical galaxies in both 24 and 70 micron luminosities remains significant along the correlation. We find that the 24 micron variation among ellipticals correlates with stellar metallicity, reflecting emission from hot dust in winds from asymptotic giant branch stars of varying metallicity. Infrared surface brightness variations in elliptical galaxies indicate that the K - 24 color profile is U-shaped for reasons that are unclear. In some elliptical galaxies cold interstellar dust emitting at 70 and 160 microns may arise from recent gas-rich mergers. However, we argue that most of the large range of 70 micron luminosity in elliptical galaxies is due to dust transported from galactic cores by feedback events in (currently IR-quiet) active galactic nuclei. Cooler dusty gas naturally accumulates in the cores of elliptical galaxies due to dust-cooled local stellar mass loss and may accrete onto the central black hole, releasing energy. AGN-heated gas can transport dust in cores 5-10 kpc out into the hot gas atmospheres where it radiates extended 70 micron emission but is eventually destroyed by sputtering. This, and some modest star formation, defines a cycle of dust creation and destruction. Elliptical galaxies evidently undergo large transient excursions in the banana plot in times comparable to the sputtering time or AGN duty cycle, 10 Myrs. Normally regarded as passive, elliptical galaxies are the most active galaxies in the IR color-color correlation.

preprint2012arXiv

Dynamics Inside the Radio and X-ray Cluster Cavities of Cygnus A and Similar FRII Sources

We describe approximate axisymmetric computations of the dynamical evolution of material inside radio lobes and X-ray cluster gas cavities in Fanaroff-Riley II sources such as Cygnus A. All energy is delivered by a jet to the lobe/cavity via a moving hotspot where jet energy dissipates in a reverse shock. Our calculations describe the evolution of hot plasma, cosmic rays (CRs) and toroidal magnetic fields flowing from the hotspot into the cavity. Many observed features are explained. Gas, CRs and field flow back along the cavity surface in a "boundary backflow" consistent with detailed FRII observations. Computed ages of backflowing CRs are consistent with observed radio-synchrotron age variations only if shear instabilities in the boundary backflow are damped and we assume this is done with viscosity of unknown origin. Magnetic fields estimated from synchrotron self-Compton (SSC) X-radiation observed near the hotspot evolve into radio lobe fields. Computed profiles of radio synchrotron lobe emission perpendicular to the jet are dramatically limb-brightened in excellent agreement with FRII observations although computed lobe fields exceed those observed. Strong winds flowing from hotspots naturally create kpc-sized spatial offsets between hotspot inverse Compton (IC-CMB) X-ray emission and radio synchrotron emission that peaks 1-2 kpc ahead where the field increases due to wind compression. In our computed version of Cygnus A, nonthermal X-ray emission increases from the hotspot (some IC-CMB, mostly SSC) toward the offset radio synchrotron peak (mostly SSC). A faint thermal jet along the symmetry axis may be responsible for redirecting the Cygnus A non-thermal jet.

preprint2012arXiv

Radiating Bondi and Cooling Site Flows

Steady accretion of a radiating gas onto a central mass point is described and compared to classic Bondi accretion. Radiation losses are essential for accretion flows to be observed. Unlike Bondi flows, radiating Bondi flows pass through a sonic point at a finite radius and become supersonic near the center. The morphology of all radiating Bondi flows is described by a single dimensionless parameter. In radiating Bondi flows the mass accretion rate varies approximately as the first power of the central mass -- this differs significantly from the quadratic dependence on the central mass in classical Bondi flows. Mass accretion rates onto galaxy or cluster-centered black holes estimated from traditional and radiating Bondi flows are significantly different. In radiating Bondi flows the gas temperature increases at large radii, as in the cores of many galaxy groups and clusters, allowing radiating Bondi flows to merge naturally with gas arriving from their cluster environments. Some radiating flows cool completely before reaching the center of the flow, and this also occurs in cooling site flows in which there is no central gravitating mass.

preprint2012arXiv

Reconciling stellar dynamical and hydrostatic X-ray mass measurements of an elliptical galaxy with gas rotation, turbulence and magnetic fields

Recent hydrostatic X-ray studies of the hot interstellar medium (ISM) in early-type galaxies underestimate the gravitating mass as compared to stellar dynamics, implying modest, but significant deviations from exact hydrostatic equilibrium. We present a method for combining X-ray measurements and stellar dynamical constraints in the context of Bayesian statistics that allows the radial distribution of the implied nonthermal pressure or bulk motions in the hot ISM to be constrained. We demonstrate the accuracy of the method with hydrodynamical simulations tailored to produce a realistic galaxy model. Applying the method to the nearby elliptical galaxy NGC4649, we find a significant but subdominant nonthermal pressure fraction (0.27+/-0.06) in the central (<5 kpc) part of the galaxy, similar to the level of deviations from hydrostatic equilibrium expected in galaxy clusters. Plausible sources of systematic error, if important, may reduce this fraction. This would imply >360 km/s random turbulence or a magnetic field B=(39+/-6)(n_e/0.1 cm^{-3})^{0.59+/-0.09} muG, whereas gas rotation alone is unlikely to explain the detailed nonthermal profile. Future observations with Astro-H will allow turbulence or gas rotation at this level to be detected.

preprint2012arXiv

The Fermi Bubbles. I. Possible Evidence for Recent AGN Jet Activity in the Galaxy

The Fermi Gamma-ray Space Telescope reveals two large gamma-ray bubbles in the Galaxy, which extend about 50 degrees (~ 10 kpc) above and below the Galactic center (GC) and are symmetric about the Galactic plane. Using axisymmetric hydrodynamic simulations with a self-consistent treatment of the dynamical cosmic ray (CR) - gas interaction, we show that the bubbles can be created with a recent active galactic nucleus (AGN) jet activity about 1 - 3 Myr ago, which was active for a duration of ~ 0.1 - 0.5 Myr. The bipolar jets were ejected into the Galactic halo along the rotation axis of the Galaxy. Near the GC, the jets must be moderately light with a typical density contrast 0.001 <~ η<~ 0.1 relative to the ambient hot gas. The jets are energetically dominated by kinetic energy, and over-pressured with either CR or thermal pressure which induces lateral jet expansion, creating fat CR bubbles as observed. The sharp edges of the bubbles imply that CR diffusion across the bubble surface is strongly suppressed. The jet activity induces a strong shock, which heats and compresses the ambient gas in the Galactic halo, potentially explaining the ROSAT X-ray shell features surrounding the bubbles. The Fermi bubbles provide plausible evidence for a recent powerful AGN jet activity in our Galaxy, shedding new insights into the origin of the halo CR population and the channel through which massive black holes in disk galaxies release feedback energy during their growth.

preprint2012arXiv

The Fermi Bubbles. II. The Potential Roles of Viscosity and Cosmic Ray Diffusion in Jet Models

The origin of the Fermi bubbles recently detected by the Fermi Gamma-ray Space Telescope in the inner Galaxy is mysterious. In the companion paper Guo & Mathews (Paper I), we use hydrodynamic simulations to show that they could be produced by a recent powerful AGN jet event. Here we further explore this scenario to study the potential roles of shear viscosity and cosmic ray (CR) diffusion on the morphology and CR distribution of the bubbles. We show that even a relatively low level of viscosity (μ_{visc} >~ 3 g cm^{-1} s^{-1}, or ~0.1% - 1% of Braginskii viscosity in this context) could effectively suppress the development of Kelvin-Helmholtz instabilities at the bubble surface, resulting in smooth bubble edges as observed. Furthermore, viscosity reduces circulating motions within the bubbles, which would otherwise mix the CR-carrying jet backflow near bubble edges with the bubble interior. Thus viscosity naturally produces an edge-favored CR distribution, an important ingredient to produce the observed flat gamma-ray surface brightness distribution. Generically, such a CR distribution often produces a limb-brightened gamma-ray intensity distribution. However, we show that by incorporating CR diffusion which is strongly suppressed across the bubble surface (as inferred from sharp bubble edges) but is close to canonical values in the bubble interior, we obtain a reasonably flat gamma-ray intensity profile. The similarity of the resulting CR bubble with the observed Fermi bubbles strengthens our previous result in Paper I that the Fermi bubbles were produced by a recent AGN jet event. Studies of the nearby Fermi bubbles may provide a unique opportunity to study the potential roles of plasma viscosity and CR diffusion on the evolution of AGN jets and bubbles.

preprint2012arXiv

Tracing the Gas to the Virial Radius (R100) in a Fossil Group

We present a Chandra, Suzaku and Rosat study of the hot Intra Group Medium (IGrM) of the relaxed fossil group/ poor cluster RXJ1159+5531. This group exhibits an advantageous combination of flat surface brightness profile, high luminosity and optimal distance, allowing the gas to be detected out to the virial radius (Rvir=R108=1100 kpc) in a single Suzaku pointing, while the complementary Chandra data reveal a round morphology and relaxed IGrM image down to kpc scales. We measure the IGrM entropy profile over 3 orders of magnitude in radius, including 3 data bins beyond 0.5R200 with azimuthal coverage (>30%). We find no evidence that the profile flattens at large scales (>R500), and when corrected for the enclosed gas fraction, the entropy profile is very close to the predictions from self-similar structure formation simulations, as seen in massive clusters. Within Rvir, we measure a baryon fraction of 0.17+/-0.02, consistent with the Cosmological value. These results are in sharp contrast to the gas behaviour at large scales recently reported in the Virgo and Perseus clusters, and indicate that substantial gas clumping cannot be ubiquitous near Rvir, at least in highly evolved (fossil) groups.

preprint2011arXiv

Cosmic ray diffusion fronts in the Virgo cluster

The pair of large radio lobes in the Virgo cluster, each about 23 kpc in radius, have curiously sharp outer edges where the radio-synchrotron continuum flux declines abruptly. However, just adjacent to this sharp transition, the radio flux increases. This radio limb-brightening is observed over at least half of the perimeter of both lobes. We describe slowly propagating steady state diffusion fronts that explain these counterintuitive features. Because of the natural buoyancy of radio lobes, the magnetic field is largely tangent to the lobe boundary, an alignment that polarizes the radio emission and dramatically reduces the diffusion coefficient of relativistic electrons. As cosmic ray electrons diffuse slowly into the cluster gas, the local magnetic field and gas density are reduced as gas flows back toward the radio lobe. Radio emission peaks can occur because the synchrotron emissivity increases with magnetic field and then decreases with the density of non-thermal electrons. A detailed comparison of steady diffusion fronts with quantitative radio observations may reveal information about the spatial variation of magnetic fields and the diffusion coefficient of relativistic electrons. On larger scales, some reduction of the gas density inside the Virgo lobes due to cosmic ray pressure must occur and may be measurable. Such X-ray observations could reveal important information about the presence of otherwise unobservable non-thermal components such as relativistic electrons of low energy or proton cosmic rays.

preprint2011arXiv

Estimate of the total mechanical feedback energy from galaxy cluster-centered black holes: implications for black hole evolution, cluster gas fraction and entropy

The total feedback energy injected into hot gas in galaxy clusters by central black holes can be estimated by comparing the potential energy of observed cluster gas profiles with the potential energy of non-radiating, feedback-free hot gas atmospheres resulting from gravitational collapse in clusters of the same total mass. Feedback energy from cluster-centered black holes expands the cluster gas, lowering the gas-to-dark matter mass ratio below the cosmic value. Feedback energy is unnecessarily delivered by radio-emitting jets to distant gas far beyond the cooling radius where the cooling time equals the cluster lifetime. For clusters of mass 4-11 X 10^14 Msun estimates of the total feedback energy, 1-3 X 10^63 ergs, far exceed feedback energies estimated from observations of X-ray cavities and shocks in the cluster gas, energies gained from supernovae, and energies lost from cluster gas by radiation. The time-averaged mean feedback luminosity is comparable to those of powerful quasars, implying that some significant fraction of this energy may arise from the spin of the black hole. The universal entropy profile in feedback-free gaseous atmospheres in NFW cluster halos can be recovered by multiplying the observed gas entropy profile of any relaxed cluster by a factor involving the gas fraction profile. While the feedback energy and associated mass outflow in the clusters we consider far exceed that necessary to stop cooling inflow, the time-averaged mass outflow at the cooling radius almost exactly balances the mass that cools within this radius, an essential condition to shut down cluster cooling flows.

preprint2011arXiv

Resolving the Bondi Accretion Flow toward the Supermassive Black Hole of NGC 3115 with Chandra

Gas undergoing Bondi accretion onto a supermassive black hole (SMBH) becomes hotter toward smaller radii. We searched for this signature with a Chandra observation of the hot gas in NGC 3115, which optical observations show has a very massive SMBH. Our analysis suggests that we are resolving, for the first time, the accretion flow within the Bondi radius of an SMBH. We show that the temperature is rising toward the galaxy center as expected in all accretion models in which the black hole is gravitationally capturing the ambient gas. There is no hard central point source that could cause such an apparent rise in temperature. The data support that the Bondi radius is at about 4 arcsec-5 arcsec (188-235 pc), suggesting an SMBH of 2 x 10^9 M_sun that is consistent with the upper end of the optical results. The density profile within the Bondi radius has a power-law index of 1.03^{+0.23}_{-0.21} which is consistent with gas in transition from the ambient medium and the accretion flow. The accretion rate at the Bondi radius is determined to be {\dot M}_B = 2.2 x 10^{-2} M_sun yr^{-1}. Thus, the accretion luminosity with 10% radiative efficiency at the Bondi radius (10^{44} erg s^{-1}) is about six orders of magnitude higher than the upper limit of the X-ray luminosity of the nucleus.

preprint2010arXiv

Cosmic Ray-Dominated AGN Jets and the Formation of X-ray Cavities in Galaxy Clusters

It is widely accepted that feedback from active galactic nuclei (AGN) plays a key role in the evolution of gas in groups and clusters of galaxies. Unequivocal evidence comes from quasi-spherical X-ray cavities observed near cluster centers having sizes ranging from a few to tens of kpc, some containing radio emission. Cavities apparently evolve from the interaction of AGN jets with the intracluster medium (ICM). However, in numerical simulations it has been difficult to create such fat cavities from narrow jets. Ultra-hot thermal jets dominated by kinetic energy typically penetrate deep into the ICM, forming radially elongated cavities at large radii unlike those observed. Here, we study very light jets dominated energetically by relativistic cosmic rays (CRs) with axisymmetric hydrodynamic simulations, investigating the jet evolution both when they are active and when they are later turned off. We find that, when the thermal gas density in a CR-dominated jet is sufficiently low, the jet has a correspondingly low inertia, and thus decelerates quickly in the ICM. Furthermore, CR pressure causes the jet to expand laterally, encounter and displace more decelerating ICM gas, naturally producing fat cavities near cluster centers similar to those observed. Our calculations of cavity formation imply that AGN jets responsible for creating fat X-ray cavities (radio bubbles) are very light, and dominated by CRs. This scenario is consistent with radio observations of Fanaroff-Riley I jets that appear to decelerate rapidly, produce strong synchrotron emission and expand typically at distances of a few kpc from the central AGN.

preprint2010arXiv

Removing Cool Cores and Central Metallicity Peaks in Galaxy Clusters with Powerful AGN Outbursts

Recent X-ray observations of galaxy clusters suggest that cluster populations are bimodally distributed according to central gas entropy and are separated into two distinct classes: cool core (CC) and non-cool core (NCC) clusters. While it is widely accepted that AGN feedback plays a key role in offsetting radiative losses and maintaining many clusters in the CC state, the origin of NCC clusters is much less clear. At the same time, a handful of extremely powerful AGN outbursts have recently been detected in clusters, with a total energy ~10^{61}-10^{62} erg. Using two dimensional hydrodynamic simulations, we show that if a large fraction of this energy is deposited near the centers of CC clusters, which is likely common due to dense cores, these AGN outbursts can completely remove CCs, transforming them to NCC clusters. Our model also has interesting implications for cluster abundance profiles, which usually show a central peak in CC systems. Our calculations indicate that during the CC to NCC transformation, AGN outbursts efficiently mix metals in cluster central regions, and may even remove central abundance peaks if they are not broad enough. For CC clusters with broad central abundance peaks, AGN outbursts decrease peak abundances, but can not effectively destroy the peaks. Our model may simultaneously explain the contradictory (possibly bimodal) results of abundance profiles in NCC clusters, some of which are nearly flat, while others have strong central peaks similar to those in CC clusters. A statistical analysis of the sizes of central abundance peaks and their redshift evolution may shed interesting insights on the origin of both types of NCC clusters and the evolution history of thermodynamics and AGN activity in clusters.

preprint2010arXiv

Self-consistent evolution of gas and cosmic rays in Cygnus A and similar FR II classic double radio sources

In Cygnus A and other classical FR II double radio sources, powerful opposing jets from the cores of halo-centered galaxies drive out into the surrounding cluster gas, forming hotspots of shocked and compressed cluster gas at the jet extremities. The moving hotspots are sandwiched between two shocks. An inner-facing shock receives momentum and cosmic rays from the jet and creates additional cosmic rays that form a radio lobe elongated along the jet axis. An outer-facing bow shock moves directly into the undisturbed group or cluster gas, creating a cocoon of shocked gas enclosing the radio lobe. We describe computations that follow the self-consistent dynamical evolution of the shocked cluster gas and the relativistic synchrotron-emitting gas inside the lobes. Relativistic and non-relativistic components exchange momentum by interacting with small magnetic fields having dynamically negligible energy densities. The evolution of Cygnus A is governed almost entirely by cosmic ray energy flowing from the hotspots. Mass flowing into hotspots from the jets is assumed to be small, greatly reducing the mass of gas flowing back along the jet, common in previous calculations, that would disrupt the spatial segregation of synchrotron-loss ages observed inside FR II radio lobes. We compute the evolution of the cocoon when the velocity and cosmic ray luminosity of the hotspots are constant and when they vary with time. If cosmic rays mix with cluster gas in hotspots before flowing into the radio lobe, the thermal gas is heated to mildly relativistic temperatures, producing an unobserved pressure inside the lobe.

preprint2010arXiv

Simulating X-ray Supercavities and Their Impact on Galaxy Clusters

Recent X-ray observations of hot gas in the galaxy cluster MS 0735.6+7421 reveal huge radio-bright, quasi-bipolar X-ray cavities having a total energy ~10^{62} ergs, the most energetic AGN outburst currently known. We investigate the evolution of this outburst with two-dimensional axisymmetric gasdynamical calculations in which the cavities are inflated by relativistic cosmic rays. Many key observational features of the cavities and associated shocks are successfully reproduced. The radial elongation of the cavities indicates that cosmic rays were injected into the cluster gas by a (jet) source moving out from the central AGN. AGN jets of this magnitude must be almost perfectly identically bipolar. The relativistic momentum of a single jet would cause a central AGN black hole of mass 10^9 M_{sun} to recoil at ~6000 km s^{-1}, exceeding kick velocities during black hole mergers, and be ejected from the cluster-center galaxy. When the cavity inflation is complete, 4PV underestimates the total energy received by the cluster gas. Deviations of the cluster gas from hydrostatic equilibrium are most pronounced during the early cavity evolution when the integrated cluster mass found from the observed gas pressure gradient can have systematic errors near the cavities of ~10-30%. The creation of the cavity with cosmic rays generates a long-lasting global cluster expansion that reduces the total gas thermal energy below that received from the cavity shock. One Gyr after this single outburst, a gas mass of ~ 6 \times 10^{11} M_{sun} is transported out beyond a cluster radius of 500 kpc. Such post-cavity outflows can naturally produce the discrepancy observed between the cluster gas mass fraction and the universal baryon fraction inferred from WMAP observations. (Abridged)

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2009arXiv

AGN Feedback in Galaxy Groups: the two interesting cases of AWM 4 and NGC 5044

We present AGN feedback in the interesting cases of two groups: AWM 4 and NGC 5044. AWM 4 is characterized by a combination of properties which seems to defy the paradigm for AGN heating in cluster cores: a flat inner temperature profile indicative of a past, major heating episode which completely erased the cool core, as testified by the high central cooling time (> 3 Gyrs) and by the high central entropy level (~ 50 keV cm^2), and yet an active central radio galaxy with extended radio lobes out to 100 kpc, revealing recent feeding of the central massive black hole. A recent Chandra observation has revealed the presence of a compact cool corona associated with the BCG, solving the puzzle of the apparent lack of low entropy gas surrounding a bright radio source, but opening the question of its origin. NGC 5044 shows in the inner 10 kpc a pair of cavities together with a set of bright filaments. The cavities are consistent with a recent AGN outburst as also indicated by the extent of dust and H_alpha emission even though the absence of extended 1.4 GHz emission remains to be explained. The soft X-ray filaments coincident with H_alpha and dust emission are cooler than those which do not correlate with optical and infrared emission, suggesting that dust-aided cooling can contribute to the overall cooling. For the first time sloshing cold fronts at the scale of a galaxy group have been observed in this object.

preprint2009arXiv

Cosmic Ray Feedback

Cosmic rays produced or deposited at sites in hot cluster gas are thought to provide the pressure that forms X-ray cavities. While cavities have a net cooling effect on cluster gas, young, expanding cavities drive shocks that increase the local entropy. Cavities also produce radial filaments of thermal gas and are sources of cluster cosmic rays that diffuse through cavity walls, as in Virgo where a radio lobe surrounds a radial thermal filament. Cosmic rays also make the hot gas locally buoyant, allowing large masses of low entropy gas to be transported out beyond the cooling radius. Successive cavities maintain a buoyant outflow that preserves the cluster gas temperature and gas fraction profiles and dramatically reduces the cooling rate onto the central black hole.

preprint2009arXiv

Spitzer Observations of Passive and Star Forming Early-type Galaxies: an Infrared Color-Color Sequence

We describe the infrared properties of a large sample of early type galaxies, comparing data from the Spitzer archive with Ks-band emission from 2MASS. While most representations of this data result in correlations with large scatter, we find a remarkably tight relation among colors formed by ratios of luminosities in Spitzer-MIPS (24, 70 and 160 um) bands and the Ks-band. Remarkably, this correlation among E and S0 galaxies follows that of nearby normal galaxies of all morphological types. In particular, the tight infrared color-color correlation for S0 galaxies alone follows that of the entire Hubble sequence of normal galaxies, roughly in order of galaxy type from ellipticals to spirals to irregulars. The specific star formation rate of S0 galaxies estimated from the 24um luminosity increases with decreasing Ks-band luminosity (or stellar mass) from essentially zero, as with most massive ellipticals, to rates typical of irregular galaxies. Moreover, the luminosities of the many infrared-luminous S0 galaxies can significantly exceed those of the most luminous (presumably post-merger) E galaxies. Star formation rates in the most infrared-luminous S0 galaxies approach 1-10 solar masses per year. Consistently with this picture we find that while most early-type galaxies populate an infrared red sequence, about 24% of the objects (mostly S0s) are in an infrared blue cloud together with late type galaxies. For those early-type galaxies also observed at radio frequencies we find that the far-infrared luminosities correlate with the mass of neutral and molecular hydrogen, but the scatter is large. This scatter suggests that the star formation may be intermittent or that similar S0 galaxies with cold gaseous disks of nearly equal mass can have varying radial column density distributions that alter the local and global SF rates.

preprint2009arXiv

X-ray Isophotes in a Rapidly Rotating Elliptical Galaxy: Evidence of Inflowing Gas

We describe two-dimensional gasdynamical computations of the X-ray emitting gas in the rotating elliptical galaxy NGC 4649 that indicate an inflow of about one solar mass per year at every radius. Such a large instantaneous inflow cannot have persisted over a Hubble time. The central constant-entropy temperature peak recently observed in the innermost 150 parsecs is explained by compressive heating as gas flows toward the central massive black hole. Since the cooling time of this gas is only a few million years, NGC 4649 provides the most acutely concentrated known example of the cooling flow problem in which the time-integrated apparent mass that has flowed into the galactic core exceeds the total mass observed there. This paradox can be resolved by intermittent outflows of energy or mass driven by accretion energy released near the black hole. Inflowing gas is also required at intermediate kpc radii to explain the ellipticity of X-ray isophotes due to spin-up by mass ejected by stars that rotate with the galaxy and to explain local density and temperature profiles. We provide evidence that many luminous elliptical galaxies undergo similar inflow spin-up. A small turbulent viscosity is required in NGC 4649 to avoid forming large X-ray luminous disks that are not observed, but the turbulent pressure is small and does not interfere with mass determinations that assume hydrostatic equilibrium.

preprint2008arXiv

Evidence of Star Formation in Local S0 Galaxies: Spitzer Observations of the SAURON Sample

We discuss infrared Spitzer observations of early type galaxies in the SAURON sample at 24, 60 and 170 microns. When compared with 2MASS Ks band luminosities, lenticular (S0) galaxies exhibit a much wider range of mid to far-infrared luminosities then elliptical (E) galaxies. Mid and far-infrared emission from E galaxies is a combination of circumstellar or interstellar emission from local mass-losing red giant stars, dust buoyantly transported from the galactic cores into distant hot interstellar gas and dust accreted from the environment. The source of mid and far-IR emission in S0 galaxies is quite different and is consistent with low levels of star formation, 0.02 - 0.2 Msol/yr, in cold, dusty gaseous disks. The infrared 24micron-70micron color is systematically lower for (mostly S0) galaxies with known molecular disks. Our observations support the conjecture that cold dusty gas in some S0 galaxies is created by stellar mass loss at approximately the same rate that it is consumed by star formation, so the mass depletion of these disks by star formation will be slow. Unlike E galaxies, the infrared luminosities of S0 galaxies correlate with both the mass of molecular gas and the stellar Hbeta spectral index, and all are related to the recent star formation rate. However, star formation rates estimated from the Hbeta emission line luminosities L_{Hbeta} in SAURON S0 galaxies are generally much smaller. Since L_{Hbeta} does not correlate with 24 microns emission from dust heated by young stars, optical emission lines appear to be a poor indicator of star formation rates in SAURON S0 galaxies. The absence of Hbeta emission may be due to a relative absence of OB stars in the initial mass function or to dust absorption of Hbeta emission lines.

preprint2008arXiv

The X-ray cavities, filaments and cold fronts in the core of the galaxy group NGC 5044

We present a two-dimensional analysis of the bright nearby galaxy group NGC 5044 using the currently available Chandra and XMM data. In the inner 10 kpc a pair of cavities are evident together with a set of bright X-ray filaments. If the cavities are interpreted as gas displaced by relativistic plasma inflated by an AGN, even in the absence of extended 1.4 GHz emission, this would be consistent with a recent outburst as also indicated by the extent of dust and H_alpha emission. The soft X-ray filaments coincident with H_alpha and dust emission are cooler than the ones which do not correlate with optical and infrared emission. We suggest that dust-aided cooling contributes to form warm (T =10^4 K) gas, emitting H_alpha radiation. At 31 kpc and 67 kpc a pair of cold fronts are present, indicative of sloshing due to a dynamical perturbation caused by accretion of a less massive group, also suggested by the peculiar velocity of the brightest galaxy NGC 5044 with respect to the mean group velocity.

preprint2007arXiv

Far Infrared Spitzer Observations of Elliptical Galaxies: Evidence for Extended Diffuse Dust

Far-infrared Spitzer observations of elliptical galaxies are inconsistent with simple steady state models of dust creation in red giant stars and destruction by grain sputtering in the hot interstellar gas at T ~ 10^7 K. The flux at 24 microns correlates with optical fluxes, suggesting that this relatively hot dust is largely circumstellar. But fluxes at 70 and 160 microns do not correlate with optical fluxes. Elliptical galaxies with similar L_B have luminosities at 70 and 160 microns (L_70 and L_160) that vary over a factor ~ 100, implying an additional source of dust unrelated to that produced by ongoing local stellar mass loss. Neither L_70/L_B nor L_160/L_B correlate with the stellar age or metallicity. Optical line fluxes from warm gas at T ~ 10^4 K correlate weakly with L_70 and L_160, suggesting that the dust may be responsible for cooling this gas. Many normal elliptical galaxies have emission at 70 microns that is extended to 5-10 kpc. Extended far-infrared emission with sputtering lifetimes of ~10^8 yrs is difficult to maintain by mergers with gas-rich galaxies. Instead, we propose that this cold dust is buoyantly transported from reservoirs of dust in the galactic cores which is supplied by mass loss from stars in the core. Intermittent energy outbursts from AGNs can drive the buoyant outflow.

preprint2007arXiv

Trouble for AGN Feedback ? The puzzle of the core of the Galaxy Cluster AWM 4

The core of the relaxed cluster AWM 4 is characterized by a unique combination of properties which defy a popular scenario for ANG heating of cluster cores. A flat inner temperature profile is indicative of a past, major heating episode which completely erased the cool core, as testified by the high central cooling time (~ 3 Gyr) and by the high central entropy level (~ 60 keV cm^2). Yet the presence of a 1.4 GHz active central radio galaxy with extended radio lobes out to 100 kpc, reveals recent feeding of the central massive black hole. A system like AWM 4 should have no radio emission at all if only feedback from the cooling hot gas regulates the AGN activity.

preprint2003arXiv

Cold Dust in Early-Type Galaxies. I. Observations

We describe far infrared observations of early-type galaxies selected from the ISO archive. This rather inhomogeneous sample includes 39 giant elliptical galaxies and 14 S0 (or later) galaxies. These galaxies were observed using a variety of different observing modes, each of which requires different and often rather elaborate photometric reduction procedures. The ISO background data agrees well with COBE-DIRBE results to which we have renormalized our calibrations. As a further check, the ISO fluxes from galaxies at 60 and 100 microns agree very well with those previously observed with IRAS at these wavelengths. The spatial resolution of ISO is several times greater than that of IRAS and the ISO observations extend out to 200 microns which views a significantly greater mass of colder dust not assessable to IRAS. Most of the galaxies are essentially point sources at ISO resolution, but a few are clearly extended at FIR wavelengths with image sizes that increase with FIR wavelength. The integrated far infrared luminosities do not correlate with optical luminosities, suggesting that the dust may have an external, merger-related origin. In general the far infrared spectral energy distributions can be modeled with dust at two temperatures, ~ 43 K and ~ 20 K, which probably represent limits of a continuous range of temperatures. The colder dust component dominates the total mass of dust, which is typically more than ten times larger than the dust masses previously estimated for the same galaxies using IRAS observations. For S0 galaxies we find that the optically normalized far-infrared luminosity $L_{FIR} / L_B$ correlates strongly with the mid-infrared luminosity $L_{15μ\rm{m}} / L_B$, but that correlation is weaker for elliptical galaxies.

preprint1995arXiv

Interstellar Grains in Elliptical Galaxies: Grain Evolution

We consider the lifecycle of dust introduced into the hot interstellar medium in isolated elliptical galaxies. Dust grains are ejected into galactic-scale cooling flows in large ellipticals by normal mass loss from evolving red giants. Newly introduced dust rapidly enters the hot gas environment and is sputtered away by thermal collisions with ions. Before the grains are completely sputtered away, they emit prodigious amounts of infrared radiation which may contribute to the large far infrared luminosities observed in ellipticals. In order to study the global properties of grains in ellipticals we construct a new series of King-type galactic models which are consistent with the fundamental plane, galactic mass to light ratios and other relevant observational correlations. We describe a new ``continuity'' procedure to construct simple time-dependent gas dynamic models for cooling flows. In all galaxy models, although grains can flow a considerable distance from their radius of origin before being sputtered away, the grain size distribution is accurately determined by assuming {\it in situ} sputtering, completely ignoring advection. The dominant source of grain heating is absorption of starlight; grain heating by collisions with energetic thermal electrons or X-ray absorption are negligible. However, we also show that when self-consistent grain sputtering is included the dust-to-gas ratio is reduced and radiative cooling, not electron-grain interactions dominates plasma cooling, even for the most massive ellipticals.