Astrophysical parameters of LS 2883 and implications for the PSR B1259-63 gamma-ray binary

Only a few binary systems with compact objects display TeV emission. The physical properties of the companion stars represent basic input for understanding the physical mechanisms behind the particle acceleration, emission, and absorption processes in these so-called gamma-ray binaries. Here we present high-resolution and high signal-to-noise optical spectra of LS 2883, the Be star forming a gamma-ray binary with the young non-accreting pulsar PSR B1259-63, showing it to rotate faster and be significantly earlier and more luminous than previously thought. Analysis of the interstellar lines suggests that the system is located at the same distance as (and thus is likely a member of) Cen OB1. Taking the distance to the association, d = 2.3 kpc, and a color excess of E(B – V) = 0.85 for LS 2883 results in MV ≈ –4.4. Because of fast rotation, LS 2883 is oblate (R eq sime 9.7 R ☉ and R pole sime 8.1 R ☉) and presents a temperature gradient (T eq≈ 27,500 K, log g eq = 3.7; T pole≈ 34,000 K, log g pole = 4.1). If the star did not rotate, it would have parameters corresponding to a late O-type star. We estimate its luminosity at log(L */L ☉) sime 4.79 and its mass at M * ≈ 30 M ☉. The mass function then implies an inclination of the binary system i orb ≈ 23°, slightly smaller than previous estimates. We discuss the implications of these new astrophysical parameters of LS 2883 for the production of high-energy and very high-energy gamma rays in the PSR B1259-63/LS 2883 gamma-ray binary system. In particular, the stellar properties are very important for prediction of the line-like bulk Comptonization component from the unshocked ultrarelativistic pulsar wind.

NGC 7419 as a template for red supergiant clusters

Context. The open cluster NGC 7419 is known to contain five red supergiants and a very high number of Be stars. However, there are conflicting reports about its age and distance that prevent a useful comparison with other clusters. Aims. We intend to obtain more accurate parameters for NGC 7419, using techniques different from those of previous authors, so that it may be used as a calibrator for more obscured clusters. Methods. We obtained Strömgren photometry of the open cluster NGC 7419, as well as classification spectroscopy of ~20 stars in the area. We then applied standard analysis and classification techniques. Results. We find a distance of 4 ± 0.4 kpc and an age of 14 ± 2 Myr for NGC 7419. The main-sequence turn-off is found at spectral type B1, in excellent agreement. We identify 179 B-type members, implying that there are more than 1200 M⊙ in B stars at present. Extrapolating this to lower masses indicates an initial cluster mass of between 7000 and 10 000 M⊙, depending on the shape of the initial mass function. We find a very high fraction (≈40%) of Be stars around the turn-off, but very few Be stars at lower masses. We also report for the first time a strong variability in the emission characteristics of Be stars. We verified that the parameters of the red supergiant members can be used to obtain accurate cluster parameters. Conclusions. NGC 7419 is sufficiently massive to serve as a testbed for theoretical predictions and as a template to compare more obscured clusters. The distribution of stars above the main-sequence turn-off is difficult to accommodate with current evolutionary tracks. Though the presence of five red supergiants is marginally consistent with theoretical expectations, the high number of Be stars and very low number of luminous evolved B stars hint at some unknown physical factor that is not considered in current synthesis models.

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.

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.

Hot thermal X-ray emission from the Be star HD 119682

We present an analysis of a series of four consecutive Chandra high-resolution transmission gratings observations, amounting to a total of 150 ks, of the Be X-ray source HD 119682 (=1WGA J1346.5–6255), a member of the new class of γ Cas analogs. The Chandra light curve shows significant brightness variations on timescales of hours. However, the spectral distribution appears rather stable within each observation and during the whole campaign. A detailed analysis is not able to detect any coherent pulsation up to a frequency of 0.05 Hz. The Chandra High Energy Transmission Gratings spectrum seems to be devoid of any strong emission line, including Fe Kα fluorescence. The continuum is well described with the addition of two collisionally ionized plasmas of temperatures kT ≈ 15 keV and 0.2 keV, respectively, by the apec model. Models using photoionized plasma components (mekal) or non-thermal components (powerlaw) give poorer fits, providing support for the pure thermal scenario. These two components are absorbed by a single column with N H = (0.20+0.15 –0.03) × 1022 cm–2 compatible with the interstellar value. We conclude that HD 119682 can be regarded as a pole-on γ Cas analog.

Chandra and Suzaku observations of the Be/X-ray star HD110432

We present an analysis of a pointed 141 ks Chandra high-resolution transmission gratings observation of the Be X-ray emitting star HD110432, a prominent member of the γ Cas analogs. This observation represents the first high-resolution spectrum taken for this source as well as the longest uninterrupted observation of any γ Cas analog. The Chandra light curve shows a high variability but its analysis fails to detect any coherent periodicity up to a frequency of 0.05 Hz. Hardness ratio versus intensity analyses demonstrate that the relative contributions of the [1.5-3] Å, [3-6] Å, and [6-16] Å energy bands to the total flux change rapidly in the short term. The analysis of the Chandra High Energy Transmission Grating (HETG) spectrum shows that, to correctly describe the spectrum, three model components are needed. Two of those components are optically thin thermal plasmas of different temperatures (kT ≈ 8-9 and 0.2-0.3 keV, respectively) described by the models vmekal or bvapec. The Fe abundance in each of these two components appears equal within the errors and is slightly subsolar with Z ≈ 0.75 Z ☉. The bvapec model better describes the Fe L transitions, although it cannot fit well the Na XI Lyα line at 10.02 Å, which appears to be overabundant. Two different models seem to describe well the third component. One possibility is a third hot optically thin thermal plasma at kT = 16-21 keV with an Fe abundance Z ≈ 0.3 Z ☉, definitely smaller than for the other two thermal components. Furthermore, the bvapec model describes well the Fe K shell transitions because it accounts for the turbulence broadening of the Fe XXV and Fe XXVI lines with a v turb ≈ 1200 km s–1. These two lines, contributed mainly by the hot thermal plasma, are significantly wider than the Fe Kα line whose FWHM < 5 mÅ is not resolved by Chandra. Alternatively, the third component can be described by a power law with a photon index of Γ = 1.56. In either case, the Chandra HETG spectrum establishes that each one of these components must be modified by distinct absorption columns. The analysis of a noncontemporaneous 25 ks Suzaku observation shows the presence of a hard tail extending up to at least 33 keV. The Suzaku spectrum is described with the sum of two components: an optically thin thermal plasma at kT ≈ 9 keV and Z ≈ 0.74 Z ☉, and a very hot second plasma with kT ≈ 33 keV or, alternatively, a power law with photon index of Γ = 1.58. In either case, each one of the two components must be affected by different absorption columns. Therefore, the kT = 8-9 keV component is definitely needed while the nature of the harder emission cannot be unambiguously established with the present data sets. The analysis of the Si XIII and S XV He-like triplets present in the Chandra spectrum points to a very dense (ne ~ 1013 cm–3) plasma located either close to the stellar surface (r < 3R *) of the Be star or, alternatively, very close (r ~ 1.5R WD) to the surface of a (hypothetical) white dwarf companion. We argue, however, that the available data support the first scenario.

A Be-type star with a black-hole companion

Stellar-mass black holes have all been discovered through X-ray emission, which arises from the accretion of gas from their binary companions (this gas is either stripped from low-mass stars or supplied as winds from massive ones). Binary evolution models also predict the existence of black holes accreting from the equatorial envelope of rapidly spinning Be-type stars1, 2, 3 (stars of the Be type are hot blue irregular variables showing characteristic spectral emission lines of hydrogen). Of the approximately 80 Be X-ray binaries known in the Galaxy, however, only pulsating neutron stars have been found as companions2, 3, 4. A black hole was formally allowed as a solution for the companion to the Be star MWC 656 (ref. 5; also known as HD 215227), although that conclusion was based on a single radial velocity curve of the Be star, a mistaken spectral classification6 and rough estimates of the inclination angle. Here we report observations of an accretion disk line mirroring the orbit of MWC 656. This, together with an improved radial velocity curve of the Be star through fitting sharp Fe ii profiles from the equatorial disk, and a refined Be classification (to that of a B1.5–B2 III star), indicates that a black hole of 3.8 to 6.9 solar masses orbits MWC 656, the candidate counterpart of the γ-ray source AGL J2241+4454 (refs 5, 6). The black hole is X-ray quiescent and fed by a radiatively inefficient accretion flow giving a luminosity less than 1.6 × 10−7 times the Eddington luminosity. This implies that Be binaries with black-hole companions are difficult to detect in conventional X-ray surveys.

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.

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.

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

On-line microwave-based preconcentration device for inductively coupled plasma atomic emission spectrometry: Application to the elemental analysis of spirit samples

A microwave-based thermal nebulizer (MWTN) has been employed for the first time as on-line preconcentration device in inductively coupled plasma atomic emission spectrometry (ICP-AES). By the appropriate selection of the experimental conditions, the MWTN could be either operated as a conventional thermal nebulizer or as on-line analyte preconcentration and nebulization device. Thus, when operating at microwave power values above 100 W and highly concentrated alcohol solutions, the amount of energy per solvent mass liquid unit (EMR) is high enough to completely evaporate the solvent inside the system and, as a consequence, the analyte is deposited (and then preconcentrated) on the inner walls of the MWTN capillary. When reducing the EMR to the appropriate value (e.g., by reducing the microwave power at a constant sample uptake rate) the retained analyte is swept along by the liquid-gas stream and an analyte-enriched aerosol is generated and next introduced into the plasma cell. Emission signals obtained with the MWTN operating in preconcentration-nebulization mode improved when increasing preconcentration time and sample uptake rate as well as when decreasing the nozzle inner diameter. When running with pure ethanol solution at its optimum experimental conditions, the MWTN in preconcentration-nebulization mode afforded limits of detection up to one order of magnitude lowers than those obtained operating the MWTN exclusively as a nebulizer. To validate the method, the multi-element analysis (i.e. Al, Ca, Cd, Cr, Cu, Fe, K, Mg, Mn, Na, Pb and Zn) of different commercial spirit samples in ICP-AES has been performed. Analyte recoveries for all the elements studied ranged between 93% and 107% and the dynamic linear range covered up to 4 orders of magnitude (i.e. from 0.1 to 1000 μg L−1). In these analysis, both MWTN operating modes afforded similar results. Nevertheless, the preconcentration-nebulization mode permits to determine a higher number of analytes due to its higher detection capabilities.