Dust formation and survival in supernova ejecta

The presence of dust at high redshift requires efficient condensation of grains in supernova (SN) ejecta, in accordance with current theoretical models. Yet observations of the few well-studied supernovae (SNe) and supernova remnants (SNRs) imply condensation efficiencies which are about two orders of magnitude smaller. Motivated by this tension, we have (i) revisited the model of Todini & Ferrara for dust formation in the ejecta of core collapse SNe, and (ii) followed, for the first time, the evolution of newly condensed grains from the time of formation to their survival – through the passage of the reverse shock – in the SNR. We find that  0.1–0.6  M_⊙  of dust form in the ejecta of 12–40 M_⊙ stellar progenitors. Depending on the density of the surrounding interstellar medium, between 2 and 20 per cent of the initial dust mass survives the passage of the reverse shock, on time-scales of about  4–8 × 10⁴  yr  from the stellar explosion. Sputtering by the hot gas induces a shift of the dust size distribution towards smaller grains. The resulting dust extinction curve shows a good agreement with that derived by observations of a reddened QSO at   z = 6.2  . Stochastic heating of small grains leads to a wide distribution of dust temperatures. This supports the idea that large amounts (∼0.1 M_⊙) of cold dust  ( T ∼ 40   K)  can be present in SNRs, without being in conflict with the observed infrared emission.     

Long Secondary Periods in variable red giants

We present a study of a sample of Large Magellanic Cloud red giants exhibiting Long Secondary Periods (LSPs). We use radial velocities obtained from VLT spectral observations and MACHO and OGLE light curves to examine properties of the stars and to evaluate models for the cause of LSPs. This sample is much larger than the combined previous studies of Hinkle et al. and Wood, Olivier & Kawaler.
Binary and pulsation models have enjoyed much support in recent years. Assuming stellar pulsation, we calculate from the velocity curves that the typical fractional radius change over an LSP cycle is greater than 30 per cent. This should lead to large changes in T _eff that are not observed. Also, the small light amplitude of these stars seems inconsistent with the radius amplitude. We conclude that pulsation is not a likely explanation for the LSPs. The main alternative, physical movement of the star – binary motion – also has severe problems. If the velocity variations are due to binary motion, the distribution of the angle of periastron in our large sample of stars has a probability of  1.4 × 10⁻³  that it comes from randomly aligned binary orbits. In addition, we calculate a typical companion mass of  0.09 M_⊙  . Less than 1 per cent of low-mass main-sequence stars have companions near this mass  (0.06–0.12 M_⊙)  whereas ∼25–50 per cent of low-mass red giants end up with LSPs. We are unable to find a suitable model for the LSPs and conclude by listing their known properties.     

The SCUBA HAlf Degree Extragalactic Survey (SHADES) – VII. Optical/IR photometry and stellar masses of submillimetre galaxies

We present estimates of the photometric redshifts, stellar masses and star formation histories of sources in the Submillimetre Common-User Bolometer Array (SCUBA) HAlf Degree Extragalactic Survey (SHADES). This paper describes the 60 SCUBA sources detected in the Lockman Hole covering an area of ∼320 arcmin². Using photometry spanning the B band to 8 μm, we find that the average SCUBA source forms a significant fraction of its stars in an early period of star formation and that most of the remainder forms in a shorter more intense burst around the redshift it is observed. This trend does not vary significantly with source redshift. However, the sources show a clear increase in stellar mass with redshift, consistent with downsizing. In terms of spectral energy distribution types, only two out of the 51 sources we have obtained photometric redshifts for are best fitted by a quasar-like spectrum, with approximately 80 per cent of the sources being best fitted with late-type spectra (Sc, Im and starburst). By including photometry at 850 μm, we conclude that the average SCUBA source is forming stars at a rate somewhere between 6 and 30 times the rate implied from the rest-frame optical in a dust obscured burst and that this burst creates 15–65 per cent of the total stellar mass. Using a simplistic calculation, we estimate from the average star formation history that between one in five and one in 15 bright  ( L _*+ 2 < L _optical < L _*− 1 mag)  galaxies in the field over the interval  0 < z < 3  will at some point in their lifetime experience a similar energetic dusty burst of star formation. Finally, we compute the evolution of the star formation rate density and find it peaks around   z ∼ 2  .     

Surface abundances of light elements for a large sample of early B-type stars – II. Basic parameters of 107 stars

Effective temperatures T _eff, surface gravities log  g and interstellar extinctions A _V are found for 107 B stars. Distances d of the stars, which are based on the derived T _eff, log  g and A _V values, show good agreement with those obtained from the Hipparcos parallaxes. Comparing the T _eff and log  g values with evolutionary computations, we infer masses, radii, luminosities, ages and relative ages of the stars. Empirical relations between the T _eff and log  g parameters, on the one hand, and the photometric indices Q , [ c ₁] and β , on the other hand, are constructed; these relations give a fast method for the T _eff and log  g estimation of early and medium B stars. Inclusion of the infrared J , H and K colours into the T _eff, log  g and A _V determination shows that (i) the T _eff and log  g parameters are altered only slightly; (ii) the A _V value is rather sensitive to these colours, so an accuracy better than 0.05 mag in the JHK data is necessary for precise A _V evaluation.     

NLTE models of line-driven stellar winds – III. Influence of X-ray radiation on wind structure of O stars

We study the influence of X-rays on the wind structure of selected O stars. For this purpose we use our non-local thermodynamic equilibrium (NLTE) wind code with inclusion of additional artificial source of X-rays, assumed to originate in the wind shocks.
We show that the influence of shock X-ray emission on wind mass-loss rate is relatively small. Wind terminal velocity may be slightly influenced by the presence of strong X-ray sources, especially for stars cooler than   T _eff≲ 35 000 K  .
We discuss the origin of the   L _X/ L ∼ 10⁻⁷  relation. For stars with thick wind this relation can be explained assuming that the cooling time depends on wind density. Stars with optically thin winds exhibiting the 'weak wind problem' display enhanced X-ray emission which may be connected with large shock cooling length. We propose that this effect can explain the 'weak wind problem'.
Inclusion of X-rays leads to a better agreement of the model ionization structure with observations. However, we do not find any significant influence of X-rays on P  v ionization fraction implying that the presence of X-rays cannot explain the P  v problem.
We study the implications of modified ionization equilibrium due to shock emission on the line transfer in the X-ray region. We conclude that the X-ray line profiles of helium-like ions may be affected by the line absorption within the cool wind.     

The EC 14026 stars – VIII. PG 1336−018: a pulsating sdB star in an HW Vir-type eclipsing binary

The sdB star PG 1336−018 is found to be a very short-period eclipsing binary system, remarkably similar to the previously unique system HW Vir. In addition, and unlike HW Vir, the sdB star in the PG 1336 system shows rapid oscillations of the type found in the recently discovered sdB pulsators, or EC 14026 stars. The orbital period, 0.101 0174 d, is one of the shortest known for a detached binary. Analysis of photoelectric and CCD photometry reveals pulsation periods near 184 and 141 s, with semi-amplitudes of ∼0.01 and ∼0.005 mag respectively. Both oscillations might have variable amplitude, and it is probable that other frequencies are present with amplitudes ∼0.003 mag or less. The 184- and 141-s pulsations are in the range of periods predicted by models for hot horizontal-branch stars. Analysis of medium-dispersion spectrograms yields T _eff=33 000±1000 K and log g =5.7±0.1 for the sdB primary star, a radial velocity semi-amplitude K ₁=78±3 km s⁻¹ and a system velocity γ=6±2 km s⁻¹. Spectrograms from the IUE Final Archive give T _eff=33 000±3000 K and E ( B − V )=0.05 for log g =6.0 models. The derived angular radius leads to a distance of 710±50 pc for the system, and an absolute magnitude for the sdB star of +4.1±0.2. A preliminary analysis of U , V and R light curves indicates the orbital inclination to be near 81° and the relative radii to be r ₁=0.19 and r ₂=0.205. Assuming the mass of the sdB primary to be 0.5 M⊙ leads to a mass ratio q =0.3 for the system, and indicates that the secondary is a late-type dwarf of type ∼M5. As with HW Vir, it is necessary to invoke small limb-darkening coefficients and high albedos for the secondary star to obtain reasonable fits to the observed light curves.     

Collisions and close encounters involving massive main-sequence stars

We study close encounters involving massive main-sequence stars and the evolution of the exotic products of these encounters as common-envelope systems or possible hypernova progenitors. We show that parabolic encounters between low- and high-mass stars and between two high-mass stars with small periastrons result in mergers on time-scales of a few tens of stellar free-fall times (a few tens of hours). We show that such mergers of unevolved low-mass stars with evolved high-mass stars result in little mass-loss  (∼0.01 M_⊙)  and can deliver sufficient fresh hydrogen to the core of the collision product to allow the collision product to burn for several million years. We find that grazing encounters enter a common-envelope phase which may expel the envelope of the merger product. The deposition of energy in the envelopes of our merger products causes them to swell by factors of ∼100. If these remnants exist in very densely populated environments  ( n ≳ 10⁷ pc⁻³)  , they will suffer further collisions which may drive off their envelopes, leaving behind hard binaries. We show that the products of collisions have cores rotating sufficiently rapidly to make them candidate hypernova/gamma-ray burst progenitors and that ∼0.1 per cent of massive stars may suffer collisions, sufficient for such events to contribute significantly to the observed rates of hypernovae and gamma-ray bursts.     

Towards a spectral technique for determining material geometry around evolved stars: application to HD 179821

HD 179821 is an evolved star of unknown progenitor mass range (either post-asymptotic giant branch or post-red supergiant) exhibiting a double-peaked spectral energy distribution (SED) with a sharp rise from ∼8 to  20 μm  . Such features have been associated with ejected dust shells or inwardly truncated circumstellar discs. In order to compare SEDs from both systems, we employ a spherically symmetric radiative transfer code and compare it to a radiative, inwardly truncated disc code. As a case study, we model the broad-band SED of HD 179821 using both codes. Shortward of 40 μm, we find that both models produce equivalent fits to the data. However, longward of 40 μm, the radial density distribution and corresponding broad range of disc temperatures produce excess emission above our spherically symmetric solutions and observations. For HD 179821, our best fit consists of a   T _eff= 7000 K  central source characterized by  τ _V ∼ 1.95  and surrounded by a radiatively driven, spherically symmetric dust shell. The extinction of the central source reddens the broad-band colours so that they resemble a   T _eff= 5750 K  photosphere. We believe that HD 179821 contains a hotter central star than previously thought. Our results provide an initial step towards a technique to distinguish geometric differences from spectral modelling.     

Infrared helium–hydrogen line ratios as a measure of stellar effective temperature

We have observed a large sample of compact planetary nebulae in the near-infrared to determine how the 2¹P–2¹S He  i line at 2.058 μm varies as a function of stellar effective temperature, T _eff. The ratio of this line with H  i Br γ at 2.166 μm has often been used as a measure of the highest T _eff present in a stellar cluster, and hence of whether there is a cut-off in the stellar initial mass function at high masses. However, recent photoionization modelling has revealed that the behaviour of this line is more complex than previously anticipated. Our work shows that in most aspects the photoionization models are correct. In particular, we confirm the weakening of the 2¹P–2¹S line as T _eff increases beyond 40 000 K. However, in many cases the model underpredicts the observed ratio when we consider the detailed physical conditions in the individual planetary nebulae. Furthermore, there is evidence that there is still significant 2¹P–2¹S He  i line emission even in the planetary nebulae with very hot     central stars. It is clear from our work that this ratio cannot be considered as a reliable measure of effective temperature on its own.     

On the possibility of observing H2 emission from primordial molecular cloud kernels

We study the prospects for observing H₂ emission during the assembly of primordial molecular cloud kernels. The primordial molecular cloud cores, which resemble those at the present epoch, can emerge around  1+ z ∼20  according to recent numerical simulations. The kernels form inside the cores, and the first stars will appear inside the kernels. A kernel typically contracts to form one of the first generation stars with an accretion rate that is as large as ∼0.01 M_⊙ yr⁻¹. This occurs owing to the primordial abundances, which result in a kernel temperature of order 1000 K, and the collapsing kernel emits H₂ line radiation at a rate ∼10³⁵ erg s⁻¹. Predominantly   J =5-3   ( v =0)  rotational emission of H₂ is expected. At redshift  1+ z ∼20  , the expected flux is ∼0.01 μJy for a single kernel. While an individual object is not observable by any facilities available in the near future, the expected assembly of primordial star clusters on subgalactic scales can result in fluxes at the sub-mJy level. This is marginally observable with ASTRO-F and ALMA. We also examine the rotational   J =2-0   ( v =0)  and vibrational   δv =1  emission lines. The former may possibly be detectable with ALMA.     

M4–18: the planetary nebula and its WC10 central star

We present a detailed analysis of the planetary nebula M4–18 (G146.7+07.6) and its WC10-type Wolf–Rayet (WR) central star, based on high‐quality optical spectroscopy (WHT/UES, INT/IDS, WIYN/DensPak) and imaging ( HST /WFPC2). From a non-LTE model atmosphere analysis of the stellar spectrum, we derive T _eff=31 kK,     v _∞=160 km s⁻¹ and abundance number ratios of H/He<0.5, C/He=0.60 and O/He=0.10. These parameters are remarkably similar to those of He 2–113 ([WC10]). Assuming an identical stellar mass to that determined by De Marco et al. for He 2–113, we obtain a distance of 6.8 kpc to M4–18 [ E ( B−V )=0.55 mag from nebular and stellar techniques]. This implies that the planetary nebula of M4–18 has a dynamical age of ∼3100 yr, in contrast to ≥270 yr for He 2–113. This is supported by the much higher electron density of the latter. These observations may be reconciled with evolutionary predictions only if [WC]-type stars exhibit a range in stellar masses.
Photoionization modelling of M4–18 is carried out using our stellar WR flux distribution, together with blackbody and Kurucz energy distributions obtained from Zanstra analyses. We conclude that the ionizing energy distribution from the WR model provides the best consistency with the observed nebular properties, although discrepancies remain.     

The far-infrared signature of dust in high-latitude regions

We present ISOPHOT observations of eight interstellar regions in the 60–200 μm wavelength range. The regions belong to mostly quiescent high-latitude clouds and have optical extinction peaks from   A_V ∼1–6 mag  . From the 150- and 200-μm emission, we derived colour temperatures for the classical big grain component which show a clear trend of decreasing temperature with increasing 200-μm emission. The 200-μm emission per unit   A_V   , however, does not drop at lower temperatures. This fact can be interpreted in terms of an increased far-infrared (FIR) emissivity of the big grains. We developed a two-component model including warm dust with the temperature of the diffuse interstellar medium (ISM) of   T = 17.5 K  , and cold dust with   T = 13.5 K  and FIR emissivity increased by a factor of >4. A mixture of the two components can reproduce the observed colour variations and the ratios   I ₂₀₀/ A_V   and  τ₂₀₀/ A_V   . The relative abundance of small grains with respect to the big grains shows significant variations from region to region at low column densities. However, in lines of sight of higher column density, our data indicate the disappearance of small grains, perhaps a signature of adsorption/coagulation of dust. The larger size and porous structure could also explain the increased FIR emissivity. Our results from eight independent regions suggest that these grains might be ubiquitous in the galactic ISM.     

Upper limit on dimming of cosmological sources by intergalactic grey dust from the soft X-ray background

Active galactic nuclei (AGN) produce a dominant fraction  ( F _AGN∼ 80 per cent)  of the soft X-ray background (SXB) at photon energies  0.5 < E < 2 keV  . If dust pervaded throughout the intergalactic medium, its scattering opacity would have produced diffuse X-ray haloes around AGN. Taking account of known galaxies and galaxy clusters, only a fraction   F _halo≲ 10 per cent  of the SXB can be in the form of diffuse X-ray haloes around AGN. We therefore limit the intergalactic opacity to optical/infrared photons from large dust grains, with radii in the range   a = 0.2–2.0 μm  , to a level  τ_GD≲ 0.15( F _halo/10 per cent)(F_AGN/80 per cent)⁻¹  to a redshift   z ∼ 1  . Our results are only weakly dependent on the grain size distribution in this size range or the redshift evolution of the intergalactic dust. Stacking X-ray images of AGN can be used to improve our constraints and diminish the importance of dust as a source of systematic uncertainty for future supernova surveys which aim to improve the precision on measuring the redshift evolution of the dark energy equation-of-state.     

Extended dust emission in NGC 7465

We present SCUBA 850-μm, JCMT  CO( J =2→1)  , B -band imaging and VLA H  i observations of the NGC 7465/4/3 group of galaxies. The 850-μm emission associated with NGC 7465 extends to at least ∼2 R ₂₅ and is well correlated with the H  i . We investigate a range of possible mechanisms by which dust beyond R ₂₅ may be heated to give the observed extended submillimetre emission. By modelling the dust heating by stars in two extreme geometries, we fail to find any reasonable star formation scenario that is consistent with both the 850-μm and optical data. Furthermore, we do not detect any  CO( J =2→1)  emission coincident with the extended dust and atomic gas as would be expected if significant star formation were occurring. We show that shock-heating of dust via cloud–cloud collisions in the stripped interstellar medium of NGC 7465 could be sufficient to explain the extended 850-μm emission and lack of optical emission in the stripped gas, and suggest that cloud–cloud collisions may be an important dust heating mechanism in gas-rich systems.     

Voids in a ΛCDM universe

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular Λ cold dark matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of −0.8 or less – roughly the equivalent of using a spherical overdensity group finder for haloes. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement when differences in particle and grid resolution are accounted for. Voids very clearly correspond to minima in the smoothed initial density field. Apart from a very weak dependence on the mass resolution, the rescaled mass profiles of voids in the different simulations agree remarkably well. We find a universal void mass profile of the form  ρ(< r )/ρ( r _eff) ∝ exp[( r / r _eff)^α]  , where r _eff is the effective radius of a void and  α∼ 2  . The mass function of haloes in voids is steeper than that of haloes that populate denser regions. In addition, the abundances of void haloes seem to evolve somewhat more strongly between redshifts ∼1 and 0 than the global abundances of haloes.     

Dust coagulation in star formation with different metallicities

Dust grains coagulate into larger aggregates in dense gas. This changes their size distribution and possibly affects the thermal evolution of star-forming clouds. We here investigate dust coagulation in collapsing pre-stellar cores with different metallicities by considering the thermal motions of grains. We show that coagulation does occur even at low metallicity  ∼10⁻⁶ Z_⊙  . However, we also find (i) that the H₂ formation rate on dust grains is reduced only after the majority of H₂ is formed and (ii) that the dust opacity is modified only after the core becomes optically thick. Therefore, we conclude that the effects of dust coagulation can safely be neglected in discussing the temperature evolution of the pre-stellar cores for any metallicity as long as the grain motions are thermal.     

The local star formation rate and radio luminosity density

We present a new determination of the local volume-averaged star formation rate from the 1.4-GHz luminosity function of star forming galaxies. Our sample, taken from the   B ≤12  Revised Shapley–Ames catalogue (231 normal spiral galaxies over an effective area of 7.1 sr) has ≃100 per cent complete radio detections and is insensitive to dust obscuration and cirrus contamination. After removal of known active galaxies, the best-fitting Schechter function has a faint-end slope of  −1.27±0.07  in agreement with the local H α luminosity function, characteristic luminosity   L ∗=(2.6±0.7)×10²² W Hz⁻¹  and density   φ ∗=(4.8±1.1)×10⁻⁴ Mpc⁻³.  The inferred local radio luminosity density of  (1.73±0.37±0.03)×10¹⁹ W Hz⁻¹ Mpc⁻³  (Poisson noise, large-scale structure fluctuations) implies a volume-averaged star formation rate ∼2 times larger than the Gallego et al. H α estimate, i.e.   ρ _1.4 GHz=(2.10±0.45±0.04)×10⁻² M_⊙ yr⁻¹ Mpc⁻³  for a Salpeter initial mass function from  0.1–125 M_⊙  and Hubble constant of 50 km s⁻¹ Mpc⁻¹. We demonstrate that the Balmer decrement is a highly unreliable extinction estimator, and argue that optical–ultraviolet (UV) star formation rates (SFRs) are easily underestimated, particularly at high redshift.     

On the time dependence of differential rotation in young late-type stars

A model for the angular momentum transfer within the convection zone of a rapidly rotating star is introduced and applied to the analysis of recent observations of temporal fluctuations of the differential rotation on the young late-type stars AB Doradus (AB Dor) and LQ Hydrae (LQ Hya). Under the hypothesis that the mean magnetic field produced by the stellar dynamo rules the angular momentum exchanges and that the angular velocity depends only on the distance s from the rotation axis and the time, the minimum azimuthal Maxwell stress  | B_sB _φ|  , averaged over the convection zone, is found to range from ∼0.04 to  ∼0.14 T²  . If the poloidal mean magnetic field   B _s   is of the order of 0.01 T, as indicated by the Zeeman–Doppler imaging maps of those stars, then the azimuthal mean field   B _φ  can reach an intensity of several teslas, which significantly exceeds equipartition with the turbulent kinetic energy. Such strong fields can account also for the orbital period modulation observed in cataclysmic variables and RS Canum Venaticorum systems with a main-sequence secondary component. Moreover, the model allows us to compute the kinetic energy dissipation rate during the maintenance of the differential rotation. Only in the case of the largest surface shear observed on LQ Hya may the dissipated power exceed the stellar luminosity, but the lack of a sufficient statistic on the occurrence of such episodes of large shear does not allow us to estimate their impact on the energy budget of the convection zone.     

The European Large Area ISO Survey: optical identifications of 15-μm and 1.4-GHz sources in N1 and N2

We present the multiwavelength properties and catalogue of the 15 μm and 1.4 GHz radio sources detected in the European Large Area ISO Survey ( ELAIS ) areas N1 and N2. Using the optical data from the Wide Field Survey we use a likelihood ratio method to search for the counterparts of the 1056 and 691 sources detected at 15 μm and 1.4 GHz, respectively, down to flux limits of   S ₁₅= 0.5 mJy  and   S _1.4 GHz= 0.135 mJy  . We find that ∼92 per cent of the 15 μm ELAIS sources have an optical counterpart down to   r '= 24  . All mid-infrared (IR) sources with fluxes   S ₁₅≥ 3 mJy  have an optical counterpart. The magnitude distribution of the sources shows a well-defined peak at relatively bright magnitudes   r '∼ 18  . The mid-IR-to-optical and radio-to-optical flux diagrams are presented and discussed in terms of actual galaxy models. About 15 per cent of the sources are bright galactic stars; of the extragalactic objects ∼65 per cent are compatible with being normal or starburst galaxies and ∼25 per cent active galactic nuclei (AGNs). Objects with mid-IR-to-optical fluxes larger than 100 are found, comprising ∼20 per cent of the sample. We suggest that that these sources are highly obscured luminous and ultraluminous starburst galaxies and AGNs.     

Dust in the core of the metal-rich globular cluster NGC 6356

We report the 3.5σ detection of dust in the core of the metal-rich ([Fe/H] =−0.4) globular cluster NGC 6356. The dust mass in the core is ∼ 4–17 × 10⁻³ M⊙, depending on the dust equilibrium temperature. We rule out sputtering as a grain destruction mechanism in NGC 6356, unless grains are very small (∼30 Å). We also place the upper limit of 6.6 × 10⁻³ on the CO-to-dust ratio in this cluster. This value is significantly smaller than the CO-to-dust ratio in the general Galactic interstellar medium, and in the outflows of evolved stars.