A new approach to the calculation of the cratering rate of the Earth over the last 125 ± 20 Myr

The recent cratering record of the surface of the Earth is re-examined using a new technique that concentrates on estimating the mean areas occupied by individual craters, together with the gradients of linear plots of crater numbers versus crater ages. This analysis indicates that the lower limit of the rate at which craters have been produced over the last 125±20 Myr is, for example, (12.0±0.7)×10⁻¹⁵ km⁻² yr⁻¹ for D 2.4 km craters, (9.5±0.6)×10⁻¹⁵ km⁻² yr⁻¹ for D 5.0 km craters, (6.5±0.5)×10⁻¹⁵ km⁻² yr⁻¹ for D 12 km craters, and (3.0±0.3)×10⁻¹⁵ km⁻² yr⁻¹ for D 22 km craters. These figures indicate that previous researchers have considerably overestimated the rate at which small (2.4< D <20 km) craters are being produced. It is also found that the relationship between crater production rate and crater diameter is not a simple power law in the 2.4< D <40 km diameter range. On the most stable areas of the Earth's continents, and over the last 125±20 Myr it seems that the rate at which craters are eroded below the detection limit does not depend on crater diameter throughout the above size range.     

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.     

Type IIP supernova SN 2004et: a multiwavelength study in X-ray, optical and radio

We present X-ray, broad-band optical and low-frequency radio observations of the bright type IIP supernova SN 2004et. The Chandra X-ray Observatory observed the supernova at three epochs, and the optical coverage spans a period of ∼470 d since explosion. The X-ray emission softens with time, and we characterize the X-ray luminosity evolution as   L _X∝ t ^−0.4  . We use the observed X-ray luminosity to estimate a mass-loss rate for the progenitor star of  ∼2 × 10⁻⁶ M_⊙ yr⁻¹  . The optical light curve shows a pronounced plateau lasting for about 110 d. Temporal evolution of photospheric radius and colour temperature during the plateau phase is determined by making blackbody fits. We estimate the ejected mass of ⁵⁶Ni to be  0.06 ± 0.03 M_⊙  . Using the expressions of Litvinova & Nadëzhin we estimate an explosion energy of  (0.98 ± 0.25) × 10⁵¹ erg  . We also present a single epoch radio observation of SN 2004et. We compare this with the predictions of the model proposed by Chevalier, Fransson & Nymark. These multiwavelength studies suggest a main-sequence progenitor mass of  ∼20 M_⊙  for SN 2004et.     

Radiative association of C and P, and Si and P atoms

The rate coefficients for the formation of carbon monophosphide (CP) and silicon monophosphide (SiP) by radiative association are estimated for temperatures ranging from 300 to 14 100 K. In this temperature range, the radiative association rate coefficients are found to vary from  1.14 × 10⁻¹⁸  to  1.62 × 10⁻¹⁸ cm³ s⁻¹  and from  3.73 × 10⁻²⁰  to  7.03 × 10⁻²⁰ cm³ s⁻¹  for CP and SiP, respectively. In both cases, rate coefficients increase slowly with the increase in temperature.     

An alternative symbiotic channel to Type Ia supernovae

By assuming an aspherical stellar wind with an equatorial disc from a red giant, we investigate the production of Type Ia supernovae (SNe Ia) via a symbiotic channel. We estimate that the Galactic birthrate of SNe Ia via the symbiotic channel is between  1.03 × 10⁻³  and  2.27 × 10⁻⁵ yr⁻¹  , while the delay time of SNe Ia has a wide range from ∼0.07 to 5 Gyr. The results are greatly affected by the outflow velocity and mass-loss rate of the equatorial disc. Using our model, we discuss the progenitors of SN 2002ic and SN 2006X.     

X-ray analysis of the cluster 3C 129

The cluster 3C 129 is classified as a rich cluster. An analysis of the properties of the cluster 3C 129 from ROSAT PSPC and HRI, Einstein IPC, and EXOSAT ME observations is presented. The mean temperature from a joint fit of the ROSAT PSPC and EXOSAT ME data is 5.5(±0.2) keV. The luminosity is 0.6×10⁴⁴ erg s⁻¹ in 0.2–2.4 keV and 2.7×10⁴⁴ erg s⁻¹ in 0.2–10 keV. We find a cooling flow with a rate of ∼84 M_⊙ yr⁻¹. The central gas density is 6×10⁻³ cm⁻³, and the ICM mass is 3.6×10¹³ M_⊙. The total cluster mass is ∼5×10¹⁴ M_⊙. The X-ray morphology shows an east–west elongation, which is evidence for a recent merger event. The radio source 3C 129.1 is located near the X-ray centre. Another cluster member galaxy (the radio galaxy 3C 129) is a prototype of head-tailed radio galaxies, and is located in the west part of the cluster. The tail points along the gradient of intracluster gas pressure. There are no significant point X-ray sources associated with the AGNs of the two radio galaxies.     

The continuing saga of Sakurai's object (V4334 Sgr): dust production and helium line emission

We report further UKIRT spectroscopic observations of Sakurai's object (V4334 Sgr) made in 1999 April/May in the 1–4.75 μm range, and find that the emission is dominated by amorphous carbon at T _d~600 K. The estimated maximum grain size is 0.6 μm, and the mass lower limit is 1.7±0.2×10⁻⁸ M_⊙ to 8.9±0.6×10⁻⁷ M_⊙ for distances of 1.1–8 kpc. For 3.8 kpc the mass is 2.0±0.1×10⁻⁷ M_⊙.
We also report strong He  i emission at 1.083 μm, in contrast to the strong absorption in this line in 1998. We conclude that the excitation is collisional, and is probably caused by a wind, consistent with the P Cygni profile observed by Eyres et al. in 1998.     

Measurement of the Sunyaev–Zel'dovich increment in massive galaxy clusters

We have detected the Sunyaev–Zel'dovich (SZ) increment at 850 μm in two galaxy clusters (Cl 0016+16 and MS 1054.4−0321) using the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. Fits to the isothermal β model yield a central Compton y parameter of  (2.2 ± 0.7) × 10⁻⁴  and a central 850-μm flux of  Δ I ₀= 2.2 ± 0.7 mJy beam⁻¹  in Cl 0016. This can be combined with decrement measurements to infer   y = (2.38 ±^0.36_0.34) × 10⁻⁴  and   v _pec= 400±¹⁹⁰⁰₁₄₀₀ km s⁻¹  . In MS 1054 we find a peak 850-μm flux of  Δ I ₀= 2.0 ± 1.0 mJy beam⁻¹  and   y = (2.0 ± 1.0) × 10⁻⁴  . To be successful such measurements require large chop throws and non-standard data analysis techniques. In particular, the 450-μm data are used to remove atmospheric variations in the 850-μm data. An explicit annular model is fit to the SCUBA difference data in order to extract the radial profile, and separately fit to the model differences to minimize the effect of correlations induced by our scanning strategy. We have demonstrated that with sufficient care, SCUBA can be used to measure the SZ increment in massive, compact galaxy clusters.     

Gravitational wave signals from the collapse of the first stars

We study the gravitational wave emission from the first stars, which are assumed to be very massive objects (VMOs). We take into account various feedback (both radiative and stellar) effects regulating the collapse of objects in the early Universe and thus derive the VMO initial mass function and formation rate. If the final fate of VMOs is to collapse, leaving very massive black hole remnants, then the gravitational waves emitted during each collapse would be seen as a stochastic background. The predicted spectral strain amplitude in a critical density cold dark matter (CDM) universe peaks in the frequency range ν ≈5×10⁻⁴–5×10⁻³ Hz, where it has a value in the range ≈10⁻²⁰–10⁻¹⁹ Hz^−1/2, and might be detected by the Laser Interferometer Space Antenna ( LISA ). The expected emission rate is roughly 4000 event yr⁻¹, resulting in a stationary discrete sequence of bursts, i.e. a shot-noise signal.     

The temperature of the warm neutral medium in the Milky Way

We report high-spectral-resolution Australia Telescope Compact Array (ATCA) H  i 21-cm observations resulting in the detection of the warm neutral medium (WNM) of the Galaxy in absorption against two extragalactic radio sources, PKS 1814−637 and PKS 0407−658. The two lines of sight were selected on the basis of the simplicity of their absorption profiles and the strength of the background sources; the high velocity resolution of the spectra then enabled us to estimate the kinetic temperatures of the absorbing gas by fitting multiple Gaussians to the absorption profiles. Four separate WNM components were detected towards the two sources, with peak optical depths  τ_max= (1.0 ± 0.08) × 10⁻², (1.4 ± 0.2) × 10⁻³, (2.2 ± 0.5) × 10⁻³  and  (3.4 ± 0.5) × 10⁻³  and kinetic temperatures   T _k= 3127 ± 300, 3694 ± 1595, 3500 ± 1354  and  2165 ± 608 K  , respectively. All four components were thus found to have temperatures in the thermally unstable range  500 < T _k < 5000 K  ; this suggests that thermal equilibrium has not been reached throughout the WNM.     

Episodic dust formation by HD 192641 (WR 137) – II

We present new infrared photometry of the WC7-type Wolf–Rayet star HD 192641 (WR 137) from 1985 to 1999. These data track the cooling of the dust cloud formed in the 1982–84 dust-formation episode from 1985 to 1991, the increase of the infrared flux from 1994.5 to a new dust-formation maximum in 1997 and its subsequent fading. From these and earlier data we derive a period of 4765±50 d (13.05±0.15 yr) for the dust-formation episodes. Between dust-emission episodes, the infrared spectral energy distribution has the form of a power law, λF _λ ∝ λ ^−1.86. The rising branch of the infrared light curve (1994–97) differs in form from that of the episodic dust-maker WR 125. Time-dependent modelling shows that this difference can be attributed to a different time dependence of dust formation in WR 137, which occurred approximately ∝ t ² until maximum, whereas that of WR 125 could be described by a step function, akin to a threshold effect. For an adopted distance of 1.6 kpc, the rate of dust formation was found to be 5.0×10⁻⁸ M_⊙ yr⁻¹ at maximum, accounting for a fraction f _C≈1.5×10⁻³ of the carbon flowing in the stellar wind. The fading branches of the light curves show evidence for secondary 'mini-eruptions' in 1987, 1988 and 1990, behaviour very different from that of the prototypical episodic dust-maker HD 193793 (WR 140), and suggesting the presence in the WR 137 stellar wind of large-scale structures that are crossed by the wind–wind collision region.     

Long-term evolution of Oort Cloud comets: capture of comets

We test different possibilities for the origin of short-period comets captured from the Oort Cloud. We use an efficient Monte Carlo simulation method that takes into account non-gravitational forces, Galactic perturbations, observational selection effects, physical evolution and tidal splittings of comets. We confirm previous results and conclude that the Jupiter family comets cannot originate in the spherically distributed Oort Cloud, since there is no physically possible model of how these comets can be captured from the Oort Cloud flux and produce the observed inclination and Tisserand constant distributions. The extended model of the Oort Cloud predicted by the planetesimal theory consisting of a non-randomly distributed inner core and a classical Oort Cloud also cannot explain the observed distributions of Jupiter family comets. The number of comets captured from the outer region of the Solar system are too high compared with the observations if the inclination distribution of Jupiter family comets is matched with the observed distribution. It is very likely that the Halley-type comets are captured mainly from the classical Oort Cloud, since the distributions in inclination and Tisserand value can be fitted to the observed distributions with very high confidence. Also the expected number of comets is in agreement with the observations when physical evolution of the comets is included. However, the solution is not unique, and other more complicated models can also explain the observed properties of Halley-type comets. The existence of Jupiter family comets can be explained only if they are captured from the extended disc of comets with semimajor axes of the comets   a <5000 au  . The original flattened distribution of comets is conserved as the cometary orbits evolve from the outer Solar system era to the observed region.     

Spectroscopic and photometric observations of the eclipsing star UV Leo

High-resolution spectroscopic observations around the Hα line and BVRI photometry of the eclipsing short-period RS CVn star UV Leo are presented. The simultaneous light-curve solution and radial velocity-curve solution led to the following values of the global parameters of the binary: temperatures   T ₁= 6000 ± 100 K  and   T ₂= 5970 ± 20 K  ; masses   M ₁= 0.976 ± 0.067 M_⊙  and   M ₂= 0.931 ± 0.052 M_⊙  ; separation   a = 3.716 ± 0.048 R_⊙  ; orbital inclination     ; radii   R ₁= 1.115 ± 0.052 R_⊙  and   R ₂= 1.078 ± 0.051 R_⊙  ; equatorial velocities   V ₁= 98.8 ± 2.3 km s⁻¹  and   V ₂= 89.6 ± 2.7 km s⁻¹  . These results lead to the conclusion that the two components of UV Leo are slightly oversized for their masses and lie within the main-sequence band on the mass–radius diagram, close to the isochrone 9 × 10¹⁰ yr.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.     

Cygnus X-3 with ISO: investigating the wind

We have observed the energetic binary Cygnus X-3 in both quiescent and flaring states between 4 and 16 μm using the ISO satellite. We find that the quiescent source shows the thermal free–free spectrum typical of a hot, fast stellar wind, such as from a massive helium star. The quiescent mass-loss rate arising from a spherically symmetric, non-accelerating wind is found to be in the range (0.4–2.9)×10⁻⁴ M_⊙ yr⁻¹, consistent with other infrared and radio observations, but considerably larger than the 10⁻⁵ M_⊙ yr⁻¹ deduced from both the orbital change and the X-ray column density. There is rapid, large-amplitude flaring at 4.5 and 11.5 μm at the same time as enhanced radio and X-ray activity, with the infrared spectrum apparently becoming flatter in the flaring state. We believe that non-thermal processes are operating, perhaps along with enhanced thermal emission.     

Does Collinder 236 host a Cepheid calibrator?

Photoelectric UBV photometry and star counts are presented for the previously unstudied open cluster Collinder 236, supplemented by observations for stars near the Cepheid WZ Car. Collinder 236 is typical of groups associated with Cepheids, with an evolutionary age of  (3.2 ± 1.1) × 10⁷ yr  , but it is  1944 ± 71  pc distant, only half the predicted distance to WZ Car. The cluster is reddened by   E ( B − V ) ≃ 0.26  , and has nuclear and coronal radii of   r _n≃ 2 arcmin (1.1 pc)  and   R _c≃ 8 arcmin (4.5 pc)  , respectively. The Cepheid is not a member of Collinder 236 on the basis of location beyond the cluster tidal radius and implied distance, but its space reddening can be established as   E ( B − V ) = 0.268 ± 0.006  s.e. from five adjacent stars. Period changes in WZ Car studied with the aid of archival data are revised. The period of WZ Car is increasing, its rate of  +8.27 ± 0.19 s yr⁻¹  being consistent with a third crossing of the instability strip.     

Magnetically collimated jets with high mass flux

We present HST /WFPC2 observations of UGC 4483, an irregular galaxy in the M81/NGC 2403 complex. Stellar photometry was carried out with HSTphot, and is complete to V ≃26.0 and I ≃24.7. We measure the red giant branch tip at I =23.56±0.10, and calculate a distance modulus of μ ₀=27.53±0.12 (corresponding to a distance of 3.2±0.2 Mpc), placing UGC 4483 within the NGC 2403 subgroup. We were able to measure properties of a previously known young star cluster in UGC 4483, finding integrated magnitudes of V =18.66±0.21 and I =18.54±0.10 for the stellar contribution (integrated light minus H α and [O  iii ] contribution), corresponding to an age of ∼10–15 Myr and an initial mass of ∼10⁴ M_⊙. This is consistent with the properties of the cluster's brightest stars, which were resolved in the data for the first time. Finally, a numerical analysis of the galaxy's stellar content yields a roughly constant star formation rate of 1.3×10⁻³ M_⊙ yr⁻¹ and mean metallicity of [Fe/H]=−1.3 dex from 15 Gyr ago to the present.     

An extremely slow nova?

The binary companion to the peculiar F supergiant HD 172481 is shown to be a Mira variable with a pulsation period of 312 d. Its characteristics are within the normal range found for solitary Miras of that period, although its pulsation amplitude and mass-loss rate ̇ ∼3×10⁻⁶ M_⊙ yr⁻¹ are higher than average. Reasons are given for suspecting that the F supergiant, which has L ∼10⁴ L_⊙, is a white dwarf burning hydrogen accreted from its companion.     

Population III stars: hidden or disappeared?

A Population III/Population II transition from massive to normal stars is predicted to occur when the metallicity of the star-forming gas crosses the critical range   Z _cr= 10^−5±1 Z_⊙  . To investigate the cosmic implications of such a process, we use numerical simulations which follow the evolution, metal enrichment and energy deposition of both Population II and Population III stars. We find that: (i) due to inefficient heavy element transport by outflows and slow 'genetic' transmission during hierarchical growth, large fluctuations around the average metallicity arise; as a result, Population III star formation continues down to   z = 2.5  , but at a low peak rate of  10⁻⁵ M_⊙ yr⁻¹ Mpc⁻³  occurring at   z ≈ 6  (about 10⁻⁴ of the Population II one); and (ii) Population III star formation proceeds in an 'inside–out' mode in which formation sites are progressively confined to the periphery of collapsed structures, where the low gas density and correspondingly long free-fall time-scales result in a very inefficient astration. These conclusions strongly encourage deep searches for pristine star formation sites at moderate  (2 < z < 5)  redshifts where metal-free stars are likely to be hidden.     

The space density of cataclysmic variables: constraints from the ROSAT North Ecliptic Pole survey

We use the ROSAT North Ecliptic Pole (NEP) survey to construct a small, but purely X-ray flux-limited sample of cataclysmic variable stars (CVs). The sample includes only four systems, two of which (RX J1715.6+6856 and RX J1831.7+6511) are new discoveries. We present time-resolved spectroscopy of the new CVs and measure orbital periods of 1.64 ± 0.02 and 4.01 ± 0.03 h for RX J1715.6+6856 and RX J1831.7+6511, respectively. We also estimate distances for all the CVs in our sample, based mainly on their apparent brightness in the infrared. The space density of the CV population represented by our small sample is  1.1^+2.3_−0.7× 10⁻⁵ pc⁻³  . We can also place upper limits on the space density of any subpopulation of CVs too faint to be included in the NEP survey. In particular, we show that if the overall space density of CVs is as high as  2 × 10⁻⁴ pc⁻³  (as has been predicted theoretically), the vast majority of CVs must be fainter than   L _X≃ 2 × 10²⁹ erg s⁻¹  .