Relativistic ionized accretion disc models of MCG–6-30-15

We present BeppoSAX observations of Nova Velorum 1999 (V382 Vel), carried out in a broad X-ray band covering 0.1–300 keV only 15 d after the discovery and again after 6 months. The nova was detected at day 15 with the BeppoSAX instruments which measured a flux F _x≃1.8×10⁻¹¹ erg cm⁻² s⁻¹ in the 0.1–10 keV range and a 2 σ upper limit F _x<6.7×10⁻¹² erg cm⁻² s⁻¹ in the 15–60 keV range. We attribute the emission to shocked nebular ejecta at a plasma temperature kT ≃6 keV . At six months no bright component emerged in the 15–60 keV range, but a bright central supersoft X-ray source appeared. The hot nebular component previously detected had cooled to a plasma temperature kT <1 keV . There was strong intrinsic absorption of the ejecta in the first observation and not in the second, because the column density of neutral hydrogen decreased from N (H)≃1.7×10²³ to N (H)≃10²¹ cm⁻² (close to the interstellar value). The unabsorbed X-ray flux also decreased from F _x=4.3×10⁻¹¹ to F _x≃10⁻¹² erg cm⁻² s⁻¹ .     

Neutral sodium atoms release from the surface of the Moon induced by meteoroid impacts

In this work, we calculate the vapour and neutral Na production rates on the Moon, as due to the impacts of meteoroids in the radius range of 10⁻⁸–0.15 m. We limit our calculations to this size range, since meteoroids with radius larger than 0.15 m have not been found to be important for the production of the exosphere in a time interval comparable with that of the observations.
We have considered a new dynamical model of the flux of meteoroids at the heliocentric distance of the Moon, regarding objects in the radius range of 10⁻²–0.15 m. Instead, the flux of smaller meteoroids (radius range 10⁻⁸–10⁻² m) has been calculated using the two distributions adopted by Cintala and Love & Brownlee.
The results of our model are that (i) the neutral Na production rate is  ∼3–4.9 × 10⁴ atoms cm⁻² s⁻¹  , slightly larger than the previous estimates  (∼2–3 × 10⁴ atoms cm⁻² s⁻¹)  , and (ii) only about 6 per cent of neutral Na is produced by the impacts of meteoroids in the size range 10⁻³–0.15 m, whereas about 94 per cent of the Na comes from the  10⁻⁵–10⁻³ m  size range.     

Radio observations of comet 9P/Tempel 1 with the Australia Telescope facilities during the Deep Impact encounter

We present radio observations of comet 9P/Tempel 1 associated with the Deep Impact spacecraft collision of 2005 July 4. Weak 18-cm OH emission was detected with the Parkes 64-m telescope, in data averaged over July 4–6, at a level of  12 ± 3 mJy km s⁻¹  , corresponding to OH production rate  2.8 × 10²⁸  molecules s⁻¹ (Despois et al. inversion model, or  1.0 × 10²⁸ s⁻¹  for the Schleicher & A'Hearn model). We did not detect the HCN 1–0 line with the Mopra 22-m telescope over the period July 2–6. The 3σ limit of 0.06 K km s⁻¹ for HCN on July 4 after the impact gives the limit to the HCN production rate of  <1.8 × 10²⁵ s⁻¹  . We did not detect the HCN 1–0 line, 6.7 GHz CH₃OH line or 3.4-mm continuum with the Australia Telescope Compact Array (ATCA) on July 4, giving further limits on any small-scale structure due to an outburst. The 3σ limit on HCN emission of 2.5 K km s⁻¹ from the ATCA around impact corresponds to limit < 4 × 10²⁹ HCN molecules released by the impact.     

Tracing gas motions in the Centaurus cluster

We apply the stochastic model of iron transport developed by Rebusco et al. to the Centaurus cluster. Using this model, we find that an effective diffusion coefficient D in the range  2 × 10²⁸–4 × 10²⁸ cm² s⁻¹  can approximately reproduce the observed abundance distribution. Reproducing the flat central profile and sharp drop around  30–70 kpc  , however, requires a diffusion coefficient that drops rapidly with radius so that   D > 4 × 10²⁸ cm² s⁻¹  only inside about  25 kpc  . Assuming that all transport is due to fully developed turbulence, which is also responsible for offsetting cooling in the cluster core, we calculate the length- and velocity-scales of energy injection. These length-scales are found to be up to a factor of ∼10 larger than expected if the turbulence is due to the inflation and rising of a bubble. We also calculate the turbulent thermal conductivity and find it is unlikely to be significant in preventing cooling.     

High-energy emission from the starburst galaxy NGC 253

Measurement sensitivity in the energetic γ-ray region has improved considerably and is about to increase further in the near future, motivating a detailed calculation of high-energy (HE; ≥100 MeV) and very high-energy (VHE; ≥100 GeV) γ-ray emission from the nearby starburst galaxy NGC 253. Adopting the convection–diffusion model for energetic electron and proton propagation, and accounting for all the relevant hadronic and leptonic processes, we determine the steady-state energy distributions of these particles by a detailed numerical treatment. The electron distribution is directly normalized by the measured synchrotron radio emission from the central starburst region; a commonly expected theoretical relation is then used to normalize the proton spectrum in this region. Doing so fully specifies the electron spectrum throughout the galactic disc and, with an assumed spatial profile of the magnetic field, the predicted radio emission from the full disc matches well the observed spectrum, confirming the validity of our treatment. The resulting radiative yields of both particles are calculated; the integrated HE and VHE fluxes from the entire disc are predicted to be   f (≥100 MeV) ≃ (1.8^+1.5_−0.8) × 10⁻⁸ cm⁻² s⁻¹  and   f (≥100 GeV) ≃ (3.6^+3.4_−1.7) × 10⁻¹² cm⁻² s⁻¹  , with a central magnetic field value   B ₀≃ 190 ± 10 μ  G. We discuss the feasibility of measuring emission at these levels with the space-borne Fermi and ground-based Cherenkov telescopes.     

Recent torque reversal of 4U 1907+09

We present timing and spectral analysis of RXTE -PCA (Proportional Counter Array) observations of the accretion powered pulsar 4U 1907+09 between 2007 June and 2008 August. 4U 1907+09 had been in a spin-down episode with a spin-down rate of  −3.54 × 10⁻¹⁴ Hz s⁻¹  before 1999. From RXTE observations after 2001 March, the source showed a ∼60 per cent decrease in spin-down magnitude, and INTEGRAL observations after 2003 March showed that source started to spin-up. We found that the source recently entered into a new spin-down episode with a spin-down rate of  −3.59 × 10⁻¹⁴ Hz s⁻¹  . This spin-down rate is pretty close to the previous long-term spin-down rate of the source measured before 1999. From the spectral analysis, we showed that hydrogen column density varies with the orbital phase.     

Deep hard X-ray source counts from a fluctuation analysis of ASCA SIS images

An analysis of the spatial fluctuations in 15 deep ASCA SIS0 images has been conducted in order to probe the 2–10 keV X-ray source counts down to a flux limit ∼ 2 × 10⁻¹⁴ erg cm⁻² s⁻¹. Special care has been taken in modelling the fluctuations in terms of the sensitivity maps of every one of the 16 regions (5.6 × 5.6 arcmin² each) into which the SIS0 has been divided, by means of ray-tracing simulations with improved optical constants in the X-ray telescope. The very extended 'sidelobes' (extending up to a couple of degrees) exhibited by these sensitivity maps make our analysis sensitive to both faint on-axis sources and brighter off-axis ones, the former being dominant. The source counts in the range (2−12) × 10⁻¹⁴ erg cm⁻² s⁻¹ are found to be close to a Euclidean form which extrapolates well to previous results from higher fluxes and are in reasonable agreement with some recent ASCA surveys. However, our results disagree with the deep survey counts by Georgantopoulos et al. The possibility that the source counts flatten to a sub-Euclidean form, as is observed at soft energies in ROSAT data, is only weakly constrained to happen at a flux < 1.8 × 10⁻¹² erg cm⁻² s⁻¹ (90 per cent confidence). Down to the sensitivity limit of our analysis, the integrated contribution of the sources the imprint of which is seen in the fluctuations amounts to ∼ 35 ± 13 per cent of the extragalactic 2–10 keV X-ray background.     

A survey of hard spectrum ROSAT sources – I. X-ray source catalogue

We present a catalogue of 147 serendipitous X-ray sources selected to have hard spectra ( α <0.5) from a survey of 188 ROSAT fields. Such sources must be the dominant contributors to the X-ray background at faint fluxes. We have used Monte Carlo simulations to verify that our technique is very efficient at selecting hard sources: the survey has 10 times as much effective area for hard sources as it has for soft sources above a 0.5–2 keV flux level of 10⁻¹⁴ erg cm⁻² s⁻¹. The distribution of best-fitting spectral slopes of the hard sources suggests that a typical ROSAT hard source in our survey has a spectral slope α ∼0. The hard sources have a steep number flux relation (d N /d S ∝ S ^{− γ} with a best-fitting value of γ =2.72±0.12) and make up about 15 per cent of all 0.5–2 keV sources with S >10⁻¹⁴ erg cm⁻² s⁻¹. If their N ( S ) continues to fainter fluxes, the hard sources will comprise ∼40 per cent of sources with 5×10⁻¹⁵< S <10⁻¹⁴ erg cm⁻² s⁻¹. The population of hard sources can therefore account for the harder average spectra of ROSAT sources with S <10⁻¹⁴ erg cm⁻² s⁻¹. They probably make a strong contribution to the X-ray background at faint fluxes and could be the solution to the X-ray background spectral paradox.     

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.     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

Detecting the first objects in the mid-infrared with the Next Generation Space Telescope

We calculate the expected mid-infrared (MIR) molecular hydrogen line emission from the first objects in the Universe. As a result of their low masses, the stellar feedback from massive stars is able to blow away their gas content and collect it into a cooling shell where H₂ rapidly forms and IR roto-vibrational (as for example the rest-frame 2.12 μm) lines carry away a large fraction (up to 10 per cent) of the explosion energy. The fluxes from these sources are in the range 10⁻²¹–10⁻¹⁷ erg s⁻¹ cm⁻² . The highest number counts are expected in the 20-μm band, where about 10⁵ sources deg⁻² are predicted at the limiting flux of 3×10⁻¹⁸ erg s⁻¹ cm⁻². Among the planned observational facilities, we find that the best detection perspectives are offered by the Next Generation Space Telescope ( NGST ), which should be able to reveal about 200 first objects in one hour observation time at its limiting flux in the above band. Therefore, mid-IR instruments appear to represent perfect tools to trace star formation and stellar feedback in the high ( z ≳5) redshift Universe.     

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.     

The accretion powered spin-up of GRO J1750−27

The timing properties of the 4.45 s pulsar in the Be X-ray binary system GRO J1750−27 are examined using hard X-ray data from INTEGRAL and Swift during a type II outburst observed during 2008. The orbital parameters of the system are measured and agree well with those found during the last known outburst of the system in 1995. Correcting the effects of the Doppler shifting of the period, due to the orbital motion of the pulsar, leads to the detection of an intrinsic spin-up that is well described by a simple model including     and     terms of  −7.5 × 10⁻¹⁰ s s⁻¹  and  1 × 10⁻¹⁶ s s⁻²  , respectively. The model is then used to compare the time-resolved variation of the X-ray flux and intrinsic spin-up against the accretion torque model of Ghosh & Lamb; this finds that GRO J1750−27 is likely located 12–22 kpc distant and that the surface magnetic field of the neutron star is  ∼2 × 10¹²  G. The shape of the pulse and the pulsed fraction shows different behaviour above and below 20 keV, indicating that the observed pulsations are the convolution of many complex components.     

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.     

Unresolved emission and ionized gas in the bulge of M31

We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 10²⁷ erg s⁻¹ L⁻¹_⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 10⁶ M_⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L _X/SFR  (star formation rate) ratio equals  ∼9 × 10³⁸ (erg s⁻¹)(M_⊙ yr⁻¹)⁻¹  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.     

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.     

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.     

The K-band luminosity function at z = 1: a powerful constraint on galaxy formation theory

We report the first detection of an inverse Compton X-ray emission, spatially correlated with a very steep spectrum radio source (VSSRS), 0038-096, without any detected optical counterpart, in cluster Abell 85. The ROSAT PSPC data and its multiscale wavelet analysis reveal a large-scale (linear diameter of the order of 500 h ⁻¹₅₀ kpc), diffuse X-ray component, in addition to the thermal bremsstrahlung, overlapping an equally large-scale VSSRS. The primeval 3 K background photons, scattering off the relativistic electrons, can produce the X-rays at the detected level. The inverse Compton flux is estimated to be (6.5 ± 0.5) × 10⁻¹³ erg s⁻¹ cm⁻² in the 0.5–2.4 keV X-ray band. A new 327-MHz radio map is presented for the cluster field. The synchrotron emission flux is estimated to be (6.6 ± 0.90) × 10⁻¹⁴ erg s⁻¹ cm⁻² in the 10–100 MHz radio band. The positive detection of both radio and X-ray emission from a common ensemble of relativistic electrons leads to an estimate of (0.95 ± 0.10) × 10⁻⁶ G for the cluster-scale magnetic field strength. The estimated field is free of the 'equipartition' conjecture, the distance, and the emission volume. Further, the radiative fluxes and the estimated magnetic field imply the presence of 'relic' (radiative lifetime ≳ 10⁹ yr) relativistic electrons with Lorentz factors γ ≈ 700–1700; this would be a significant source of radio emission in the hitherto unexplored frequency range ν ≈ 2–10 MHz.     

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⁻¹  .