Prolate Jaffe models for galaxies

We introduce a class of prolate Jaffe models for elliptical galaxies, which are a further extension of Jaffe's spherical models of axisymmetric elliptical systems, and study the properties of their densities, circular velocities, velocity dispersions and two-integral even distribution functions. The form of the potential allows the density to be expressed simply as a function of the potential and radial coordinate R . The models have finite total mass and their densities at large distances decay radially as r ⁻⁴, except on the major axis, where the densities decay as r ⁻³. It is known from Hunter's formulae that the velocity dispersions for prolate models can be expressed in terms of elementary functions of R and z , unlike those for the oblate Jaffe models recently given by Jiang, and that the prolate models have anisotropic velocity distributions. Thus the prolate models are easier to study than the oblate models. It is also found that the two-integral even distribution functions on the physical boundary of the galaxies increase monotonically with the relative energy, for the prolate models. Furthermore, numerical calculation shows that the two-integral even distribution functions generated from their densities are non-negative, even for very 'squeezed' prolate Jaffe models. However, the edge-on projected surface densities for these prolate models cannot be expressed as simply as for the oblate models.     

Application of a clumpy core model to NGC 2071

The multi-transitional observations of CS molecules towards the NGC 2071 core have been re-analysed by using a tri-dimensional Monte Carlo radiative transfer code. Better agreement with the observations is made by an introduction of clumpiness to this model than by smoothly varying density to the 1D microturbulent one. The best-fitting model shows that, when a unique density is assumed for clumps, the volume filling factor of the clumps varies as r ⁻² with an average of ∼5 per cent over the entire core, and that the H₂ number density and the CS abundance of the clump relative to H₂ are ∼ 2 × 10⁶ cm⁻³ and ∼ 6 × 10⁻¹⁰, respectively. The radial density gradient ∝ r ⁻² obtained from our clumpy core model is steeper than that (∝ r ^−1.3) obtained from the microturbulent model. Since all clumps are subject to random bulk motions in this 3D clumpy macroturbulent model, synthesized line profiles do not show self-absorption dips even for opaque transitions and the resulting linewidth is in good accordance with the observations.     

N-body simulations for testing the stability of triaxial galaxies in MOND

We perform a stability test of triaxial models in Modified Newtonian Dynamics (MOND) using N -body simulations. The triaxial models considered here have densities that vary with   r ⁻¹  in the centre and   r ⁻⁴  at large radii. The total mass of the model varies from 10⁸ to  10¹⁰ M_⊙  , representing the mass scale of dwarfs to medium-mass elliptical galaxies, respectively, from deep MOND to quasi-Newtonian gravity. We build triaxial galaxy models using the Schwarzschild technique, and evolve the systems for 200 Keplerian dynamical times (at the typical length-scale of 1.0 kpc). We find that the systems are virial overheating, and in quasi-equilibrium with the relaxation taking approximately 5 Keplerian dynamical times (1.0 kpc). For all systems, the change of the inertial (kinetic) energy is less than 10 per cent (20 per cent) after relaxation. However, the central profile of the model is flattened during the relaxation and the (overall) axis ratios change by roughly 10 per cent within 200 Keplerian dynamical times (at 1.0 kpc) in our simulations. We further find that the systems are stable once they reach the equilibrium state.     

The origin and formation of cuspy density profiles through violent relaxation of stellar systems

It is shown that the cuspy density distributions observed in the cores of elliptical galaxies can be realized by dissipationless gravitational collapse. The initial models consist of power-law density spheres such as ρ ∝ r ⁻¹ with anisotropic velocity dispersions. Collapse simulations are carried out by integrating the collisionless Boltzmann equation directly, on the assumption of spherical symmetry. From the results obtained, the extent of constant density cores, formed through violent relaxation, decreases as the velocity anisotropy increases radially, and practically disappears for extremely radially anisotropic models. As a result, the relaxed density distributions become more cuspy with increasing radial velocity anisotropy. It is thus concluded that the velocity anisotropy could be a key ingredient for the formation of density cusps in a dissipationless collapse picture. The velocity dispersions increase with radius in the cores according to the nearly power-law density distributions. The power-law index, n , of the density profiles, defined as ρ ∝ r ^{− n} , changes from n ≈2.1 at intermediate radii to a shallower power than n ≈2.1 toward the centre. This density bend can be explained from our postulated local phase-space constraint that the phase-space density accessible to the relaxed state is determined at each radius by the maximum phase-space density of the initial state.     

Low-frequency gravitational waves from cosmological compact binaries

We consider gravitational waves emitted by various populations of compact binaries at cosmological distances. We use population synthesis models to characterize the properties of double neutron stars, double black holes and double white dwarf binaries, and white dwarf–neutron star, white dwarf–black hole and black hole–neutron star systems.
We use the observationally determined cosmic star formation history to reconstruct the redshift distribution of these sources and their merging rate evolution.
The gravitational signals emitted by each source during its early spiralling in phase add randomly to produce a stochastic background in the low-frequency band with spectral strain amplitude between ~10⁻¹⁸ and ~5×10⁻¹⁷ Hz^−1/2 at frequencies in the interval ~5×10⁻⁶–5×10⁻⁵ Hz.
The overall signal, which at frequencies above 10⁻⁴ Hz is largely dominated by double white dwarf systems, might be detectable with LISA in the frequency range 1–10 mHz and acts like a confusion-limited noise component, which might limit the LISA sensitivity at frequencies above 1 mHz.     

Dynamics of the boxy elliptical galaxy NGC 1600

We use three-integral models to infer the distribution function (DF) of the boxy E3–E4 galaxy NGC 1600 from surface brightness and line-profile data on the minor and major axes. We assume axisymmetry and that the mass-to-light ratio is constant in the central ∼1 R _e. Stars in the resulting gravitational potential move mainly on regular orbits. We use an approximate third integral K from perturbation theory and write the DF as a sum of basis functions in the three integrals E , L _z and K . We then fit the projected moments of these basis functions to the kinematic observables and deprojected density, using a non-parametric algorithm. The deduced dynamical structure is radially anisotropic, with σ _θ σ _r ≈ σ _φ σ _r ≈0.7 on the major axis. Both on the minor axis and near the centre the velocity distribution is more isotropic; thus the model is flattened by equatorial radial orbits. The kinematic data are fitted without the need for a central black hole; the central mass determined previously from ground-based data therefore overestimates the actual black-hole mass. The mass-to-light ratio of the stars is M L _V =6  h ₅₀. The anisotropy structure of NGC 1600 with a radially anisotropic main body and more nearly isotropic centre is similar to that found recently in NGC 1399, 2434, 3379 and 6703, suggesting that this pattern may be common amongst massive elliptical galaxies. We discuss a possible merger origin of NGC 1600 in the light of these results.     

The size and shape of local voids

We study the size and shape of low-density regions in the local Universe, which we identify in the smoothed density field of the PSCz flux-limited IRAS galaxy catalogue. After quantifying the systematic biases that enter the detection of voids using our data set and method, we identify, using a smoothing length of 5  h ⁻¹ Mpc, 14 voids within 80  h ⁻¹ Mpc (having volumes 10³  h ⁻³ Mpc³) and, using a smoothing length of 10  h ⁻¹ Mpc, eight voids within 130  h ⁻¹ Mpc (having volumes  8×10³ h⁻³ Mpc³)  . We study the void size distribution and morphologies and find that there is roughly an equal number of prolate and oblate-like spheroidal voids. We compare the measured PSCz void shape and size distributions with those expected in six different cold dark matter (CDM) models and find that only the size distribution can discriminate between models. The models preferred by the PSCz data are those with intermediate values of   σ ₈(≃0.83)  , independent of cosmology.     

Large-Scale Coherent Dipole Anisotropy?

We have reanalysed and compared the dipoles of the 1.2-Jy and 0.6-Jy (QDOT) IRAS galaxy samples. We find strong indications from both samples for (a) significant contributions to the gravitational field that shapes the Local Group motion from depths up to ∼ 170 h ⁻¹ Mpc and (b) a large-scale coherence of the dipole anisotropy, indications provided mainly by the fact that the differential dipoles of large equal-volume shells are aligned with the CMB dipole and exhibit significant dipole signals. The two IRAS dipoles are indistinguishable within 50 h ⁻¹ Mpc and beyond ∼ 130 h ⁻¹ Mpc while the QDOT dipole, having a lower flux limit, continues growing with respect to the 1.2-Jy sample up to ∼ 100 h ⁻¹ Mpc in agreement with Rowan-Robinson et al.     

Helical fields and filamentary molecular clouds – I

We study the equilibrium of pressure truncated, filamentary molecular clouds that are threaded by rather general helical magnetic fields. We first apply the virial theorem to filamentary molecular clouds, including the effects of non-thermal motions and the turbulent pressure of the surrounding ISM. When compared with the data, we find that many filamentary clouds have a mass per unit length that is significantly reduced by the effects of external pressure, and that toroidal fields play a significant role in squeezing such clouds.
We also develop exact numerical MHD models of filamentary molecular clouds with more general helical field configurations than have previously been considered. We examine the effects of the equation of state by comparing 'isothermal' filaments, with constant total (thermal plus turbulent) velocity dispersion, with equilibria constructed using a logatropic equation of state.
Our theoretical models involve three parameters: two to describe the mass loading of the toroidal and poloidal fields, and a third that describes the radial concentration of the filament. We thoroughly explore our parameter space to determine which choices of parameters result in models that agree with the available observational constraints. We find that both equations of state result in equilibria that agree with the observational results. Moreover, we find that models with helical fields have more realistic density profiles than either unmagnetized models or those with purely poloidal fields; we find that most isothermal models have density distributions that fall off as r ^−1.8 to r ⁻², while logatropes have density profiles that range from r ⁻¹ to r ^−1.8. We find that purely poloidal fields produce filaments with steep radial density gradients that are not allowed by the observations.     

Calculated spectra for HeH+ and its effect on the opacity of cool metal-poor stars

The wavelength and Einstein A coefficient are calculated for all rotation–vibration transitions of  ⁴He¹H⁺, ³He¹H⁺, ⁴He²H⁺  and  ³ He²H⁺  , giving a complete line list and the partition function for  ⁴HeH⁺  and its isotopologues. This opacity is included in the calculation of the total opacity of low-metallicity stars and its effect is analysed for different conditions of temperature, density and hydrogen number fraction. For a low helium number fraction (as in the Sun), it is found that HeH⁺ has a visible but small effect for very low densities  (ρ≤ 10⁻¹⁰ g cm⁻³)  , at temperatures around 3500 K. However, for high helium number fraction, the effect of HeH⁺ becomes important for higher densities  (ρ≤ 10⁻⁶ g cm⁻³)  , its effect being most important for a temperature around 3500 K. Synthetic spectra for a variety of different conditions are presented.     

Clustering of galaxies at 3.6 μm in the Spitzer Wide-area Infrared Extragalactic legacy survey

We investigate the clustering of galaxies selected in the 3.6 μm band of the Spitzer Wide-area Infrared Extragalactic (SWIRE) legacy survey. The angular two-point correlation function is calculated for 11 samples with flux limits of S _3.6≥ 4–400 μJy, over an 8 deg² field. The angular clustering strength is measured at >5σ significance at all flux limits, with amplitudes of A = (0.49–29) × 10⁻³ at 1°, for a power-law model, A θ^−0.8. We estimate the redshift distributions of the samples using phenomological models, simulations and photometric redshifts, and so derive the spatial correlation lengths. We compare our results with the Galaxies In Cosmological Simulations (GalICS) models of galaxy evolution and with parametrized models of clustering evolution. The GalICS simulations are consistent with our angular correlation functions, but fail to match the spatial clustering inferred from the phenomological models or the photometric redshifts. We find that the uncertainties in the redshift distributions of our samples dominate the statistical errors in our estimates of the spatial clustering. At low redshifts (median z ≤ 0.5), the comoving correlation length is approximately constant,   r ₀= 6.1 ± 0.5  h ⁻¹  Mpc, and then decreases with increasing redshift to a value of 2.9 ± 0.3  h ⁻¹ Mpc for the faintest sample, for which the median redshift is z ∼ 1. We suggest that this trend can be attributed to a decrease in the average galaxy and halo mass in the fainter flux-limited samples, corresponding to changes in the relative numbers of early- and late-type galaxies. However, we cannot rule out strong evolution of the correlation length over  0.5 < z < 1  .     

The dark matter distribution in disc galaxies

We use high-quality optical rotation curves of nine low-luminosity disc galaxies to obtain the velocity profiles of the surrounding dark matter haloes. We find that they increase linearly with radius at least out to the edge of the stellar disc, implying that, over the entire stellar region, the density of the dark halo is about constant.
The properties of the mass structure of these haloes are similar to those found for a number of dwarf and low surface brightness galaxies, but provide a more substantial evidence of the discrepancy between the halo mass distribution predicted in the cold dark matter scenario and those actually detected around galaxies. We find that the density law proposed by Burkert reproduces the halo rotation curves, with halo central densities ( ρ ₀∼1–4×10⁻²⁴ g cm⁻³) and core radii ( r ₀∼5–15 kpc) scaling as ρ ₀∝ r ₀^−2/3.     

The 2dF Galaxy Redshift Survey: the nature of the relative bias between galaxies of different spectral type

We present an analysis of the relative bias between early- and late-type galaxies in the Two-degree Field Galaxy Redshift Survey (2dFGRS) – as defined by the η parameter of Madgwick et al., which quantifies the spectral type of galaxies in the survey. We calculate counts in cells for flux-limited samples of early- and late-type galaxies, using approximately cubical cells with sides ranging from 7 to  42 h ⁻¹ Mpc  . We measure the variance of the counts in cells using the method of Efstathiou et al., which we find requires a correction for a finite volume effect equivalent to the integral constraint bias of the autocorrelation function. Using a maximum-likelihood technique we fit lognormal models to the one-point density distribution, and develop methods of dealing with biases in the recovered variances resulting from this technique. We then examine the joint density distribution function,   f (δ_E, δ_L)  , and directly fit deterministic bias models to the joint counts in cells. We measure a linear relative bias of ≈1.3, which does not vary significantly with ℓ. A deterministic linear bias model is, however, a poor approximation to the data, especially on small scales  (ℓ≤ 28  h ⁻¹ Mpc)  where deterministic linear bias is excluded at high significance. A power-law bias model with index   b ₁≈ 0.75  is a significantly better fit to the data on all scales, although linear bias becomes consistent with the data for  ℓ≳ 40  h ⁻¹ Mpc  .     

Lens models with density cusps

Lens models appropriate for representing cusped galaxies and clusters are developed. The analogue of the odd-number theorem for cusped density distributions is given. Density cusps are classified into strong, isothermal or weak, according to their lensing properties. Strong cusps cause multiple imaging for any source position, whereas isothermal and weak cusps give rise to only one image for distant sources. Isothermal cusps always possess a pseudo-caustic. When the source crosses the pseudo-caustic, the number of images changes by unity.   Two families of cusped galaxy and cluster models are examined in detail. The double power-law family has an inner cusp, followed by a transition region and an outer envelope. One member of this family — the isothermal double power-law model — possesses an exceedingly scarce property, namely the lens equation is exactly solvable for any source position. This means that the magnifications, the time delay and the lensing cross-sections are all readily available. The model has a three-dimensional density that is cusped like r ⁻² at small radii and falls off like r ⁻⁴ asymptotically. Thus, it provides a very useful representation of the lensing properties of a galaxy or cluster of finite total mass with a flat rotation curve. The second set of models studied is the single power-law family. These are single density cusps of infinite extent. The properties of the critical curves and caustics and the behaviour of the lenses in the presence of external shear are all discussed in some detail.     

The polarimetric nature of HD 108

V -band polarimetric observations of HD 108 were obtained during 1994 August and September. This is the first time that the temporal polarimetric variability of this star has been investigated. Its percentage polarization and position angle vary randomly, and no evidence of its proposed binary nature is detected. The random variability is consistent with that seen for previously observed Wolf–Rayet stars and OB supergiants. Nightly variations do show some systematic behaviour consistent with the blob ejection model of Underhill & Fehey (1984). From the data it is estimated that the mass-loss rate resulting from blobs is ∼1-10⁻⁷ M⊙ yr⁻¹. The total mass-loss rate of the star is estimated to be ∼1-10⁻⁵ M⊙ yr⁻¹. The stellar rotation rate is estimated to be V _rot∼400 km s⁻¹ with an inclination of i <20c. The data are also interpreted in terms of perturbations occurring in an equatorial disc plus bipolar jets viewed equatorially (Underhill 1994). It is found that the perturbation mass-loss rate is about ∼3-10⁻⁷ M⊙ yr⁻¹, giving a total mass-loss rate for the star of 13-10⁻⁵ M⊙ yr⁻¹. Both models are consistent with HD 108 being an OB supergiant or of the class Ofpe/WR9.     

Three-integral models for axisymmetric galactic discs

We present new equilibrium component distribution functions that depend on three analytic integrals in a Stäckel potential, and that can be used to model stellar discs of galaxies. These components are generalizations of two-integral ones and can thus provide thin discs in the two-integral approximation. Their most important properties are the partly analytical expression for their moments, the disc-like features of their configuration space densities (exponential decline in the galactic plane and finite extent in the vertical direction) and the anisotropy of their velocity dispersions. We further show that a linear combination of such components can fit a van der Kruit disc.     

New biorthogonal potential–density basis functions

We use the weighted integral form of spherical Bessel functions and introduce a new analytical set of complete and biorthogonal potential–density basis functions. The potential and density functions of the new set have finite central values and they fall off, respectively, similar to   r ^{−(1+ l )}  and   r ^{−(4+ l )}  at large radii, where l is the latitudinal quantum number of spherical harmonics. The lowest order term associated with   l = 0  is the perfect sphere of de Zeeuw. Our basis functions are intrinsically suitable for the modelling of three-dimensional, soft-centred stellar systems and they complement the basis sets of Clutton-Brock, Hernquist & Ostriker and Zhao. We test the performance of our functions by expanding the density and potential profiles of some spherical and oblate galaxy models.     

The redshift evolution of clustering in the Hubble Deep Field

We present a correlation function analysis for the catalogue of photometric redshifts obtained from the Hubble Deep Field image by Fernandez-Soto, Lanzetta & Yahil. By dividing the catalogue into redshift bins of width Δ z =0.4 we measured the angular correlation function w ( θ ) as a function of redshift up to z ∼4.8. From these measurements we derive the trend of the correlation length r ₀. We find that r ₀( z ) is roughly constant with look-back time up to z ≃2, and then increases to higher values at z ≳2.4. We estimate the values of r ₀, assuming ξ ( r , z )=[ r r ₀( z )]^{− γ} , γ =1.8 and various geometries. For Ω₀=1 we find r ₀( z =3)≃7.00±4.87  h ⁻¹ Mpc, in good agreement with the values obtained from analysis of the Lyman break galaxies.     

Alignment between luminous red galaxies and surrounding structures at z∼ 0.5

We analyse a high-redshift sample (0.4 < z < 0.5) of luminous red galaxies (LRGs) extracted from the Sloan Digital Sky Survey data release 4 and their surrounding structures to explore the presence of alignment effects of these bright galaxies with neighbour objects. In order to avoid projection effects, we compute photometric redshifts for galaxies within 3  h ⁻¹ Mpc in projection of LRGs and calculate the relative angle between the LRG major axis and the direction to neighbours within 1000 km s⁻¹. We find a clear signal of alignment between LRG orientations and the distribution of galaxies within 1.5  h ⁻¹ Mpc. The alignment effects are present only for the red population of tracers; LRG orientation is uncorrelated to the blue population of neighbour galaxies. These results add evidence to the alignment effects between primaries and satellites detected at low redshifts. We conclude that such alignments were already present at z ∼ 0.5.