The supergiant B[e] star LHA 115-S 18 - binary and/or luminous blue variable?

Context. The mechanism by which supergiant (sg)B[e] stars support cool, dense dusty discs/tori and their physical relationship with other evolved, massive stars such as luminous blue variables is uncertain. Aims. In order to investigate both issues we have analysed the long term behaviour of the canonical sgB[e] star LHA 115-S 18. Methods. We employed the OGLE II-IV lightcurve to search for (a-)periodic variability and supplemented these data with new and historic spectroscopy. Results. In contrast to historical expectations for sgB[e] stars, S18 is both photometrically and spectroscopically highly variable. The lightcurve is characterised by rapid aperiodic ` aring' throughout the 16 years of observations. Changes in the high excitation emission line component of the spectrum imply evolution in the stellar temperature - as expected for luminous blue variables - although somewhat surprisingly, spectroscopic and photometric variability appears not to be correlated. Characterised by emission in low excitation metallic species, the cool circumstellar torus appears largely unaffected by this behaviour. Finally, in conjunction with intense, highly variable He ii emission, X-ray emission implies the presence of an unseen binary companion. Conclusions. S18 provides observational support for the putative physical association of (a subset of) sgB[e] stars and luminous blue variables. Given the nature of the circumstellar environment of S18 and that luminous blue variables have been suggested as SN progenitors, it is tempting to draw a parallel to the progenitors of SN1987A and SN2009ip. Moreover the likely binary nature of S18 strengthens the possibility that the dusty discs/tori that characterise sgB[e] stars are the result of binary-driven mass-loss; consequently such stars may provide a window on the short lived phase of mass-transfer in massive compact binaries.

The circumstellar environment and evolutionary status of the supergiant B[e] star Wd1-9

Context. Historically, supergiant (sg)B[e] stars have been difficult to include in theoretical schemes for the evolution of massive OB stars. Aims. The location of Wd1-9 within the coeval starburst cluster Westerlund 1 means that it may be placed into a proper evolutionary context and we therefore aim to utilise a comprehensive multiwavelength dataset to determine its physical properties and consequently its relation to other sgB[e] stars and the global population of massive evolved stars within Wd1. Methods. Multi-epoch R- and I-band VLT/UVES and VLT/FORS2 spectra are used to constrain the properties of the circumstellar gas, while an ISO-SWS spectrum covering 2.45−45μm is used to investigate the distribution, geometry and composition of the dust via a semi-analytic irradiated disk model. Radio emission enables a long term mass-loss history to be determined, while X-ray observations reveal the physical nature of high energy processes within the system. Results. Wd1-9 exhibits the rich optical emission line spectrum that is characteristic of sgB[e] stars. Likewise its mid-IR spectrum resembles those of the LMC sgB[e] stars R66 and 126, revealing the presence of equatorially concentrated silicate dust, with a mass of ~10−4M⊙. Extreme historical and ongoing mass loss (≳ 10−4M⊙yr−1) is inferred from the radio observations. The X-ray properties of Wd1-9 imply the presence of high temperature plasma within the system and are directly comparable to a number of confirmed short-period colliding wind binaries within Wd1. Conclusions. The most complete explanation for the observational properties of Wd1-9 is that it is a massive interacting binary currently undergoing, or recently exited from, rapid Roche-lobe overflow, supporting the hypothesis that binarity mediates the formation of (a subset of) sgB[e] stars. The mass loss rate of Wd1-9 is consistent with such an assertion, while viable progenitor and descendent systems are present within Wd1 and comparable sgB[e] binaries have been identified in the Galaxy. Moreover, the rarity of sgB[e] stars - only two examples are identified from a census of ~ 68 young massive Galactic clusters and associations containing ~ 600 post-Main Sequence stars - is explicable given the rapidity (~ 104yr) expected for this phase of massive binary evolution.

Accreting magnetars: a new type of high-mass X-ray binaries?

The discovery of very slow pulsations (Pspin =5560 s) has solved the long-standing question of the nature of the compact object in the high-mass X-ray binary 4U 2206+54 but has posed new ones. According to spin evolutionary models in close binary systems, such slow pulsations require a neutron star magnetic field strength larger than the quantum critical value of 4.4 × 1013 G, suggesting the presence of a magnetar. We present the first XMM–Newton observations of 4U 2206+54 and investigate its spin evolution. We find that the observed spin-down rate agrees with the magnetar scenario. We analyse Integral Spacecraft Gamma-Ray Imager (ISGRI)/INTErnational Gamma-RAy Laboratory (INTEGRAL) observations of 4U 2206+54 to search for the previously suggested cyclotron resonance scattering feature at ∼30 keV. We do not find a clear indication of the presence of the line, although certain spectra display shallow dips, not always at 30 keV. The association of these dips with a cyclotron line is very dubious because of its apparent transient nature. We also investigate the energy spectrum of 4U 2206+54 in the energy range 0.3–10 keV with unprecedented detail and report for the first time the detection of very weak 6.5 keV fluorescence iron lines. The photoelectric absorption is consistent with the interstellar value, indicating very small amount of local matter, which would explain the weakness of the florescence lines. The lack of matter locally to the source may be the consequence of the relatively large orbital separation of the two components of the binary. The wind would be too tenuous in the vicinity of the neutron star.

IRAS 18357-0604 – an analogue of the galactic yellow hypergiant IRC +10420?

Context. Yellow hypergiants represent a short-lived evolutionary episode experienced by massive stars as they transit to and from a red supergiant phase. As such, their properties provide a critical test of stellar evolutionary theory, while recent observations unexpectedly suggest that a subset may explode as Type II supernovae. Aims. The galactic yellow hypergiant IRC +10420 is a cornerstone system for understanding this phase since it is the strongest post-RSG candidate known, has demonstrated real-time evolution across the Hertzsprung-Russell diagram and been subject to extensive mass loss. In this paper we report on the discovery of a twin of IRC +10420 - IRAS 18357-0604. Methods. Optical and near-IR spectroscopy are used to investigate the physical properties of IRAS 18357-0604 and also provide an estimate of its systemic velocity, while near- to mid-IR photometry probes the nature of its circumstellar environment. Results. These observations reveal pronounced spectral similarities between IRAS 18357-0604 and IRC +10420, suggesting comparable temperatures and wind geometries. IR photometric data reveals a similarly dusty circumstellar environment, although historical mass loss appears to have been heavier in IRC +10420. The systemic velocity implies a distance compatible with the red supergiant-dominated complex at the base of the Scutum Crux arm; the resultant luminosity determination is consistent with a physical association but suggests a lower initial mass than inferred for IRC +10420 (≲20 M⊙ versus ~40 M⊙). Evolutionary predictions for the physical properties of supernova progenitors derived from ~18–20 M⊙ stars – or ~12–15 M⊙ stars that have experienced enhanced mass loss as red supergiants – compare favourably with those of IRAS 18357-0604, which in turn appears to be similar to the the progenitor of SN2011dh; it may therefore provide an important insight into the nature of the apparently H-depleted yellow hypergiant progenitors of some Type IIb SNe.

An XMM-Newton view of FeKα in high-mass X-ray binaries

We present a comprehensive analysis of the whole sample of available XMM-Newton observations of high-mass X-ray binaries (HMXBs) until August 2013, focusing on the FeKα emission line. This line is key to better understanding the physical properties of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We collected observations from 46 HMXBs and detected FeKα in 21 of them. We used the standard classification of HMXBs to divide the sample into different groups. We find that (1) different classes of HMXBs display different qualitative behaviours in the FeKα spectral region. This is visible especially in SGXBs (showing ubiquitous Fe fluorescence but not recombination Fe lines) and in γ Cass analogues (showing both fluorescent and recombination Fe lines). (2) FeKα is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncertainties, this value is compatible with ionization states that are lower than Fe xviii. (3) The flux of the continuum is well correlated with the flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding the X-ray source. (4) We observe an inverse correlation between the X-ray luminosity and the equivalent width of FeKα (EW). This phenomenon is known as the X-ray Baldwin effect. (5) FeKα is narrow (σline< 0.15 keV), reflecting that the reprocessing material does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line blending, Compton scattering, and Doppler shifts (with velocities of the reprocessing material V ~ 1000 km s-1). (6) The equivalent hydrogen column (NH) directly correlates to the EW of FeKα, displaying clear similarities to numerical simulations. It highlights the strong link between the absorbing and the fluorescent matter. (7) The observed NH in supergiant X-ray binaries (SGXBs) is in general higher than in supergiant fast X-ray transients (SFXTs). We suggest two possible explanations: different orbital configurations or a different interaction compact object – wind. (8) Finally, we analysed the sources IGR J16320-4751 and 4U 1700-37 in more detail, covering several orbital phases. The observed variation in NH between phases is compatible with the absorption produced by the wind of their optical companions. The results clearly point to a very important contribution of the donor’s wind in the FeKα emission and the absorption when the donor is a supergiant massive star.

Orbital phase-resolved spectroscopy of 4U 1538−52 with MAXI

Context. 4U 1538−52, an absorbed high mass X-ray binary with an orbital period of ~3.73 days, shows moderate orbital intensity modulations with a low level of counts during the eclipse. Several models have been proposed to explain the accretion at different orbital phases by a spherically symmetric stellar wind from the companion. Aims. The aim of this work is to study both the light curve and orbital phase spectroscopy of this source in the long term. In particular, we study the folded light curve and the changes in the spectral parameters with orbital phase to analyse the stellar wind of QV Nor, the mass donor of this binary system. Methods. We used all the observations made from the Gas Slit Camera on board MAXI of 4U 1538−52 covering many orbits continuously. We obtained the good interval times for all orbital phase ranges, which were the input for extracting our data. We estimated the orbital period of the system and then folded the light curves, and we fitted the X-ray spectra with the same model for every orbital phase spectrum. We also extracted the averaged spectrum of all the MAXI data available. Results. The MAXI spectra in the 2–20 keV energy range were fitted with an absorbed Comptonisation of cool photons on hot electrons. We found a strong orbital dependence of the absorption column density but neither the fluorescence iron emission line nor low energy excess were needed to fit the MAXI spectra. The variation in the spectral parameters over the binary orbit were used to examine the mode of accretion onto the neutron star in 4U 1538−52. We deduce a best value of Ṁ/v∞ = 0.65 × 10-9M⊙ yr-1/ (km s-1) for QV Nor.

The HIPASS catalogue - II. Completeness reliability and parameter accuracy

The H I Parkes All Sky Survey (HIPASS) is a blind extragalactic H I 21-cm emission-line survey covering the whole southern sky from declination -90degrees to +25degrees. The HIPASS catalogue (HICAT), containing 4315 H I-selected galaxies from the region south of declination +2degrees, is presented in Meyer et al. (Paper I). This paper describes in detail the completeness and reliability of HICAT, which are calculated from the recovery rate of synthetic sources and follow-up observations, respectively. HICAT is found to be 99 per cent complete at a peak flux of 84 mJy and an integrated flux of 9.4 Jy km. s(-1). The overall reliability is 95 per cent, but rises to 99 per cent for sources with peak fluxes >58 mJy or integrated flux >8.2 Jy km s(-1). Expressions are derived for the uncertainties on the most important HICAT parameters: peak flux, integrated flux, velocity width and recessional velocity. The errors on HICAT parameters are dominated by the noise in the HIPASS data, rather than by the parametrization procedure.

The 2dF-SDSS LRG and QSO survey: the star formation histories of luminous red galaxies

We present a detailed investigation into the recent star formation histories of 5697 luminous red galaxies (LRGs) based on the H delta (4101 angstrom), and [O II] (3727 angstrom) lines and the D4000 index. LRGs are luminous (L > 3L*) galaxies which have been selected to have photometric properties consistent with an old, passively evolving stellar population. For this study, we utilize LRGs from the recently completed 2dF-SDSS LRG and QSO Survey (2SLAQ). Equivalent widths of the H delta and [O II] lines are measured and used to define three spectral types, those with only strong H delta absorption (k+a), those with strong [O II] in emission (em) and those with both (em+a). All other LRGs are considered to have passive star formation histories. The vast majority of LRGs are found to be passive (similar to 80 per cent); however, significant numbers of k+a (2.7 per cent), em+a (1.2 per cent) and em LRGs (8.6 per cent) are identified. An investigation into the redshift dependence of the fractions is also performed. A sample of SDSS MAIN galaxies with colours and luminosities consistent with the 2SLAQ LRGs is selected to provide a low-redshift comparison. While the em and em+a fractions are consistent with the low-redshift SDSS sample, the fraction of k+a LRGs is found to increase significantly with redshift. This result is interpreted as an indication of an increasing amount of recent star formation activity in LRGs with redshift. By considering the expected lifetime of the k+a phase, the number of LRGs which will undergo a k+a phase can be estimated. A crude comparison of this estimate with the predictions from semi-analytic models of galaxy formation shows that the predicted level of k+a and em+a activities is not sufficient to reconcile the predicted mass growth for massive early types in a hierarchical merging scenario.

SHARDS: survey for high-Z absorption red & dead sources

SHARDS, an ESO/GTC Large Program, is an ultra-deep (26.5 mag) spectro-photometric survey with GTC/OSIRIS designed to select and study massive passively evolving galaxies at z=1.0-2.3 in the GOODS-N field using a set of 24 medium-band filters (FWHM ∼ 17 nm) covering the 500-950 nm spectral range. Our observing strategy has been planned to detect, for z>1 sources, the prominent Mg absorption feature (at rest-frame ∼ 280 nm), a distinctive, necessary, and sufficient feature of evolved stellar populations (older than 0.5 Gyr). These observations are being used to: (1) derive for the first time an unbiased sample of high-z quiescent galaxies, which extends to fainter magnitudes the samples selected with color techniques and spectroscopic surveys; (2) derive accurate ages and stellar masses based on robust measurements of spectral features such as the Mg_UV or D(4000) indices; (3) measure their redshift with an accuracy Δz/(1+z)<0.02; and (4) study emission-line galaxies (starbursts and AGN) up to very high redshifts. The well-sampled optical SEDs provided by SHARDS for all sources in the GOODS-N field are a valuable complement for current and future surveys carried out with other telescopes (e.g., Spitzer, HST, and Herschel).

Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 1043 erg s-1. Here, we introduce a simple treatment of clumping, and find that a filling factor of ˜0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C IV 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.

The VLT-FLAMES Tarantula Survey. XV. VFTS 822: A candidate Herbig B[e] star at low metallicity

We report the discovery of the B[e] star VFTS 822 in the 30 Doradus star-forming region of the Large Magellanic Cloud, classified by optical spectroscopy from the VLT-FLAMES Tarantula Survey and complementary infrared photometry. VFTS 822 is a relatively low-luminosity (log L = 4.04 ± 0.25 L_·) B8[e] star. In this Letter, we evaluate the evolutionary status of VFTS 822 and discuss its candidacy as a Herbig B[e] star. If the object is indeed in the pre-main sequence phase, it would present an exciting opportunity to spectroscopically measure mass accretion rates at low metallicity, to probe the effect of metallicity on accretion rates. 

The VLT-FLAMES Tarantula Survey:V. the peculiar B[e]-like supergiant, VFTS698, in 30 Doradus

Aims. We present an analysis of a peculiar supergiant B-type star (VFTS698/Melnick 2/Parker 1797) in the 30 Doradus region of the Large Magellanic Cloud which exhibits characteristics similar to the broad class of B[e] stars. Methods. We analyse optical spectra from the VLT-FLAMES survey, together with archival optical and infrared photometry and X-ray imaging to characterise the system. Results. We find radial velocity variations of around 400 km s ^-1 in the high excitation Si iv, N iii and He ii spectra, and photometric variability of ∼0.6 mag with a period of 12.7 d. In addition, we detect long-term photometric variations of ∼0.25 mag, which may be due to a longer-term variability with a period of ∼400 d. Conclusions. We conclude that VFTS698 is likely an interacting binary comprising an early B-type star secondary orbiting a veiled, more massive companion. Spectral evidence suggests a mid-to-late B-type primary, but this may originate from an optically-thick accretion disc directly surrounding the primary. © 2012 ESO.

LOTUS: A low cost, ultraviolet spectrograph

We describe the design, construction and commissioning of LOTUS; a simple, low-cost long-slit spectrograph for the Liverpool Telescope. The design is optimized for near-UV and visible wavelengths and uses all transmitting optics. It exploits the instrument focal plane field curvature to partially correct axial chromatic aberration. A stepped slit provides narrow (2.5x95 arcsec) and wide (5x25 arcsec) options that are optimized for spectral resolution and flux calibration respectively. On sky testing shows a wavelength range of 3200-6300 Angstroms with a peak system throughput (including detector quantum efficiency) of 15 per cent and wavelength dependant spectral resolution of R=225-430. By repeated observations of the symbiotic emission line star AG Peg we demonstrate the wavelength stability of the system is less than 2 Angstroms rms and is limited by the positioning of the object in the slit. The spectrograph is now in routine operation monitoring the activity of comet 67P/Churyumov-Gerasimenko during its current post-perihelion apparition.

BARIUM CONCENTRATIONS IN ROCK SALT BY LASER INDUCED BREAKDOWN SPECTROSCOPY

Time-REsolved Laser Induced Breakdown Spectroscopy (TRELIBS) was used to determine the elemental concentration of barium in Texas Dome rock salt. TRELIBS allows for an efficient and in situ concentration analysis technique that detects a wide range of elements with no sample preparation. TRELIBS measurements were made in the 545nm to 594nm wavelength range. The proximity of a strong barium emission line (553.5481 nm) to the sodium doublet (588.9950 nm and 589.5924 nm) allowed for measurement within a single frame of the spectrograph. This barium emission line was compared to the sodium doublet for relative intensity. A homemade calibration sample containing known amounts of barium and sodium was used to determine the relative concentrations. By approximating the sodium content in the rock salt as 50%, we determined the absolute concentration of barium in the salt to be (0.13±0.03)%.

Space density distribution of galaxies in the absolute magnitude - rotation velocity plane: a volume-complete Tully-Fisher relation from CALIFA stellar kinematics

We measured the distribution in absolute magnitude - circular velocity space for a well-defined sample of 199 rotating galaxies of the Calar Alto Legacy Integral Field Area Survey (CALIFA) using their stellar kinematics. Our aim in this analysis is to avoid subjective selection criteria and to take volume and large-scale structure factors into account. Using stellar velocity fields instead of gas emission line kinematics allows including rapidly rotating early-type galaxies. Our initial sample contains 277 galaxies with available stellar velocity fields and growth curve r-band photometry. After rejecting 51 velocity fields that could not be modelled because of the low number of bins, foreground contamination, or significant interaction, we performed Markov chain Monte Carlo modelling of the velocity fields, from which we obtained the rotation curve and kinematic parameters and their realistic uncertainties. We performed an extinction correction and calculated the circular velocity v_circ accounting for the pressure support of a given galaxy. The resulting galaxy distribution on the M-r - v(circ) plane was then modelled as a mixture of two distinct populations, allowing robust and reproducible rejection of outliers, a significant fraction of which are slow rotators. The selection effects are understood well enough that we were able to correct for the incompleteness of the sample. The 199 galaxies were weighted by volume and large-scale structure factors, which enabled us to fit a volume-corrected Tully-Fisher relation (TFR). More importantly, we also provide the volume-corrected distribution of galaxies in the M_r - v_circ plane, which can be compared with cosmological simulations. The joint distribution of the luminosity and circular velocity space densities, representative over the range of -20 > M_r > -22 mag, can place more stringent constraints on the galaxy formation and evolution scenarios than linear TFR fit parameters or the luminosity function alone.