Survival and disruption of subsystems during a cold collapse

Cold collapse of a cluster composed of small identical clumps, each of which is in virial equilibrium, is considered. Since the clumps have no relative motion with respect to each other initially, the cluster collapses under its own gravity. At the first collapse of the cluster, most of the clumps are destroyed, but some survive. In order to find the condition for the clumps to survive, we made a systematic study in two-parameter space: the number of the clumps N _c and the size of the clump r _v . We obtained the condition N _c ≫ 1 and n _k  ≥ 1, where n _k is related to r _v and the initial radius of the cluster R _ini through the relation R _ini/ r _v  = 2 N ^{( n k +5)/6}_c. A simple analytical argument supports the numerical result. This n _k corresponds to the index of the power spectrum of the density fluctuation in the cosmological hierarchical clustering, and thus our result may suggest that in the systems smaller than 2/Ω h ²)Mpc, the first violent collapse is strong enough to sweep away all the substructures that exist before the collapse.     

A silhouette envelope around GGD30IR detected by Spitzer

This paper presents a study of the envelope of the young stellar object (YSO) GGD30IR. What distinguishes this from most other YSOs is the elongated absorption feature seen in silhouette against the background emission in the Spitzer Galactic Legacy Infrared Midplane Survey Extraordinaire 8 μm Infrared Array Camera image of the region. The size and the symmetrical placement of GGD30IR in the centre of this feature suggest that it is an extended envelope, perhaps the remnant of the collapse of the GGD30 core. We have used the extinction in the envelope measured from (i) the reduction in the 8 μm background intensity and (ii) field star colour excesses, to estimate the envelope mass, obtaining values of 0.6 ± 0.2 and  0.5 ± 0.3 M_⊙  , respectively. To investigate the envelope further, we have obtained Australia Telescope Compact Array 3 mm continuum and HCO+ line observations of the region. The continuum emission at 3 mm arises from both a compact (unresolved; ≤730 au) core embedded in an extended envelope ∼18 000 au × 38 000 au in extent. We estimate the core mass to be 0.11 ± 0.02  M_⊙  . The HCO+ emission is extended in a direction perpendicular to the long axis of the envelope, suggesting it comes from an outflow. The spectral energy distribution (SED) provides a 2–24 μm spectral index, α= 1.0, which places GGD30IR in the Class I YSO category. Integrating the SED provides a luminosity of   L _*≃ 25 ± 5 L_⊙  .     

Spectral analysis of the high-gravity extreme helium star LS IV+6°2

The optical spectrum of the early B hydrogen-deficient star LS IV+6°2 has been analysed. With T _eff = 31000 K, log  g  = 4.05, n _H ∼ 10⁻⁴ and n _c ∼ 0.1 (relative abundances by number), it is the hottest high-gravity extreme helium star (EHe) yet studied. The He  i spectrum shows all predicted permitted and forbidden transitions in absorption. LS IV+6°2 is a comparatively metal-rich EHe star; abundances of C, N, O, Ne, Mg, Al and P are typical of other EHes, whilst Si and S are somewhat deficient. With the surface parameters given, LS IV+6°2 lies close to the boundary of the helium star pulsation instability finger near T _eff ∼ 27000 K. Available data indicate that the radial velocity is variable, but give no indication of amplitude or period.     

Rotation and activity in the solar-type stars of NGC 2547

We present high-resolution spectroscopy of a sample of 24 solar-type stars in the young (15–40 Myr), open cluster NGC 2547. We use our spectra to confirm cluster membership in 23 of these stars, to determine projected equatorial velocities and chromospheric activity, and to search for the presence of accretion discs. We find examples of both fast ( v _e sin  i >50 km s⁻¹) and slow ( v _e sin  i <10 km s⁻¹) rotators, but no evidence for active accretion in any of the sample. The distribution of projected rotation velocities is indistinguishable from the slightly older IC 2391 and IC 2602 clusters, implying similar initial angular momentum distributions and circumstellar disc lifetimes. The presence of very slow rotators indicates either that long (10–40 Myr) disc lifetimes or internal differential rotation are needed, or that NGC 2547 (and IC 2391/2602) were born with more slowly rotating stars than are presently seen in even younger clusters and associations. The solar-type stars in NGC 2547 follow a similar rotation–activity relationship to that seen in older clusters. X-ray activity increases until a saturation level is reached for v _e sin  i >15–20 km s⁻¹. We are unable to explain why this saturation level, of log( L _x L _bol)≃−3.3, is a factor of 2 lower than in other clusters, but rule out anomalously slow rotation rates or uncertainties in X-ray flux calculations.     

An explosive end to intermediate-mass zero-metallicity stars and early Universe nucleosynthesis

We use the Cambridge stellar evolution code stars to model the evolution of 5 and  7 M_⊙  zero-metallicity stars. With enhanced resolution at the hydrogen- and helium-burning shell in the asymptotic giant branch (AGB) phases, we are able to model the entire thermally pulsing AGB (TP-AGB) phase. The helium luminosities of the thermal pulses are significantly lower than in higher metallicity stars so there is no third dredge-up. The envelope is enriched in nitrogen by hot-bottom burning of carbon that was previously mixed in during second dredge-up. There is no s -process enrichment owing to the lack of third dredge-up. The thermal pulses grow weaker as the core mass increases and they eventually cease. From then on the star enters a quiescent burning phase which lasts until carbon ignites at the centre of the star when the CO core mass is  1.36 M_⊙  . With such a high degeneracy and a core mass so close to the Chandrasekhar mass, we expect these stars to explode as type 1.5 supernovae, very similar to type Ia supernovae but inside a hydrogen-rich envelope.     

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.     

Structural parameters of Mayall II = G1 in M31

Mayall II = G1 is one of the most luminous globular clusters (GCs) known in M31. New deep, high-resolution observations with the Advanced Camera for Surveys on the Hubble Space Telescope are used to provide accurate photometric data to the smallest radii yet. In particular, we present the precise variation of ellipticity and position angle, and of surface brightness for the core of the object. Based on these accurate photometric data, we redetermine the structural parameters of G1 by fitting a single-mass isotropic King model. We derive a core radius,   r _c= 0.21 ± 0.01  arcsec (= 0.78 ± 0.04  pc)  , a tidal radius,   r _t= 21.8 ± 1.1  arcsec (= 80.7 ± 3.9  pc)  , and a concentration index   c = log ( r _t/ r _c) = 2.01 ± 0.02  . The central surface brightness is 13.510 mag arcsec⁻². We also calculate the half-light radius, at   r _h= 1.73 ± 0.07  arcsec (= 6.5 ± 0.3  pc)  . The results show that, within 10 core radii, a King model fits the surface brightness distribution well. We find that this object falls in the same region of the   M_V   versus  log   R _h  diagram as ω Centauri, M54 and NGC 2419 in the Milky Way. All three of these objects have been claimed to be the stripped cores of now defunct dwarf galaxies. We discuss in detail whether GCs, stripped cores of dwarf spheroidals and normal dwarf galaxies form a continuous distribution in the   M_V   versus  log   R _h  plane, or if GCs and dwarf spheroidals constitute distinct classes of objects; we present arguments in favour of this latter view.     

Streaming motions of galaxy clusters within 12 000 km s−1– II. New photometric data for the Fundamental Plane

We present new R -band photometric data for 447 galaxies, gathered for the 'Streaming Motions of Abell Clusters' (SMAC) project. The data comprise 629 individual measurements of the Fundamental Plane (FP) parameters effective radius ( R _e) and surface brightness (〈 μ 〉_e), derived from r ^1/4-law profile fitting. More than a third of the galaxies were observed more than once. The photometric precision is ∼0.02 mag as judged from comparisons of aperture photometry between repeat observations of galaxies. The combination     which enters into the Fundamental Plane relation, has internal uncertainties of ∼0.008, corresponding to < 2 per cent in estimated distance. Taken individually, the (correlated) internal errors in R _e and 〈 μ 〉_e are ∼8 per cent and ∼0.12 mag respectively. Comparisons with literature data constrain the external random errors to ≲5 per cent in distance (per observation), which is small in comparison to the ∼20 per cent scatter in the FP. The data will form part of a merged catalogue of FP parameters, presented in a companion paper.     

A constant dark matter halo surface density in galaxies

We confirm and extend the recent finding that the central surface density  μ_0D≡ r ₀ρ₀  of galaxy dark matter haloes, where r ₀ and  ρ₀  are the halo core radius and central density, is nearly constant and independent of galaxy luminosity. Based on the co-added rotation curves (RCs) of ∼1000 spiral galaxies, the mass models of individual dwarf irregular and spiral galaxies of late and early types with high-quality RCs, and the galaxy–galaxy weak-lensing signals from a sample of spiral and elliptical galaxies, we find that  log μ_0D= 2.15 ± 0.2  in units of  log(M_⊙ pc⁻²)  . We also show that the observed kinematics of Local Group dwarf spheroidal galaxies are consistent with this value. Our results are obtained for galactic systems spanning over 14 mag, belonging to different Hubble types and whose mass profiles have been determined by several independent methods. In the same objects, the approximate constancy of  μ_0D  is in sharp contrast to the systematical variations, by several orders of magnitude, of galaxy properties, including  ρ₀  and central stellar surface density.     

Helium recombination spectra as temperature diagnostics for planetary nebulae

Electron temperatures derived from the He  i recombination line ratios, designated T _e(He  i ), are presented for 48 planetary nebulae (PNe). We study the effect that temperature fluctuations inside nebulae have on the T _e(He  i ) value. We show that a comparison between T _e(He  i ) and the electron temperature derived from the Balmer jump of the H  i recombination spectrum, designated T _e(H  i ), provides an opportunity to discriminate between the paradigms of a chemically homogeneous plasma with temperature and density variations, and a two-abundance nebular model with hydrogen-deficient material embedded in diffuse gas of a 'normal' chemical composition (i.e. ∼solar), as the possible causes of the dichotomy between the abundances that are deduced from collisionally excited lines and those deduced from recombination lines. We find that T _e(He  i ) values are significantly lower than T _e(H  i ) values, with an average difference of  〈 T _e(H  i ) − T _e(He  i )〉= 4000 K  . The result is consistent with the expectation of the two-abundance nebular model but is opposite to the prediction of the scenarios of temperature fluctuations and/or density inhomogeneities. From the observed difference between T _e(He  i ) and T _e(H  i ), we estimate that the filling factor of hydrogen-deficient components has a typical value of 10⁻⁴. In spite of its small mass, the existence of hydrogen-deficient inclusions may potentially have a profound effect in enhancing the intensities of He  i recombination lines and thereby lead to apparently overestimated helium abundances for PNe.     

Electron density diagnostic potential of the Si x ion and its application in Procyon

Theoretical electron density sensitive line ratios   R ₁– R ₆  of Si  x soft X-ray emission lines are presented. We found that these line ratios are sensitive to electron density n _e, and the ratio R ₁ is insensitive to electron temperature T _e. For reliable determination of the electron density of laboratory and astrophysical plasmas, atomic data, such as electron impact excitation rates, are very important. Our results reveal that the discrepancy of the line ratios from different atomic data calculated with the distorted wave (DW) approximation and the R-matrix method is up to 19 per cent at   n _e= 2 × 10⁸ cm⁻³  . We applied the theoretical intensity ratio R ₁ to the Low Energy Transmission Grating Spectrometer (LETGS) spectrum of the solar-like star Procyon. By comparing the observed value (1.29) with the theoretical calculation, the derived electron density n _e is  2.6 × 10⁸ cm⁻³  , which is consistent with that derived from  (C  v < 8.3 × 10⁸ cm⁻³)  . When the temperature structure of the Procyon corona is taken into account, the derived electron density increases from   n _e= 2.6 × 10⁸  to  2.8 × 10⁸ cm⁻³  .     

Measuring Ω0 using cluster evolution

The evolution of the abundance of galaxy clusters depends sensitively on the value of the cosmological density parameter, Ω₀. Recent ASCA data are used to quantify this evolution as measured by the cluster X-ray temperature function. A χ² minimization fit to the cumulative temperature function, as well as a maximum-likelihood estimate (which requires additional assumptions about cluster luminosities), leads to the estimate Ω₀ ≈ 0.45 ± 0.25 (1σ statistical error). Various systematic uncertainties are considered, none of which significantly enhances the probability that Ω₀ = 1. These conclusions hold for models with or without a cosmological constant, i.e., with Λ₀ = 0 or Λ₀ = 1 − Ω₀. The statistical uncertainties are at least as large as any of the individual systematic errors that have been considered here, suggesting that additional temperature measurements of distant clusters will allow an improvement in this estimate. An alternative method that uses the highest redshift clusters to place an upper limit on Ω₀ is also presented and tentatively applied, with the result that Ω₀  1 can be ruled out at the 98 per cent confidence level. Whilst this method does not require a well-defined statistical sample of distant clusters, there are still modelling uncertainties that preclude a firmer conclusion at this time.     

Old high-redshift galaxies and primordial density fluctuation spectra

We have discovered a population of extremely red galaxies at z  ≃ 1.5 which have apparent stellar ages of ≳ 3 Gyr, based on detailed spectroscopy in the rest-frame ultraviolet. In order for galaxies to have existed at the high collapse redshifts indicated by these ages, there must be a minimum level of power in the density fluctuation spectrum on galaxy scales. This paper compares the required power with that inferred from other high-redshift populations: damped Lyα absorbers and Lyman-limit galaxies at z  ≃ 3.2. If the collapse redshifts for the old red galaxies are in the range z _c ≃ 6–8, there is general agreement between the various tracers on the required inhomogeneity on 1-Mpc scales. This level of small-scale power requires the Lyman-limit galaxies to be approximately ν ≃ 3.0 fluctuations, implying a very large bias parameter b  ≃ 6. If the collapse redshifts of the red galaxies are indeed in the range z _c = 6–8 required for power spectrum consistency, their implied ages at z  ≃ 1.5 are between 3 and 3.8 Gyr for essentially any model universe of current age 14 Gyr. The age of these objects as deduced from gravitational collapse thus provides independent support for the ages estimated from their stellar populations. Such early-forming galaxies are rare, and their contribution to the cosmological stellar density is consistent with an extrapolation to higher redshifts of the star formation rate measured at z  < 5; there is no evidence for a general era of spheroid formation at extreme redshifts.     

The Gunn–Peterson effect and the Lyman alpha forest

We show that spatial correlations in a stochastic large-scale velocity field in an otherwise smooth intergalactic medium (homogeneous comoving density) superposed on the general Hubble flow may cause a 'line-like' structure in QSO spectra similar to the population of unsaturated Lyα forest lines which usually are attributed to individual clouds with 10¹¹ ≲ N _{H i}  5 × 10¹³ cm⁻². Therefore there is no clear observational distinction between a diffuse intergalactic medium and discrete intergalactic clouds. It follows that the H  i density in the diffuse intergalactic medium might be substantially underestimated if it is determined from the observed intensity distribution near the apparent continuum in high-resolution spectra of QSOs. Our tentative estimate implies a diffuse neutral hydrogen opacity τ_GP ∼ 0.3 at z  ∼ 3 and a current baryon density Ω_IGM ≃ 0.08, assuming a Hubble constant H ₀ = 70 km s⁻¹ Mpc⁻¹.     

Gamma-ray burst afterglows: effects of radiative corrections and non-uniformity of the surrounding medium

The afterglow of a gamma-ray burst (GRB) is commonly thought to be the result of continuous deceleration of a relativistically expanding fireball in the surrounding medium. Assuming that the expansion of the fireball is adiabatic and that the density of the medium is a power-law function of shock radius, i.e. n _ext ∝  R ^{− k} , we study the effects of the first-order radiative correction and the non-uniformity of the medium on a GRB afterglow analytically. We first derive a new relation among the observed time, the shock radius and the Lorentz factor of the fireball: t _⊕ =  R /4(4− k ) γ²c, and also derive a new relation among the comoving time, the shock radius and the Lorentz factor of the fireball: t _co = 2 R /(5− k ) γc. We next study the evolution of the fireball by using the analytic solution of Blandford &38; McKee. The radiation losses may not significantly influence this evolution. We further derive new scaling laws both between the X-ray flux and observed time and between the optical flux and observed time. We use these scaling laws to discuss the afterglows of GRB 970228 and GRB 970616, and find that if the spectral index of the electron distribution is p  = 2.5, implied from the spectra of GRBs, the X-ray afterglow of GRB 970616 is well fitted by assuming k  = 2.     

Detection of OH and wide H i absorption toward B0218+357

Large-scale polarization of the cosmic microwave background measured by the WMAP satellite requires a mean optical depth to Thomson scattering,  τ_e∼ 0.17  . The reionization of the Universe must therefore have begun at relatively high redshift. We have studied the reionization process using supercomputer simulations of a large and representative region of a universe which has cosmological parameters consistent with the WMAP results (  Ω_m= 0.3, Ω_Λ= 0.7, h = 0.7, Ω_b= 0.04, n = 1  and  σ₈= 0.9  ). Our simulations follow both the radiative transfer of ionizing photons and the formation and evolution of the galaxy population which produces them. A previously published model with ionizing photon production as expected for zero-metallicity stars distributed according to a standard stellar initial mass function (IMF) (10⁶¹ photons per unit solar mass of formed stars) and with a moderate photon escape fraction from galaxies (5 per cent), produces  τ_e= 0.104  , which is within 1.0 to  1.5σ  of the 'best' WMAP value. Values of up to 0.16 can be produced by taking larger escape fractions or a top-heavy IMF. The data do not require a separate populations of 'miniquasars' or of stars forming in objects with total masses below  10⁹ M_⊙  . Reconciling such early reionization with the observed Gunn–Peterson troughs in   z > 6  quasars may be challenging. Possible resolutions of this problem are discussed.     

Photoionized Hβ emission in NGC 5548: it breathes!

Emission-line regions in active galactic nuclei (AGNs) and other photoionized nebulae should become larger in size when the ionizing luminosity increases. This 'breathing' effect is observed for the Hβ emission in NGC 5548 by using Hβ and optical continuum light curves from the 13-yr (1989–2001) AGN Watch monitoring campaign. To model the breathing, we use two methods to fit the observed light curves in detail: (i) parametrized models and, (ii) the memecho reverberation-mapping code. Our models assume that optical continuum variations track the ionizing radiation, and that the Hβ variations respond with time-delays τ due to light travel-time. By fitting the data using a delay-map  Ψ(τ, F _c)  that is allowed to change with continuum flux F _c, we find that the strength of the Hβ response decreases and the time-delay increases with ionizing luminosity. The parametrized breathing models allow the time-delay and the Hβ flux to depend on the continuum flux so that,  τ∝ F ^β_c  and   F _Hβ∝ F ^α_c  . Our fits give  0.1 < β < 0.46  and  0.57 < α < 0.66. α  is consistent with previous work by Gilbert and Peterson, and Goad, Korista and Knigge. Although we find β to be flatter than previously determined by Peterson et al. using cross-correlation methods, it is closer to the predicted values from recent theoretical work by Korista and Goad.     

A spectroscopic study of IRAS F10214 + 4724

The z  = 2.286  IRAS galaxy F10214 + 4724 remains one of the most luminous galaxies in the Universe, despite its gravitational lens magnification. We present optical and near-infrared spectra of F10214 + 4724, with clear evidence for three distinct components: lines of width ∼ 1000 km s⁻¹ from a Seyfert 2 nucleus; ≲ 200 km s⁻¹ lines which are likely to be associated with star formation; and a broad (∼ 4000 km s⁻¹) C  III ] 1909-Å emission line which is blueshifted by ∼ 1000 km s⁻¹ with respect to the Seyfert 2 lines. Our study of the Seyfert 2 component leads to several new results. (i) From the double-peaked structure in the Lyα line, and the lack of Lyβ, we argue that the Lyα photons have emerged through a neutral column of N _H ∼ 2.5 × 10²⁵ m⁻², possibly located within the AGN narrow-line region, as proposed for several high-redshift radio galaxies. (ii) The resonant O  VI 1032, 1036-Å doublet (previously identified as Lyβ) is in an optically thick (1:1) ratio. At face value this implies an extreme density ( n _e ∼ 10¹⁷ m⁻³) more typical of broad-line region clouds. However, we attribute this instead to the damping wings of Lyβ from the resonant absorption. (iii) A tentative detection of He  II 1086 suggests little extinction in the rest frame ultraviolet.     

An iterative filter to reconstruct planetary transit signals in the presence of stellar variability

We report the identification, from a photometric, astrometric and spectroscopic study, of a massive white dwarf member of the nearby, approximately solar metallicity, Coma Berenices open star cluster (Melotte 111). We find the optical to near-infrared energy distribution of WD 1216+260 to be entirely consistent with that of an isolated DA and determine the effective temperature and surface gravity of this object to be   T _eff= 15 739⁺¹⁹⁷₋₁₉₆ K  and  log  g = 8.46^+0.03_−0.02  . We set tight limits on the mass of a putative cool companion,   M ≳ 0.036 M_⊙  (spatially unresolved) and   M ≳ 0.034 M_⊙  (spatially resolved and   a ≲ 2500 au  ). Based on the predictions of CO core, thick H layer evolutionary models we determine the mass and cooling time of WD 1216+260 to be   M _WD= 0.90 ± 0.04 M_⊙  and  τ_cool= 363⁺⁴⁶₋₄₁ Myr  , respectively. For an adopted cluster age of  τ= 500 ± 100 Myr  we infer the mass of its progenitor star to be   M _init= 4.77^+5.37_−0.97 M_⊙  . We briefly discuss this result in the context of the form of the stellar initial mass–final mass relation.     

Spitzer observations of acetylene bands in carbon-rich asymptotic giant branch stars in the Large Magellanic Cloud

We investigate the molecular bands in carbon-rich asymptotic giant branch (AGB) stars in the Large Magellanic Cloud (LMC), using the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope ( SST ) over the 5–38 μm range. All 26 low-resolution spectra show acetylene (C₂H₂) bands at 7 and 14 μm. The hydrogen cyanide (HCN) bands at these wavelengths are very weak or absent. This is consistent with low nitrogen abundances in the LMC. The observed 14 μm C₂H₂  band is reasonably reproduced by an excitation temperature of 500 K. There is no clear dilution of the 14 μm C₂H₂  band by circumstellar dust emission. This 14-μm band originates from molecular gas in the circumstellar envelope in these high mass-loss rate stars, in agreement with previous findings for Galactic stars. The C₂H₂ column density, derived from the 13.7 μm band, shows a gas mass-loss rate in the range 3 × 10⁻⁶ to 5 × 10⁻⁵ M_⊙ yr⁻¹. This is comparable with the total mass-loss rate of these stars estimated from the spectral energy distribution. Additionally, we compare the line strengths of the 13.7 μm C₂H₂  band of our LMC sample with those of a Galactic sample. Despite the low metallicity of the LMC, there is no clear difference in the C₂H₂  abundance among LMC and Galactic stars. This reflects the effect of the third dredge-up bringing self-produced carbon to the surface, leading to high carbon-to-oxygen ratio at low metallicity.