New Method for the Estimation of Viscosity of Pure and Mixtures of Ionic Liquids Based on the UNIFAC–VISCO Model

A modified UNIFAC–VISCO group contribution method was developed for the correlation and prediction of viscosity of ionic liquids as a function of temperature at 0.1 MPa. In this original approach, cations and anions were regarded as peculiar molecular groups. The significance of this approach comes from the ability to calculate the viscosity of mixtures of ionic liquids as well as pure ionic liquids. Binary interaction parameters for selected cations and anions were determined by fitting the experimental viscosity data available in literature for selected ionic liquids. The temperature dependence on the viscosity of the cations and anions were fitted to a Vogel–Fulcher–Tamman behavior. Binary interaction parameters and VFT type fitting parameters were then used to determine the viscosity of pure and mixtures of ionic liquids with different combinations of cations and anions to ensure the validity of the prediction method. Consequently, the viscosities of binary ionic liquid mixtures were then calculated by using this prediction method. In this work, the viscosity data of pure ionic liquids and of binary mixtures of ionic liquids are successfully calculated from 293.15 K to 363.15 K at 0.1 MPa. All calculated viscosity data showed excellent agreement with experimental data with a relative absolute average deviation lower than 1.7%.

New Method Based on the UNIFAC–VISCO Model for the Estimation of Ionic Liquids Viscosity Using the Experimental Data Recommended by Mathematical Gnostics

The viscosity of ionic liquids (ILs) has been modeled as a function of temperature and at atmospheric pressure using a new method based on the UNIFAC–VISCO method. This model extends the calculations previously reported by our group (see Zhao et al. J. Chem. Eng. Data 2016, 61, 2160–2169) which used 154 experimental viscosity data points of 25 ionic liquids for regression of a set of binary interaction parameters and ion Vogel–Fulcher–Tammann (VFT) parameters. Discrepancies in the experimental data of the same IL affect the quality of the correlation and thus the development of the predictive method. In this work, mathematical gnostics was used to analyze the experimental data from different sources and recommend one set of reliable data for each IL. These recommended data (totally 819 data points) for 70 ILs were correlated using this model to obtain an extended set of binary interaction parameters and ion VFT parameters, with a regression accuracy of 1.4%. In addition, 966 experimental viscosity data points for 11 binary mixtures of ILs were collected from literature to establish this model. All the binary data consist of 128 training data points used for the optimization of binary interaction parameters and 838 test data points used for the comparison of the pure evaluated values. The relative average absolute deviation (RAAD) for training and test is 2.9% and 3.9%, respectively.

The VLT-FLAMES Tarantula Survey: XIX. B-type supergiants: Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

Context: Model atmosphere analyses have been previously undertaken for both Galactic and extragalactic B-type supergiants. By contrast, little attention has been given to a comparison of the properties of single supergiants and those that are members of multiple systems. 

Aims: Atmospheric parameters and nitrogen abundances have been estimated for all the B-type supergiants identified in the VLT-FLAMES Tarantula survey. These include both single targets and binary candidates. The results have been analysed to investigate the role of binarity in the evolutionary history of supergiants. 

Methods: tlusty non-local thermodynamic equilibrium (LTE) model atmosphere calculations have been used to determine atmospheric parameters and nitrogen abundances for 34 single and 18 binary supergiants. Effective temperatures were deduced using the silicon balance technique, complemented by the helium ionisation in the hotter spectra. Surface gravities were estimated using Balmer line profiles and microturbulent velocities deduced using the silicon spectrum. Nitrogen abundances or upper limits were estimated from the Nii spectrum. The effects of a flux contribution from an unseen secondary were considered for the binary sample. Results. We present the first systematic study of the incidence of binarity for a sample of B-type supergiants across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the B-type supergiants it may be necessary to extend the TAMS to lower temperatures. This is also consistent with the derived distribution of mass discrepancies, projected rotational velocities and nitrogen abundances, provided that stars cooler than this temperature are post-red supergiant objects. For all the supergiants in the Tarantula and in a previous FLAMES survey, the majority have small projected rotational velocities. The distribution peaks at about 50 km s^-1 with 65% in the range 30 km s^-1 ≤ ν_e sin i ≤ 60 km s^-1. About ten per cent have larger ve sin i (≥100 km s^-1), but surprisingly these show little or no nitrogen enhancement. All the cooler supergiants have low projected rotational velocities of ≤70 km s^-1 and high nitrogen abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. Additionally, there is a lack of cooler binaries, possibly reflecting the small sample sizes. Single-star evolutionary models, which include rotation, can account for all of the nitrogen enhancement in both the single and binary samples. The detailed distribution of nitrogen abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. 

Conclusions: The first comparative study of single and binary B-type supergiants has revealed that the main sequence may be significantly wider than previously assumed, extending to T_eff = 20 000 K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolutionary status of blue supergiants.

The VLT-FLAMES Tarantula Survey: XVIII. Classifications and radial velocities of the B-type stars

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 ± 12.2 km s^-1 from 273 presumed single stars. Employing a 3σ criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (≥345 km s^-1) or small (≤65 km s^-1) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 ± 16.2 km s^-1), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12 + log (N/H) ≳ 8.5. Combined with a large rotational velocity (υ_e sin i = 345 ± 22 km s^-1), this is suggestive of past binary interaction for this star.

The VLT-FLAMES Tarantula Survey XXI. Stellar spin rates of O-type spectroscopic binaries

Context: The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are, however, found in multiple systems. 

Aims: By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary subpopulations with one another and with that of presumed-single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the spin rates of massive stars. 

Methods: We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (ν_esini) for components of 114 spectroscopic binaries in 30 Doradus. The ν_esini values are derived from the full width at half maximum (FWHM) of a set of spectral lines, using a FWHM vs. ν_esini calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. 

Results: The overall ν_esini distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at ν_esini<200 kms^-1) and a shoulder at intermediate velocities (200 <ν_esini<300 kms^-1). The distributions of binaries and single stars, however, differ in two ways. First, the main peak at ν_esini ~ 100kms^-1 is broader and slightly shifted towards higher spin rates in the binary distribution than that of the presumed-single stars. This shift is mostly due to short-period binaries (P_orb~<10 d). Second, the ν_esini distribution of primaries lacks a significant population of stars spinning faster than 300 kms^-1, while such a population is clearly present in the single-star sample. The ν_esini distribution of binaries with amplitudes of radial velocity variation in the range of 20 to 200 kms^-1 (mostly binaries with P_orb ~ 10-1000 d and/or with q<0.5) is similar to that of single O stars below ν_esini_~<170kms^-1. 

Conclusions: Our results are compatible with the assumption that binary components formed with the same spin distribution as single stars, and that this distribution contains few or no fast-spinning stars. The higher average spin rate of stars in short-period binaries may either be explained by spin-up through tides in such tight binary systems, or by spin-down of a fraction of the presumed-single stars and long-period binaries through magnetic braking (or by a combination of both mechanisms). Most primaries and secondaries of SB2 systems with P_orb~<10 d appear to have similar rotational velocities. This is in agreement with tidal locking in close binaries where the components have similar radii. The lack of very rapidly spinning stars among binary systems supports the idea that most stars with ν_esini_~> 300kms^-1 in the single-star sample are actually spun-up post-binary interaction products. Finally, the overall similarities (low-velocity peak and intermediate-velocity shoulder) of the spin distribution of binary and single stars argue for a massive star formation process in which the initial spin is set independently of whether stars are formed as single stars or as components of a binary system.

The VLT-FLAMES Tarantula survey: XX. The nature of the X-ray bright emission-line star VFTS 399

Context: The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. 

Aims: One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity-in this study we aim to determine its physical nature and the cause of this behaviour. 

Methods: To accomplish this we performed a time-resolved analysis of optical, infrared and X-ray observations. 

Results: We found VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select He i lines-taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. 

Conclusions: VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor, sharing many observational characteristics with classical Be X-ray binaries. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass 25 M_⊙ for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating "single" O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.

The VLT-FLAMES Tarantula Survey: XXIII. Two massive double-lined binaries in 30 Doradus

Aims: We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems in the Large Magellanic Cloud, VFTS 450 (O9.7 II-Ib + O7::) and VFTS 652 (B1 Ib + O9: III:). 

Methods: We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the "primary" as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. 

Results: We estimate T_eff = 27 kK, log g = 2.9 (cgs) for the VFTS 450 primary spectrum (34 kK, 3.6: for the secondary spectrum); and T_eff = 22 kK, log g = 2.8 for the VFTS 652 primary spectrum (35 kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89 d and 8.59 d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range ∼20-50 M_⊙. 

Conclusions: The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks.

The VLT-FLAMES tarantula survey:XII. Rotational velocities of the single O-type stars

Context. The 30 Doradus (30 Dor) region of the Large Magellanic Cloud, also known as the Tarantula nebula, is the nearest starburst region. It contains the richest population of massive stars in the Local Group, and it is thus the best possible laboratory to investigate open questions on the formation and evolution of massive stars. Aims. Using ground-based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Dor. The sample is large enough to obtain statistically significant information and to search for variations among subpopulations - in terms of spectral type, luminosity class, and spatial location - in the field of view. Methods. We measured projected rotational velocities, 3e sin i, by means of a Fourier transform method and a profile fitting method applied to a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, P(3e), of the equatorial rotational velocity, 3e. Results. The distribution of 3e sin i shows a two-component structure: a peak around 80 km s1 and a high-velocity tail extending up to 600 km s^-1 This structure is also present in the inferred distribution P(3e) with around 80% of the sample having 0 <3e ≤ 300 km s^-1 and the other 20% distributed in the high-velocity region. The presence of the low-velocity peak is consistent with what has been found in other studies for late O- and early B-type stars. Conclusions. Most of the stars in our sample rotate with a rate less than 20% of their break-up velocity. For the bulk of the sample, mass loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. If massive-star formation results in stars rotating at birth with a large portion of their break-up velocities, an alternative braking mechanism, possibly magnetic fields, is thus required to explain the present-day rotational properties of the O-type stars in 30 Dor. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Even though we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. That the highvelocity tail may be populated primarily (and perhaps exclusively) by post-binary interaction products has important implications for the evolutionary origin of systems that produce gamma-ray bursts. © 2013 Author(s).

The crowded magnetosphere of the post-common-envelope binary QS Virginis

We present high-speed photometry and high-resolution spectroscopy of the eclipsing post-common-envelope binary QS Virginis (QS Vir). Our Ultraviolet and Visual Echelle Spectrograph (UVES) spectra span multiple orbits over more than a year and reveal the presence of several large prominences passing in front of both the M star and its white dwarf companion, allowing us to triangulate their positions. Despite showing small variations on a time-scale of days, they persist for more than a year and may last decades. One large prominence extends almost three stellar radii from the M star. Roche tomography reveals that the M star is heavily spotted and that these spots are long-lived and in relatively fixed locations, preferentially found on the hemisphere facing the white dwarf. We also determine precise binary and physical parameters for the system. We find that the 14 220 ± 350 K white dwarf is relatively massive, 0.782 ± 0.013 M⊙, and has a radius of 0.010 68 ± 0.000 07 R⊙, consistent with evolutionary models. The tidally distorted M star has a mass of 0.382 ± 0.006 M⊙ and a radius of 0.381 ± 0.003 R⊙, also consistent with evolutionary models. We find that the magnesium absorption line from the white dwarf is broader than expected. This could be due to rotation (implying a spin period of only ˜700 s), or due to a weak (˜100 kG) magnetic field, we favour the latter interpretation. Since the M star's radius is still within its Roche lobe and there is no evidence that it is overinflated, we conclude that QS Vir is most likely a pre-cataclysmic binary just about to become semidetached.

Rotational velocities of single and binary O-type stars in the Tarantula Nebula

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, υ e sin i, of ~330 O-type objects, i.e. ~210 spectroscopic single stars and ~110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30 Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the υ e sin i distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100 km s^-1. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300 km s^-1, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part - and could potentially be completely due - to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.