Time evolution of galactic warps in prolate haloes

A recent observation with the Hipparcos satellite and some numerical simulations imply that the interaction between an oblate halo and a disc is inappropriate for the persistence of galactic warps. Following on from this , we have compared the time evolution of galactic warps in a prolate halo with that in an oblate halo. The haloes were approximated as fixed potentials, while the discs were represented by N -body particles. We have found that the warping in the oblate halo continues to wind up, and finally disappears. On the other hand, for the prolate halo model, the precession rate of the outer disc increases when the precession of the outer disc recedes from that of the inner disc, and vice versa. Consequently, the warping in the prolate halo persisted to the end of the simulation by retaining the alignment of the line of nodes of the warped disc. Therefore, our results suggest that prolate haloes could sustain galactic warps. The physical mechanism of the persistence of warp is discussed on the basis of the torque between a halo and a disc and that between the inner and outer regions of the disc.     

Large-scale structure, the cosmic microwave background and primordial non-Gaussianity

Cosmic microwave background and large-scale structure data will shortly improve dramatically with the Microwave Anisotropy Probe and Planck Surveyor , and the Anglo-Australian 2-Degree Field and Sloan Digital Sky Survey. It is therefore timely to ask which of the microwave background and large-scale structure will provide a better probe of primordial non-Gaussianity. In this paper we consider this question, using the bispectrum as a discriminating statistic. We consider several non-Gaussian models and find that in each case the microwave background will provide a better probe of primordial non-Gaussianity. Our results suggest that if microwave background maps appear Gaussian, then apparent deviations from Gaussian initial conditions in galaxy surveys can be attributed with confidence to the effects of biasing. We demonstrate this precisely for the spatial bispectrum induced by local non-linear biasing.     

Approximate third integrals for axisymmetric potentials using local Stäckel fits

We use a set of Stäckel potentials to obtain a local approximation for an effective third integral in axisymmetric systems. We present a study on the feasibility and effectiveness of this approach. We apply it to three trial potentials of various flattenings, corresponding to nearly ellipsoidal, discy and boxy density isophotes. In all three cases, a good fit to the potential requires only a small set of Stäckel potentials, and the associated Stäckel third integral provides a very satisfactory, yet analytically simple, approximation to the trial potentials effective third integral.     

Galaxy destruction and diffuse light in clusters

Deep images of the Centaurus and Coma clusters reveal two spectacular arcs of diffuse light that stretch for over 100 kpc, yet are just a few kiloparsecs wide. At a surface brightness of m _b ∼27–28 mag arcsec⁻² , the Centaurus arc is the most striking example known of structure in the diffuse light component of a rich galaxy cluster. We use numerical simulations to show that the Centaurus feature can be reproduced by the tidal debris of a spiral galaxy that has been tidally disrupted by the gravitational potential of NGC 4709. The surface brightness and narrow dimensions of the diffuse light suggest that the disc was corotating with its orbital path past pericentre. Features this prominent in clusters will be relatively rare, although at fainter surface brightness levels the diffuse light will reveal a wealth of structure. Deeper imaging surveys may be able to trace this feature for several times its presently observed extent, and somewhere along the tidal debris, a fraction of the original stellar component of the disc will remain bound, but transformed into a faint spheroidal galaxy. It should be possible to confirm the galactic origin of the Centaurus arc by observing planetary nebulae along its length with redshifts close to that of NGC 4709.     

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.     

Non-linear hydrodynamical evolution of rotating relativistic stars: numerical methods and code tests

We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, non-linear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes (based on approximate Riemann solvers) and compare their accuracy in preserving uniformly rotating stationary initial configurations in long-term evolutions. Among these four schemes, we find that the third-order piecewise parabolic method scheme is superior in maintaining the initial rotation law in long-term evolutions, especially near the surface of the star. It is further shown that HRSC schemes are suitable for the evolution of perturbed neutron stars and for the accurate identification (via Fourier transforms) of normal modes of oscillation. This is demonstrated for radial and quadrupolar pulsations in the non-rotating limit, where we find good agreement with frequencies obtained with a linear perturbation code. The code can be used for studying small-amplitude or non-linear pulsations of differentially rotating neutron stars, while our present results serve as testbed computations for three-dimensional general-relativistic evolution codes.     

Interstellar CO towards X Persei (HD 24534) – II. Two-component model

The observations made by the Goddard High Resolution Spectrograph (GHRS) aboard the Hubble Space Telescope ( HST ) of molecular CO in absorbing gas towards X Persei are reported. The two-component statistical equilibrium model incorporating radiative excitation of CO by line emission at the same velocity that originates in nearby molecular clouds has been used to reproduce high-resolution GHRS spectra. Earlier analysis indicates that the cloud has a complex structure and at least a two-component model should be used to obtain accurate results. The spectra obtained from the International Ultraviolet Explorer ( IUE ) were used to complement GHRS data and constrain the space of possible solutions. The new oscillator strengths recommended by Eidelsberg et al. for A–X bands have been used. The results show that one of the components may be attributed to the Perseus OB2 molecular cloud, and the other component to an extension of the Taurus dark cloud. The total CO column density N (CO)=(1.0±0.2)×10¹⁶ cm⁻² has been determined. According to the results about 85 per cent of the observed CO belongs to an extension of the Taurus dark cloud. The CO radiation that originates in nearby molecular clouds may be the dominant excitation mechanism of the observed CO. The early results of ¹³CO line analysis indicate a ¹³CO/¹²CO ratio of about 40.     

Ionizing radiation in smoothed particle hydrodynamics

A new method for the inclusion of ionizing radiation from uniform radiation fields into 3D smoothed particle hydrodynamics (SPHI) simulations is presented. We calculate the optical depth for the Lyman continuum radiation from the source towards the SPHI particles by ray-tracing integration. The time-dependent ionization rate equation is then solved locally for the particles within the ionizing radiation field. Using test calculations, we explore the numerical behaviour of the code with respect to the implementation of the time-dependent ionization rate equation. We also test the coupling of the heating caused by the ionization to the hydrodynamical part of the SPHI code.     

Isotropic wavelets: a powerful tool to extract point sources from cosmic microwave background maps

It is the aim of this paper to introduce the use of isotropic wavelets to detect and determine the flux of point sources appearing in cosmic microwave background (CMB) maps. The most suitable wavelet to detect point sources filtered with a Gaussian beam is the 'Mexican Hat'. An analytical expression of the wavelet coefficient obtained in the presence of a point source is provided and used in the detection and flux estimation methods presented. For illustration the method is applied to two simulations (assuming Planck mission characteristics) dominated by CMB (100 GHz) and dust (857 GHz), as these will be the two signals dominating at low and high frequencies respectively in the Planck channels. We are able to detect bright sources above 1.58 Jy at 857 GHz (82 per cent of all sources) and above 0.36 Jy at 100 GHz (100 per cent of all), with errors in the flux estimation below 25 per cent. The main advantage of this method is that nothing has to be assumed about the underlying field, i.e. about the nature and properties of the signal plus noise present in the maps. This is not the case in the detection method presented by Tegmark & Oliveira-Costa. Both methods are compared, producing similar results.