The little-studied cluster Berkeley 90 I. LS III +46 11: a very massive O3.5 If* + O3.5 If* binary

Context. It appears that most (if not all) massive stars are born in multiple systems. At the same time, the most massive binaries are hard to find owing to their low numbers throughout the Galaxy and the implied large distances and extinctions. Aims. We want to study LS III +46 11, identified in this paper as a very massive binary; another nearby massive system, LS III +46 12; and the surrounding stellar cluster, Berkeley 90. Methods. Most of the data used in this paper are multi-epoch high S/N optical spectra, although we also use Lucky Imaging and archival photometry. The spectra are reduced with dedicated pipelines and processed with our own software, such as a spectroscopic-orbit code, CHORIZOS, and MGB. Results. LS III +46 11 is identified as a new very early O-type spectroscopic binary [O3.5 If* + O3.5 If*] and LS III +46 12 as another early O-type system [O4.5 V((f))]. We measure a 97.2-day period for LS III +46 11 and derive minimum masses of 38.80 ± 0.83 M⊙ and 35.60 ± 0.77 M⊙ for its two stars. We measure the extinction to both stars, estimate the distance, search for optical companions, and study the surrounding cluster. In doing so, a variable extinction is found as well as discrepant results for the distance. We discuss possible explanations and suggest that LS III +46 12 may be a hidden binary system where the companion is currently undetected.

Earth’s Rotation: A Challenging Problem in Mathematics and Physics

A suitable knowledge of the orientation and motion of the Earth in space is a common need in various fields. That knowledge has been ever necessary to carry out astronomical observations, but with the advent of the space age, it became essential for making observations of satellites and predicting and determining their orbits, and for observing the Earth from space as well. Given the relevant role it plays in Space Geodesy, Earth rotation is considered as one of the three pillars of Geodesy, the other two being geometry and gravity. Besides, research on Earth rotation has fostered advances in many fields, such as Mathematics, Astronomy and Geophysics, for centuries. One remarkable feature of the problem is in the extreme requirements of accuracy that must be fulfilled in the near future, about a millimetre on the tangent plane to the planet surface, roughly speaking. That challenges all of the theories that have been devised and used to-date; the paper makes a short review of some of the most relevant methods, which can be envisaged as milestones in Earth rotation research, emphasizing the Hamiltonian approach developed by the authors. Some contemporary problems are presented, as well as the main lines of future research prospected by the International Astronomical Union/International Association of Geodesy Joint Working Group on Theory of Earth Rotation, created in 2013.

Pulsar timing irregularities and the imprint of magnetic field evolution

Context. The rotational evolution of isolated neutron stars is dominated by the magnetic field anchored to the solid crust of the star. Assuming that the core field evolves on much longer timescales, the crustal field evolves mainly though Ohmic dissipation and the Hall drift, and it may be subject to relatively rapid changes with remarkable effects on the observed timing properties. Aims. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable consequences in the braking index n. This is the most sensitive quantity to reflect small variations of the timing properties that are caused by magnetic field rearrangements. Methods. We performed axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art microphysical inputs to calculate the evolution of the braking index. Relatively rapid magnetic field modifications can be expected only in the crust of neutron stars, where we focus our study. Results. We find that the effect of the magnetic field evolution on the braking index can be divided into three qualitatively different stages depending on the age and the internal temperature: a first stage that may be different for standard pulsars (with n ~ 3) or low field neutron stars that accreted fallback matter during the supernova explosion (systematically n < 3); in a second stage, the evolution is governed by almost pure Ohmic field decay, and a braking index n > 3 is expected; in the third stage, at late times, when the interior temperature has dropped to very low values, Hall oscillatory modes in the neutron star crust result in braking indices of a high absolute value and both positive and negative signs. Conclusions. Current magneto-thermal evolution models predict a large contribution to the timing noise and, in particular, to the braking index, from temporal variations of the magnetic field. Models with strong (≳ 1014 G) multipolar or toroidal components, even with a weak (~1012 G) dipolar field are consistent with the observed trend of the timing properties.

A VLT/FLAMES survey for massive binaries in Westerlund 1. III. The WC9d binary W239 and implications for massive stellar evolution

Context. There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars – frequency, mass ratio & orbital separation – are still poorly constrained. Aims. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-IR photometric observations in order to search for binary signatures. Comparison of these data to synthetic non-LTE model atmosphere spectra were used to derive the fundamental properties of the WC9 primary. Results. We found W239 to have an orbital period of only ~5.05 days, making it one of the most compact WC binaries yet identified. Analysis of the long term near-IR lightcurve reveals a significant flare between 2004-6. We interpret this as evidence for a third massive stellar component in the system in a long period (>6 yr), eccentric orbit, with dust production occuring at periastron leading to the flare. The presence of a near-IR excess characteristic of hot (~1300 K) dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short (<1 yr) period binaries, although confirmation will require further observations. Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent with other Galactic examples. Conclusions. The simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss – likely via case A mass transfer or a contact configuration – or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants.

The masses of the neutron and donor star in the high-mass X-ray binary IGR J18027-2016

Aims. We report near-infrared observations of the supergiant donor to the eclipsing high mass X-ray binary pulsar IGR J18027-2016. We aim to determine its spectral type and measure its radial velocity curve and hence determine the stellar masses of the components. Methods. ESO/VLT observations of the donor utilising the NIR spectrograph ISAAC were obtained in the H and K bands. The multi-epoch H band spectra were cross-correlated with RV templates in order to determine a radial solution for the system. Results. The spectral type of the donor was confirmed as B0-1 I. The radial velocity curve constructed has a semi-amplitude of 23.8 ± 3.1 km s-1. Combined with other measured system parameters, a dynamically determined neutron star mass of 1.4  ±  0.2–1.6  ±  0.3 M⊙ is found. The mass range of the B0-B1 I donor was 18.6  ±  0.8–21.8  ±  2.4 M⊙. These lower and upper limits were obtained under the assumption that the system is viewed edge-on (i = 90° with β = 0.89) for the lower limit and the donor fills its Roche lobe (β = 1 with i = 73.1°) for the upper limit respectively.

Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way

Aims. We study the optical and near-infrared colour excesses produced by circumstellar emission in a sample of Be/X-ray binaries. Our main goals are exploring whether previously published relations, valid for isolated Be stars, are applicable to Be/X-ray binaries and computing the distance to these systems after correcting for the effects of the circumstellar contamination. Methods. Simultaneous UBVRI photometry and spectra in the 3500−7000 Å spectral range were obtained for 11 optical counterparts to Be/X-ray binaries in the LMC, 5 in the SMC and 12 in the Milky Way. As a measure of the amount of circumstellar emission we used the Hα equivalent width corrected for photospheric absorption. Results. We find a linear relationship between the strength of the Hα emission line and the component of E(B − V) originating from the circumstellar disk. This relationship is valid for stars with emission lines weaker than EW ≈ −15   Å. Beyond this point, the circumstellar contribution to E(B − V) saturates at a value ≈0.17   mag. A similar relationship is found for the (V − I) near infrared colour excess, albeit with a steeper slope and saturation level. The circumstellar excess in (B − V) is found to be about five times higher for Be/X-ray binaries than for isolated Be stars with the same equivalent width EW(Hα), implying significant differences in the physical properties of their circumstellar envelopes. The distance to Be/X-ray binaries (with non-shell Be star companions) can only be correctly estimated by taking into account the excess emission in the V band produced by free-free and free-bound transitions in the circumstellar envelope. We provide a simple method to determine the distances that includes this effect.

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.

No anticorrelation between cyclotron line energy and X-ray flux in 4U 0115+634

We report on an outburst of the high mass X-ray binary 4U 0115+634 with a pulse period of 3.6 s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star’s luminosity varied by a factor of 10 in the 3–50 keV band. In agreement with earlier work we find evidence of five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7 keV. Previous work had found an anticorrelation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anticorrelation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10 keV, the correlation between the line energy and the flux vanishes, and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anticorrelation is an artifact of the assumptions adopted in the modeling of the continuum.

Astrophysical parameters and orbital solution of the peculiar X-ray transient IGR J00370+6122

Context. BD + 60° 73 is the optical counterpart of the X-ray source IGR J00370+6122, a probable accretion-powered X-ray pulsar. The X-ray light curve of this binary system shows clear periodicity at 15.7 d, which has been interpreted as repeated outbursts around the periastron of an eccentric orbit. Aims. We aim to characterise the binary system IGR J00370+6122 by deriving its orbital and physical parameters. Methods. We obtained high-resolution spectra of BD + 60° 73 at different epochs. We used the fastwind code to generate a stellar atmosphere model to fit the observed spectrum and obtain physical magnitudes. The synthetic spectrum was used as a template for cross-correlation with the observed spectra to measure radial velocities. The radial velocity curve provided an orbital solution for the system. We also analysed the RXTE/ASM and Swift/BAT light curves to confirm the stability of the periodicity. Results. BD + 60° 73 is a BN0.7 Ib low-luminosity supergiant located at a distance ~3.1 kpc, in the Cas OB4 association. We derive Teff = 24 000 K and log gc = 3.0, and chemical abundances consistent with a moderately high level of evolution. The spectroscopic and evolutionary masses are consistent at the 1-σ level with a mass M∗ ≈ 15 M⊙. The recurrence time of the X-ray flares is the orbital period of the system. The neutron star is in a high-eccentricity (e = 0.56 ± 0.07) orbit, and the X-ray emission is strongly peaked around orbital phase φ = 0.2, though the observations are consistent with some level of X-ray activity happening at all orbital phases. Conclusions. The X-ray behaviour of IGR J00370+6122 is reminiscent of “intermediate” supergiant X-ray transients, though its peak luminosity is rather low. The orbit is somewhat wider than those of classical persistent supergiant X-ray binaries, which when combined with the low luminosity of the mass donor, explains the low X-ray luminosity. IGR J00370+6122 will very likely evolve towards a persistent supergiant system, highlighting the evolutionary connection between different classes of wind-accreting X-ray sources.

New galactic star clusters discovered in the VVV survey. Candidates projected on the inner disk and bulge

Context. VISTA Variables in the Vía Láctea (VVV) is one of six ESO Public Surveys using the 4 meter Visible and Infrared Survey Telescope for Astronomy (VISTA). The VVV survey covers the Milky Way bulge and an adjacent section of the disk, and one of the principal objectives is to search for new star clusters within previously unreachable obscured parts of the Galaxy. Aims. The primary motivation behind this work is to discover and analyze obscured star clusters in the direction of the inner Galactic disk and bulge. Methods. Regions of the inner disk and bulge covered by the VVV survey were visually inspected using composite JHKS color images to select new cluster candidates on the basis of apparent overdensities. DR1, DR2, CASU, and point spread function photometry of 10 × 10 arcmin fields centered on each candidate cluster were used to construct color–magnitude and color–color diagrams. Follow-up spectroscopy of the brightest members of several cluster candidates was obtained in order to clarify their nature. Results. We report the discovery of 58 new infrared cluster candidates. Fundamental parameters such as age, distance, and metallicity were determined for 20 of the most populous clusters.

A VLT/FLAMES survey for massive binaries in Westerlund 1. IV. Wd1-5 – binary product and a pre-supernova companion for the magnetar CXOU J1647-45?

Context. The first soft gamma-ray repeater was discovered over three decades ago, and was subsequently identified as a magnetar, a class of highly magnetised neutron star. It has been hypothesised that these stars power some of the brightest supernovae known, and that they may form the central engines of some long duration gamma-ray bursts. However there is currently no consenus on the formation channel(s) of these objects. Aims. The presence of a magnetar in the starburst cluster Westerlund 1 implies a progenitor with a mass ≥40 M⊙, which favours its formation in a binary that was disrupted at supernova. To test this hypothesis we conducted a search for the putative pre-SN companion. Methods. This was accomplished via a radial velocity survey to identify high-velocity runaways, with subsequent non-LTE model atmosphere analysis of the resultant candidate, Wd1-5. Results. Wd1-5 closely resembles the primaries in the short-period binaries, Wd1-13 and 44, suggesting a similar evolutionary history, although it currently appears single. It is overluminous for its spectroscopic mass and we find evidence of He- and N-enrichement, O-depletion, and critically C-enrichment, a combination of properties that is difficult to explain under single star evolutionary paradigms. We infer a pre-SN history for Wd1-5 which supposes an initial close binary comprising two stars of comparable (~ 41 M⊙ + 35 M⊙) masses. Efficient mass transfer from the initially more massive component leads to the mass-gainer evolving more rapidly, initiating luminous blue variable/common envelope evolution. Reverse, wind-driven mass transfer during its subsequent WC Wolf-Rayet phase leads to the carbon pollution of Wd1-5, before a type Ibc supernova disrupts the binary system. Under the assumption of a physical association between Wd1-5 and J1647-45, the secondary is identified as the magnetar progenitor; its common envelope evolutionary phase prevents spin-down of its core prior to SN and the seed magnetic field for the magnetar forms either in this phase or during the earlier episode of mass transfer in which it was spun-up. Conclusions. Our results suggest that binarity is a key ingredient in the formation of at least a subset of magnetars by preventing spin-down via core-coupling and potentially generating a seed magnetic field. The apparent formation of a magnetar in a Type Ibc supernova is consistent with recent suggestions that superluminous Type Ibc supernovae are powered by the rapid spin-down of these objects.

The Gaia-ESO Survey: radial metallicity gradients and age-metallicity relation of stars in the Milky Way disk

We study the relationship between age, metallicity, and α-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpcto 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages >9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more α-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

The Gaia-ESO Survey: Empirical determination of the precision of stellar radial velocities and projected rotation velocities

Context. The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. Aims. A key aim is to provide precise radial velocities (RVs) and projected equatorial velocities (vsini) for representative samples of Galactic stars, which will complement information obtained by the Gaia astrometry satellite. Methods. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and vsini using spectra from repeated exposures of the same stars. Results. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S/N) and vsini, that the uncertainties become larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty distributions have extended tails that are better represented by Student’s t-distributions than by normal distributions. Conclusions. Parametrised results are provided, which enable estimates of the RV precision for almost all GES measurements, and estimates of the vsini precision for stars in young clusters, as a function of S/N, vsini and stellar temperature. The precision of individual high S/N GES RV measurements is 0.22–0.26 km s-1, dependent on instrumental configuration.

Open cluster Dolidze 25: Stellar parameters and the metallicity in the Galactic anticentre

Context. The young open cluster Dolidze 25, in the direction of the Galactic anticentre, has been attributed a very low metallicity, with typical abundances between −0.5 and −0.7 dex below solar. Aims. We intend to derive accurate cluster parameters and accurate stellar abundances for some of its members. Methods. We have obtained a large sample of intermediate- and high-resolution spectra for stars in and around Dolidze 25. We used the fastwind code to generate stellar atmosphere models to fit the observed spectra. We derive stellar parameters for a large number of OB stars in the area, and abundances of oxygen and silicon for a number of stars with spectral types around B0. Results. We measure low abundances in stars of Dolidze 25. For the three stars with spectral types around B0, we find 0.3 dex (Si) and 0.5 dex (O) below the values typical in the solar neighbourhood. These values, even though not as low as those given previously, confirm Dolidze 25 and the surrounding H ii region Sh2-284 as the most metal-poor star-forming environment known in the Milky Way. We derive a distance 4.5 ± 0.3 kpc to the cluster (rG ≈ 12.3 kpc). The cluster cannot be older than ~3 Myr, and likely is not much younger. One star in its immediate vicinity, sharing the same distance, has Si and O abundances at most 0.15 dex below solar. Conclusions. The low abundances measured in Dolidze 25 are compatible with currently accepted values for the slope of the Galactic metallicity gradient, if we take into account that variations of at least ±0.15 dex are observed at a given radius. The area traditionally identified as Dolidze 25 is only a small part of a much larger star-forming region that comprises the whole dust shell associated with Sh2-284 and very likely several other smaller H ii regions in its vicinity.