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.

SN2014J gamma rays from the 56Ni decay chain

Context. The detection and measurement of gamma-ray lines from the decaychain of 56Ni provides unique information about the explosionin supernovae. SN2014J at 3.3 Mpc is a sufficiently-nearby supernova oftype Ia so that such measurements have been feasible with the gamma-rayspectrometer SPI on ESA's INTEGRAL gamma-ray observatory.

Aims:The 56Ni freshly produced in the supernova is understood topower the optical light curve, because it emits gamma rays upon itsradioactive decay first to 56Co and then to 56Fe.Gamma-ray lines from 56Co decay are expected to becomedirectly visible through the white dwarf material several weeks afterthe explosion, as they progressively penetrate the overlying material ofthe supernova envelope, which is diluted as it expands. The lines areexpected to be Doppler-shifted or broadened from the kinematics of the56Ni ejecta. We aim to exploit high-resolution gamma-rayspectroscopy with the SPI spectrometer on INTEGRAL toward constrainingthe 56Ni distribution and kinematics in this supernova.

Methods: We use the observations with the SPI spectrometer onINTEGRAL, together with an improved instrumental background method.

Results: We detect the two main lines from 56Co decay at847 and 1238 keV, which are significantly Doppler-broadened, and atintensities (3.65 ± 1.21) × 10-4 and (2.27± 0.69) × 10-4 ph cm-2s-1, respectively, at their brightness maximum. We measuretheir rise toward a maximum after about 60-100 days and a declinethereafter. The intensity ratio of the two lines is found to beconsistent with expectations from 56Co decay (0.62 ±0.28 at brightness maximum, the expected ratio is 0.68). We find thatthe broad lines seen in the late, gamma-ray transparent phase are notrepresentative of the early gamma-ray emission, and notice instead thatthe emission spectrum is complex and irregular until the supernova isfully transparent to gamma rays, with progressive uncovering of the bulkof 56Ni. We infer that the explosion morphology is notspherically symmetric, both in the distribution of 56Ni andin the unburnt material which occults the 56Co emission.After we compare light curves from different plausible models, theresulting 56Ni mass is determined to be 0.49 ± 0.09M⊙.

Locations of peculiar supernovae as a diagnostic of their origins

We put constraints on the properties of the progenitors of peculiar calcium-rich transients using the distribution of locations within their host galaxies. We confirm that this class of transients do not follow the galaxy stellar mass profile and are more likely to be found in remote locations of their apparent hosts. We test the hypothesis that these transients are from low-metallicity progenitors by comparing their spatial distributions with the predictions of self-consistent cosmological simulations that include star formation and chemical enrichment. We find that while metal-poor stars and our transient sample show a consistent preference for large offsets, metallicity alone cannot explain the extreme cases. Invoking a lower age limit on the progenitor helps to improve the match, indicating these events may result from a very old metal-poor population. We also investigate the radial distribution of globular cluster systems, and show that they too are consistent with the class of calcium-rich transients. Because photometric upper limits exist for globular clusters for some members of the class, a production mechanism related to the dense environment of globular clusters is not favoured for the calcium-rich events. However, the methods developed in this paper may be used in the future to constrain the effects of low metallicity on radially distant core-collapse events or help establish a correlation with globular clusters for other classes of peculiar explosions.

Self-regulated Shocks in Massive Star Binary Systems

In an early-type, massive star binary system, X-ray bright shocks result from the powerful collision of stellar winds driven by radiation pressure on spectral line transitions. We examine the influence of the X-rays from the wind-wind collision shocks on the radiative driving of the stellar winds using steady-state models that include a parameterized line force with X-ray ionization dependence. Our primary result is that X-ray radiation from the shocks inhibits wind acceleration and can lead to a lower pre-shock velocity, and a correspondingly lower shocked plasma temperature, yet the intrinsic X-ray luminosity of the shocks, L X, remains largely unaltered, with the exception of a modest increase at small binary separations. Due to the feedback loop between the ionizing X-rays from the shocks and the wind driving, we term this scenario as self-regulated shocks. This effect is found to greatly increase the range of binary separations at which a wind-photosphere collision is likely to occur in systems where the momenta of the two winds are significantly different. Furthermore, the excessive levels of X-ray ionization close to the shocks completely suppress the line force, and we suggest that this may render radiative braking less effective. Comparisons of model results against observations reveal reasonable agreement in terms of log (L X/L bol). The inclusion of self-regulated shocks improves the match for kT values in roughly equal wind momenta systems, but there is a systematic offset for systems with unequal wind momenta (if considered to be a wind-photosphere collision).

EARLY RADIO AND X-RAY OBSERVATIONS OF THE YOUNGEST NEARBY TYPE Ia SUPERNOVA PTF 11kly (SN 2011fe)

On 2011 August 24 (UT) the Palomar Transient Factory (PTF) discovered PTF11kly (SN 2011fe), the youngest and most nearby Type Ia supernova (SN Ia) in decades. We followed this event up in the radio (centimeter and millimeter bands) and X-ray bands, starting about a day after the estimated explosion time. We present our analysis of the radio and X-ray observations, yielding the tightest constraints yet placed on the pre-explosion mass-loss rate from the progenitor system of this supernova. We find a robust limit of from sensitive X-ray non-detections, as well as a similar limit from radio data, which depends, however, on assumptions about microphysical parameters. We discuss our results in the context of single-degenerate models for SNe Ia and find that our observations modestly disfavor symbiotic progenitor models involving a red giant donor, but cannot constrain systems accreting from main-sequence or sub-giant stars, including the popular supersoft channel. In view of the proximity of PTF11kly and the sensitivity of our prompt observations, we would have to wait for a long time (a decade or longer) in order to more meaningfully probe the circumstellar matter of SNe Ia.

SN 2011dh: DISCOVERY OF A TYPE IIb SUPERNOVA FROM A COMPACT PROGENITOR IN THE NEARBY GALAXY M51

On 2011 May 31 UT a supernova (SN) exploded in the nearby galaxy M51 (the Whirlpool Galaxy). We discovered this event using small telescopes equipped with CCD cameras and also detected it with the Palomar Transient Factory survey, rapidly confirming it to be a Type II SN. Here, we present multi-color ultraviolet through infrared photometry which is used to calculate the bolometric luminosity and a series of spectra. Our early-time observations indicate that SN 2011dh resulted from the explosion of a relatively compact progenitor star. Rapid shock-breakout cooling leads to relatively low temperatures in early-time spectra, compared to explosions of red supergiant stars, as well as a rapid early light curve decline. Optical spectra of SN 2011dh are dominated by H lines out to day 10 after explosion, after which He I lines develop. This SN is likely a member of the cIIb (compact IIb) class, with progenitor radius larger than that of SN 2008ax and smaller than the eIIb (extended IIb) SN 1993J progenitor. Our data imply that the object identified in pre-explosion Hubble Space Telescope images at the SN location is possibly a companion to the progenitor or a blended source, and not the progenitor star itself, as its radius (~1013 cm) would be highly inconsistent with constraints from our post-explosion spectra.

Young Stellar Object VARiability (YSOVAR): Long Timescale Variations in the Mid-infrared

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 and 4.5 μm) time series photometry of the Orion Nebula Cluster plus smaller footprints in 11 other star-forming cores (AFGL 490, NGC 1333, Mon R2, GGD 12-15, NGC 2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC 1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs roughly twice per day for ~40 days. We also define a "standard sample of members" consisting of all the IR-selected members and X-ray-selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data—the Stetson index, a χ2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of six to seven years by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data set; out of members and field objects combined, at most 0.02% may have transient IR excesses.

YSOVAR: Mid-infrared Variability in the Star-forming Region Lynds 1688

The emission from young stellar objects (YSOs) in the mid-infrared (mid-IR) is dominated by the inner rim of their circumstellar disks. We present IR data from the Young Stellar Object VARiability (YSOVAR) survey of ~800 objects in the direction of the Lynds 1688 (L1688) star-forming region over four visibility windows spanning 1.6 yr using the Spitzer Space Telescope in its warm mission phase. Among all light curves, 57 sources are cluster members identified based on their spectral energy distribution and X-ray emission. Almost all cluster members show significant variability. The amplitude of the variability is larger in more embedded YSOs. Ten out of 57 cluster members have periodic variations in the light curves with periods typically between three and seven days, but even for those sources, significant variability in addition to the periodic signal can be seen. No period is stable over 1.6 yr. Nonperiodic light curves often still show a preferred timescale of variability that is longer for more embedded sources. About half of all sources exhibit redder colors in a fainter state. This is compatible with time-variable absorption toward the YSO. The other half becomes bluer when fainter. These colors can only be explained with significant changes in the structure of the inner disk. No relation between mid-IR variability and stellar effective temperature or X-ray spectrum is found.

Star-planet interactions and selection effects from planet detection methods

Planets may have effects on their host stars by tidal or magnetic interaction. Such star-planet interactions are thought to enhance the activity level of the host star. However, stellar activity also affects the sensitivity of planet detection methods. Samples of planet-hosting stars which are investigated for such star-planet interactions are therefore subject to strong selection effects which need to be taken into account.

YSOVAR: Mid-infrared Variability of Young Stellar Objects and Their Disks in the Cluster IRAS 20050+2720

We present a time-variability study of young stellar objects (YSOs) in the cluster IRAS 20050+2720, performed at 3.6 and 4.5 μm with the Spitzer Space Telescope; this study is part of the Young Stellar Object VARiability (YSOVAR) project. We have collected light curves for 181 cluster members over 60 days. We find a high variability fraction among embedded cluster members of ca. 70%, whereas young stars without a detectable disk display variability less often (in ca. 50% of the cases) and with lower amplitudes. We detect periodic variability for 33 sources with periods primarily in the range of 2–6 days. Practically all embedded periodic sources display additional variability on top of their periodicity. Furthermore, we analyze the slopes of the tracks that our sources span in the color–magnitude diagram (CMD). We find that sources with long variability time scales tend to display CMD slopes that are at least partially influenced by accretion processes, while sources with short variability timescales tend to display extinction-dominated slopes. We find a tentative trend of X-ray detected cluster members to vary on longer timescales than the X-ray undetected members.

YSOVAR: Mid-infrared Variability Among YSOs in the Star Formation Region GGD12-15

We present an IR-monitoring survey with the Spitzer Space Telescope of the star-forming region GGD 12-15. More than 1000 objects were monitored, including about 350 objects within the central 5′, which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability project. The region was also the subject of a contemporaneous 67 ks Chandra observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5′ are young stellar objects (YSOs) based on a combination of their spectral energy distribution, IR variability, and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-infrared (mid-IR) variability. There are 39 periodic sources, and all but one is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded YSOs. Most of the Class I/ II objects exhibit redder colors in a fainter state, which is compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR light curves. However, we find that X-ray detected Class II sources have longer timescales for change in the MIR than a similar set of non-X-ray detected Class IIs.

YSOVAR: Mid-infrared Variability in NGC 1333

As part of the Young Stellar Object VARiability (YSOVAR) program, wemonitored NGC 1333 for ∼35 days at 3.6 and 4.5 μm using theSpitzer Space Telescope. We report here on the mid-infrared variabilityof the point sources in the ∼10‧ × ∼20‧ areacentered on 03:29:06, +31:19:30 (J2000). Out of 701 light curves ineither channel, we find 78 variables over the YSOVAR campaign. Abouthalf of the members are variable. The variable fraction for the mostembedded spectral energy distributions (SEDs) (Class I, flat) is higherthan that for less embedded SEDs (Class II), which is in turn higherthan the star-like SEDs (Class III). A few objects have amplitudes(10–90th percentile brightness) in [3.6] or [4.5] > 0.2 mag; amore typical amplitude is 0.1–0.15 mag. The largest color changeis >0.2 mag. There are 24 periodic objects, with 40% of them beingflat SED class. This may mean that the periodic signal is primarily fromthe disk, not the photosphere, in those cases. We find 9 variableslikely to be “dippers,” where texture in the disk occultsthe central star, and 11 likely to be “bursters,” whereaccretion instabilities create brightness bursts. There are 39 objectsthat have significant trends in [3.6]–[4.5] color over thecampaign, about evenly divided between redder-when-fainter (consistentwith extinction variations) and bluer-when-fainter. About a third of the17 Class 0 and/or jet-driving sources from the literature are variableover the YSOVAR campaign, and a larger fraction (∼half) are variablebetween the YSOVAR campaign and the cryogenic-era Spitzer observations(6–7 years), perhaps because it takes time for the envelope torespond to changes in the central source. The NGC 1333 brown dwarfs donot stand out from the stellar light curves in any way except there is amuch larger fraction of periodic objects (∼60% of variable browndwarfs are periodic, compared to ∼30% of the variables overall).

Where are the mid-sized flares of ultracool M dwarfs?

We propose to observe the M8.5 dwarf SCR J1845-6357 with XMM-Newton EPIC for 60 ks. Very low-mass M dwarfs show a distinct drop in X-ray luminosity compared to slightly more massive M dwarfs. Surprisingly, this does not happen at the mass threshold where M dwarfs become fully convective (M4), but at significantly lower masses (M8). These very low mass stars seem to have a flaring behaviour different from earlier type stars: they display either occasional large flares or a very low-level "flickering" in their X-ray light curves, but not the canonical power-law flare-energy distribution observed for the Sun and other cool stars. Our aim is to collect a long-duration light curve for one of the most nearby ultracool dwarfs to quantify how its flare-energy distribution differs from earlier type stars.

X-ray driven evaporation of exoplanet atmospheres: discovery of a uniquely suited test system

The evaporation of exoplanetary atmospheres is thought to be driven by high-energy irradiation. However, the actual mass loss rates are not well constrained. Co-I Kipping has recently discovered that the star KOI-314, an M1V dwarf at 65 pc distance, is orbited by two earth-sized planets, the inner one of them rocky and the outer one gaseous (P_orb = 14d and 23d). Other recent works have shown an abundance of small rocky planets in very close orbits around their host stars, suggesting that the stellar high-energy irradiation evaporates away gaseous envelopes. KOI-314 is the first nearby system in which earth-sized planets of both types are detected, allowing us to constrain the efficiency of planetary evaporation if the stellar X-ray irradiation is measured. We therefore propose a 10 ks Chandra ACIS-S pointing to determine the stellar X-ray luminosity and hardness ratio. The accuracy of the orbital solution decreases quickly due to Transit-Timing Variations, which is why we ask for DDT.