The metal content of bulge field stars from FLAMES-GIRAFFE spectra - I. Stellar parameters and iron abundances

Aims. We determine the iron distribution function (IDF) for bulge field stars, in three different fields along the Galactic minor axis and at latitudes b = -4 degrees, b = -6 degrees, and b = -12 degrees. A fourth field including NGC 6553 is also included in the discussion. Methods. About 800 bulge field K giants were observed with the GIRAFFE spectrograph of FLAMES@VLT at spectral resolution R similar to 20 000. Several of them were observed again with UVES at R similar to 45 000 to insure the accuracy of the measurements. The LTE abundance analysis yielded stellar parameters and iron abundances that allowed us to construct an IDF for the bulge that, for the first time, is based on high-resolution spectroscopy for each individual star. Results. The IDF derived here is centered on solar metallicity, and extends from [Fe/H] similar to -1.5 to [Fe/H] similar to + 0.5. The distribution is asymmetric, with a sharper cutoff on the high-metallicity side, and it is narrower than previously measured. A variation in the mean metallicity along the bulge minor axis is clearly between b = -4 degrees and b = -6 degrees ([Fe/H] decreasing similar to by 0.6 dex per kpc). The field at b = -12 degrees. is consistent with the presence of a gradient, but its quantification is complicated by the higher disk/bulge fraction in this field. Conclusions. Our findings support a scenario in which both infall and outflow were important during the bulge formation, and then suggest the presence of a radial gradient, which poses some challenges to the scenario in which the bulge would result solely from the vertical heating of the bar.

OXYGEN ABUNDANCES IN LOW- AND HIGH-alpha FIELD HALO STARS AND THE DISCOVERY OF TWO FIELD STARS BORN IN GLOBULAR CLUSTERS

Oxygen abundances of 67 dwarf stars in the metallicity range -1.6 < [Fe/H] < -0.4 are derived from a non-LTE analysis of the 777 nm O I triplet lines. These stars have precise atmospheric parameters measured by Nissen and Schuster, who find that they separate into three groups based on their kinematics and alpha-element (Mg, Si, Ca, Ti) abundances: thick disk, high-alpha halo, and low-alpha halo. We find the oxygen abundance trends of thick-disk and high-alpha halo stars very similar. The low-alpha stars show a larger star-to-star scatter in [O/Fe] at a given [Fe/H] and have systematically lower oxygen abundances compared to the other two groups. Thus, we find the behavior of oxygen abundances in these groups of stars similar to that of the a elements. We use previously published oxygen abundance data of disk and very metal-poor halo stars to present an overall view (-2.3 < [Fe/H] < +0.3) of oxygen abundance trends of stars in the solar neighborhood. Two field halo dwarf stars stand out in their O and Na abundances. Both G53-41 and G150-40 have very low oxygen and very high sodium abundances, which are key signatures of the abundance anomalies observed in globular cluster (GC) stars. Therefore, they are likely field halo stars born in GCs. If true, we estimate that at least 3% +/- 2% of the local field metal-poor star population was born in GCs.

Library of medium-resolution Fiber Optic Echelle spectra of F, G, K, and M field dwarfs to giant stars

We present a library of Penn State Fiber Optic Echelle (FOE) observations of a sample of field stars with spectral types F to M and luminosity classes V to I. The spectral coverage is from 3800 to 10000 Å with a nominal resolving power of 12,000. These spectra include many of the spectral lines most widely used as optical and near-infrared indicators of chromospheric activity such as the Balmer lines (Hα to H epsilon), Ca II H & K, the Mg I b triplet, Na I D_1, D_2, He I D_3, and Ca II IRT lines. There are also a large number of photospheric lines, which can also be affected by chromospheric activity, and temperature-sensitive photospheric features such as TiO bands. The spectra have been compiled with the goal of providing a set of standards observed at medium resolution. We have extensively used such data for the study of active chromosphere stars by applying a spectral subtraction technique. However, the data set presented here can also be utilized in a wide variety of ways ranging from radial velocity templates to study of variable stars and stellar population synthesis. This library can also be used for spectral classification purposes and determination of atmospheric parameters (T_eff, log g, [Fe/H]). A digital version of all the fully reduced spectra is available via ftp and the World Wide Web (WWW) in FITS format.

The imprint of the crustal magnetic field on the thermal spectra and pulse profiles of isolated neutron stars

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the-art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double-peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50–60 per cent level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated ‘observed’ spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ∼1–2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements.

Li-rich red giant branch stars in the Galactic bulge

Aims. We present lithium abundance determination for a sample of K giant stars in the Galactic bulge. The stars presented here are the only 13 stars with a detectable lithium line (6767.18 angstrom) among similar to 400 stars for which we have spectra in this wavelength range, half of them in Baade's Window (b = -4 degrees) and half in a field at b = -6 degrees. Methods. The stars were observed with the GIRAFFE spectrograph of FLAMES mounted on VLT, with a spectral resolution of R similar to 20 000. Abundances were derived from spectral synthesis and the results are compared with those of stars with similar parameters, but no detectable Li line. Results. We find 13 stars with a detectable Li line, among which 2 have abundances A(Li) > 2.7. No clear correlations were found between the Li abundance and those of other elements. With the exception of the two most Li rich stars, the others follow a fairly tight A(Li) - T(eff) correlation. Conclusions. There is strong indication of a Li production phase during the red giant branch (RGB), acting either on a very short timescale, or selectively only in some stars. That the proposed Li production phase is associated with the RGB bump cannot be excluded, although our targets are significantly brighter than the predicted RGB bump magnitude for a population at 8 kpc.

Detection of moving groups among early type stars

Our procedure to detect moving groups in the solar neighbourhood (Chen et al., 1997) in the four-dimensional space of the stellar velocity components and age has been improved. The method, which takes advantadge of non-parametric estimators of density distribution to avoid any a priori knowledge of the kinematic properties of these stellar groups, now includes the effect of observational errors on the process to select moving group stars, uses a better estimation of the density distribution of the total sample and field stars, and classifies moving group stars using all the available information. It is applied here to an accurately selected sample of early-type stars with known radial velocities and Strömgren photometry. Astrometric data are taken from the HIPPARCOS catalogue (ESA, 1997), which results in an important decrease in the observational errors with respect to ground-based data, and ensures the uniformity of the observed data. Both the improvement of our method and the use of precise astrometric data have allowed us not only to confirm the existence of classical moving groups, but also to detect finer structures that in several cases can be related to kinematic properties of nearby open clusters or associations.

A spectroscopy study of nearby late-type stars, possible members of stellar kinematic groups

Context. Nearby late-type stars are excellent targets for seeking young objects in stellar associations and moving groups. The origin of these structures is still misunderstood, and lists of moving group members often change with time and also from author to author. Most members of these groups have been identified by means of kinematic criteria, leading to an important contamination of previous lists by old field stars. Aims. We attempt to identify unambiguous moving group members among a sample of nearby-late type stars by studying their kinematics, lithium abundance, chromospheric activity, and other age-related properties. Methods. High-resolution echelle spectra (R ~ 57 000) of a sample of nearby late-type stars are used to derive accurate radial velocities that are combined with the precise Hipparcos parallaxes and proper motions to compute galactic-spatial velocity components. Stars are classified as possible members of the classical moving groups according to their kinematics. The spectra are also used to study several age-related properties for young late-type stars, i.e., the equivalent width of the lithium Li i 6707.8 Å line or the R'_HK index. Additional information like X-ray fluxes from the ROSAT All-Sky Survey or the presence of debris discs is also taken into account. The different age estimators are compared and the moving group membership of the kinematically selected candidates are discussed. Results. From a total list of 405 nearby stars, 102 have been classified as moving group candidates according to their kinematics. i.e., only ~25.2% of the sample. The number reduces when age estimates are considered, and only 26 moving group candidates (25.5% of the 102 candidates) have ages in agreement with the star having the same age as an MG member.

Quantifying the contamination by old main-sequence stars in young moving groups: the case of the Local Association

Context. The associations and moving groups of young stars are excellent laboratories for investigating stellar formation in the solar neighborhood. Previous results have confirmed that a non-negligible fraction of old main-sequence stars is present in the lists of possible members of young stellar kinematic groups. A detailed study of the properties of these samples is needed to separate the young stars from old main-sequence stars with similar space motion, and identify the origin of these structures. Aims. Our intention is to characterize members of the young moving groups, determine their age distribution, and quantify the contamination by old main-sequence stars, in particular, for the Local Association. Methods. We used stars possible members of the young (~10-650 Myr) moving groups from the literature. To determine the age of the stars, we used several suitable age indicators for young main sequence stars, i.e., X-ray fluxes from the Rosat All-sky Survey database, photometric data from the Tycho-2, Hipparcos, and 2MASS database. We also used spectroscopic data, in particular the equivalent width of the lithium line Li I λ6707.8 Å and H_α, to constrain the range of ages of the stars. Results. By combining photometric and spectroscopic data, we were able to separate the young stars (10-650 Myr) from the old (> 1 Gyr) field ones. We found, in particular, that the Local Association is contaminated by old field stars at the level of ~30%. This value must be considered as the contamination for our particular sample, and not of the entire Local Association. For other young moving groups, it is more difficult to estimate the fraction of old stars among possible members. However, the level of X-ray emission can, at least, help to separate two age populations: stars with <200 Myr and stars older than this. Conclusions. Among the candidate members of the classical moving groups, there is a non-negligible fraction of old field stars that should be taken into account when studying the stellar birthrate in the solar neighborhood. Our results are consistent with a scenario in which the moving groups contain both groups of young stars formed in a recent star-formation episode and old field stars with similar space motion. Only by combining X-ray and optical spectroscopic data is it possible to distinguish between these two age populations.

Force-free twisted magnetospheres of neutron stars

Context. The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. Aims. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. Methods. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. Results. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. Conclusions. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.

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.

Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radio-pulsars, thus confirming the idea that magnetic fields are involved in their X-ray emission. Here we present the results of 2D simulations of the fully coupled evolution of temperature and magnetic field in neutron stars, including the state-of-the-art kinetic coefficients and, for the first time, the important effect of the Hall term. After gathering and thoroughly re-analysing in a consistent way all the best available data on isolated, thermally emitting neutron stars, we compare our theoretical models to a data sample of 40 sources. We find that our evolutionary models can explain the phenomenological diversity of magnetars, high-B radio-pulsars, and isolated nearby neutron stars by only varying their initial magnetic field, mass and envelope composition. Nearly all sources appear to follow the expectations of the standard theoretical models. Finally, we discuss the expected outburst rates and the evolutionary links between different classes. Our results constitute a major step towards the grand unification of the isolated neutron star zoo.

The outburst decay of the low magnetic field magnetar SGR 0418+5729

We report on the long-term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729 using all the available X-ray data obtained with RXTE, Swift, Chandra, and XMM-Newton observations from the discovery of the source in 2009 June up to 2012 August. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: ˙ P = 4(1) × 10−15 s s−1, significant at a ∼3.5σ confidence level. This leads to a surface dipolar magnetic field of Bdip 6 × 1012 G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ∼1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ∼0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about three orders of magnitude: from about 1.4 × 10−11 to 1.2 × 10−14 erg s−1 cm−2. Deep radio, millimeter, optical, and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ∼550 kyr, and a dipolar magnetic field at birth of ∼1014 G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such a result in several other astrophysical contexts.

Central compact objects and the hidden magnetic field scenario

Central compact objects (CCOs) are X-ray sources lying close to the centre of supernova remnants, with inferred values of the surface magnetic fields significantly lower (≲1011 G) than those of standard pulsars. In this paper, we revise the hidden magnetic field scenario, presenting the first 2D simulations of the submergence and re-emergence of the magnetic field in the crust of a neutron star. A post-supernova accretion stage of about 10−4–10−3 M⊙ over a vast region of the surface is required to bury the magnetic field into the inner crust. When accretion stops, the field re-emerges on a typical time-scale of 1–100 kyr, depending on the submergence conditions. After this stage, the surface magnetic field is restored close to its birth values. A possible observable consequence of the hidden magnetic field is the anisotropy of the surface temperature distribution, in agreement with observations of several of these sources. We conclude that the hidden magnetic field model is viable as an alternative to the antimagnetar scenario, and it could provide the missing link between CCOs and the other classes of isolated neutron stars.

Population synthesis of isolated neutron stars with magneto-rotational evolution

We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in parameter space that are more consistent with the observational data. In particular, we find that any broad distribution of birth spin periods with P0 ≲ 0.5 s can fit the data, and that if the alignment angle is allowed to vary consistently with the torque model, realistic magnetospheric models are favoured compared to models with classical magneto-dipolar radiation losses. Assuming that the initial magnetic field is given by a lognormal distribution, our optimal model has mean strength 〈log B0[G]〉 ≈ 13.0–13.2 with width σ(log B0) = 0.6–0.7. However, there are strong correlations between parameters. This degeneracy in the parameter space can be broken by an independent estimate of the pulsar birth rate or by future studies correlating this information with the population in other observational bands (X-rays and γ-rays).

Population synthesis of isolated neutron stars with magneto-rotational evolution – II. From radio-pulsars to magnetars

Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include X-ray thermal emitting pulsars within the same evolutionary model as radio-pulsars. We find that the cumulative distribution of the number of X-ray pulsars can be well reproduced by several models that, simultaneously, reproduce the characteristics of the radio-pulsar distribution. However, even considering the most favourable magneto-thermal evolution models with fast field decay, lognormal distributions of the initial magnetic field overpredict the number of visible sources with periods longer than 12 s. We then show that the problem can be solved with different distributions of magnetic field, such as a truncated lognormal distribution, or a binormal distribution with two distinct populations. We use the observational lack of isolated neutron stars (NSs) with spin periods P > 12 s to establish an upper limit to the fraction of magnetars born with B > 1015 G (less than 1 per cent). As future detections keep increasing the magnetar and high-B pulsar statistics, our approach can be used to establish a severe constraint on the maximum magnetic field at birth of NSs.