A Herbig–Haro object associated with GGD30 and its exciting source

GGD30 has been suggested to be either a small reflection nebulosity or a Herbig–Haro (HH) object formed in the outflow from a nearby obscured star. Observations to date have not been able to distinguish between these two scenarios. In addition, there are conflicting proposals for the location of the exciting source for GGD30. To resolve these questions, we have carried out optical spectroscopy and near-infrared ( J , K and 3.6-μm) imaging of GGD30. Taken together, these observations reveal that the bright optical knot in GGD30 must be a HH object, excited by the outflow from an optically obscured pre-main-sequence (PMS) star located ∼3 arcsec to the southwest. Based on mid-infrared fluxes from the Mid-course Space Experiment ( MSX ) satellite, we estimate the luminosity of this PMS star to be  ∼12.5 L_⊙  which suggests it is an intermediate-mass object rather than low-mass as previously proposed. The optical spectroscopy indicates projected velocities of  ∼−270 km s⁻¹  associated with the HH object. The fact that these velocities are blueshifted and relatively high compared to the velocities typical of HH flows suggests that the outflow from the PMS star must be almost aligned with the line of sight. There is an additional low-velocity  (∼−70 km s⁻¹) Hα  component but its origin is not clear.     

The effect of non-gravitational forces on the median inclination of short-period comets

The numbered Jupiter family comets (orbital periods   P < 20 yr  ) have a median orbital inclination of about     . In this paper, we integrate the orbits of these comets into the future, under the influence of both typical non-gravitational forces and planetary perturbation, using a Bulirsch–Stoer integrator. In the case where non-gravitational forces were not acting, the median inclination of those comets that remained on   P < 20 yr  orbits increased at the rate of  (1.92 ± 0.12) × 10⁻³ deg yr⁻¹  for the first 3600 yr of the integration. During this time the population of the original family decreases, such that the half-life is about 13 200 ± 800 yr. The introduction of non-gravitational forces slows down the rate of increase in inclination to a value of around  (1.23 ± 0.16) × 10⁻³ deg yr⁻¹  . This rate of increase in inclination was found to be only weakly dependent on the non-gravitational parameters used during the integration. After a few thousand years, the rate of change in inclination decreases, and after 20 000 yr the inclinations of those initial Jupiter family members that still have orbits with   P < 20 yr  become constant at about     , independent of whether non-gravitational forces are acting or not. The presently known Jupiter family of comets is losing members at the rate of one in every 67 yr. To maintain the family in equilibrium, Jupiter has to capture comets at a similar rate, and these captured comets have to be of low inclination to compensate for the pumping up of inclinations by gravitational perturbation.     

Circumstellar structure of RU Lupi down to au scales

We have used the technique of spectro-astrometry to study the milliarcsecond scale structure of the emission lines in the T Tauri star RU Lupi. The wings of the H α emission are found to be displaced from the star towards the south-west (blue wing) and north-east (red wing) with angular scales of 20–30 mas. This structure is consistent with a bipolar outflow from the star. From a study of the variability of the intensity and position spectra, we argue that a combination of magnetically driven bipolar outflow and accreting gas contributes to the H α emission. On the other hand, the [O  i ] and [S  ii ] emission are displaced from the star to the south-west but at much larger distances than the H α , hundreds of milliarcseconds for the high-velocity component (HVC) and down to 30 mas for the low-velocity components (LVCs). The presence of both redshifted and blueshifted outflows in H α but only a blueshifted outflow in the forbidden lines can be explained if the disc obscures the redshifted forbidden line outflow, but a disc gap with outer radius 3–4 au allows the redshifted H α to be seen. This gap could be induced by an unseen companion.     

The emission and kinematic structures of the bipolar planetary nebula M 1-8

We have undertaken echelle spectroscopy and narrow-band line imaging of the bipolar planetary nebula M 1-8. This has permitted us to map the outflow in [N  ii ]λλ 6548+6583 Å, Hα, and in the v = 1–0 S(1) transition of H₂ at λ 2.122 μm. It has also permitted us to acquire high-resolution spectra for [N  ii ]λ 6583 Å, Hα and He  ii λ 6560 Å. Our observations support the results of a previous 2MASS analysis by two of the authors (J. P. Phillips and G. Ramos-Larios), and confirm that there is strong H₂ emission outside of the ionized zone, as well as along the major axis of the outflow. Finally, we have investigated the spatial structure of the outflow in low and high excitation lines, and noted evidence for strong ionization stratification within the envelope of the source. We also note that major axis spectra show asymmetries attributable to outflow along the lobes, oriented at an angle i ∼ 35°–40° to the line of sight. Asymmetries along the minor axis, by contrast, appear to be associated with the central collimating disc, and may be interpretable in terms of asymmetries in disc structure, or rotation at an angular velocity of Ω∼ 1.4 10⁻¹² rad s⁻¹. If the disc arises due to common-envelope evolution, then it seems that angular momentum constraints must be relatively tight, and can only be satisfied given fairly extreme physical assumptions (such as low disc mass, high primary star mass, a low distance to the source and so forth).     

Maximum-likelihood method for estimating the mass and period distributions of extrasolar planets

We investigate the distribution of mass M and orbital period P of extrasolar planets, taking account of selection effects caused by the limited velocity precision and duration of existing surveys. We fit the data on 72 planets to a power-law distribution of the form  d n = CM ^−α P ^−β(d M / M )(d P / P )  , and find  α= 0.11 ± 0.10  ,  β=−0.27 ± 0.06  for   M ≲ 10 M _J  , where   M _J  is the mass of Jupiter. The correlation coefficient between these two exponents is −0.31, indicating that uncertainties in the two distributions are coupled. We estimate that 4 per cent of solar-type stars have companions in the range  1 M _J < M < 10 M _J  ,  2_d < P < 10 yr  .     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Inflow and outflow from the accretion disc of the microquasar SS 433: UKIRT spectroscopy

A succession of near-infrared (near-IR) spectroscopic observations, taken nightly throughout an entire cycle of SS 433's orbit, reveal (i) the persistent signature of SS 433's accretion disc, having a rotation speed of  ∼500 km s⁻¹  , (ii) the presence of circumbinary disc recently discovered at optical wavelengths by Blundell, Bowler & Schmidtobreick (2008) and (iii) a much faster outflow than has previously been measured for the disc wind, with a terminal velocity of  ∼1500 km s⁻¹  . The increased wind terminal velocity results in a mass-loss rate of  ∼10⁻⁴ M_⊙ yr⁻¹  . These, together with the newly (upwardly) determined masses for the components of the SS 433 system, result in an accurate diagnosis of the extent to which SS 433 has super-Eddington flows. Our observations imply that the size of the companion star is comparable with the semiminor axis of the orbit which is given by     , where e is the eccentricity. Our relatively spectral resolution at these near-IR wavelengths has enabled us to deconstruct the different components that comprise the Brackett-γ (Brγ) line in this binary system, and their physical origins. With this line being dominated throughout our series of observations by the disc wind, and the accretion disc itself being only a minority (∼15 per cent) contribution, we caution against use of the unresolved Brγ line intensity as an 'accretion signature' in X-ray binaries or microquasars in any quantitative way.     

Peculiar nature of hard X-ray eclipse in SS433 from INTEGRAL observations

The analysis of hard X-ray INTEGRAL observations (2003–2008) of superaccreting Galactic microquasar SS433 at precessional phases of the source with the maximum disc opening angle is carried out. It is found that the shape and width of the primary X-ray eclipse are strongly variable, suggesting additional absorption in dense stellar wind and gas outflows from the optical A7I component and the wind–wind collision region. The independence of the observed hard X-ray spectrum on the accretion disc precessional phase suggests that hard X-ray emission (20–100 keV) is formed in an extended, hot, quasi-isothermal corona, probably heated by interaction of relativistic jet with inhomogeneous wind outflow from the precessing supercritical accretion disc. A joint modelling of X-ray eclipsing and precessional hard X-ray variability of SS433 revealed by INTEGRAL by a geometrical model suggests the binary mass ratio   q = m_x / m _v ≃  0.25–0.5. The absolute minimum of joint orbital and precessional  χ²  residuals is reached at   q ≃ 0.3  . The found binary mass ratio range allows us to explain the substantial precessional variability of the minimum brightness at the middle of the primary optical eclipse. For the mass function of the optical star   f _v = 0.268 M_⊙  as derived from Hillwig & Gies data, the obtained value of   q ≃ 0.3  yields the masses of the components   m_x ≃ 5.3 M_⊙, m _v ≃ 17.7 M_⊙  , confirming the black hole nature of the compact object in SS433.     

Detection of superhumps in XTE J1118+480 approaching quiescence

We present the results of our monitoring of the halo black hole soft X-ray transient (SXT) XTE J1118+480 during its decline to quiescence. The system has decayed 0.5 mag from 2000 December to its present near-quiescent level at   R ≃18.65  (2001 June). The ellipsoidal light curve is distorted by an additional modulation that we interpret as a superhump of   P _sh=0.17049(1) d  i.e. 0.3 per cent longer than the orbital period. This implies a disc precession period   P _prec∼52 d  . After correcting the average phase-folded light curve for veiling, the amplitude difference between the minima suggests that the binary inclination angle lies in the range   i =71–82°  . However, we urge caution in the interpretation of these values because of residual systematic contamination of the ellipsoidal light curve by the complex form of the superhump modulation. The orbital-mean H α profiles exhibit clear velocity variations with ∼500 km s⁻¹ amplitude. We interpret this as the first spectroscopic evidence of an eccentric precessing disc.     

Quark-nova remnants: application to RRATs

This is the third paper of a series of papers where we explore the evolution of iron-rich ejecta from quark-novae. In the first paper, we explored the case where a quark-nova ejecta forms a degenerate shell, supported by the star's magnetic field, with applications to SGRs. In the second paper, we considered the case where the ejecta would have sufficient angular momentum to form a degenerate Keplerian torus and applied such a system to two AXPs, namely 1E2259+586 and 4U0142+615. Here, we explore the late evolution of the degenerate torus and find that it can remain unchanged for  ∼10⁶ yr  before it becomes non-degenerate. This transition from a degenerate torus (accretion-dominated) to a non-degenerate disc (no accretion) occurs about 10⁶ yr following the quark-nova, and exhibits features that are reminiscent of observed properties of Rotation RAdio Transients (RRATs). Using this model, we can account for the duration of both the radio bursts and the quiet phase, as well as the observed radio flux from RRATs. We discuss a connection between XDINs and RRATs and argue that some XDINs may be 'dead RRATs' that have already consumed their non-degenerate disc.     

The nature of TW Pictoris

The results of X-ray and optical observations of the candidate intermediate polar TW Pic are presented in an attempt to understand its nature. We find no sign of the previously proposed ∼2 h white-dwarf spin period and ∼6 h orbital period of TW Pic in its X-ray light curve. There is therefore no convincing evidence in support of its previous classification. The lack of X-ray pulsation could be the result of a low inclination angle, but in that case there would be no reason why an optical pulsation should have been seen previously. Negative results from polarimetry also preclude TW Pic from being a polar. One possibility may be that the shorter of the two periods is in fact the orbital period, whilst the longer one is a harmonic of a disc precession period. Alternatively, both the high accretion rate and period structure of TW Pic indicate that it may be a system that displays persistent negative superhumps. In this case the true orbital period of the binary may be around 6.36 h and the shorter of the two previously identified periods, 1.996 h, represents the (shifted) second harmonic of a negative superhump period of 6.06 h. Under this interpretation, it would be the longest period system to display such a phenomenon. Finally there is also evidence that TW Pic may be a VY Scl star, in which case it would be the longest period member of that subclass too.     

Timing of the accreting millisecond pulsar XTE J1814−338

We present a precise timing analysis of the accreting millisecond pulsar XTE J1814−338 during its 2003 outburst, observed by RXTE . A full orbital solution is given for the first time; Doppler effects induced by the motion of the source in the binary system were corrected, leading to a refined estimate of the orbital period,   P _orb= 15 388.7229(2)  s, and of the projected semimajor axis,   a sin  i / c = 0.390633(9)  light-second. We could then investigate the spin behaviour of the accreting compact object during the outburst. We report here a refined value of the spin frequency  (ν= 314.356 108 79(1) Hz)  and the first estimate of the spin frequency derivative of this source while accreting     . This spin-down behaviour arises when both the fundamental frequency and the second harmonic are taken into consideration. We discuss this in the context of the interaction between the disc and the quickly rotating magnetosphere, at accretion rates sufficiently low to allow a threading of the accretion disc in regions where the Keplerian velocity is slower than the magnetosphere velocity. We also present indications of a jitter of the pulse phases around the mean trend, which we argue results from movements of the accreting hotspots in response to variations of the accretion rate.     

Mass and spin co-evolution during the alignment of a black hole in a warped accretion disc

In this paper, we explore the gravitomagnetic interaction of a black hole (BH) with a misaligned accretion disc to study BH spin precession and alignment jointly with BH mass M _BH and spin parameter a evolution, under the assumption that the disc is continually fed, in its outer region, by matter with angular momentum fixed on a given direction     . We develop an iterative scheme based on the adiabatic approximation to study the BH–disc co-evolution: in this approach, the accretion disc transits through a sequence of quasi-steady warped states (Bardeen–Petterson effect) and interacts with the BH until the spin   J _BH  aligns with     . For a BH aligning with a corotating disc, the fractional increase in mass is typically less than a few per cent, while the spin modulus can increase up to a few tens of per cent. The alignment time-scale     is of  ∼10⁵–10⁶ yr  for a maximally rotating BH accreting at the Eddington rate. BH–disc alignment from an initially counter-rotating disc tends to be more efficient compared to the specular corotating case due to the asymmetry seeded in the Kerr metric: counter-rotating matter carries a larger and opposite angular momentum when crossing the innermost stable orbit, so that the spin modulus decreases faster and so the relative inclination angle.     

The excitation and kinematical properties of H2 and [Fe ii] in the HH 46/47 bipolar outflow

Long-slit spectra of the molecular outflow Herbig–Haro (HH) 46/47 have been taken in the J and K near-infrared bands. The observed H₂ line emission confirms the existence of a bright and extended redshifted counter-jet outflow south-west of HH 46. In contrast with the optical appearance of this object, we show that this outflow seems to be composed of two different emission regions characterized by distinct heliocentric velocities. This implies an acceleration of the counter-jet.
The observed [Fe  ii ] emission suggests an average extinction of 7–9 visual magnitudes for the region associated with the counter-jet.
Through position–velocity diagrams, we show the existence of different morphologies for the H₂ and [Fe  ii ] emission regions in the northern part of the HH 46/47 outflow. We have detected for the first time high-velocity (−250 km s⁻¹) [Fe  ii ] emission in the region bridging HH 46 to HH 47A. The two strong peaks detected can be identified with the optical positions B8 and HH 47B.
The H₂ excitation diagrams for the counter-jet shock suggest an excitation temperature for the gas of T _ex≈2600 K . The lack of emission from the higher energy H₂ lines, such as the 4–3 S(3) transition, suggests a thermal excitation scenario for the origin of the observed emission. Comparison of the H₂ line ratios with various shock models yielded useful constraints about the geometry and type of these shocks. Planar shocks can be ruled out whereas curved or bow shocks (both J- and C-type) can be parametrized to fit our data.     

A shallow though extensive H2 2.122-μm imaging survey of Taurus–Auriga–Perseus – I. NGC 1333, L1455, L1448 and B1

We discuss wide-field near-infrared (near-IR) imaging of the NGC 1333, L1448, L1455 and B1 star-forming regions in Perseus. The observations have been extracted from a much larger narrow-band imaging survey of the Taurus–Auriga–Perseus complex. These H₂ 2.122-μm observations are complemented by broad-band K imaging, mid-IR imaging and photometry from the Spitzer Space Telescope , and published submillimetre CO   J = 3–2  maps of high-velocity molecular outflows. We detect and label 85 H₂ features and associate these with 26 molecular outflows. Three are parsec-scale flows, with a mean flow lobe length exceeding 11.5 arcmin. 37 (44 per cent) of the detected H₂ features are associated with a known Herbig–Haro object, while 72 (46 per cent) of catalogued HH objects are detected in H₂ emission. Embedded Spitzer sources are identified for all but two of the 26 molecular outflows. These candidate outflow sources all have high near-to-mid-IR spectral indices (mean value of  α∼ 1.4  ) as well as red IRAC 3.6–4.5 μm and IRAC/MIPS 4.5–24.0 μm colours: 80 per cent have [3.6]–[4.5] > 1.0 and [4.5]–[24] > 1.5. These criteria – high α and red [4.5]–[24] and [3.6]–[4.5] colours – are powerful discriminants when searching for molecular outflow sources. However, we find no correlation between α and flow length or opening angle, and the outflows appear randomly orientated in each region. The more massive clouds are associated with a greater number of outflows, which suggests that the star formation efficiency is roughly the same in each region.     

Probing the microarcsecond structure of the Circinus megamasers through diffractive interstellar scintillation

We have used the Long Baseline Array to observe the Circinus Galaxy H₂O megamasers on two epochs separated by 49 d. The masers were observed to be highly variable, consistent with previous studies. We observe rapid variations in the linewidth and central velocity of the maser's spectral profile. The characteristics of these spectral variations suggest that the masers have an anisotropic structure, with axial ratios ≲5, and display velocity gradients in their microarcsecond structure. We also find a strong anticorrelation between the linewidth and the modulation index of the maser features and a relationship between the modulation index and characteristic time-scale of the variations. These various phenomena can all be simply explained through a model where the rapid variations in the masers are produced by quenched diffractive scintillation in the presence of modest velocity gradients in the maser emission.
The very long baseline interferometry (VLBI) images are consistent with previous high-resolution observations and show maser emission arising both in a circumnuclear accretion disc, and in a low-velocity outflow. We find no systematic differences in the variability properties of the disc and outflow masers. The VLBI images reveal that there is considerable spectral blending in single-dish spectra of the Circinus megamasers.
In addition to the VLBI imaging we also used these observations to look for a time delay in the intensity variations observed at the different telescopes in the array, as expected for scintillation-induced variability. No time delay was detected, which allows us to place a lower limit on the velocity of the interstellar medium during these observations 22 km s⁻¹ for the second epoch. We conclude that any future attempts to measure a time delay in this source should utilize intercontinental baselines to increase the chances of success.     

Emission-line oscillations in the dwarf nova V2051 Ophiuchi

We have detected coherent oscillations, at multiple frequencies, in the line and continuum emission of the eclipsing dwarf nova V2051 Ophiuchi using the 10-m Keck II telescope. Our own novel data acquisition system allowed us to obtain very fast spectroscopy using a continuous readout of the CCD on the LRIS spectrograph. This is the first time that dwarf nova oscillations have been detected and resolved in the emission lines. The accretion disc is highly asymmetric with a stronger contribution from the blueshifted side of the disc during our observations. The disc extends from close to the white dwarf out to the outer regions of the primary Roche lobe.
Continuum oscillations at 56.12 s and its first harmonic at 28.06 s are most likely to originate on the surface of a spinning white dwarf with the fundamental period corresponding to the spin period. Balmer and helium emission lines oscillate with a period of 29.77 s at a mean amplitude of 1.9 per cent. The line kinematics and the eclipse constraints indicate an origin in the accretion disc at a radius of 12±2 R _wd. The amplitude of the emission-line oscillation modulates (0–4 per cent) at a period of 488 s, corresponding to the Kepler period at R =12 R _wd. This modulation is caused by the beating between the white dwarf spin and the orbital motion in the disc.
The observed emission-line oscillations cannot be explained by a truncated disc as in the intermediate polars. The observations suggest a non-axisymmetric bulge in the disc, orbiting at 12 R _wd, is required. The close correspondence between the location of the oscillations and the circularization radius of the system suggests that stream overflow effects may be of relevance.     

Disc precession in the M31 dipping X-ray binary Bo 158?

We present results from three XMM–Newton observations of the M31 low mass X-ray binary (LMXB) XMMU J004314.4+410726.3 (Bo 158), spaced over 3 d in 2004 July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-h period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anticorrelated with source intensity. In light of the 2004 XMM–Newton observations of Bo 158, we suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio ≲0.3 (period ≲4 h), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of  ∼11 P _orb  , and prograde disc precession on a period of  29 ± 1 P _orb  . This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on an  81 ± 3 h  cycle.     

A new determination of the orbit and masses of the Be binary system δ Scorpii

The binary star δ Sco (HD143275) underwent remarkable brightening in the visible in 2000, and continues to be irregularly variable. The system was observed with the Sydney University Stellar Interferometer (SUSI) in 1999, 2000, 2001, 2006 and 2007. The 1999 observations were consistent with predictions based on the previously published orbital elements. The subsequent observations can only be explained by assuming that an optically bright emission region with an angular size of  ≳2 ± 1 mas  formed around the primary in 2000. By 2006/2007 the size of this region grew to an estimated ≳4 mas.
We have determined a consistent set of orbital elements by simultaneously fitting all the published interferometric and spectroscopic data as well as the SUSI data reported here. The resulting elements and the brightness ratio for the system measured prior to the outburst in 2000 have been used to estimate the masses of the components. We find   M_A = 15 ± 7 M_⊙  and   M_B = 8.0 ± 3.6 M_⊙  . The dynamical parallax is estimated to be  7.03 ± 0.15 mas  , which is in good agreement with the revised Hipparcos parallax.     

Cyclical period changes in the dwarf novae V2051 Oph and V4140 Sgr

We report on the identification of cyclical changes in the orbital period of the eclipsing dwarf novae V2051 Ophiuchi and V4140 Sagittarii. We used sets of white dwarf mid-eclipse timings to construct observed-minus-calculated diagrams covering, respectively, 25 and 16 yr of observations. The V2051 Oph data present cyclical variations that can be fitted by a linear plus sinusoidal function with period of  22 ± 2 yr  and amplitude of  17 ± 3 s  . The statistical significance of this period by an F-test is larger than 99.9 per cent. The V4140 Sgr data present cyclical variations of similar amplitude and period of  6.9 ± 0.3 yr  which are statistically significant at the 99.7 per cent level. We derive upper limits for secular period changes of     and     for V2051 Oph and V4140 Sgr, respectively.
We have combined our results with those in the literature to construct a diagram of the amplitude versus period of the modulation for a sample of 11 eclipsing cataclysmic variables (CVs). If the cyclical period changes are the consequence of a solar-type magnetic activity cycle in the secondary star, then magnetic activity is a widespread phenomenon in CVs, being equally common among long- and short-period systems. This gives independent evidence that the magnetic field (and activity) of the secondary stars of CVs do not disappear when they become fully convective. We also find that the fractional cycle period changes of the short-period CVs are systematically smaller than those of the long-period CVs.