Self-consistent modelling of line-driven hot-star winds with Monte Carlo radiation hydrodynamics

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multiline effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach are demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally, we discuss some of the current limitations of our method, some possible extensions and potential future applications.

CSI 2264: Simultaneous Optical and Infrared Light Curves of Young Disk-bearing Stars in NGC 2264 with CoRoT and Spitzer—Evidence for Multiple Origins of Variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30 day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" with discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with Hα emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk. Based on data from the Spitzer and CoRoT missions. The CoRoT space mission was developed and is operated by the French space agency CNES, with participation of ESA's RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain.

Hα Line Profile Asymmetries and the Chromospheric Flare Velocity Field

The asymmetries observed in the line profiles of solar flares can provide important diagnostics of the properties and dynamics of the flaring atmosphere. In this paper the evolution of the Hα and Ca ii λ8542 lines are studied using high spatial, temporal, and spectral resolution ground-based observations of an M1.1 flare obtained with the Swedish 1 m Solar Telescope. The temporal evolution of the Hα line profiles from the flare kernel shows excess emission in the red wing (red asymmetry) before flare maximum and excess in the blue wing (blue asymmetry) after maximum. However, the Ca ii λ8542 line does not follow the same pattern, showing only a weak red asymmetry during the flare. RADYN simulations are used to synthesize spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. We show that the red asymmetry observed in Hα is not necessarily associated with plasma downflows, and the blue asymmetry may not be related to plasma upflows. Indeed, we conclude that the steep velocity gradients in the flaring chromosphere modify the wavelength of the central reversal in the Hα line profile. The shift in the wavelength of maximum opacity to shorter and longer wavelengths generates the red and blue asymmetries, respectively.

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

Searching for swept-up hydrogen and helium in the late-time spectra of 11 nearby Type Ia supernovae

The direct detection of a stellar system that explodes as a Type Ia supernova (SN Ia) has not yet been successful. Various indirect methods have been used to investigate SN Ia progenitor systems but none have produced conclusive results. A prediction of single-degenerate models is that H- (or He-) rich material from the envelope of the companion star should be swept up by the SN ejecta in the explosion. Seven SNe Ia have been analysed to date looking for signs of H-rich material in their late-time spectra and none were detected. We present results from new late-time spectra of 11 SNe Ia obtained at the Very Large Telescope using XShooter and FORS2. We present the tentative detection of Hα emission for SN 2013ct, corresponding to ∼0.007 M_⊙ of stripped/ablated companion star material (under the assumptions of the spectral modelling). This mass is significantly lower than expected for single-degenerate scenarios, suggesting that >0.1 M_⊙ of H-rich is present but not observed. We do not detect Hα emission in the other 10 SNe Ia. This brings the total sample of normal SNe Ia with non-detections (<0.001–0.058 M_⊙) of H-rich material to 17 events. The simplest explanation for these non-detections is that these objects did not result from the explosion of a CO white dwarf accreting matter from a H-rich companion star via Roche lobe overflow or symbiotic channels. However, further spectral modelling is needed to confirm this. We also find no evidence of He-emission features, but models with He-rich companion stars are not available to place mass limits.

ULTRAVIOLET EMISSION LINES of Si II in QUASARS - INVESTIGATING the Si II DISASTER

The observed line intensity ratios of the Si ii λ1263 and λ1307 multiplets to that of Si ii λ1814 in the broad-line region (BLR) of quasars are both an order of magnitude larger than the theoretical values. This was first pointed out by Baldwin et al., who termed it the "Si ii disaster," and it has remained unresolved. We investigate the problem in the light of newly published atomic data for Si ii. Specifically, we perform BLR calculations using several different atomic data sets within the CLOUDY modeling code under optically thick quasar cloud conditions. In addition, we test for selective pumping by the source photons or intrinsic galactic reddening as possible causes for the discrepancy, and we also consider blending with other species. However, we find that none of the options investigated resolve the Si ii disaster, with the potential exception of microturbulent velocity broadening and line blending. We find that a larger microturbulent velocity () may solve the Si ii disaster through continuum pumping and other effects. The CLOUDY models indicate strong blending of the Si ii λ1307 multiplet with emission lines of O i, although the predicted degree of blending is incompatible with the observed λ1263/λ1307 intensity ratios. Clearly, more work is required on the quasar modeling of not just the Si ii lines but also nearby transitions (in particular those of O i) to fully investigate whether blending may be responsible for the Si ii disaster.

Chemistry in AGB stars: successes and challenges

Emission and absorption line observations of molecules in late-type stars are a vital component in our understanding of stellar evolution, dust formation and mass loss in these objects. The molecular composition of the gas in the circumstellar envelopes of AGB stars reflects chemical processes in gas whose properties are strong functions of radius with density and temperature varying by more than ten and two orders of magnitude, respectively. In addition, the interstellar UV field plays a critical role in determining not only molecular abundances but also their radial distributions. In this article, I shall briefly review some recent successful approaches to describing chemistry in both the inner and outer envelopes and outline areas of challenge for the future.

Photoionization of Co+ and electron-impact excitation of Co2 + using the Dirac R-matrix method

Modelling of massive stars and supernovae (SNe) plays a crucial role in understanding galaxies. From this modelling we can derive fundamental constraints on stellar evolution, mass-loss processes, mixing, and the products of nucleosynthesis. Proper account must be taken of all important processes that populate and depopulate the levels (collisional excitation, de-excitation, ionization, recombination, photoionization, bound–bound processes). For the analysis of Type Ia SNe and core collapse SNe (Types Ib, Ic and II) Fe group elements are particularly important. Unfortunately little data is currently available and most noticeably absent are the photoionization cross-sections for the Fe-peaks which have high abundances in SNe. Important interactions for both photoionization and electron-impact excitation are calculated using the relativistic Dirac atomic R-matrix codes (DARC) for low-ionization stages of Cobalt. All results are calculated up to photon energies of 45 eV and electron energies up to 20 eV. The wavefunction representation of Co III has been generated using GRASP0 by including the dominant 3d⁷, 3d⁶[4s, 4p], 3p⁴3d⁹ and 3p⁶3d⁹ configurations, resulting in 292 fine structure levels. Electron-impact collision strengths and Maxwellian averaged effective collision strengths across a wide range of astrophysically relevant temperatures are computed for Co III. In addition, statistically weighted level-resolved ground and metastable photoionization cross-sections are presented for Co II and compared directly with existing work.

Molecular Hydrogen Emission from Protoplanetary Disks

We have modeled self-consistently the density and temperature profiles of gas and dust in protoplanetary disks, taking into account irradiation from a central star. Making use of this physical structure, we have calculated the level populations of molecular hydrogen and the line emission from the disks. As a result, we can reproduce the observed strong line spectra of molecular hydrogen from protoplanetary disks, both in the ultraviolet (UV) and the near-infrared, but only if the central star has a strong UV excess radiation.

A marine reservoir correction database and on-line interface

Calibration is essential for interpretation of radiocarbon dates, especially when the 14C dates are compared to historical or climatic records with a different chronological basis. 14C ages of samples from the marine environment, such as shells or fish bones, or samples with a marine component, such as human bone in coastal regions, require an additional consideration because of the reservoir age of the ocean. While the pre-industrial global mean reservoir correction, R(t), is about 400 years, local variations (?R) can be several hundred years or more. ?R compilations on a global scale have been undertaken previously (Stuiver et al. 1986; Stuiver and Braziunas 1993), but have not been updated recently. Here we describe an on-line reservoir correction database accessed via mapping software. Rather than publishing a static ?R compilation, new data will be incorporated when it becomes available. The on-line marine reservoir correction database can be accessed at the website http://www.calib.org/.

High-resolution Keck I spectroscopy of Galactic halo post-asymptotic giant branch stars

Absolute and differential abundance analyses have been performed from high-resolution, high signal-to-noise ratio optical (Keck I) spectra for three evolved Galactic halo stars, namely PG 1704 + 222, HD 341617 and LSIV -0401. Their derived atmospheric parameters indicate that all three objects are undergoing a post-asymptotic giant branch (post-AGB) phase of evolution. A differential abundance analysis reveals HD 341617 as having a mild carbon deficiency of 0.74 dex, possibly due to the star having evolved off the AGB before the onset of the third dredge-up. Although such carbon underabundances are typical of hot post-AGB objects, the same trend is not observed in PG 1704 + 222, where the carbon abundance is found to be consistent with those derived for nitrogen and oxygen. Hence, a dredge-up scenario need not be invoked to explain the chemical composition of PG 1704 + 222. For LSIV -0401 no iron deficiency is apparent relative to magnesium and silicon, and hence a gas- dust separation event in the AGB progenitor need not be invoked for this star.

CaII K interstellar observations towards early-type disc and halo stars

We present high-resolution (R = lambda/Deltalambda similar to 40 000) Ca II K interstellar observations (lambda(air) = 3933.66Angstrom) towards 88 mainly B-type stars, of which 74 are taken from the Edinburgh-Cape or Palomar-Green surveys, and 81 have > 25degrees. The majority of the data come from previously existing spectroscopy, although also included are 18 new observations of stars with echelle spectra taken with UVES on the Very Large Telescope UT2 (Kueyen). Some 49 of the sample stars have distance estimates above the Galactic plane (z) greater than or equal to 1 kpc, and are thus good probes of the halo interstellar medium. Of the 362 interstellar Ca K components that we detect, 75 (21 per cent) have absolute values of their LSR velocity values exceeding 40 km s(-1). In terms of the deviation velocity for the sightlines with distance estimates, 46/273 (17 per cent) of components have velocity values exceeding those predicted by standard Galactic rotation by more than 40 km s(-1). Combining this data set with previous observations, we find that the median value of the reduced equivalent width (REW) of stars with z greater than or equal to 1 kpc (EW x sin ) is similar to 115 mAngstrom (n = 80), similar to that observed in extragalactic sightlines by Bowen. Using data of all z distances, the REW at infinity is found to be similar to 130 mAngstrom, with the scaleheight (1) of the Ca II K column density distribution being;z 800 pc (n = 196) and reduced column density at infinity of log[N(Ca II K) cm(-2)] similar to 12.24. This implies that similar to30 per cent of Ca II K absorption occurs at distances exceeding similar to1 kpc. For nine sightlines, with distance exceeding 1 kpc and with a companion object within 5degrees, we find that all but two have values of Ca II reduced equivalent width the same to within similar to20 per cent, when the REW of the nearest object is extrapolated to the distance of the further of the pair, and assuming 1 = 800 pc. For 29 of our sightlines with z greater than or equal to 1 kpc and a H I detection from the Leiden-Dwingeloo survey (beamsize of 0.5degrees), we find log(N(Ca II K)IN(H I)) ranging from -7.4 to - 8.4. Values of the Ca II K abundance relative to neutral hydrogen (log[N(Ca II K) cm(-2)] - log[N(H I) cm(-2)]) are found to be more than similar to0.5 dex higher in stars with distances exceeding approximate to100 pc, when compared with the (log[N(Ca II K) cm(-2)] -log[N(H-tot) cm(-2)]) values found in nearby sightlines such as those in Wakker & Mathis (2000). Finally, stellar Ca II K equivalent widths of the sample are determined for 26 objects.

The VLT-FLAMES survey of massive stars: Observations in the Galactic clusters NGC 3293, NGC 4755 and NGC 6611

We introduce a new survey of massive stars in the Galaxy and the Magellanic Clouds using the Fibre Large Array Multi- Element Spectrograph ( FLAMES) instrument at the Very Large Telescope ( VLT). Here we present observations of 269 Galactic stars with the FLAMES- Giraffe Spectrograph ( R similar or equal to 25 000), in fields centered on the open clusters NGC3293, NGC4755 and NGC6611. These data are supplemented by a further 50 targets observed with the Fibre- Fed Extended Range Optical Spectrograph ( FEROS, R = 48 000). Following a description of our scientific motivations and target selection criteria, the data reduction methods are described; of critical importance the FLAMES reduction pipeline is found to yield spectra that are in excellent agreement with less automated methods. Spectral classifications and radial velocity measurements are presented for each star, with particular attention paid to morphological peculiarities and evidence of binarity. These observations represent a significant increase in the known spectral content of NGC3293 and NGC4755, and will serve as standards against which our subsequent FLAMES observations in the Magellanic Clouds will be compared.

The periodic variations of a white-light flare observed with ULTRACAM

High time resolution observations of a white-light flare on the active star EQ PegB show evidence of intensity variations with a period of ≈10 s. The period drifts to longer values during the decay phase of the flare. If the oscillation is interpreted as an impulsively-excited, standing-acoustic wave in a flare loop, the period implies a loop length of ≈3.4 Mm and ≈6.8 Mm for the case of the fundamental mode and the second harmonic, respectively. However, the small loop lengths imply a very high modulation depth making the acoustic interpretation unlikely. A more realistic interpretation may be that of a fast-MHD wave, with the modulation of the emission being due to the magnetic field. Alternatively, the variations could be due to a series of reconnection events. The periodic signature may then arise as a result of the lateral separation of individual flare loops or current sheets with oscillatory dynamics (i.e., periodic reconnection).

Proof of polar ejection from the close-binary core of the planetary nebula Abell 63

We present the first detailed kinematical analysis of the planetary nebula Abell 63, which is known to contain the eclipsing close-binary nucleus UU Sge. Abell 63 provides an important test case in investigating the role of close-binary central stars on the evolution of planetary nebulae. Longslit observations were obtained using the Manchester echelle spectrometer combined with the 2.1-m San Pedro Martir Telescope. The spectra reveal that the central bright rim of Abell 63 has a tube-like structure. A deep image shows collimated lobes extending from the nebula, which are shown to be high-velocity outflows. The kinematic ages of the nebular rim and the extended lobes are calculated to be 8400 +/- 500 and 12900 +/- 2800 yr, respectively, which suggests that the lobes were formed at an earlier stage than the nebular rim. This is consistent with expectations that disc-generated jets form immediately after the common envelope phase. A morphological-kinematical model of the central nebula is presented and the best-fitting model is found to have the same inclination as the orbital plane of the central binary system; this is the first proof that a close-binary system directly affects the shaping of its nebula. A Hubble-type flow is well-established in the morphological-kinematical modelling of the observed line profiles and imagery. Two possible formation models for the elongated lobes of Abell 63 are considered, (i) a low-density, pressure-driven jet excavates a cavity in the remnant asymptotic giant branch (AGB) envelope; (ii) high-density bullets form the lobes in a single ballistic ejection event.