A hydrodynamical model of the circumstellar bubble created by two massive stars

Numerical models of the wind-blown bubble of massive stars usually only account for the wind of a single star. However, since massive stars are usually formed in clusters, it would be more realistic to follow the evolution of a bubble created by several stars. We develope a two-dimensional (2D) model of the circumstellar bubble created by two massive stars, a 40 solar mass star and a 25 solar mass star, and follow its evolution. The stars are separated by approximately 16 pc and surrounded by a cold medium with a density of 20 particles per cubic cm. We use the MPI-AMRVAC hydrodynamics code to solve the conservation equations of hydrodynamics on a 2D cylindrical grid using time-dependent models for the wind parameters of the two stars. At the end of the stellar evolution (4.5 and 7.0 million years for the 40 and 25 solar mass stars, respectively), we simulate the supernova explosion of each star. Each star initially creates its own bubble. However, as the bubbles expand they merge, creating a combined, aspherical bubble. The combined bubble evolves over time, influenced by the stellar winds and supernova explosions. The evolution of a wind-blown bubble created by two stars deviates from that of the bubbles around single stars. In particular, once one of the stars has exploded, the bubble is too large for the wind of the remaining star to maintain and the outer shell starts to disintegrate. The lack of thermal pressure inside the bubble also changes the behavior of circumstellar features close to the remaining star. The supernovae are contained inside the bubble, which reflects part of the energy back into the circumstellar medium.Comment: Accepted for publication in A&A. Six .avi files to be published online (uploaded to ArXiv DC and available as ancillary files) (updated after language corrections

Continuum driven winds from super-Eddington stars. A tale of two limits

Continuum driving is an effective method to drive a strong stellar wind. It is governed by two limits: the Eddington limit and the photon-tiring limit. A star must exceed the effective Eddington limit for continuum driving to overcome the stellar gravity. The photon-tiring limit places an upper limit on the mass loss rate that can be driven to infinity, given the energy available in the radiation field of the star. Since continuum driving does not require the presence of metals in the stellar atmosphere it is particularly suited to removing mass from low- and zero-metallicity stars and can play a crucial part in their evolution. Using a porosity length formalism we compute numerical simulations of super-Eddington, continuum driven winds to explore their behaviour for stars both below and above the photon-tiring limit. We find that below the photon tiring limit, continuum driving can produce a large, steady mass loss rate at velocities on the order of the escape velocity. If the star exceeds the photon-tiring limit, a steady solution is no longer possible. While the effective mass loss rate is still very large, the wind velocity is much smallerComment: to be published in the conference proceedings of: First Stars III, Santa Fe, 200

Eyes in the sky: Interactions between AGB winds and the interstellar magnetic field

We aim to examine the role of the interstellar magnetic field in shaping the extended morphologies of slow dusty winds of Asymptotic Giant-branch (AGB) stars in an effort to pin-point the origin of so-called eye shaped CSE of three carbon-rich AGB stars. In addition, we seek to understand if this pre-planetary nebula (PN) shaping can be responsible for asymmetries observed in PNe. Hydrodynamical simulations are used to study the effect of typical interstellar magnetic fields on the free-expanding spherical stellar winds as they sweep up the local interstellar medium (ISM). The simulations show that typical Galactic interstellar magnetic fields of 5 to 10 muG, are sufficient to alter the spherical expanding shells of AGB stars to appear as the characteristic eye shape revealed by far-infrared observations. The typical sizes of the simulated eyes are in accordance with the observed physical sizes. However, the eye shapes are of transient nature. Depending on the stellar and interstellar conditions they develop after 20,000 to 200,000yrs and last for about 50,000 to 500,000 yrs, assuming that the star is at rest relative to the local interstellar medium. Once formed the eye shape will develop lateral outflows parallel to the magnetic field. The "explosion" of a PN in the center of the eye-shaped dust shell gives rise to an asymmetrical nebula with prominent inward pointing Rayleigh-Taylor instabilities. Interstellar magnetic fields can clearly affect the shaping of wind-ISM interaction shells. The occurrence of the eyes is most strongly influenced by stellar space motion and ISM density. Observability of this transient phase is favoured for lines-of-sight perpendicular to the interstellar magnetic field direction. The simulations indicate that shaping of the pre-PN envelope can strongly affect the shape and size of PNe.Comment: Accepted for publication in A&A. Final version will contain animated result

Constraints on gamma-ray burst and supernova progenitors through circumstellar absorption lines. (II): Post-LBV Wolf-Rayet stars

Van Marle et al. (2005) showed that circumstellar absorption lines in early Type Ib/c supernova and gamma-ray burst afterglow spectra may reveal the progenitor evolution of the exploding Wolf-Rayet star. While the quoted paper deals with Wolf-Rayet stars which evolved through a red supergiant stage, we investigate here the initially more massive Wolf-Rayet stars which are thought to evolve through a Luminous Blue Variable (LBV) stage. We perform hydrodynamic simulations of the evolution of the circumstellar medium around a 60 Msol star, from the main sequence through the LBV and Wolf-Rayet stages, up to core collapse. We then compute the column density of the circumstellar matter as a function of radial velocity, time and angle. This allows a comparison with the number and blue-shifts, of absorption components in the spectra of LBVs, Wolf-Rayet stars, Type Ib/c supernovae and gamma-ray burst afterglows. Our simulation for the post-LBV stage shows the formation of various absorption components, which are, however, rather short lived; they dissipate on time scales shorter than 50,000yr. As the LBV stage is thought to occur at the beginning of core helium burning, the remaining Wolf-Rayet life time is expected to be one order of magnitude larger. When interpreting the absorption components in the afterglow spectrum of GRB-021004 as circumstellar, it can be concluded that the progenitor of this source did most likely not evolve through an LBV stage. However, a close binary with late common-envelope phase (Case C) may produce a circumstellar medium that closely resembles the LBV to Wolf-Rayet evolution, but with a much shorter Wolf-Rayet period.Comment: accepted for publication by A&

On the observability of bow shocks of Galactic runaway OB stars

Massive stars that have been ejected from their parent cluster and supersonically sailing away through the interstellar medium (ISM) are classified as exiled. They generate circumstellar bow shock nebulae that can be observed. We present two-dimensional, axisymmetric hydrodynamical simulations of a representative sample of stellar wind bow shocks from Galactic OB stars in an ambient medium of densities ranging from n_ISM=0.01 up to 10.0/cm3. Independently of their location in the Galaxy, we confirm that the infrared is the most appropriated waveband to search for bow shocks from massive stars. Their spectral energy distribution is the convenient tool to analyze them since their emission does not depend on the temporary effects which could affect unstable, thin-shelled bow shocks. Our numerical models of Galactic bow shocks generated by high-mass (~40 Mo) runaway stars yield H$\alpha$ fluxes which could be observed by facilities such as the SuperCOSMOS H-Alpha Survey. The brightest bow shock nebulae are produced in the denser regions of the ISM. We predict that bow shocks in the field observed at Ha by means of Rayleigh-sensitive facilities are formed around stars of initial mass larger than about 20 Mo. Our models of bow shocks from OB stars have the emission maximum in the wavelength range 3 <= lambda <= 50 micrometer which can be up to several orders of magnitude brighter than the runaway stars themselves, particularly for stars of initial mass larger than 20 Mo.Comment: 13 pages, 12 figures. Accepted to MNRAS (2016

Thin shell morphology in the circumstellar medium of massive binaries

We investigate the morphology of the collision front between the stellar winds of binary components in two long-period binary systems, one consisting of a hydrogen rich Wolf-Rayet star (WNL) and an O-star and the other of a Luminous Blue Variable (LBV) and an O-star. Specifically, we follow the development and evolution of instabilities that form in such a shell, if it is sufficiently compressed, due to both the wind interaction and the orbital motion. We use MPI-AMRVAC to time-integrate the equations of hydrodynamics, combined with optically thin radiative cooling, on an adaptive mesh 3D grid. Using parameters for generic binary systems, we simulate the interaction between the winds of the two stars. The WNL+O star binary shows a typical example of an adiabatic wind collision. The resulting shell is thick and smooth, showing no instabilities. On the other hand, the shell created by the collision of the O star wind with the LBV wind, combined with the orbital motion of the binary components, is susceptible to thin shell instabilities, which create a highly structured morphology. We identify the nature of the instabilities as both linear and non-linear thin-shell instabilities, with distinct differences between the leading and the trailing parts of the collision front. We also find that for binaries containing a star with a (relatively) slow wind, the global shape of the shell is determined more by the slow wind velocity and the orbital motion of the binary, than the ram pressure balance between the two winds. The interaction between massive binary winds needs further parametric exploration, to identify the role and dynamical importance of multiple instabilities at the collision front, as shown here for an LBV+O star system.Comment: 10 pages, 13 figures. Accepted for publication in A&

Cryo-Scanning Electron Microscopy Investigation of the Texture of Cooked Potatoes

The texture of steam-cooked potatoes was investigated by examining the fracture planes of four different cultivars, using cryo-scanning electron microscopy (cryo-SEM), which yielded a good preservation of the hydrated structures in potato tissue . For all cultivars , fracturing after steam -cooking took place between cells preferentially alongside the cell walls. However, textural difference appeared from the degree of intercellular contact, the cell shape and the appearance of cell surfaces. Cells in the fracture planes of firm potatoes had large intercellular contacts. In this case, most of the cells were flat and cell surfaces showed folds and cracks . For mealy potatoes , it appeared that cells in the fracture planes had little intercellular contacts . The cells were round and turgid and had smooth surfaces . In conclusion, the structure of cell wall and middle lamella , and the starch content , appear to be important parameters to distingui sh firm and mealy cultivars

The origin of blue-shifted absorption lines in a gamma-ray burst afterglow

The afterglow spectrum of GRB 021004 shows a system of blueshifted absorption lines, indicating the presence of matter moving towards us at discrete velocities in the range from 0 to more than 2000 km/s. We propose that these lines are the result of absorption by circumstellar matter, which was ejected by gamma-ray burst progenitor (a massive star) during its evolution. We have simulated the evolution of the circumstellar medium around such a star and find that the evolutionary sequence: main-sequence, Red Supergiant, Wolf-Rayet star can qualitatively reproduce the various absorption lines systems

Forming a constant density medium close to long gamma-ray bursts

The progenitor stars of long Gamma-Ray Bursts (GRBs) are thought to be Wolf-Rayet stars, which generate a massive and energetic wind. Nevertheless, about 25 percent of all GRB afterglows light curves indicate a constant density medium close to the exploding star. We explore various ways to produce this, by creating situations where the wind termination shock arrives very close to the star, as the shocked wind material has a nearly constant density. Typically, the distance between a Wolf-Rayet star and the wind termination shock is too large to allow afterglow formation in the shocked wind material. Here, we investigate possible causes allowing for a smaller distance: A high density or a high pressure in the surrounding interstellar medium (ISM), a weak Wolf-Rayet star wind, the presence of a binary companion, and fast motion of the Wolf-Rayet star relative to the ISM. We find that all four scenarios are possible in a limited parameter space, but that none of them is by itself likely to explain the large fraction of constant density afterglows. A low GRB progenitor metallicity, and a high GRB energy make the occurrence of a GRB afterglow in a constant density medium more likely. This may be consistent with constant densities beingpreferentially found for energetic, high redshift GRBs.Comment: 13 pages, 13 figures, new version: as accepted by Astronomy & Astrophysic

The Environments around Long-Duration Gamma-Ray Burst Progenitors

Gamma-ray burst (GRB) afterglow observations have allowed us to significantly constrain the engines producing these energetic explosions. Te redshift and position information provided by these afterglows have already allowed us to limit the progenitors of GRBs to only a few models. The afterglows may also provide another observation that can place further constraints on the GRB progenitor: measurements telling us about the environments surrounding GRBs. Current analyses of GRB afterglows suggest that roughly half of long-duration gamma-ray bursts occur in surroundings with density profiles that are uniform. We study the constraints placed by this observation on both the classic ``collapsar'' massive star progenitor and its relative, the ``helium-merger'' progenitor. We study several aspects of wind mass-loss and find that our modifications to the standard Wolf-Rayet mass-loss paradigm are not sufficient to produce constant density profiles. Although this does not rule out the standard ``collapsar'' progenitor, it does suggest a deficiency with this model. We then focus on the He-merger models and find that such progenitors can fit this particular constraint well. We show how detailed observations can not only determine the correct progenitor for GRBs, but also allow us to study binary evolution physics.Comment: 44 pages including 11 figure