We present new ALMA CO(2--1) observations of two well studied group-centered elliptical galaxies: NGC~4636 and NGC~5846. In addition, we include a revised analysis of Cycle 0 ALMA observations of the central galaxy in the NGC~5044 group that has been previously published. We find evidence that molecular gas, in the form of off-center orbiting clouds, is a common presence in bright group-centered galaxies (BGG). CO line widths are $\gtrsim 10$ times broader than Galactic molecular clouds, and using the reference Milky Way $X_{CO}$, the total molecular mass ranges from as low as $2.6\times 10^5 M_\odot$ in NGC~4636 to $6.1\times 10^7 M_\odot$ in NGC~5044. With these parameters the virial parameters of the molecular structures is $\gg 1$. Complementary observations of NGC~5846 and NGC~4636 using the ALMA Compact Array (ACA) do not exhibit any detection of a CO diffuse component at the sensitivity level achieved by current exposures. The origin of the detected molecular features is still uncertain, but these ALMA observations suggest that they are the end product of the hot gas cooling process and not the result of merger events. Some of the molecular clouds are associated with dust features as revealed by HST dust extinction maps suggesting that these clouds formed from dust-enhanced cooling. The global nonlinear condensation may be triggered via the chaotic turbulent field or buoyant uplift. The large virial parameter of the molecular structures and correlation with the warm ($10^3 - 10^5 K$)/hot ($\ge10^6$) phase velocity dispersion provide evidence that they are unbound giant molecular associations drifting in the turbulent field, consistently with numerical predictions of the chaotic cold accretion process. Alternatively, the observed large CO line widths may be generated by molecular gas flowing out from cloud surfaces due to heating by the local hot gas atmosphere.

We use 2.0 Msec of Chandra observations to investigate the cocoon shocks of Cygnus A and some implications for its lobes and jet. Measured shock Mach numbers vary in the range 1.18-1.66 around the cocoon. We estimate a total outburst energy of $\simeq 4.7\times10^{60}\rm\ erg$, with an age of $\simeq 2 \times 10^{7}\rm\ yr$. The average postshock pressure is found to be $8.6 \pm 0.3 \times 10^{-10}\rm\ erg\ cm^{-3}$, which agrees with the average pressure of the thin rim of compressed gas between the radio lobes and shocks, as determined from X-ray spectra. However, average rim pressures are found to be lower in the western lobe than in the eastern lobe by $\simeq 20\%$. Pressure estimates for hotspots A and D from synchrotron self-Compton models imply that each jet exerts a ram pressure $\gtrsim$ 3 times its static pressure, consistent with the positions of the hotspots moving about on the cocoon shock over time. A steady, one-dimensional flow model is used to estimate jet properties, finding mildly relativistic flow speeds within the allowed parameter range. Models in which the jet carries a negligible flux of rest mass are consistent with with the observed properties of the jets and hotspots. This favors the jets being light, implying that the kinetic power and momentum flux are carried primarily by the internal energy of the jet plasma rather than by its rest mass.

We use new and archival Chandra observations of Cygnus A, totalling $\sim$1.9 Ms, to investigate the distribution and temperature structure of gas lying within the projected extent of the cocoon shock and exhibiting a rib-like structure. We confirm that the X-rays are dominated by thermal emission with an average temperature of around 4 keV, and have discovered an asymmetry in the temperature gradient, with the southwestern part of the gas cooler than the rest by up to 2 keV. Pressure estimates suggest that the gas is a coherent structure of single origin located inside the cocoon, with a mass of roughly $2\times10^{10} M_{\odot}$. We conclude that the gas is debris resulting from disintegration of the cool core of the Cygnus A cluster after the passage of the jet during the early stages of the current epoch of activity. The 4 keV gas now lies on the central inside surface of the hotter cocoon rim. The temperature gradient could result from an offset between the centre of the cluster core and the Cygnus A host galaxy at the switch-on of current radio activity.

OVI absorption in quasar spectra caused by intervening circumgalactic atmospheres suggests a downturn in the atmospheric column density in sightlines passing beyond about 100 kpc from central star-forming galaxies. This turnover supports the hypothesis that the oxygen originates in the central galaxies. When converted into oxygen space density using an Abel integral inversion, the OVI columns require greater than $\approx 10^9 M_\odot$ of oxygen concentrated near 100 kpc. Circumgalactic gas within this radius cools in less than 1 Gyr and radiates $\sim 10^{42.2}$ erg s$^{-1}$ overall. The feedback power necessary to maintain such oxygen-rich atmospheres for many Gyrs cannot be easily supplied by galactic supernovae. However, massive central black holes in star-forming galaxies may generate sufficient accretion power and intermittent shock waves at $r \sim 100$ kpc to balance circumgalactic radiation losses in late-type $L^\star$ galaxies. The relative absence of OVI absorption observed in early-type, passive $L^{\star}$ galaxies may arise from enhanced AGN feedback from their more massive central black holes.

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.

A short 30 minute ALMA observation of the early-type galaxy NGC 5044, which resides at the center of an X-ray bright group with a moderate cooling flow, has detected 24 molecular structures within the central 2.5 kpc. The masses of the molecular structures vary from 3e5 to 1e7 Mo3 and the CO(2-1) linewidths vary from 15 to 65 km/s. Given the large CO(2-1) linewidths, the observed structures are likely giant molecular associations (GMAs) and not individual molecular clouds (GMCs). Only a few of the GMAs are spatially resolved with the cycle 0 ALMA beam and the average density of these GMAs yields a GMC volume filling factor of about 15%. The observed masses of the resolved GMAs are insufficient for them to be gravitationally bound, however, the most massive GMA does contain a less massive component with a linewidth of 5.5 km/s (typical of an individual virialized GMC). We also show that the GMAs cannot be pressure confined by the hot gas. Given the observed CO(2-1) linewidths of the GMAs (i.e., the velocity dispersion of the embedded GMCs) they will likely disperse on a timescale of about 12 Myr, which is less than the central cooling time of the hot gas, so the embedded GMCs within a GMA must condense out of the hot gas at the same time and arise from local concentrations of thermally unstable parcels of hot gas. There are no indications of any disk-like molecular structures and all indications suggest that the molecular gas follows ballistic trajectories after condensing out of the thermally unstable hot gas. The 230 GHz luminosity of the central continuum source is 500 times greater than its low frequency radio luminosity and probably reflects a recent accretion event by the central supermassive black hole. The spectrum of the central continuum source also exhibits an absorption feature.

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.

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.

We present a combined X-ray, optical, and radio analysis of the galaxy group IC 1860 using the currently available Chandra and XMM data, literature multi-object spectroscopy data and GMRT data. The Chandra and XMM imaging and spectroscopy reveal two surface brightness discontinuities at 45 and 76 kpc shown to be consistent with a pair of cold fronts. These features are interpreted as due to sloshing of the central gas induced by an off-axis minor merger with a perturber. This scenario is further supported by the presence of a peculiar velocity of the central galaxy IC 1860 and the identification of a possible perturber in the optically disturbed spiral galaxy IC 1859. The identification of the perturber is consistent with the comparison with numerical simulations of sloshing. The GMRT observation at 325 MHz shows faint, extended radio emission contained within the inner cold front, as seen in some galaxy clusters hosting diffuse radio mini-halos. However, unlike mini-halos, no particle reacceleration is needed to explain the extended radio emission, which is consistent with aged radio plasma redistributed by the sloshing. There is strong analogy of the X-ray and optical phenomenology of the IC 1860 group with two other groups, NGC 5044 and NGC 5846, showing cold fronts. The evidence presented in this paper is among the strongest supporting the currently favored model of cold-front formation in relaxed objects and establishes the group scale as a chief environment to study this phenomenon.

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.

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.

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 <~ \eta <~ 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.

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 (\mu_{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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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)

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.

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.

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.

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.

We present new mass measurements for the supermassive black holes (SMBHs) in the centres of three early-type galaxies. The gas pressure in the surrounding, hot interstellar medium (ISM) is measured through spatially resolved spectroscopy with the Chandra X-ray observatory, allowing the SMBH mass (Mbh) to be inferred directly under the hydrostatic approximation. This technique does not require calibration against other SMBH measurement methods and its accuracy depends only on the ISM being close to hydrostatic, which is supported by the smooth X-ray isophotes of the galaxies. Combined with results from our recent study of the elliptical galaxy NGC4649, this brings to four the number of galaxies with SMBHs measured in this way. Of these, three already have mass determinations from the kinematics of either the stars or a central gas disc, and hence join only a handful of galaxies with Mbh measured by more than one technique. We find good agreement between the different methods, providing support for the assumptions implicit in both the hydrostatic and the dynamical models. The stellar mass-to-light ratios for each galaxy inferred by our technique are in agreement with the predictions of stellar population synthesis models assuming a Kroupa initial mass function (IMF). This concurrence implies that no more than ~10-20% of the ISM pressure is nonthermal, unless there is a conspiracy between the shape of the IMF and nonthermal pressure. Finally, we compute Bondi accretion rates, finding that the two galaxies with the highest rates exhibit little evidence of X-ray cavities, suggesting that the correlation with the AGN jet power takes time to be established.

Multi-Gyr two-dimensional calculations describe the gasdynamical evolution of hot gas in the Virgo cluster resulting from intermittent cavities formed with cosmic rays. Without cosmic rays, the gas evolves into a cooling flow, depositing about 85 solar masses per year of cold gas in the cluster core -- such uninhibited cooling conflicts with X-ray spectra and many other observations. When cosmic rays are produced or deposited 10 kpc from the cluster center in bursts of about 10^{59} ergs lasting 20 Myrs and spaced at intervals of 200 Myrs, the central cooling rate is greatly reduced to 0.1 - 1 solar masses per year, consistent with observations. After cosmic rays diffuse through the cavity walls, the ambient gas density is reduced and is buoyantly transported 30-70 kpc out into the cluster. Cosmic rays do not directly heat the gas and the modest shock heating around young cavities is offset by global cooling as the cluster gas expands. After several Gyrs the hot gas density and temperature profiles remain similar to those observed, provided the time-averaged cosmic ray luminosity is about 10^{43} erg/s, approximately equal to the bolometric X-ray luminosity within only 56 kpc. If an appreciable fraction of the relativistic cosmic rays are protons, gamma rays produced by pion decay following inelastic p-p collisions may be detected with the Fermi Gamma Ray Telescope.

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.

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.

We describe the formation and evolution of X-ray cavities in the hot gas of galaxy clusters. The cavities are formed only with relativistic cosmic rays that eventually diffuse into the surrounding gas. We explore the evolution of cavities formed with a wide range of cosmic ray diffusion rates. In previous numerical simulations cavities are formed by injecting ultra-hot but non-relativistic gas which increases the global thermal energy, offsetting radiative losses in the gas and helping to solve the cooling flow problem. Contrary to these results, we find that X-ray cavities formed solely by cosmic rays have a global cooling effect. As the cluster gas is displaced by cosmic rays, a global expansion of the cluster gas occurs with associated cooling that exceeds the heating by shock waves as the cavity forms. Most cosmic rays in our cavity evolutions do not move beyond the cooling radius even after 1 Gyr. The gas density is depressed by cosmic rays, becomes buoyant, and undergoes a significant outward mass transfer within the cooling radius, carrying cosmic rays and relatively low entropy gas to distant regions in the cluster where it remains for times exceeding the local cooling time in the hot gas. This post-cavity mass outflow due to cosmic ray buoyancy may contribute toward solving the cooling flow problem. We describe the energetics, size, stability and buoyant rise of X-ray cavities in detail, showing how each depends on the rate of cosmic ray diffusion.

We present a Chandra study of the hot ISM in the giant elliptical galaxy NGC4649. In common with other group-centred ellipticals, its temperature profile rises with radius in the outer parts of the galaxy, from ~0.7keV at 2kpc to ~0.9keV by 20kpc. However, within the central ~2kpc the trend reverses and the temperature peaks at ~1.1keV within the innermost 200pc. Under the assumption of hydrostatic equilibrium, we demonstrate that the central temperature spike arises due to the gravitational influence of a quiescent central super-massive black hole. We constrain the black hole mass (MBH) to $(3.35^{+0.67}_{-0.95})\times 10^9$Msun (90% confidence), in good agreement with stellar kinematics measurements. This is the first direct measurement of MBH based on studies of hydrostatic X-ray emitting gas, which are sensitive to the most massive black holes, and is a crucial validation of both mass-determination techniques. This agreement clearly demonstrates the gas must be close to hydrostatic, even in the very centre of the galaxy, which is consistent with the lack of morphological disturbances in the X-ray image. NGC4649 is now one of only a handful of galaxies for which MBH has been measured by more than one method. At larger radii, we were able to decompose the gravitating mass profile into stellar and dark matter (DM) components. Unless one accounts for the DM, a standard Virial analysis of the stars dramatically over-estimates the stellar mass of the galaxy. We find the measured J-band stellar mass-to-light ratio, 1.37+/-0.10 Msun/Lsun, is in good agreement with simple stellar population model calculations for this object.

Young cavities in the X-ray emitting hot gas in galaxy clusters are often filled with radio synchrotron emission and it is widely thought that the cavities are inflated by these cosmic rays. At a later stage of its evolution, when the cavity becomes buoyant, the converging flow of gas beneath the cavity results in a filament of thermal gas, a cavity jet, that moves radially outward at large subsonic velocities. As the cavity jet forms, the cosmic ray electrons may diffuse through the cavity walls, filling a large volume surrounding the cavity jet, as observed in M87/Virgo and elsewhere, sometimes referred to as relic radio sources. We compute the combined evolution of cosmic rays, cavities and cavity jets. The observed pattern in M87/Virgo can be reached in 100 Myrs, matching the synchrotron age of the extended radio source. A 20-30 kpc long cavity jet is surrounded by a quasi-spherical radio lobe 40 kpc in diameter, but the initial cavity has disappeared. At later times the cavity jet will fall back to the origin, leaving only the extended radio source. The combined jet-lobe evolution in M87/Virgo requires a total cosmic ray energy that is more than 10 times larger than that usually assumed, 4PV.

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.

We describe photometry at mid-infrared passbands (1.2 - 24 microns) for a sample of 18 elliptical galaxies. All surface brightness distributions resemble de Vaucouleurs profiles, indicating that most of the emission arises from the photospheres or circumstellar regions of red giant stars. The spectral energy distribution peaks near 1.6 microns, but the half-light or effective radius has a pronounced minimum near the K band (2.15 microns). Apart from the 24 micron passband, all sample-averaged radial color profiles have measurable slopes within about twice the (K band) effective radius. Evidently this variation arises because of an increase in stellar metallicity toward the galactic cores. For example, the sampled-averaged color profile (K - 5.8 microns) has a positive slope although no obvious absorption feature is observed in spectra of elliptical galaxies near 5.8 microns. This, and the minimum in the effective radius, suggests that the K band may be anomalously luminous in metal-rich stars in galaxy cores. Unusual radial color profiles involving the 24 micron passband may suggest that some 24 micron emission comes from interstellar not circumstellar dust grains.

Spitzer observations of extended dust in two optically normal elliptical galaxies provide a new confirmation of buoyant feedback outflow in the hot gas atmospheres around these galaxies. AGN feedback energy is required to prevent wholesale cooling and star formation in these group-centered galaxies. In NGC 5044 we observe interstellar (presumably PAH) emission at 8 microns out to about 5 kpc. Both NGC 5044 and 4636 have extended 70 microns emission from cold dust exceeding that expected from stellar mass loss. The sputtering lifetime of this extended dust in the ~1keV interstellar gas, ~10^7 yrs, establishes the time when the dust first entered the hot gas. Evidently the extended dust originated in dusty disks or clouds, commonly observed in elliptical galaxy cores, that were disrupted, heated and buoyantly transported outward. The surviving central dust in NGC 5044 and 4636 has been disrupted into many small filaments. It is remarkable that the asymmetrically extended 8 micron emission in NGC 5044 is spatially coincident with Halpha+[NII] emission from warm gas. A calculation shows that dust-assisted cooling in buoyant hot gas moving out from the galactic core can cool within a few kpc in about ~10^7 yrs, explaining the optical line emission observed. The X-ray images of both galaxies are disturbed. All timescales for transient activity - restoration of equilibrium and buoyant transport in the hot gas, dynamics of surviving dust fragments, and dust sputtering - are consistent with a central release of feedback energy in both galaxies about 10^7 yrs ago.

We present radial mass profiles within 0.3 r_vir for 16 relaxed galaxy groups-poor clusters (kT range 1-3 keV) selected for optimal mass constraints from the Chandra and XMM data archives. After accounting for the mass of hot gas, the resulting mass profiles are described well by a two-component model consisting of dark matter (DM), represented by an NFW model, and stars from the central galaxy. The stellar component is required only for 8 systems, for which reasonable stellar mass-to-light ratios (M/L_K) are obtained, assuming a Kroupa IMF. Modifying the NFW dark matter halo by adiabatic contraction does not improve the fit and yields systematically lower M/L_K. In contrast to previous results for massive clusters, we find that the NFW concentration parameter (c_vir) for groups decreases with increasing M_vir and is inconsistent with no variation at the 3 sigma level. The normalization and slope of the c_vir-M_vir relation are consistent with the standard LambdaCDM cosmological model with sigma_8 = 0.9. The small intrinsic scatter measured about the c_vir-M_vir relation implies the groups represent preferentially relaxed, early forming systems. The mean gas fraction (f =0.05 +/- 0.01) of the groups measured within an overdensity Delta=2500 is lower than for hot, massive clusters, but the fractional scatter (sigma_f/f=0.2) for groups is larger, implying a greater impact of feedback processes on groups, as expected.

Based on results from cold dark matter N-body simulations we develop a dynamical model for the evolution of subhaloes within host haloes of galaxy groups. Only subhaloes more massive than 5 times 10^8 M_{sol} at the time of accretion are examined because they are massive enough to possibly host luminous galaxies. As they orbit within a growing host potential the subhaloes are subject to tidal stripping and dynamical friction. We consider groups of equal mass (M_{vir} = 3.9 times 10^{13} M_{sol}) at redshift z=0 but with different concentrations associated with different formation times. We investigate the variation of subhaloe (or satellite galaxy) velocity dispersion with host concentration and/or formation time. In agreement with the Jeans equation the velocity dispersion of subhaloes increases with the host concentration. Between concentrations ~5 and ~20 the subhaloe velocity dispersions increase by ~25 per cent. By applying a simplified tidal disruption criterion, i.e. rejection of all subhaloes with a tidal truncation radius below 3 kpc at z=0, the central velocity dispersion of 'surviving' subhaloes increases substantially for all concentrations. The enhanced central velocity dispersion among surviving subhaloes is caused by a lack of slow tangential motions. Additionally, we present a fitting formula for the velocity anisotropy parameter \beta(r) which does not depend on concentration if the group-centric distances are scaled by r_s, the characteristic radius of the NFW-profile.

We present the concentration (c)-virial mass (M) relation of 39 galaxy systems ranging in mass from individual early-type galaxies up to the most massive galaxy clusters, (0.06-20) x 10^{14} M_sun. We selected for analysis the most relaxed systems possessing the highest quality data currently available in the Chandra and XMM public data archives. A power-law model fitted to the X-ray c-M relation requires at high significance (6.6 sigma) that c decreases with increasing M, which is a general feature of CDM models. The median and scatter of the c-M relation produced by the flat, concordance LCDM model (Omega_m=0.3, sigma_8=0.9) agrees with the X-ray data provided the sample is comprised of the most relaxed, early forming systems, which is consistent with our selection criteria. Holding the rest of the cosmological parameters fixed to those in the concordance model the c-M relation requires 0.76< sigma_8 <1.07 (99% conf.), assuming a 10% upward bias in the concentrations for early forming systems. The tilted, low-sigma_8 model suggested by a new WMAP analysis is rejected at 99.99% confidence, but a model with the same tilt and normalization can be reconciled with the X-ray data by increasing the dark energy equation of state parameter to w ~ -0.8. When imposing the additional constraint of the tight relation between sigma_8 and Omega_m from studies of cluster abundances, the X-ray c-M relation excludes (>99% conf.) both open CDM models and flat CDM models with Omega_m ~1. This result provides novel evidence for a flat, low-Omega_m universe with dark energy using observations only in the local (z << 1) universe. Possible systematic errors in the X-ray mass measurements of a magnitude ~10% suggested by CDM simulations do not change our conclusions.

We report the discovery of a galaxy cluster serendipitously detected as an extended X-ray source in an offset observation of the group NGC 5044. The cluster redshift, z=0.281, determined from the optical spectrum of the brightest cluster galaxy, agrees with that inferred from the X-ray spectrum using the Fe K alpha complex of the hot ICM (z=0.27 +/- 0.01). Based on the 50 ks XMM observation, we find that within a radius of 383 kpc the cluster has an unabsorbed X-ray flux, f_X (0.5-2 keV) = 3.34 (+0.08, -0.13) x 10^{-13} erg/cm^2/s, a bolometric X-ray luminosity, L_X = 2.21 (+0.34, -0.19) x 10^{44} erg/s, kT = 3.57 +/- 0.12 keV, and metallicity, 0.60 +/- 0.09 solar. The cluster obeys the scaling relations for L_X and T observed at intermediate redshift. The mass derived from an isothermal NFW model fit is, M_vir = 3.89 +/- 0.35 x 10^{14} solar masses, with a concentration parameter, c = 6.7 +/- 0.4, consistent with the range of values expected in the concordance cosmological model for relaxed clusters. The optical properties suggest this could be a ``fossil cluster''.

We describe how AGN-produced cosmic rays form large X-ray cavities and radio lobes in the hot diffuse gas in galaxy groups and clusters. Cosmic rays are assumed to be produced in a small shocked region near the cavity center, such as at the working surface of a radio jet. The coupled equations for gasdynamics and cosmic ray diffusion are solved with various assumptions about the diffusion coefficient. To form large, long-lived cavities similar to those observed, the diffusion coefficient must not exceed kappa = 10^28 cm^2/s in the hot gas, very similar to values required in models of cosmic ray diffusion in the Milky Way. When kappa does not exceed 10^28, cosmic rays are confined within the cavities for times comparable to the cavity buoyancy time, as implied by observations of X-ray cavities and their radio synchrotron emission. Collisions of proton cosmic rays with thermal plasma nuclei followed by pion decay can result in enhanced gamma ray emission from the cavity walls.

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.

Within the cosmic framework clusters of galaxies are relatively young objects. Many of them have recently experienced major mergers. Here we investigate an equal mass merging event at z = ~0.6 resulting in a dark matter haloe of ~2.2 times 10^{14} Msol/h at z=0. The merging process is covered by 270 outputs of a high resolution cosmological N-body simulation performed with the ART (adaptive refinement tree) code. Some 2 Gyrs elapse between the first peri-centre passage of the progenitor cores and their final coalescence. During that phase the cores experience six peri-centre passages with minimal distances declining from ~30 to ~2 kpc/h. The time intervals between the peri-centre passages continuously decrease from 9 to 1 times 10^8 yrs. We follow the mean density, the velocity dispersion and the entropy of the two progenitors within a set of fixed proper radii (25, 50, 100, 250, 500, 1000 kpc/h). During the peri-centre passages we find sharp peaks of the mean densities within these radii, which exceed the sum of the corresponding progenitor densities. In addition to the intermixing of the merging haloes, the densities increase due to contraction caused by the momentary deepening of the potential well. At the end of the oscillatory relaxation phase the material originating from the less concentrated of the two equal mass progenitors is deposited at larger radii and shows a slightly more radially anisotropic velocity dispersion compared to the material coming from the more concentrated progenitor. Every peri-centre passage is accompanied by a substantial drop of the central potential well. We briefly discuss the possibility that AGN outbursts are triggered by the periodically changing potential.

The enormous range of X-ray luminosities among elliptical galaxies with similar optical luminosities can be attributed to large scatter and systematic variations in the dark halo mass determined from X-ray observations. The mean halo mass decreases sharply with both the X-ray and K-band luminosities. Smaller halos contain less diffuse hot gas and have lower X-ray luminosities. In addition the hot gas in low mass halos is more easily depleted by supernova-driven outflows, further lowering the X-ray luminosity.

We present an X-ray analysis of the radial mass profile of the radio-quiet galaxy cluster A2589 between 0.015-0.25 r_vir using an XMM-Newton observation. Except for a ~16 kpc shift of the X-ray center of the R=45-60 kpc annulus, A2589 possesses a remarkably symmetrical X-ray image and is therefore an exceptional candidate for precision studies of its mass profile by applying hydrostatic equilibrium. The total gravitating matter profile is well described by the NFW model (fractional residuals <~10%) with c_vir=6.1 +/- 0.3 and M_vir = 3.3 +/- 0.3 x 10^{14} M_sun (r_vir = 1.74 +/- 0.05 Mpc) in excellent agreement with LCDM. When the mass of the hot ICM is subtracted from the gravitating matter profile, the NFW model fitted to the resulting dark matter (DM) profile produces essentially the same result. However, if a component accounting for the stellar mass (M_*) of the cD galaxy is included, then the NFW fit to the DM profile is substantially degraded in the central r ~50 kpc for reasonable M_*/L_V. Modifying the NFW DM halo by adiabatic contraction arising from the early condensation of stellar baryons in the cD galaxy further degrades the fit. The fit is improved substantially with a Sersic-like model recently suggested by high resolution N-body simulations but with an inverse Sersic index, alpha ~0.5, a factor of ~3 higher than predicted. We argue that neither random turbulent motions nor magnetic fields can provide sufficient non-thermal pressure support to reconcile the XMM mass profile with adiabatic contraction of a CDM halo assuming reasonable M_*/L_V. Our results support the scenario where, at least for galaxy clusters, processes during halo formation counteract adiabatic contraction so that the total gravitating mass in the core approximately follows the NFW profile.

We present a Chandra study of mass profiles in 7 elliptical galaxies, of which 3 have galaxy-scale and 4 group-scale halos demarcated at 1E13Msun. These represent the best available data for nearby objects with comparable X-ray luminosities. We measure ~flat mass-to-light (M/L) profiles within an optical half-light radius (Reff), rising by an order of magnitude at ~10Reff, which confirms the presence of dark matter (DM). The data indicate hydrostatic equilibrium, which is also supported by agreement with studies of stellar kinematics in elliptical galaxies. The data are well-fitted by a model comprising an NFW DM profile and a baryonic component following the optical light. The distribution of DM halo concentration parameters (c) versus Mvir agrees with LCDM predictions and our observations of bright groups. Concentrations are slightly higher than expected, which is most likely a selection effect. Omitting the stellar mass drastically increases c, possibly explaining large concentrations found by some past observers. The stellar M/LK agree with population synthesis models, assuming a Kroupa IMF. Allowing adiabatic compression (AC) of the DM halo by baryons made M/L more discrepant, casting some doubt on AC. Our best-fitting models imply total baryon fractions \~0.04--0.09, consistent with models of galaxy formation incorporating strong feedback. The groups exhibit positive temperature gradients, consistent with the "Universal" profiles found in other groups and clusters, whereas the galaxies have negative gradients, suggesting a change in the evolutionary history of the systems around Mvir=1E13 Msun.

We describe 2D gasdynamical models of jets that carry mass as well as energy to the hot gas in galaxy clusters. These flows have many attractive attributes for solving the galaxy cluster cooling flow problem: Why the hot gas temperature and density profiles resemble cooling flows but show no spectral evidence of cooling to low temperatures. Using an approximate model for the cluster A1795, we show that mass-carrying jets can reduce the overall cooling rate to or below the low values implied by X-ray spectra. Biconical subrelativistic jets, described with several ad hoc parameters, are assumed to be activated when gas flows toward or cools near a central supermassive black hole. As the jets proceed out from the center they entrain more and more ambient gas. The jets lose internal pressure by expansion and are compressed by the ambient cluster gas, becoming rather difficult to observe. For a wide variety of initial jet parameters and several feedback scenarios the global cooling can be suppressed for many Gyrs while maintaining cluster temperature profiles similar to those observed. The intermittancy of the feedback generates multiple generations of X-ray cavities similar to those observed in the Perseus Cluster and elsewhere.

Mid-infrared observations (3.6 - 24 microns) of normal giant elliptical galaxies with the Spitzer space telescope are consistent with pure populations of very old stars with no evidence of younger stars. Most of the stars in giant elliptical galaxies are old but the mean stellar age determined from Balmer absorption in optical spectra can appear much younger due to a small admixture of younger stars. The mean stellar age can also be determined from the spectral energy distribution in the mid-infrared which decreases with time relative to the optical emission and shifts to shorter wavelengths. The observed flux ratios F_8um/F_3.6um and F_24um/F_3.6um for elliptical galaxies with the oldest Balmer line ages are lower than predicted by recent models of single stellar populations. For ellipticals with the youngest Balmer line ages in our sample, 3-5 Gyrs, the flux ratios F_24um/F_3.6um are identical to those of the oldest stars. When theoretical mid-IR spectra of old (12 Gyr) and young stellar populations are combined, errors in the F_24um/F_3.6um observations are formally inconsistent with a mass fraction of young stars that exceeds ~1%. This is less than the fraction of young stars expected in discussions of recent surveys of elliptical galaxies at higher redshifts. However, this inconsistancy between Balmer line ages and those inferred from mid-IR observations must be regarded as provisional until more accurate observations and theoretical spectra become available. Finally, there is no evidence to date that central disks or patches of dust commonly visible in optical images of elliptical galaxies contribute sensibly to the mid-IR spectrum.

Elliptical galaxies and their groups having the largest L_x/L_B lie close to the locus in the L_x,L_B diagram expected for closed systems with baryon fractions equal to the cosmic mean value, f_b = 0.16. The estimated baryon fractions for several of these galaxies/groups are also close to 0.16 when the gas density is extrapolated to the virial radius. Evidently they are the least massive baryonically closed systems. Gas retention in these groups implies that non-gravitational heating cannot exceed about 1 keV per particle, consistent with the heating required to produce the deviation of groups from the L_x - T correlation for more massive clusters. Isolated galaxies/groups with X-ray luminosities significantly lower than baryonically closed groups may have undermassive dark halos, overactive central AGNs, or higher star formation efficiencies. The virial mass and hot gas temperatures of nearly or completely closed groups correlate with the group X-ray luminosities and the optical luminosities of the group-centered elliptical galaxy, an expected consequence of their merging history. The ratio of halo mass to the mass of the central galaxy for X-ray luminous galaxy/groups is about 80.

The discovery of extended, approximately spherical weak shock waves in the hot intercluster gas in Perseus and Virgo has precipitated the notion that these waves may be the primary heating process that explains why so little gas cools to low temperatures. This type of heating has received additional support from recent gasdynamical models. We show here that outward propagating, dissipating waves deposit most of their energy near the center of the cluster atmosphere. Consequently, if the gas is heated by (intermittent) weak shocks for several Gyrs, the gas within 30-50 kpc is heated to temperatures that far exceed observed values. This heating can be avoided if dissipating shocks are sufficiently infrequent or weak so as not to be the primary source of global heating. Local PV and viscous heating associated with newly formed X-ray cavities are likely to be small, which is consistent with the low gas temperatures generally observed near the centers of groups and clusters where the cavities are located.

The Chandra and XMM data for a sample of 19 relaxed groups/poor clusters, covering the temperature range 1-3 keV and selected to have the best available data, reveal a remarkable similarity in their temperature profile: cool cores outside of which the temperatures reach a peak for radii less than 0.1 of the virial radius and then decline. We fitted the derived mass profiles using an NFW model, which provides a good fit to the data when accounting for the central galaxy in the inner region. The concentration parameters and virial masses are in the range c = 5-22 and M = 2 x 10^13 - 4 x 10^14 solar masses, in general agreement with the concentrations found in numerical simulations.

Cosmological simulations of galaxy formation predict a universal form for the mass profile of dark matter (DM) halos from cluster to galaxy scales. Remarkably few interesting constraints exist, however, on DM halos in early-type galaxies. Using Chandra we present the temperature, density and mass profiles of a small sample of early-type galaxies, revealing significant DM in each case. When a component is included to account for stellar mass and the DM halo is allowed to respond adiabatically to the baryonic condensation into stars, the mass profiles are well-fitted by the universal profile, with Virial masses and concentrations in agreement with simulations. However, only ~half, or less, of the mass within R_e seems attributable to the stars, implying stellar M_*/L_B ~1-5.

We have analyzed an XMM-Newton observation of the cluster Abell 2589. Apart from a low-level asymmetry in the central region, the cluster appears very relaxed and does not show presence of a central AGN. We derived constraints for the radial temperature, density and, assuming hydrostatic equilibrium, mass profiles. We find that the best fit to the dark matter profile is given by the Sersic-like profile proposed by Navarro et al. (2004). The NFW model does not provide a good fit. We also tested whether the central stellar component could affect the profile through the adiabatic contraction model but were unable to distinguish it from a simple "dark matter + stars" modeling.

The NGC 5044 galaxy group is dominated by a luminous elliptical galaxy which is surrounded by ~160 dwarf satellites. The projected number density profile of this dwarf population deviates within ~1/3 of the virial radius from a projected NFW-profile, which is assumed to approximate the underlying total matter distribution. By means of a semi-analytic model we demonstrate that the interplay between gravitation, dynamical friction and tidal mass loss and destruction can explain the observed number density profile. We use only two parameters in our models: the total to stellar mass fraction of the satellite halos and the disruption efficiency. The disruption efficiency is expressed by a minimum radius. If the tidal radius of the galaxy (halo) falls below this radius the galaxies is assumed to become unobservable. The preferred parameters are an initial total to stellar mass fraction of ~20 and a disruption radius of 4 kpc. In that model about 20% of all the satellites are totally disrupted on their orbits within the group environment. Dynamical friction is less important in shaping the inner slope of the number density profile since the reduction in mass by tidal forces lowers the impact of the friction term. The main destruction mechanism is tide. In the preferred model the total B-band luminosity of all the disrupted galaxies is about twice the observed luminosity of the central elliptical galaxy, indicating that a significant fraction of stars are scattered into the intragroup medium.(Abridged)

In explaining the relative metal abundances observed in galaxy groups and clusters, it is generally assumed that all metals created by supernovae are present either in visible stars or the hot gas. We discuss here the possibility that some of the iron expelled into the hot gas by Type Ia supernovae may have radiatively cooled, avoiding detection by X-ray and optical observers. Hydrodynamic models of Type Ia explosions in the hot gas inside elliptical galaxies result in a gas of nearly pure iron that is several times hotter than the local interstellar gas. We describe the subsequent thermal evolution of the iron-rich gas as it radiates and thermally mixes with the surrounding gas. There is a critical time by which the iron ions must mix into the ambient gas to avoid rapid radiative cooling. We find that successful mixing is possible if the iron ions diffuse with large mean free paths, as in an unmagnetized plasma. However, in microgauss fields the Larmor radii of the iron ions are exceptionally small, so the field geometry must be highly tangled or radial to allow the iron to mix by diffusion faster than it cools by radiative losses. The possibility that some of the supernova iron cools cannot be easily discounted.

The mid-infrared (10-20 um) luminosity of elliptical galaxies is dominated by the integrated emission from circumstellar dust in red giant stars. As a single stellar population evolves, the rate of dusty mass loss from red giant stars decreases with time, so the mid-infrared luminosity should also decline with stellar age. To seek such a correlation, we have used archival ISO observations to determine surface brightness profiles and central fluxes at 15 um in 17 early-type galaxies for which stellar ages have been determined from optical spectral indices. The radial surface brightness distributions at 15 um generally follow the stellar de Vaucouleurs profile as expected. We find that the surface brightness ratio mu_{15um}/mu_{I-band} is systematically higher in elliptical galaxies with ages < 5 Gyrs and in galaxies that exhibit evidence of recent mergers. Within the accuracy of our observations, mu_{15um}/mu_{I-band} shows no age dependence for ages > 5 Gyrs. The corresponding flux ratios F_{15um}/F_{I-band} within apertures scaled to the effective radius (R_e/8) are proportional to the mu_{15um}/mu_{I-band} ratios at larger galactic radii, indicating that no 15 um emission is detected from central dust clouds visible in optical images in some of our sample galaxies. Emission at 15 um is observed in non-central massive clouds of dust and cold gas in NGC1316, an elliptical galaxy that is thought to have had a recent merger. Recent {\it Spitzer Space Telescope} data also indicate the presence of PAH emission at 8 um. Several ellipticals have extended regions of 15 um emission that have no obvious counterparts at other frequencies.

The galaxy group NGC 5044 consists of a luminous giant elliptical galaxy surrounded by a cluster of approximately 160 low luminosity and dwarf galaxies, mostly of early type. The cumulative projected radial distribution of dwarf galaxies in the NGC 5044 group, unlike distributions of more luminous galaxies in rich clusters, does not follow a projected dark matter (NFW) profile. Instead, a deficiency or absence of low luminosity galaxies is apparent in NGC 5044 within about 350 kpc, about one third of the current group virial radius. Most of the dwarf galaxies identified in NGC 5044 entered the virial radius at redshifts less than about 2 or 3, and very few entered at higher redshifts. The peculiar, non-NFW shape of the projected cumulative dwarf galaxy distribution in NGC 5044 within 350 kpc resembles the characteristic shallow cumulative distribution of dark subhalos that are also known to be relatively young. Dynamical friction is unlikely to explain the apparent lack of group member galaxies at small radii in NGC 5044.

We describe a new type of dynamical model for hot gas in galaxy groups and clusters in which gas moves simultaneously in both radial directions. Circulation flows are consistent with (1) the failure to observe cooling gas in X-ray spectra, (2) multiphase gas observed near the centers of these flows and (3) the accumulation of iron in the hot gas from Type Ia supernovae in the central galaxy. Dense inflowing gas cools, producing a positive central temperature gradient, as in normal cooling flows. Bubbles of hot, buoyant gas flow outward. Circulation flows eventually cool catastrophically if the outward flowing gas transports mass but no heat; to maintain the circulation both mass and energy must be supplied to the inflowing gas over a large volume, extending to the cooling radius. The rapid radial recirculation of gas produces a flat central core in the gas iron abundance, similar to many observations. We believe the circulation flows described here are the first gasdynamic, long-term evolutionary models that are in good agreement with all essential features observed in the hot gas: little or no gas cools as required by XMM spectra, the gas temperature increases outward near the center, and the gaseous iron abundance is about solar near the center and decreases outward.

A new XMM observation of the outer regions of the galaxy group NGC 5044 indicates hot gas iron abundances of only ~0.15 solar between r=0.2 - 0.4 r_vir. While the total baryon mass within the virial radius may be close to the cosmic mean value observed in rich clusters, the ratio of total iron mass to optical light in NGC 5044 is about 3 times lower than that in rich clusters. The remarkably low iron abundance over a large volume of the intergroup gas in the outer regions of NGC 5044 cannot be easily understood in terms of the outflow of enriched gas in a group wind during its early history or by the long term enrichment by the group member galaxies that currently occupy this region. It is possible that the stars in NGC 5044 did not produce iron with the same efficiency as in clusters, or the iron resides in non-luminous clouds or stars, or the entropy of the iron-enriched gas created in early galactic starburst winds was too high to penetrate the group gas of lower entropy.

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\mu\rm{m}} / L_B$, but that correlation is weaker for elliptical galaxies.

We review the origin, evolution and physical nature of hot gas in elliptical galaxies and associated galaxy groups. Unanticipated recent X-ray observations with Chandra and XMM indicate much less cooling than previously expected. Consequently, many long-held assumptions need to be reexamined or discarded and new approaches must be explored. Chief among these are the role of heating by active galactic nuclei, the influence of radio lobes on the hot gas, details of the cooling process, possible relation between the hot and colder gas in elliptical galaxies, and the complexities of stellar enrichment of the hot gas.

We draw attention to the close relationship between the anisotropy parameter beta(r) for stellar orbits in elliptical galaxies and the temperature profile T(r) of the hot interstellar gas. For nearly spherical galaxies the gas density can be accurately determined from X-ray observations and the stellar luminosity density can be accurately found from the optical surface brightness. The Jeans equation and hydrostatic equilibrium establish a connection between beta(r) and T(r) that must be consistent with the observed stellar velocity dispersion. Purely optical observations of the bright elliptical galaxy NGC 4472 indicate beta(r) < 0.35 within the effective radius. However, the X-ray gas temperature profile T(r) for NGC 4472 requires significantly larger anisotropy, beta = 0.6 - 0.7, about twice the optical value. This strong preference for radial stellar orbits must be understood in terms of the formation history of massive elliptical galaxies. Conversely, if the smaller, optically determined anisotropy is indeed correct, we are led to the important conclusion that the temperature profile T(r) of the hot interstellar gas in NGC 4472 must differ from that indicated by X-ray observations, or that the hot gas is not in hydrostatic equilibrium.

The failure of the XMM-Newton and Chandra X-ray telescopes to detect cooling gas in elliptical galaxies and clusters of galaxies has led many to adopt the position that the gas is not cooling at all and that heating by an active nucleus in the central E or cD galaxy is sufficient to offset radiative cooling. In this paper we explore an idealized limiting example of this point of view in which hot, buoyant bubbles formed near the center return the inflowing, radiatively cooling gas to distant regions in the flow. We show that idealized steady state, centrally heated non-cooling flows can indeed be constructed. In addition, the emission-weighted temperature profiles in these circulating flows resemble those of normal cooling flows. However, these solutions are valid only (1) for a range of bubble parameters for which there is no independent justification, (2) for a limited spatial region in the cooling flow and (3) for a limited period of time after which cooling seems inevitable. Our exploration of non-cooling flows is set in the context of galaxy/group flows.

Using new XMM and Chandra observations we present an analysis of the metal abundances of the hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. Motivated by the inconsistent abundance and temperature determinations obtained by different observers for X-ray groups, we provide a detailed investigation of the systematic errors on the derived abundances considering the effects of the temperature distribution, calibration, plasma codes, bandwidth, Galactic Nh, and background rate. The iron abundance (Fe) drops from Fe ~1 solar within R ~50 kpc to Fe ~0.4 solar near R=100 kpc. This radial decline in Fe is highly significant: Fe=1.09 +/- 0.04 solar (stat) +/- 0.05 solar + 0.18 solar (sys) within R=48 kpc (5') compared to Fe=0.44 +/- 0.02 solar (stat) +/- 0.10 solar + 0.13 solar (sys) over R=48-96 kpc (5'-10'). The data rule out with high confidence a very sub-solar value for Fe within R=48 kpc confirming that previous claims of very sub-solar central Fe values in NGC 5044 were primarily the result of the Fe Bias: i.e., the incorrect assumption of spatially isothermal and single-phase gas when in fact temperature variations exist. Within R=48 kpc we obtain Si/Fe = 0.83 +/- 0.02 (stat) +/- 0.02 + 0.07 (sys) and S/Fe = 0.54 +/- 0.02 (stat) +/- 0.01 + 0.01 (sys) in solar units. These ratios are consistent with their values at larger radii and imply that SNIa have contributed ~80% of the iron mass within a 100 kpc radius of NGC 5044. This SNIa fraction is similar to the Sun and suggests an IMF similar to that of the Milky Way. At the very center (R ~2 kpc) the XMM and Chandra CCDs and the XMM RGS show that the Fe drops to ~50% of its value at immediately larger radius analogously to that seen in some galaxy clusters. (Abridged)

We propose a stellar origin for the central dust clouds observed in most giant elliptical galaxies. Dusty gas ejected from evolving red giant stars in E or cD galaxies can cool rapidly even after entering the hot, X-ray emitting gas. Cooling by thermal collisions with dust grains can be faster than either the dynamical time in the galactic potential or the grain sputtering time. Some grains survive in the cooled gas. Dusty stellar outflows cool more efficiently in the central regions where the stellar metallicity is higher. Mergers with gas and dust-rich dwarf galaxies may occasionally occur but are not required to explain the observed dust clouds.

Using new XMM and Chandra observations we present an analysis of the temperature structure of the hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. A spectral deprojection analysis of data extracted from circular annuli reveals that a two-temperature model (2T) of the hot gas is favored over single-phase or cooling flow (M_dot = 4.5 +/- 0.2 M_{sun}/yr) models within the central ~30 kpc. Alternatively, the data can be fit equally well if the temperature within each spherical shell varies continuously from ~T_h to T_c ~ T_h/2, but no lower. The high spatial resolution of the Chandra data allows us to determine that the temperature excursion T_h --> T_c required in each shell exceeds the temperature range between the boundaries of the same shell in the best-fitting single-phase model. This is strong evidence for a multi-phase gas having a limited temperature range. The cooler component of the 2T model has a temperature (T_c \~ 0.7 keV) similar to the kinetic temperature of the stars. The hot phase has a temperature (T_h ~ 1.4 keV) characteristic of the virial temperature of the \~10^{13} M_{sun} halo expected in the NGC 5044 group. However, in view of the morphological disturbances and X-ray holes visible in the Chandra image within R ~10 kpc, bubbles of gas heated to ~T_h in this region may be formed by intermittent AGN feedback. Some additional heating at larger radii may be associated with the evolution of the cold front near R ~50 kpc, as suggested by the sharp edge in the EPIC images.

We present far-IR ISO observations of three early-type galaxies in the Virgo cluster. The data were recorded using the ISOPHOT instrument in both the P32 oversampled maps and the P39/37 sparse maps. The maps reach the limiting sensitivity of the ISOPHOT instrument at 60, 90, and 180$\mu$m. Two of the most prominent elliptical galaxies in Virgo -- NGC 4472 and NGC 4649 -- clearly show no emission at far-IR wavelengths at a level of few tens of mJy, but NGC 4636 is detected at all three wavelengths. We have computed the far-IR emission from dust for NGC4472 and NGC4636 under the assumption that dusty outflows from evolving red giant stars are continuously supplying dust to the interstellar medium and that the grains, once diffused into the interstellar medium, are sputtered away by collisions with ions. While the calculated fluxes are consistent with the observed upper limits for NGC4472, the dust emission detected in NGC4636 support the hypothesis that additional dust has been accreted in a very recent ($\lta$ few $10^8$ yrs) merger with a dusty, gas-rich galaxy.

Cluster cooling flow models that include both thermal conduction and AGN heating have lower overall mass cooling rates and simultaneously sustain density and temperature profiles similar to those observed with no ad hoc mass dropout. To achieve this agreement, the thermal conductivity must be about 0.35 +/- 0.10 of the Spitzer value, similar to that advocated by Narayan & Medvedev. However, when applied to galaxy/group scales the synergistic combination of AGN heating and conduction is less satisfactory. When the computed density profile and the global cooling rate are lowered by AGN heating to match observations of these smaller scale flows, the gas temperatures within about 10 kpc are too large. In addition, best-fitting flows in galaxy/groups with AGN heating and thermal conduction require larger conductivities, about 0.5 - 1.0 of the Spitzer value. When the gas density in galaxy/group flows is lowered by heating, the iron enrichment by Type Ia supernovae is more effective and the hot gas iron abundance greatly exceeds observed values. Consequently, galactic/group flows with central heating and conduction require an additional process that removes the iron: failure of Type Ia ejecta to go into the hot phase, selective cooling, etc.

We present ISOPHOT P32 oversampled maps and P37/39 sparse maps, of three bright elliptical galaxies in the Virgo Cluster. The maps reach the limiting sensitivity of the ISOPHOT instrument at 60, 100, 170 and 200 microns. Two elliptical galaxies show no emission at all far-IR ISOPHOT wavelengths at a level of few tens of mJy. The null detection provides a test of the evolution of dust in elliptical galaxies and its size distribution and composition. As previous studies have shown, in many elliptical galaxies both IRAS and ISO have detected mid-IR excess 6-15 micron emission relative to the stellar continuum indicating emission from circumstellar dust. Under the assumption that these dusty outflows from evolving red giant stars and planetary nebulae are continuously supplying dust to the interstellar medium, we have computed the infrared luminosity at the ISOPHOT bands appropriate for NGC4472. The null far-IR ISOPHOT observations exceed the far-IR flux expected from dust expelled from a normal old stellar population.

The density irregularities and holes visible in many Chandra X-ray images of cluster and galactic cooling flows can be produced by symmetrically heated gas near the central galactic black hole. As the heated gas rises away from the galactic center, a relatively small number of large plumes and bubbles are formed in qualitative agreement with the observed features. The expanding centrally heated gas drives a shock into the surrounding gas, displacing it radially. Both computational and analytic results show that the ambient gas near the bubble is cooled by expansion, accounting for the cool rims commonly observed around X-ray holes in cooling flows.

In conventional models of galactic and cluster cooling flows widespread cooling (mass dropout) is assumed to avoid accumulation of unacceptably large central masses. However, recent XMM observations have failed to find spectral evidence for locally cooling gas. This has revived the notion that cooling flows are heated by some process such as an intermittent, low-level AGN involving supermassive black holes in the central galaxy. To explore this hypothesis, we consider the gasdynamical consequences of galactic cooling flows heated by many different scenarios without specifying the detailed physics of the heating process. We are unable to find a single acceptable heated flow in reasonable agreement with well observed hot gas temperature and density profiles, even using finely tuned parameters. Idealized flows in which radiative cooling is perfectly balanced by global heating are grossly incompatible with observations. Flows heated by episodic central feedback generate quasi-cyclic changes in the hot gas density profile which are not supported by current observations. Paradoxically, centrally heated (or pressurized) cooling flows experience spontaneous non-linear compressions that result in spatially widespread cooling instabilities. Therefore, spectral evidence for cooling gas is difficult to avoid by central heating.

X-ray observations of M87 in the Virgo cluster and NGC 4874 in Coma reveal that the gas temperature beyond about 50 kpc from these cD galaxies is comparable to the virial temperature of the cluster, 3 or 9 keV respectively, but within the optical galaxy the temperature drops to the galactic virial temperature, about 1 keV. We show that these steep thermal (and density) gradients follow naturally from the usual cooling inflow assumptions without recourse to thermal conductivity. To avoid unobserved central mass concentrations, most of the gas must radiatively cool (``dropout'') before it flows to the galactic core, i.e. the gas must be multiphase. Recent XMM spectra of M87 indicate single phase flow at every radius with no apparent radiative cooling to low temperatures. However, the X-ray spectral evidence for multiphase cooling beginning at lower temperatures near 1 keV (within about 10 kpc) may be less apparent and may have escaped detection. Finally, we show that the standard decomposition method used by X-ray observers to determine the mass flow dM/dt(r) may fail rather badly due to small deviations from perfectly steady state flow. When applied to our slightly non-steady computed flows, this decomposition (incorrectly) gives the usual result, dM/dt proportional to r, and overestimates the mass dropout dM/dt(r) at large radii.

Was the diffuse gas in galaxy groups and clusters heated at high redshift before it entered a massive halo, or was the heating produced inside collapsed objects by SNII accompanying normal star formation? We compare here two radically different models corresponding to the two scenarios described above. Our results indicate that internal heating by SNII works better than the extreme version of external heating that we adopt in reproducing the observed L-T and entropy-T relations.

Recent XMM-Newton observations of the cooling flow gas in M87 indicate sharply decreasing oxygen, iron and silicon abundances within $\sim 5$ kpc of the galactic center. This result is unexpected since stellar mass loss and Type Ia supernovae are expected to produce pronounced central abundance maxima for all three elements. However, it has been suggested that many of the strong X-ray lines are optically thick and diffuse to larger radii in the cooling flow before escaping, falsifying the central abundances. We verify with radiation transfer calculations that this effect does indeed occur in the M87 cooling flow, but that it is insufficient to account for the M87 observations. We suggest that some source of continuous opacity is required to reduce the central X-ray line emission, perhaps by warm gas at $T \sim 10^5 - 10^6$ K. The radial surface brightness profiles of X-ray resonance lines are also sensitive to turbulence in cooling flows which reduces the line center optical depths considerably. Turbulence may provide sufficient energy to continuously heat the warm absorbing gas.

X-ray observations of hot gas in galaxy groups indicate higher entropies than can be achieved in the accretion shocks as it fell into the dark halos. It has been proposed that this entropy excess results from some universal external heating process in the early universe - AGN, population III stars, etc. - that occurred before most of the gas entered the dark halos. Alternatively, the hot gas in groups may be heated internally by Type II supernovae when the galactic stars formed. We investigate here the gas-dynamical consequences of various amounts of external, high redshift heating and compare with internal supernova heating. Radiative losses and distributed mass dropout are included. We find that externally heated flows are unsatisfactory; when the heating is high enough to bring the X-ray luminosities into agreement with observations, the gas entropy is too high. Flows internally heated by Type II supernovae are preferred. Our results are insensitive to the epoch and spatial distribution of the supernova heating process. However, with a Salpeter IMF a large fraction of the energy produced by supernovae must heat the hot gas. Alternatively, only about 20 percent of the Type II supernova energy would be required to heat the gas if the initial mass function is flatter than Salpeter, as suggested by stellar mass to light ratios.

We review the origin and structure of hot (cooling flow) gas in elliptical galaxies. X-ray observations can be used to determine the stellar mass to light ratio, the mass profiles of dark matter halos, and the interstellar magnetic field. Interstellar gas cools over a large volume, forming stars with a bottom-heavy IMF. For consistency with the thin fundamental plane, young stars must be optically luminous. Circular X-ray isophotes in rotating elliptical galaxies indicate distributed radiative cooling or strong interstellar turbulence.

We explore the possibly important influence of cooled interstellar gas on the fundamental plane of elliptical galaxies. Interstellar cooling is described by a parameterized sink term in the equation of continuity. Parameters that give the best fits to the X-ray observations of NGC 4472 are used as a template for the radial distribution of interstellar cooling in structurally homologous elliptical galaxies of lower mass. Gas that cools within an effective radius can contribute an additional 10 - 30 percent to the mass of the old stellar population. If the cooled gas forms into stars of very low mass, $\ll M_{\odot}$, as is commonly assumed, the cooled mass is optically dark. As a result, the mass to light ratios determined from stellar velocities systematically overestimate that of the old stellar population. Moreover, the total mass and spatial distribution of the optically dark young stellar population does not scale homologously with galactic luminosity or radius and the total stellar mass to light ratio varies with galactic radius. We investigate the non-homologous perturbations of cooled gas on the mass to light ratio for several idealized homologous elliptical galaxies and show that they appear to be incompatible with the observed thinness of the fundamental plane. If optically luminous young stars formed from the cooled gas, the disturbance of the fundamental plane would be lessened.

If mass and angular momentum were conserved in cooling flows associated with luminous, slowly rotating elliptical galaxies, the inflowing hot gas would spin up, resulting in disks of cold gas and X-ray images that are highly flattened along the equatorial plane out to several effective radii. Such X-ray flattening is not observed at the spatial resolution currently available to X-ray observations. Evidently mass and angular momentum are not in fact conserved. If cooling flows are depleted by localized radiative cooling at numerous sites distributed throughout the flows, then disks of cooled gas do not form and the X-ray images appear nearly circular. However, the distribution of young stars formed from the cooled gas is still somewhat flattened relative to the stellar light. X-ray images of galactic cooling flows can also be circularized by the turbulent diffusion of angular momentum away from the axis of rotation, but the effective viscosity of known processes -- stellar mass loss, supernovae, cooling site evolution, etc. -- is insufficient to appreciably circularize the X-ray images. Radial gradients in the interstellar iron abundance are unaffected by the expected level of interstellar turbulence since these gradients are continuously re-established by Type Ia supernovae.

Although only about 5 percent of the total baryonic mass in luminous elliptical galaxies is in the form of cooled interstellar gas, it is concentrated within the optical effective radius r_e where it influences the local dynamical mass. The mass of cooled gas must be spatially distributed since it greatly exceeds the masses of central black holes. We explore here the proposition that a population of low mass, optically dark stars is created from the cooled gas. We consider a wide variety of radial distributions for the interstellar cooling, but only a few are consistent with observed X-ray surface brightness profiles. In a region of concentrated interstellar cooling, the X-ray emission can exceed that observed, suggesting the presence of additional support by magnetic stresses or non-thermal pressure. In general we find that the mass of cooled gas contributes significantly to stellar dynamical mass to light ratios which vary with galactic radius. If the stars formed from cooled interstellar gas are optically luminous, their influence on the the mass to light ratio would be reduced. The mass of cooled gas inside r_e is sensitive to the rate that old stars lose mass, which is nearly independent of the initial mass function of the old stellar population.

Anticipating forthcoming observations with the Chandra X-ray telescope, we describe the continuation of interstellar cooling flows deep into the cores of elliptical galaxies. Interstellar gas within about r = 50 parsecs from the massive black hole is heated to T > 1 keV and should be visible unless thermal heating is diluted by non-thermal pressure. Since our flows are subsonic near the massive black holes, distributed cooling continues within 300 pc from the center. Dark, low mass stars formed in this region may be responsible for some of the mass attributed to central black holes.

Observational constraints and theoretical arguments indicate that cooled interstellar gas in bright elliptical galaxies forms into a young stellar population throughout the region within the half-light radius. The young population has a bottom-heavy, but optically luminous IMF extending to 1 - 2 M_sun. When the colors and spectral features of this young population are combined with those of the underlying old stellar population, the apparent ages are significantly reduced, similar to the relatively young apparent ages observed in many ellipticals. Galactic mergers are not required to resupply young stars. The sensitivity of continuous star formation to L_B and L_x/L_B is likely to account for the observed spread in apparent ages among elliptical galaxies. Local star formation is accompanied by enhanced stellar H_beta equivalent widths, stronger optical emission lines, enhanced thermal X-ray emission and lower apparent temperatures in the hot gas. The young stars should cause M/L to vary with galactic radius, perturbing the fundamental plane occupied by the old stars.

X-ray emission from hot (T = 10^7 K) interstellar gas in massive elliptical galaxies indicates that 10^{10} M_sun has cooled over a Hubble time, but optical and radio evidence for this cold gas is lacking. We provide detailed theoretical support for the hypothesis that this gas has formed into low luminosity stars. Within several kpc of the galactic center, interstellar gas first cools to T = 10^4 K where it is heated by stellar UV and emits the observed diffuse optical line emission. This cooling occurs at a large number (10^6) of isolated sites. After less than a solar mass of gas has accumulated (10^{-6} M_sun/yr) at a typical cooling site, a neutral (HI or H_2) core develops in the HII cloud where gas temperatures drop to T = 15 K and the ionization level (from thermal X-rays) is very low (x = 10^{-6}). We show that the maximum mass of cores that become gravitationally unstable is only about 2 M_sun. No star can exceed this mass. Fragmentation of collapsing cores produces a population of low mass stars with a bottom-heavy IMF and radial orbits. Gravitational collapse and ambipolar diffusion are rapid. The total mass of star-forming (dust-free) HI or H_2 cores in a typical bright elliptical is only 10^6 M_sun, below current observational thresholds.

Hot, X-ray emitting gaseous halos around massive elliptical galaxies are a result of both stellar mass loss and inflow toward the overdensity from which giant ellipticals and their associated galaxy groups formed. The metal abundance in this gas contains important information about early star formation and past supernova activity. We find that Type II supernovae based on a Salpeter IMF, plus a small number of additional Type Ia supernovae, can explain the the density, temperature and abundance profiles currently observed in gaseous halos around massive ellipticals. Within the central, optically bright regions of luminous ellipticals, approximately half of the interstellar iron is produced by Type Ia supernovae and half by mass lost from evolving stars which were originally enriched by Type II supernovae. However, iron and silicon abundances in the intracluster gas within rich clusters suggest enrichment by a larger number of supernovae per unit optical light than we require for massive ellipticals, either more Type Ia or more Type II from a flatter IMF. Since the enrichment histories of massive ellipticals and rich clusters are fundamentally different, E and SO galaxies may not be the only sources of metal enrichment in rich cluster gas.

We study the complete dynamical evolution of hot interstellar gas in massive elliptical galaxies born into a simple flat universe beginning with an overdense perturbation. Within the turn-around radius dark matter flows in a self-similar fashion into a stationary Navarro-Frenk-White halo and the baryonic gas shocks. After a few gigayears, when enough gas accumulates within the accretion shock, the de Vaucouleurs stellar system is constructed and the energy from Type II supernovae is released. The stars and dark halo are matched to NGC 4472. Gas continues to enter the galaxy by secondary infall and by stellar mass loss based on a Salpeter IMF. After about 13 Gyrs the temperature and density distribution in the hot gas agree quite well with the hot interstellar gas observed in NGC 4472. As a result of supernova-driven outflow, the present day baryonic fraction has a deep minimum in the outer galactic halo. When relatively gas-rich, X-ray luminous models are spatially truncated at early times, simulating tidal events that may have occurred during galaxy group dynamics, the current locus of truncated models lies just along the $L_x$, X-ray size correlation among well-observed ellipticals, providing another striking confirmation of our simple model of elliptical evolution.

The variation of temperature and density in the hot, X-ray emitting gas around massive, group dominant elliptical galaxies can be understood as a combination of gas ejected from evolving galactic stars and gas that accumulates in the outer halo by secondary cosmic infall. Beginning with an overdensity perturbation in a simple flat cosmology, we can duplicate observed properties of the hot gas. At some early time we form the stellar galaxy and release supernova energy, conserving dark and baryonic matter. Using a potential appropriate to the large elliptical NGC 4472, we follow the evolution of intergalactic and interstellar gas to the present time when the computed gas density and temperature agree with X-ray observations of NGC 4472. If the hot gas and dark matter halos are subject to differential tidal truncations or mass exchanges between group members, then the correlation between $L_x/L_B$ and the relative sizes of galactic X-ray images can be generated. The physical properties of hot interstellar gas observed in bright ellipticals today are sensitive to the cosmic baryon fraction, the time of maximum star formation and the amount of ``feedback'' energy delivered to the gas by Type II supernovae at the epoch of galaxy formation.

Optical parameters of elliptical galaxies are tightly correlated, but their x-ray parameters vary widely. The x-ray luminosity L_x ranges over more than an order of magnitude for ellipticals having similar optical luminosity L_B. The source of this scatter has been elusive. We show here that the dispersion in L_x for fixed optical luminosity L_B correlates strongly with the dimensionless ratio of the sizes of the x-ray and optical images, r_ex/r_e. Specifically, we find that (L_x/L_B) is proportional to (r_{ex}/r_e)^{0.60 \pm 0.30}, a version of the correlation that is independent of distance. This correlation may be a natural result of mergings and tidal truncations that are expected during the formation and early evolution of ellipticals in groups of galaxies. The radial structure of x-ray images also varies: some are compact (e.g. NGC 4649, 7626, 5044), others diffuse (e.g. NGC 4636, 1399).

We review recent progress in understanding the evolution of hot interstellar gas in the halos of elliptical galaxies. Non-homologous variations in the physical size of the hot gas account for the large variations of x-ray luminosity among ellipticals of similar L_B. By combining ROSAT and Einstein data we derive the distribution of total mass in NGC 4472. From 0.1r_e to 1r_e the total mass is identical to the expected stellar mass. Therefore stellar mass to light ratios can be determined from x-ray observations! Also the widely used ``mass dropout'' assumption must be incorrect in this important part of the cooling flow. Recent ROSAT observations indicate gas temperatures in excess of the virial stellar temperature, totally unlike standard cooling flow models. However, these new results can be understood if an additional massive component of ``circumgalactic'' gas is assumed to fill the outer galactic halos beyond most of the stars. This old hot gas, first heated during the epoch of galaxy formation, continues to flow into the stellar parts of ellipticals today, combining with gas expelled from evolving stars. This dual origin of hot interstellar gas further complicates recent discussions of abundances in the hot interstellar gas.

X-ray emitting gaseous halos, such as that in elliptical galaxies like NGC 4472, cannot have been produced solely from gas expelled from galactic stars. In traditional models for the evolution of hot interstellar gas (cooling flows) in ellipticals, the galaxies are assumed to have been cleared of gas by SNII-driven winds at some early time then gas is subsequently replenished by mass loss from an evolving population of old stars. To test this, we accurately determine the stellar and dark halo mass of NGC 4472 using hydrostatic equilibrium, then solve the standard time-dependent cooling flow equations to recover the observed hot gas temperature and density distributions when evolved to the present time. This procedure fails: the computed gas density gradient is too steep, the total gas mass is too low, and the gas temperatures are much too low. All variants on this basic procedure also fail: increasing the SNIa rate, using the mass dropout assumption, arbitrarily adjusting uncertain coefficients, etc. However, agreement is achieved if the galaxy is supplied with additional, spatially-extended hot gas early in its evolution. This old ``circumgalactic'' gas can be retained to the present time and may be related to cosmological ``secondary infall''.

We compare optical and high quality x-ray data for three bright ellipticals in the Virgo cluster, NGC 4472, 4649, and 4636. The distribution of total mass in NGC 4472 and 4649 determined from x-ray data is sensitive to the stellar mass over a considerable range in galactic radius extending to r = r_e, the effective radius. The agreement of x-ray and optically-determined stellar masses provides a unique verification of the stellar mass to light ratio which is essentially constant over 0.1 < r/r_e < 1. However, dark matter in NGC 4636 is important at all radii > 0.35r_e. Evidently the dark to stellar mass ratio varies in quite different ways in ellipticals of comparable optical luminosity, implying the radial structure of dark halos may not be universal. There is some evidence in NGC 4636 for additional support of the hot interstellar gas at r < 0.35 r_e; either a field B = 10^{-4} G or a small (but mechanically unstable) central region of high gas temperature (T = 10^7 K) is required. The global temperature structure in the hot interstellar medium of many recently observed ellipticals is very similar, reaching a maximum near 3 - 4 r_e. This feature, which may suggest a new structural scale, is inconsistent with current theoretical gas dynamical models.

We describe the evolution of hot interstellar gas in a family of low luminosity elliptical galaxies all having $M_B = -20$ but with different degrees of flattening (E0, E2, and E6) caused by rotation. The interstellar gas is furnished by stellar mass loss. The soft X-ray luminosity Lx of the hot interstellar gas after 15 Gyrs decreases dramatically with galactic rotation. The hot gas cools onto a large disk having a size that depends inversely on the assumed supernova rate. Although Lx of the hot, rotating interstellar gas may be difficult to observe, the photoionized part of the cooled gaseous disk gas can be observed in optical line emission. The mass of HII gas ($\sim 10^8$ $M_{\odot}$) may be much greater than previously realized since large regions of rotationally supported, low density HII are optically faint. We interpret the stellar disks that are common in low luminosity ellipticals as stars that have formed in the non-ionized part of the cold disk gas. The total mass and distribution of cold disk gas available for star formation are similar to those of stellar disks observed. The high stellar H$\beta$ photometric index observed in disky ellipticals can be understood by combining the light of young disk stellar populations with that of the old bulge population.

Interstellar magnetic fields in elliptical galaxies are assumed to have their origin in stellar fields that accompany normal mass loss from an evolving population of old stars. The seed fields are amplified by interstellar turbulence driven by stellar mass loss and supernova events. These disordered fields are further amplified by time-dependent compression in the inward moving galactic cooling flow and are expected to dominate near the galactic core. Under favorable circumstances, fields similar in strength to those observed $B \sim 1-10~(r/10~kpc)^{-1.2}\mu$G can be generated solely from these natural galactic processes. In general the interstellar field throughout elliptical galaxies is determined by the outermost regions in the interstellar gas where the turbulent dynamo process can occur. Because of the long hydrodynamic flow times in galactic cooling flows, currently observed magnetic fields may result from periods of intense turbulent field amplification that occurred in the outer galaxy in the distant past. Particularly strong fields in ellipticals may result from ancient galactic mergers or shear turbulence introduced at the boundary between the interstellar gas and ambient cluster gas.

The temperature of hot interstellar gas at large radii in elliptical galaxies can be lower than the mean galactic virial temperature. If so, a nonlinear cooling wave can form in the hot interstellar gas and propagate slowly toward the galactic core. If the cooling wave survives hydrodynamic instabilities, it can intermittently deposit cold gas within about 15 effective radii. For a bright elliptical the total mass deposited in this manner can approach 10^10 solar masses. The cold gas that drips out at large galactic radii may account for the young stellar populations and extended gas at $\sim 10^4$ K observed in many ellipticals, features that are often attributed to galactic mergers. Galactic drips are expected in relatively isolated (field) ellipticals provided (i) the galactic stellar velocity ellipsoids are radially oriented at large galactic radii and (ii) the current Type Ia supernova rate is sufficiently small to be consistent with interstellar iron abundances found in recent X-ray studies. Galactic drips are surpressed in ellipticals located within clusters of galaxies; when the pressure in the ambient cluster gas exceeds that in the outer parts of the galactic interstellar medium, some cluster gas flows into the galaxy which surpresses the drips.

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.