Parallel Tree–SPH: A Tool for Galaxy Formation

We describe a new implementation of a parallel Tree-SPH code with the aim of simulating galaxy formation and evolution. The code has been parallelized using SHMEM, a Cray proprietary library to handle communications between the 256 processors of the Silicon Graphics T3E massively parallel supercomputer hosted by the Cineca Super-computing Center (Bologna, Italy). The code combines the smoothed particle hydrodynamics (SPH) method to solve hydrodynamical equations with the popular Barnes and Hut (1986) tree-code to perform gravity calculation with a N × log N scaling, and it is based on the scalar Tree-SPH code developed by Carraro et al. (1998). Parallelization is achieved by distributing particles along processors according to a workload criterion. Benchmarks of the code, in terms of load balance and scalability, are analysed and critically discussed against the adiabatic collapse of an isothermal gas sphere test using 2 × 10⁴ particles on eight processors. The code turns out to be balanced at more than 95% level. If the number of processors is increased, the load balance worsens slightly. The deviation from perfect scalability at increasing number of processors is negligible up to 64 processors. Additionally we have incorporated radiative cooling, star formation, feedback and an algorithm to follow the chemical enrichment of the interstellar medium. This revised version was published online in July 2006 with corrections to the Cover Date.     

ASCA spectroscopy of the luminous infrared galaxy NGC 6240: X-ray emission from a starburst and a buried active nucleus

We present an X-ray spectroscopic study of the prototype far-infrared galaxy NGC 6240 from ASCA . The soft X-ray spectrum (below 2 keV) shows clear signatures of thermal emission well described by a multitemperature optically thin plasma, which probably originates in a powerful starburst. Strong hard X-ray emission is also detected with ASCA and its spectrum above 3 keV is extremely flat with a prominent iron K line complex, very similar to that seen in the Seyfert 2 galaxy NGC 1068 but about an order of magnitude more luminous ( L _3−10keV ≈ 1.4 × 10⁴² erg s⁻¹). The hard X-ray spectrum indicates that only reflected X-rays of an active galactic nucleus (AGN) buried in a heavy obscuration ( N _H > 2 × 10²⁴ cm⁻²) are visible. This is evidence for an AGN in NGC 6240, emitting possibly at a quasar luminosity (∼ 10⁴⁵ erg s⁻¹), and suggests its significant contribution to the far-infrared luminosity.     

X-ray analysis of the cluster 3C 129

The cluster 3C 129 is classified as a rich cluster. An analysis of the properties of the cluster 3C 129 from ROSAT PSPC and HRI, Einstein IPC, and EXOSAT ME observations is presented. The mean temperature from a joint fit of the ROSAT PSPC and EXOSAT ME data is 5.5(±0.2) keV. The luminosity is 0.6×10⁴⁴ erg s⁻¹ in 0.2–2.4 keV and 2.7×10⁴⁴ erg s⁻¹ in 0.2–10 keV. We find a cooling flow with a rate of ∼84 M_⊙ yr⁻¹. The central gas density is 6×10⁻³ cm⁻³, and the ICM mass is 3.6×10¹³ M_⊙. The total cluster mass is ∼5×10¹⁴ M_⊙. The X-ray morphology shows an east–west elongation, which is evidence for a recent merger event. The radio source 3C 129.1 is located near the X-ray centre. Another cluster member galaxy (the radio galaxy 3C 129) is a prototype of head-tailed radio galaxies, and is located in the west part of the cluster. The tail points along the gradient of intracluster gas pressure. There are no significant point X-ray sources associated with the AGNs of the two radio galaxies.     

A Chandra detection of the radio hotspot of 3C 123

Chandra X-ray Observatory observations of the powerful, peculiar radio galaxy 3C 123 have resulted in an X-ray detection of the bright eastern hotspot, with a 1-keV flux density of ∼5 nJy. The X-ray flux and spectrum of the hotspot are consistent with the X-rays being inverse-Compton scattering of radio synchrotron photons by the population of electrons responsible for the radio emission ('synchrotron self-Compton emission') if the magnetic fields in the hotspot are close to their equipartition values. 3C 123 is thus the third radio galaxy to show X-ray emission from a hotspot which is consistent with being in equipartition. Chandra also detects emission from a moderately rich cluster surrounding 3C 123, with L _X(2–10 keV)=2×10⁴⁴ erg s⁻¹ and kT ∼5 keV, and absorbed emission from the active nucleus, with an inferred intrinsic column density of 1.7×10²² cm⁻² and an intrinsic 2–10 keV luminosity of 10⁴⁴ erg s⁻¹.     

The anomalous properties of Markarian 1460

We present and discuss optical, near-infrared and H  i measurements of the galaxy Markarian 1460 at a distance of 19 Mpc in the Ursa Major Cluster. This low-luminosity ( M _B =−14) galaxy is unusual because (i) it is blue ( B − R =0.8) and has the spectrum of an H  ii galaxy, (ii) it has a light profile that is smooth and well fitted by an r ^1/4 and not an exponential function at all radii larger than the seeing, and (iii) it has an observed central brightness of about μ _B =20 mag arcsec⁻² , intermediate between those of elliptical galaxies (on the bright μ _B side) and normal low-luminosity dwarf irregular (on the low μ _B side) galaxies. No other known galaxy exhibits all these properties in conjunction. On morphological grounds this galaxy looks like a normal distant luminous elliptical galaxy, since the Fundamental Plane tells us that higher luminosity normal elliptical galaxies tend to have lower surface-brightnesses. Markarian 1460 has 2×10⁷ M_⊙ of H  i and a ratio M (H  i )/ L _B of 0.2, which is low compared to the typical values for star-forming dwarf galaxies. From the high surface-brightness and r ^1/4 profile, we infer that the baryonic component of Markarian 1460 has become self-gravitating through dissipative processes. From the colours, radio continuum, H  i and optical emission line properties, and yet smooth texture, we infer that Markarian 1460 has had significant star formation as recently as ∼1 Gyr ago but not today.     

Australia Telescope Compact Array H i observations of the NGC 6845 galaxy group

We present the results of Australia Telescope Compact Array (ATCA) H  i line and 20-cm radio continuum observations of the galaxy quartet NGC 6845. The H  i emission extends over all four galaxies but can only be associated clearly with the two spiral galaxies, NGC 6845A and B, which show signs of strong tidal interaction. We derive a total H  i mass of at least  1.8 × 10¹⁰ M_⊙  , most of which is associated with NGC 6845A, the largest galaxy of the group. We investigate the tidal interaction between NGC 6845A and B by studying the kinematics of distinct H  i components and their relation to the known H  ii regions. No H  i emission is detected from the two lenticular galaxies, NGC 6845C and D. A previously uncatalogued dwarf galaxy, ATCA  J2001−4659  , was detected 4.4 arcmin NE from NGC 6845B and has an H  i mass of  ∼5 × 10⁸ M_⊙  . No H  i bridge is visible between the group and its newly detected companion. Extended 20-cm radio continuum emission is detected in NGC 6845A and B as well as in the tidal bridge between the two galaxies. We derive star formation rates of  15–40 M_⊙ yr⁻¹  .     

ASCA observations of deep ROSAT fields — III. The discovery of an obscured Type 2 AGN at z = 0.67

We report on the discovery of a narrow-emission-line object at z  = 0.672 detected in a deep ASCA survey. The object, AXJ 0341.4–4453, has a flux in the 2–10 keV band of 1.1 ± 0.27 × 10⁻¹³ erg s⁻¹ cm⁻², corresponding to a luminosity of 1.8 × 10⁴⁴ erg s⁻¹ ( q ₀ = 0.5, H ₀ = 50 km s⁻¹ Mpc⁻¹). It is also marginally detected in the ROSAT 0.5–2 keV band with a flux 5.8 × 10⁻¹⁵ erg s⁻¹ cm⁻². Both the ASCA data alone and the combined ROSAT/ASCA data show a very hard X-ray spectrum, consistent with either a flat power law (α < 0.1) or photoelectric absorption with a column of n _H > 4 × 10²² cm⁻² (α = 1). The optical spectrum shows the high-ionization, narrow emission lines typical of a Seyfert 2 galaxy. We suggest that this object may be typical of the hard sources required to explain the remainder of the X-ray background at hard energies.     

Self-similar shocked accretion of collisional gas with radiative cooling

We describe similarity solutions that characterize the collapse of collisional gas on to scale-free perturbations in an Einstein–de Sitter universe. We consider the effects of radiative cooling and derive self-similar solutions under the assumption that the cooling function is a power law of density and temperature, Λ( T , ρ )∝ ρ ^3/2 T . We use these results to test the ability of smooth particle hydrodynamics (SPH) techniques to follow the collapse and accretion of shocked, rapidly cooling gas in a cosmological context. Our SPH code reproduces the analytical results very well in cases that include or exclude radiative cooling. No substantial deviations from the predicted central mass accretion rates or from the temperature, density and velocity profiles are observed in well-resolved regions inside the shock radius. This test problem lends support to the reliability of SPH techniques to model the complex process of galaxy formation.     

The cooling function of HD molecule revisited

We report new calculations of the cooling rate of primordial gas by the HD molecule, taking into account its ro-vibrational structure. The HD cooling function is calculated including radiative and collisional transitions for   J ≤ 8  rotational levels, and for the vibrational levels v = 0, 1, 2 and 3. The ro-vibrational level population is calculated from the balance equation assuming steady state. The cooling function is evaluated in the ranges of the kinetic temperatures, T _k, from 10² to  2 × 10⁴ K  and the number densities, n _H, from 1 to  10⁸ cm⁻³  . We find that the inclusion of collisional ro-vibrational transitions increases significantly the HD cooling efficiency, in particular for high densities and temperatures. For   n _H≳ 10⁵  and   T _k∼ 10⁴ K  the cooling function becomes more than an order of magnitude higher than previously reported. We give also the HD cooling rate in the presence of the cosmic microwave radiation field for radiation temperatures of 30, 85 and 276 K (redshifts of 10, 30 and 100). The tabulated cooling functions are available at http://www.cifus.uson.mx/Personal_Pages/anton/DATA/HD_cooling/HD_cool.html . We discuss the relevance to explore the effects of including our results into models and simulations of galaxy formation, especially in the regime when gas cools down from temperatures above ∼3000 K.     

Free–free absorption in the Seyfert nucleus of NGC 4151

We present subarcsecond MERLIN 0.4-GHz (73 cm) and 1.6-GHz (18 cm) radio measurements of the nuclear region of the Seyfert galaxy NGC 4151. By comparison with higher frequency observations, we deduce that one component (C4) shows a low-frequency turnover which we interpret as evidence for free–free absorption by ionized gas with an emission measure between 3 × 10⁵ and 10⁶ pc cm⁻⁶. The free–free absorption appears to be localized to a region ∼50 pc in extent, and we consider models in which the ionized gas may be closely associated with a neutral molecular torus.     

A Chandra observation of the distant radio galaxy B2 0902+343: a powerful obscured active galaxy

We report a Chandra observation of the   z =3.395  radio galaxy B2 0902+343. The unresolved X-ray source is centred on the active nucleus. The spectrum is well fitted by a flat power law of photon index of  Γ∼1.1  with intrinsic absorption of  8×10²² cm^-2  , and an intrinsic  2–10 keV  luminosity of  3.3×10⁴⁵ erg s^-1  . More complex models that allow for a steeper spectral index cause the column density and intrinsic luminosity to increase. The data limit any thermal luminosity of the hot magnetized medium, assumed responsible for high Faraday rotation measures seen in the radio source, to less than ∼10⁴⁵ erg s⁻¹.     

A ultraluminous X-ray source associated with a cloud collision in M 99

The Sc galaxy M 99 in the Virgo Cluster has been strongly affected by tidal interactions and recent close encounters, responsible for an asymmetric spiral pattern and a high star formation rate. Our XMM–Newton study shows that the inner disc is dominated by hot plasma at kT ≈ 0.30 keV, with a total X-ray luminosity of ≈10⁴¹ erg s⁻¹ in the 0.3–12 keV band. At the outskirts of the galaxy, away from the main star-forming regions, there is an ultraluminous X-ray source (ULX) with an X-ray luminosity of ≈2 × 10⁴⁰ erg s⁻¹ and a hard spectrum well fitted by a power law of photon index Γ≈ 1.7. This source is close to the location where a massive H  i cloud appears to be falling on to the M 99 disc at a relative speed of >100 km s⁻¹. We suggest that there may be a direct physical link between fast cloud collisions and the formation of bright ULXs, which may be powered by accreting black holes with masses ∼100 M_⊙. External collisions may trigger large-scale dynamical collapses of protoclusters, leading to the formation of very massive (≳200 M_⊙) stellar progenitors; we argue that such stars may later collapse into massive black holes if their metal abundance is sufficiently low.     

Chandra unveils a binary active galactic nucleus in Mrk 463

We analyse Chandra , XMM–Newton and Hubble Space Telescope ( HST ) data of the double-nucleus Ultraluminous Infrared Galaxy (ULIRG), Mrk 463. The Chandra detection of two luminous  ( L _2–10 keV= 1.5 × 10⁴³  and  3.8 × 10⁴² erg cm⁻² s⁻¹)  , unresolved nuclei in Mrk 463 indicates that this galaxy hosts a binary active galactic nucleus (AGN), with a projected separation of ≃3.8 kpc (  3.83 ± 0.01  arcsec). While the East nucleus was already known to be a type 2 Seyfert (and this is further confirmed by our Chandra detection of a neutral iron line), this is the first unambiguous evidence in favour of the AGN nature of the West nucleus. Mrk 463 is therefore the clearest case so far for a binary AGN, after NGC 6240.     

The Durham/UKST Galaxy Redshift Survey – III. Large-scale structure via the two-point correlation function

We have investigated the statistical clustering properties of galaxies by calculating the two-point galaxy correlation function from the optically selected Durham/UKST Galaxy Redshift Survey. This survey is magnitude-limited to b _J∼17, contains ∼2500 galaxies sampled at a rate of one-in-three and surveys a ∼4×10⁶ ( h ⁻¹ Mpc)³ volume of space. We have empirically determined the optimal method of estimating the two-point correlation function from just such a magnitude-limited survey. Applying our methods to this survey, we find that our redshift-space results agree well with those from previous optical surveys. In particular, we confirm the previously claimed detections of large-scale power out to ∼40 h ⁻¹ Mpc scales. We compare with two common models of cosmological structure formation and find that our two-point correlation function has power significantly in excess of the standard cold dark matter model in the 10–30 h ⁻¹ Mpc region. We therefore support the observational results of the APM galaxy survey. Given that only the redshift-space clustering can be measured directly, we use standard modelling methods and indirectly estimate the real-space two-point correlation function from the projected two-point correlation function. We then invert this projected correlation function to obtain an estimate of the spatial two-point correlation function in real space. This correlation function in real space has a lower amplitude than that in redshift space, but a steeper slope.     

Understanding the halo-mass and galaxy-mass cross-correlation functions

We use the Millennium Simulation (MS) to measure the cross-correlation between halo centres and mass (or equivalently the average density profiles of dark haloes) in a Lambda cold dark matter (ΛCDM) cosmology. We present results for radii in the range  10  h ⁻¹ kpc < r < 30  h ⁻¹ Mpc  and for halo masses in the range  4 × 10¹⁰ < M ₂₀₀ < 4 × 10¹⁴  h ⁻¹ M_⊙  . Both at   z = 0  and at   z = 0.76  these cross-correlations are surprisingly well fitted if the inner region is approximated by a density profile of NFW or Einasto form, the outer region by a biased version of the linear mass autocorrelation function, and the maximum of the two is adopted where they are comparable. We use a simulation of galaxy formation within the MS to explore how these results are reflected in cross-correlations between galaxies and mass. These are directly observable through galaxy–galaxy lensing. Here also we find that simple models can represent the simulation results remarkably well, typically to ≲10 per cent. Such models can be used to extend our results to other redshifts, to cosmologies with other parameters, and to other assumptions about how galaxies populate dark haloes. Our galaxy formation simulation already reproduces current galaxy–galaxy lensing data quite well. The characteristic features predicted in the galaxy–galaxy lensing signal should provide a strong test of the ΛCDM cosmology as well as a route to understanding how galaxies form within it.     

Diffuse X-ray emission from late-type galaxy haloes

Current theories of galaxy formation predict that spiral galaxies are embedded in a reservoir of hot gas. This gas is able to cool on to the galaxy, replenishing cold gas that is consumed by star formation. Estimates of the X-ray luminosity emitted in the cooling region suggest a bolometric luminosity of the order of 10×10⁴¹ erg s⁻¹ in massive systems. We have used ROSAT PSPC data to search for extended X-ray emission from the haloes of three nearby, massive, late-type galaxies: NGC 2841, 4594 and 5529. We infer 95 per cent upper limits on the bolometric X-ray luminosities of the haloes of NGC 2841, 4594 and 5529 of 0.4, 1.2 and 3.8×10⁴¹ erg s⁻¹ respectively. Thus, the true luminosity lies well below the straightforward theoretical prediction. We discuss this discrepancy and suggest a number of ways in which the theoretical model might be brought into agreement with the observational results. A possible solution is that the gravitational potentials of the dark matter haloes of these galaxies are weaker than assumed in the current model. Alternatively, the present-day accretion may be substantially less than is required on average to build the disc over the Hubble time. Our results are, however, based on only three galaxies, none of which is ideal for this kind of study. A larger data set is required to explore this important problem further.     

Galactic winds and circulation of the interstellar medium in dwarf galaxies

We study, through 2D hydrodynamical simulations, the feedback of a starburst on the ISM of typical gas-rich dwarf galaxies. The main goal is to address the circulation of the ISM and metals following the starburst. We assume a single-phase rotating ISM in equilibrium in the galactic potential generated by a stellar disc and a spherical dark halo. The starburst is assumed to occur in a small volume in the centre of the galaxy, and it generates a mechanical power of 3.8×10³⁹ or 3.8×10⁴⁰ erg s⁻¹ for 30 Myr. We find, in accordance with previous investigations, that the galactic wind is not very effective in removing the ISM. The metal-rich stellar ejecta, however, can be efficiently expelled from the galaxy and dispersed in the intergalactic medium.
Moreover, we find that the central region of the galaxy is always replenished with cold and dense gas a few 100 million years after the starburst, achieving the requisite for a new star formation event in ≈0.5–1 Gyr. The hydrodynamical evolution of galactic winds is thus consistent with the episodic star formation regime suggested by many chemical evolution studies.
We also discuss the X-ray emission of these galaxies and find that the observable (emission-averaged) abundance of the hot gas underestimates the real one if thermal conduction is effective. This could explain the very low hot-gas metallicities estimated in starburst galaxies.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.     

Neutral hydrogen absorption observations of the central region of NGC 5929

We present 0.15-arcsec (25-pc) resolution MERLIN observations of neutral hydrogen absorption detected towards the nuclear region of the type 2 Seyfert galaxy NGC 5929. Absorption is detected only towards the north-eastern radio component with a column density of (6.5 ± 0.6) × 10²¹ cm⁻². Based on comparison with an HST WFPC2 continuum image, we propose that the absorption is caused by a 1.5-arcsec structure of neutral gas and dust offset 0.3 arcsec south-east of the nucleus and running NE–SW. A separate cloud of dust is apparent 1.5 arcsec to the south-west of the nucleus in the HST image. A comparison of the centroid velocity (2358 ± 5 km s⁻¹) and full width at half-maximum (43 ± 6 km s⁻¹) of the absorbing gas with previous [O  III ] observations suggests that both the neutral and ionized gas are undergoing galactic rotation towards the observer in the north-east and away from the observer in the south-west. The main structure is consistent with an inclined ring of gas and dust encircling the active galactic nucleus (AGN); alternatively it may be a bar or inner spiral arm. We do not detect neutral hydrogen absorption or dust obscuration against the radio nucleus (column density < 3.1 × 10²¹ cm⁻²) expected by a torus of neutral gas and dust in unified models of AGNs for a type 2 Seyfert galaxy.     

Simulating galactic outflows with kinetic supernova feedback

Feedback from star formation is thought to play a key role in the formation and evolution of galaxies, but its implementation in cosmological simulations is currently hampered by a lack of numerical resolution. We present and test a subgrid recipe to model feedback from massive stars in cosmological smoothed particle hydrodynamics simulations. The energy is distributed in kinetic form among the gas particles surrounding recently formed stars. The impact of the feedback is studied using a suite of high-resolution simulations of isolated disc galaxies embedded in dark haloes with total mass 10¹⁰ and  10¹²  h ⁻¹ M_⊙  . We focus, in particular, on the effect of pressure forces on wind particles within the disc, which we turn off temporarily in some of our runs to mimic a recipe that has been widely used in the literature. We find that this popular recipe gives dramatically different results because (ram) pressure forces on expanding superbubbles determine both the structure of the disc and the development of large-scale outflows. Pressure forces exerted by expanding superbubbles puff up the disc, giving the dwarf galaxy an irregular morphology and creating a galactic fountain in the massive galaxy. Hydrodynamic drag within the disc results in a strong increase in the effective mass loading of the wind for the dwarf galaxy, but quenches much of the outflow in the case of the high-mass galaxy.