Self-similar spherical collapse with non-radial motions

We derive the asymptotic mass profile near the collapse centre of an initial spherical density perturbation, δ ∝ M ^{− ε} , of collisionless particles with non-radial motions. We show that angular momenta introduced at the initial time do not affect the mass profile. Alternatively, we consider a scheme in which a particle moves on a radial orbit until it reaches its turnaround radius, r ∗. At turnaround the particle acquires an angular momentum L =ℒ√ GM _* r _* per unit mass, where M ∗ is the mass interior to r ∗. In this scheme, the mass profile is M ∝ r ^{3/(1+3 ε )} for all ε >0 , in the region r / r _t ≪ℒ , where r _t is the current turnaround radius. If ℒ≪1 then the profile in the region ℒ≪ r / r _t ≪1 is M ∝ r for ε <2/3 , and remains M ∝ r ^{3/(1+3 ε )} for ε ≥2/3 . The derivation relies on a general property of non-radial orbits which is that the ratio of the pericentre to apocentre is constant in a force field k ( t ) r ⁿ with k ( t ) varying adiabatically.     

Star formation triggered by supernova explosions in young galaxies

We study the evolution of supernova remnants in a low-metallicity medium   Z /Z_⊙= 10⁻⁴ to 10⁻²  in the early universe, using one-dimensional hydrodynamics with non-equilibrium chemistry. Once a post-shock layer is able to cool radiatively, a dense shell forms behind the shock. If this shell becomes gravitationally unstable and fragments into pieces, next-generation stars are expected to form from these fragments. To explore the possibility of this triggered star formation, we apply a linear perturbation analysis of an expanding shell to our results and constrain the parameter range of ambient density, explosion energy and metallicity where fragmentation of the shell occurs. For the explosion energy of  10⁵¹ erg (10⁵² erg)  , the shell fragmentation occurs for ambient densities higher than  ≳10² cm⁻³ (10 cm⁻³  ), respectively. This condition depends little on the metallicity in the ranges we examined. We find that the mode of star formation triggered occurs only in massive  (≳10⁸ M_⊙)  haloes.     

New limits on the possible variation of physical constants

Recent detections of high-redshift absorption by both atomic hydrogen and molecular gas in the radio spectra of quasars have provided a powerful tool for measuring possible temporal and spatial variations of physical 'constants' in the universe.
We compare the frequency of high-redshift hydrogen 21-cm absorption with that of associated molecular absorption in two quasars to place new (1σ) upper limits on any variation in y≡g_pα² (where α is the fine-structure constant, and g_p is the proton g -factor of
   
at redshifts z = 0.25 and 0.68. These quasars are separated by a comoving distance of 3000 Mpc ( H ₀= 75 km s⁻¹ Mpc⁻¹ and q ₀). We also derive limits on the time rates of change
   
    between the present epoch and z = 0.68. These limits are more than an order of magnitude smaller than previous results derived from high-redshift measurements.     

Microwave polarization in the direction of galaxy clusters induced by the CMB quadrupole anisotropy

Electron scattering induces a polarization in the cosmic microwave background (CMB) signal measured in the direction of a galaxy cluster owing to the presence of a quadrupole component in the CMB temperature distribution. Measuring the polarization towards distant clusters provides the unique opportunity to observe the evolution of the CMB quadrupole at moderate redshifts, z ∼0.5–3. We demonstrate that for the local cluster population the polarization degree will depend on the cluster celestial position. There are two extended regions in the sky, which are opposite to each other, where the polarization is maximal, ∼0.1( τ /0.02) μK in the Rayleigh–Jeans part of the CMB spectrum ( τ being the Thomson optical depth across the cluster). This value exceeds the polarization introduced by the cluster transverse peculiar motion if v _t<1300 km s⁻¹. One can hope to detect this small signal by measuring a large number of clusters, thereby effectively removing the systematic contribution from other polarization components produced in clusters. These polarization effects, which are of the order of ( v _t c )² τ , ( v _t c ) τ ² and ( kT _e m _e c ²) τ ², as well as the polarization owing to the CMB quadrupole, were previously given by Sunyaev and Zel'dovich for the Rayleigh–Jeans part of the spectrum. We fully confirm their earlier results and present exact frequency dependences for all these effects. The polarization degree is considerably higher in the Wien region.     

Global textures and the formation of self-bound gravitational systems

Using a Newtonian approximation we developed a quantitative criterion for the collapse of a spherical distribution of matter under an isolated texture field. In particular, we found the evolution of an overdense region is strongly determined by two parameters: the energy scale of symmetry breaking ( η ) and the initial radius of the system. Applying our collapse criterion to typical galaxy scales we verified the formation of 10¹¹‐M_⊙ objects at z ≲9 and 10¹²‐M_⊙ objects at z ≲5.     

Microlensing by cosmic strings

We consider the signature and detectability of gravitational microlensing of distant quasars by cosmic strings. Because of the simple image configuration such events will have a characteristic lightcurve, in which a source would appear to brighten by exactly a factor of 2, before reverting to its original apparent brightness. We calculate the optical depth and event rate, and conclude that current predictions and limits on the total length of strings on the sky imply optical depths of ≲ 10⁻⁸ and event rates of fewer than one event per 10⁹ sources per year. Disregarding those predictions but replacing them with limits on the density of cosmic strings from the cosmic microwave background fluctuation spectrum, leaves only a small region of parameter space (in which the sky contains about 3 × 10⁵ strings with deficit angle of the order of 0.3 milli-seconds) for which a microlensing survey of exposure 10⁷ source years, spanning a 20–40-year period, might reveal the presence of cosmic strings.     

Prompt optical emission and synchrotron self-absorption constraints on emission site of GRBs

We constrain the distance of the gamma-ray burst (GRB) prompt emission site from the explosion centre R , by determining the location of the electron's self-absorption frequency in the GRB prompt optical-to-X/γ-ray spectral energy distribution, assuming that the optical and the γ-ray emissions are among the same synchrotron radiation continuum of a group of hot electrons. All possible spectral regimes are considered in our analysis. The method has only two assumed parameters, namely the bulk Lorentz factor of the emitting source Γ and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of four bursts that satisfy the following three criteria: (1) they all have simultaneous optical and γ-ray detections in multiple observational time intervals, (2) they all show temporal correlations between the optical and γ-ray light curves and (3) the optical emission is consistent with belonging to the same spectral component as the γ-ray emission. For all the time intervals of these four bursts, it is inferred that   R ≥ 10¹⁴  (Γ/300)^3/4 ( B /10⁵ G)^1/4  cm. For a small fraction of the sample, the constraint can be pinned down to   R ≈ 10¹⁴–10¹⁵ cm  for  Γ∼ 300  . For a second sample of bursts with prompt optical non-detections, only upper limits on R can be obtained. We find no inconsistency between the R -constraints for this non-detection sample and those for the detection sample.     

Can a slowly rotating neutron star be a radio pulsar?

It is shown that the radius of curvature of magnetic field lines in the polar region of a rotating magnetized neutron star can be significantly less than the usual radius of curvature of the dipole magnetic field. The magnetic field in the polar cap is distorted by toroidal electric currents flowing in the neutron star crust. These currents close up the magnetospheric currents driven by the electron–positron plasma generation process in the pulsar magnetosphere. Owing to the decrease in the radius of curvature, electron–positron plasma generation becomes possible even for slowly rotating neutron stars, with   PB ^−2/3₁₂ < 10 s  , where P is the period of star rotation and   B ₁₂= B /10¹² G  is the magnitude of the magnetic field on the star surface.     

A circular statistical method for extracting rotation measures

We propose a new method for the extraction of rotation measures from spectral polarization data. The method is based on maximum likelihood analysis and takes into account the circular nature of the polarization data. The method is unbiased and statistically more efficient than the standard χ ² procedure. We also find that the method is computationally much more convenient than the standard χ ² procedure if the number of data points is very large. We find that for most sources the method gives results very close to the standard χ ² minimization procedure. We give results for all the cases for which the method gives significantly different results from χ ² minimization. We also make a λ ³ fit to the data in order to extract non-Faraday rotation behaviour for those sources for which a large number of data points are available.     

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.     

Dark matter subhaloes in numerical simulations

We use cosmological Λ cold dark matter (CDM) numerical simulations to model the evolution of the substructure population in 16 dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10¹¹ to  10¹⁴ M_⊙  at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation  d n /d v = (1/8)( v _cmax/ v _cmax,host)⁻⁴  . Halo to halo variance in the VDF is a factor of roughly 2 to 4. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for   r ≳ 0.3 r _vir  , where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease towards the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias towards more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power-law cosmologies of   P ∝ kⁿ   , and demonstrate that a steeper spectral index, n , results in significantly less substructure.     

Principal-component analysis of the cosmic microwave background anisotropies: revealing the tensor degeneracy

A principal-component analysis of cosmic microwave background (CMB) anisotropy measurements is used to investigate degeneracies among cosmological parameters. The results show that a degeneracy with tensor modes – the 'tensor degeneracy'– dominates uncertainties in estimates of the baryon and cold dark matter densities,   ω _b=Ω_b  h ²  ,   ω _c=Ω_c  h ²  , ¹ from an analysis of CMB anisotropies alone. The principal-component analysis agrees well with a maximum-likelihood analysis of the observations, identifying the main degeneracy directions and providing an impression of the effective dimensionality of the parameter space.     

Multi-Conjugate Adaptive Optics with laser guide stars: performance in the infrared and visible

The assumption of the Gaussianity of primordial perturbations plays an important role in modern cosmology. The most direct test of this hypothesis consists of testing the Gaussianity of cosmic microwave background (CMB) maps. Counting the pixels with the temperatures in given ranges and thus estimating the one-point probability function of the field is the simplest of all the tests. Other usually more complex tests of Gaussianity generally use a great deal of the information already contained in the probability function. However, the most interesting outcome of such a test would be the signal of non-Gaussianity independent of the probability function. It is shown that the independent information has purely morphological character i.e. it depends on the geometry and topology of the level contours only. As an example we discuss in detail the quadratic model   v = u + α ( u ²-1)  ( u is a Gaussian field with   u¯ =0  and  〈 u ²〉=1  , α is a parameter) that may arise in slow-roll or two-field inflation models. We show that in the limit of small amplitude α the full information about the non-Gaussianity is contained in the probability function. If other tests are performed on this model they simply recycle the same information. A simple procedure allowing us to assess the sensitivity of any statistics to the morphological information is suggested. We provide an analytic estimate of the statistical limit for detecting the quadratic non-Gaussianity α _c as a function of the map size in the ideal situation when the scale of the field is resolved. This estimate is in a good agreement with the results of the Monte Carlo simulations of 256² and 1024² maps. The effect of resolution on the detection quadratic non-Gaussianity is also briefly discussed.     

IRAM observations of JVAS/CLASS gravitational lenses

We have searched for molecular absorption lines at millimetre wavelengths in 11 gravitational lens systems discovered in the JVAS/CLASS surveys of flat spectrum radio sources. Spectra of only one source 1030+074 were obtained in the 3-, 2- and 1.3-mm bands at the frequencies corresponding to common molecular transitions of CO and HCO⁺ as continuum emission was not found in any of the other sources. We calculated upper limits to the column density in molecular absorption for 1030+074, using an excitation temperature of 15 K, to be N _CO<6.3×10¹³ cm⁻² and N _HCO+<1.3×10¹¹ cm⁻² , equivalent to hydrogen column density of the order N _H<10¹⁸ cm⁻² , assuming standard molecular abundances. We also present the best upper limits of the continuum at the lower frequency for the other 10 gravitational lenses.     

Molecular hydrogen emission in Cygnus A

We present J , H and K -band spectroscopy of Cygnus A, spanning 1.0–2.4 μm in the rest-frame and hence several rovibrational H₂, H recombination and [Fe  ii ] emission lines. The lines are spatially extended by up to 6 kpc from the nucleus, but their distinct kinematics indicate that the three groups (H, H₂ and [Fe  ii ]) are not wholly produced in the same gas. The broadest line, [Fe  ii ] λ 1.644, exhibits a non-Gaussian profile with a broad base (FWHM≃1040 km s⁻¹), perhaps because of the interaction with the radio source. Extinctions to the line-emitting regions substantially exceed earlier measurements based on optical H recombination lines.
Hard X-rays from the quasar nucleus are likely to dominate the excitation of the H₂ emission. The results of Maloney, Hollenbach & Tielens are thus used to infer the total mass of gas in H₂ v=1–0 S(1)-emitting clouds as a function of radius, for gas densities of 10³ and 10⁵ cm⁻³, and stopping column densities N _H=10²²–10²⁴ cm⁻². Assuming azimuthal symmetry, at least 2.3×10⁸ M_⊙ of such material is present within 5 kpc of the nucleus, if the line-emitting clouds see an unobscured quasar spectrum. Alternatively, if the bulk of the X-ray absorption to the nucleus inferred by Ueno et al. actually arises in a circumnuclear torus, the implied gas mass rises to ∼10¹⁰ M_⊙. The latter plausibly accounts for 10⁹ yr of mass deposition from the cluster cooling flow, for which within this radius.     

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.     

The effective acceleration of plasma outflow in the paraboloidal magnetic field

The problem of the efficiency of particle acceleration for a paraboloidal poloidal magnetic field is considered within the approach of steady axisymmetric magnetohydrodynamic (MHD) flow. For the large Michel magnetization parameter σ it is possible to linearize the stream equation near the force-free solution and to solve the problem self-consistently as was done by Beskin, Kuznetsova & Rafikov for a monopole magnetic field. It is shown that, on the fast magnetosonic surface (FMS), the particle Lorentz factor γ does not exceed the standard value  σ^1/3  . On the other hand, in the supersonic region, the Lorentz factor grows with the distance z from the equatorial plane as  γ≈ ( z / R _L)^1/2  up to the distance   z ≈σ² R _L  , where   R _L= c /Ω_F  is the radius of the light cylinder. Thus, the maximal Lorentz factor is  γ_max≈σ  , which corresponds to almost the full conversion of the Poynting energy flux into the particle kinetic one.     

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.     

Primordial non-Gaussianity: local curvature method and statistical significance of constraints on fNL from WMAP data

We test the consistency of estimates of the non-linear coupling constant f _NL using non-Gaussian cosmic microwave background (CMB) maps generated by the method described in the work of Liguori, Matarrese & Moscardini. This procedure to obtain non-Gaussian maps differs significantly from the method used in previous works on the estimation of f _NL. Nevertheless, using spherical wavelets, we find results in very good agreement with Mukherjee & Wang, showing that the two ways of generating primordial non-Gaussian maps give equivalent results. Moreover, we introduce a new method for estimating the non-linear coupling constant from CMB observations by using the local curvature of the temperature fluctuation field. We present both Bayesian credible regions (assuming a flat prior) and proper (frequentist) confidence intervals on f _NL, and discuss the relation between the two approaches. The Bayesian approach tends to yield lower error bars than the frequentist approach, suggesting that a careful analysis of the different interpretations is needed. Using this method, we estimate   f _NL=−10⁺²⁷⁰₋₂₆₀  at the 2σ level (Bayesian) and   f _NL=−10⁺³¹⁰₋₂₇₀  (frequentist). Moreover, we find that the wavelet and the local curvature approaches, which provide similar error bars, yield approximately uncorrelated estimates of f _NL and therefore, as advocated in the work of Cabella et al., the estimates may be combined to reduce the error bars. In this way, we obtain   f _NL=−5 ± 85  and   f _NL=−5 ± 175  at the 1σ and 2σ level respectively using the frequentist approach.     

Detection of an H92α recombination line in the starburst galaxy NGC 660

We have used the Very Large Array (VLA) to search for the H92α radio recombination line (RRL) in four starburst galaxies. In NGC 660, the line was detected over a 17Å‐8 arcsec² region near its starburst nucleus. The line and continuum emission indicate that the RRL-emitting gas is most likely in the form of a cluster of H ii regions with a small filling factor. Using a simple model we find that the total ionized mass in the nuclear region is in the range 2–8Å‐10⁴ M⊙ and the rate of production of UV photons N _Lyc∼1–3Å‐10⁵³ s⁻¹. The ratio of H92α and Brγ line intensities in NGC 660 indicates that extinction is significant even at λ=2 μm. The velocity field of the ionized gas is consistent with a rotating disc with an average velocity gradient of ∼15 km s⁻¹ arcsec⁻¹. The dynamical mass within the central 500 pc is ∼4Å‐10⁸ M⊙ and may be about ∼6Å‐10⁷ M⊙ within the central 120 pc. No line was detected in the other galaxies (NGC 520, NGC 1614 and NGC 6946) to a 3σ limit of 300 μJy. In the starburst galaxies in which RRLs have been detected, we find that there is a rough correlation between the integrated H92α line flux density and both the total far-infrared flux density and the radio continuum emission from the central region.