Observational limits to source confusion in the millimetre/submillimetre waveband

The first observations to detect a population of distant galaxies directly in the submillimetre waveband have recently been made using the new Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). The results indicate that a large number of distant galaxies are radiating strongly in this waveband. Here we discuss their significance for source confusion in future millimetre/submillimetre-wave observations of both distant galaxies and cosmic microwave background radiation (CMBR) anisotropies. Earlier estimates of such confusion involved significant extrapolation of the results of observations of galaxies at low redshifts; our new estimates do not, as they are derived from direct observations of distant galaxies in the submillimetre waveband. The results have important consequences for the design and operation of existing and proposed millimetre/submillimetre-wave telescopes: the Planck Surveyor survey will be confusion-limited at frequencies greater than 350 GHz, even in the absence of Galactic dust emission; a 1σ confusion noise limit of about 0.44 mJy beam⁻¹ is expected for the JCMT/SCUBA at a wavelength of 850 μm; and the subarcsecond resolution of large millimetre/submillimetre-wave interferometer arrays will be required in order to execute very deep galaxy surveys.     

Clustering of galaxies at high redshifts

Recent observations show a large concentration of galaxies at high redshift. At first sight, strong clustering of galaxies at high redshifts seems to be in contradiction with the models of structure formation. In this paper we show that such structures are a manifestation of the strong clustering of rare peaks in the density field. We compute the frequency of occurrence of such large concentrations of galaxies in some models of structure formation.     

On the physical connections between galaxies of different types

Galaxies can be classified in two broad sequences which are likely to reflect their formation mechanism. The 'main sequence', consisting of spirals, irregulars and all dwarf galaxies, is probably produced by gas settling within dark matter haloes. We show that the sizes and surface densities along this sequence are primarily determined by the distributions of the angular momentum and formation time of dark haloes. They are well reproduced by current cosmogonies provided that galaxies form late, at z  ≲ 2. In this scenario, dwarf ellipticals were small 'discs' at z  ∼ 1 and become 'ellipticals' after they fall into cluster environments. The strong clustering of dwarf ellipticals is then a natural by-product of the merging and transformation process. The number of dwarf galaxies predicted in a cluster such as Virgo is in good agreement with the observed number. On the other hand, the 'giant branch', consisting of giant ellipticals and bulges, is probably produced by the merging of disc galaxies. Based on the observed phase-space densities of galaxies, we show that the main bodies of all giant ellipticals can be produced by dissipationless mergers of high-redshift discs. However, high-redshift discs, although denser than present-day ones, are still not compact enough to produce the high central phase-space density of some low-luminosity ellipticals. Dissipation must have occurred in the central parts of these galaxies during the merger which formed them.     

On a time-symmetric Hermite integrator for planetary N-body simulation

We describe a P(EC) ⁿ Hermite scheme for planetary N -body simulation. The fourth-order implicit Hermite scheme is a time-symmetric integrator that has no secular energy error for the integration of periodic orbits with time-symmetric time-steps. In general N -body problems, however, this advantage is of little practical significance, since it is difficult to achieve time-symmetry with individual variable time-steps. However, we can easily enjoy the benefit of the time-symmetric Hermite integrator in planetary N -body systems, where all bodies spend most of the time on nearly circular orbits. These orbits are integrated with almost constant time-steps even if we adopt the individual time-step scheme. The P(EC) ⁿ Hermite scheme and almost constant time-steps reduce the integration error greatly. For example, the energy error of the P(EC)² Hermite scheme is two orders of magnitude smaller than that of the standard PEC Hermite scheme in the case of an N  = 100,  m  = 10²⁵ g planetesimal system with the rms eccentricity 〈 e ²〉^1/2 ≲0.03.     

Equivalent width, shape and proper motion of the iron fluorescent line emission from molecular clouds as an indicator of the illuminating source X-ray flux history

Observations of the diffuse emission in the 8–22 keV energy range, elongated parallel to the Galactic plane, and detection of the strong 6.4-keV fluorescent line with ∼ 1 keV equivalent width from some giant molecular clouds (e.g. Sgr B2) in the Galactic Centre region suggest that the neutral matter of these clouds is (or was) illuminated by powerful X-ray radiation, which gave rise to the reprocessed radiation. The source of this radiation remains unknown. A transient source close to the Sgr B2 cloud, or a short outburst of the X-ray emission from a supermassive black hole at the Galactic Centre are the two prime candidates under consideration. We argue that a new generation of X-ray telescopes combining very high sensitivity and excellent energy and angular resolutions would be able to discriminate between these two possibilities when studying time-dependent changes of the morphology of the surface brightness distribution, the equivalent width and the shape of the fluorescent line in Sgr B2 and other molecular clouds in the region. We note also that detection of broad and complex structures near the 6.4-keV line in the spectra of distant AGNs, which are X-ray weak now, may prove the presence of violent activity in the central engines of these objects in the past. Accurate measurements of the line shape may provide information on the time elapsed since the outburst. Proper motion (super- or subluminal) of the fluorescent radiation wave front can give additional information on the location of the source. Observations of the described effects can provide unique information on the matter distribution inside Sgr B2 and other giant molecular clouds.     

Formation and evolution of elliptical galaxies and QSO activity

We present the results of a numerical code that combines multi-zone chemical evolution with 1D hydrodynamics to follow in detail the evolution and radial behaviour of gas and stars during the formation of elliptical galaxies. We use the model to explore the links between the evolution and formation of elliptical galaxies and QSO activity. The knowledge of the radial gas flows in the galaxy allows us to trace metallicity gradients, and, in particular, the formation of a high-metallicity core in ellipticals. The high-metallicity core is formed soon enough to explain the metal abundances inferred in high-redshift quasars. The star formation rate and the subsequent feedback regulate the episodes of wind, outflow and cooling flow, thus affecting the recycling of the gas and the chemical enrichment of the intergalactic medium. The evolution of the galaxy shows several stages, some of which are characterized by a complex flow pattern, with inflow in some regions and outflow in other regions. All models, however, exhibit during their late evolution a galactic wind at the outer boundary and, during their early evolution, an inflow towards the galactic nucleus. The characteristics of the inner inflow could explain the bolometric luminosity of a quasar lodged at the galactic centre as well as the evolution of the optical luminosity of quasars.     

Collapse and fragmentation of rotating, prolate clouds

We present the results of collapse calculations for uniformly rotating, prolate clouds performed using the numerical method: smoothed particle hydrodynamics (SPH). The clouds considered are isothermal, prolate spheroids with different axial ratios ( a/b ), and with different values of β, the ratio of the rotational to gravitational energy. Small density perturbations are added to the clouds, and different initial perturbation spectra are studied. All of the clouds considered are strongly unstable to gravitational contraction, and so collapse to form a spindle configuration. Such a linear structure is unstable to fragmentation, so that the clouds break up into a number of subcondensations. The long-term evolution of the system is then determined by the angular momentum possessed by these fragments.
It is found that a number of the calculations performed result in the formation of orbitally stable binary systems, composed of two rotationally supported discs in orbit about their common centre of mass. Tidal interactions during closest approach, close three-body interactions and the continued accretion of material with high specific angular momentum are all found to increase the orbital separation during these calculations, ensuring that the systems do not merge at later times. The calculations are therefore relevant to the problem of binary star formation, though the systems produced tend to have large orbital separations and periods. One of the strong points of the models presented, however, is their ability to produce systems with a range of mass ratios and orbital eccentricities, without the explicit inclusion of biases in the initial conditions.     

Galaxy formation and evolution – I. The Padua tree-sph code (pd-sph)

In this paper we report on pd-sph , the new tree-sph code developed in Padua. The main features of the code are described and the results of a new and independent series of 1D and 3D tests are shown. The paper is mainly dedicated to the presentation of the code and to the critical discussion of its performance. In particular, great attention is devoted to the convergency analysis. The code is highly adaptive in space and time by means of individual smoothing lengths and individual time-steps. At present it contains both dark and baryonic matter, this latter in the form of gas and stars, cooling, thermal conduction, star formation, feedback from Type I and II supernovae, stellar winds, and ultraviolet flux from massive stars, and finally chemical enrichment. New cooling rates that depend on the metal abundance of the interstellar medium are employed, and the differences with respect to the standard ones are outlined. Finally, we show the simulation of the dynamical and chemical evolution of a disc-like galaxy with and without feedback. The code is suitably designed to study in a global fashion the problem of formation and evolution of elliptical galaxies, and in particular to feed a spectrophotometric code from which the integrated spectra, magnitudes and colours (together with their spatial gradients) can be derived.     

A numerical algorithm for dissipative Keplerian particle discs

We present a numerical algorithm designed to study the evolution of a distribution of solid particles orbiting around the Sun that could be applied to similar systems in which f _c ≤ Ω/2π, where f _c is the frequency of collisions, and Ω is the orbital angular speed. A number of sample particles are used to represent the whole system. Binary collisions are treated in a novel way using a Monte Carlo method that works as follows. Orbits are locally sampled to compute the velocity dispersion. Then the velocity vectors of the sample particles are modified according to random collisions with virtual particles which have velocities taken from a normal distribution computed using the previously found local velocity dispersion. The energy and momentum taken up by the virtual particles are redistributed among the neighbours of the sample particle undergoing the collision, so that conservation laws are satisfied. Simulations using this model give an estimation of the final distribution of inclinations and the associated evolutionary time-scale.     

Gas cooling by dust during dynamical fragmentation

We suggest that the abrupt switch, from hierarchical clustering on scales ≳ 0.04 pc, to binary (and occasionally higher multiple) systems on smaller scales, which Larson has deduced from his analysis of the grouping of pre-main-sequence stars in Taurus, arises because pre-protostellar gas becomes thermally coupled to dust at sufficiently high densities. The resulting change — from gas cooling by molecular lines at low densities to gas cooling by dust at high densities — enables the matter to radiate much more efficiently, and hence to undergo dynamical fragmentation. We derive the domain in which gas cooling by dust facilitates dynamical fragmentation. Low-mass (∼ M⊙) clumps — those supported mainly by thermal pressure — can probably access this domain spontaneously, albeit rather quasi-statically, provided that they exist in a region in which external perturbations are few and far between. More massive clumps probably require an impulsive external perturbation, for instance a supersonic collision with another clump, in order for the gas to reach sufficiently high density to couple thermally to the dust. Impulsive external perturbations should promote fragmentation, by generating highly non-linear substructures which can then be amplified by gravity during the subsequent collapse.