Radiative hydrodynamic modelling and observations of the X-class solar flare on 2011 March 9

Aims. We investigated the response of the solar atmosphere to non-thermal electron beam heating using the radiative transfer and hydrodynamics modelling code RADYN. The temporal evolution of the parameters that describe the non-thermal electron energy distribution were derived from hard X-ray observations of a particular flare, and we compared the modelled and observed parameters.

Methods. The evolution of the non-thermal electron beam parameters during the X1.5 solar flare on 2011 March 9 were obtained from analysis of RHESSI X-ray spectra. The RADYN flare model was allowed to evolve for 110 s, after which the electron beam heating was ended, and was then allowed to continue evolving for a further 300 s. The modelled flare parameters were compared to the observed parameters determined from extreme-ultraviolet spectroscopy.

Results. The model produced a hotter and denser flare loop than that observed and also cooled more rapidly, suggesting that additional energy input in the decay phase of the flare is required. In the explosive evaporation phase a region of high-density cool material propagated upward through the corona. This material underwent a rapid increase in temperature as it was unable to radiate away all of the energy deposited across it by the non-thermal electron beam and via thermal conduction. A narrow and high-density (ne ≤ 1015 cm-3) region at the base of the flare transition region was the source of optical line emission in the model atmosphere. The collision-stopping depth of electrons was calculated throughout the evolution of the flare, and it was found that the compression of the lower atmosphere may permit electrons to penetrate farther into a flaring atmosphere compared to a quiet Sun atmosphere.

Energy levels, radiative rates and electron impact excitation rates for transitions in AlX

Energy levels, radiative rates and lifetimes are calculated among the lowest 98 levels of the n ≤4 configurations of Be-like Al X. The GRASP (General-purpose Relativistic Atomic Structure Package) is adopted and data are provided for all E1, E2, M1 and M2 transitions. Similar data are also obtained with the FAC (Flexible Atomic Code) to assess the accuracy of the calculations. Based on comparisons between calculations with the two codes as well as with available measurements, our listed energy levels are assessed to be accurate to better than 0.3 per cent. However, the accuracy for radiative rates and lifetimes is estimated to be about 20 per cent. Collision strengths are also calculated for which the DARC (Dirac Atomic R-matrix Code) is used. A wide energy range (up to 380 Ryd) is considered and resonances resolved in a fine energy mesh in the thresholds region. The collision strengths are subsequently averaged over a Maxwellian velocity distribution to determine effective collision strengths up to a temperature of 1.6 × 10⁷ K. Our results are compared with the previous (limited) atomic data and significant differences (up to a factor of 4) are noted for several transitions, particularly those which are not allowed in jj coupling. 

Early-type stars observed in the ESO UVES Paranal observatory project - IV. studies of CN, CH+ and CH in the interstellar medium

High spectral resolution (~80 000) and signal-to-noise observations from the Ultraviolet and Visual Echelle Spectrograph Paranal Observatory Project (UVES-POP) are used to study the interstellarmolecular lines CN (3874 Å), CH⁺ (3957, 4232 Å) and CH (3886, 4300 Å) towards 74 O- and B-type stellar sightlines. Additionally, archive data are presented for 140 ELODIE early-type stellar sightlines at R = 42 000, plus 25 FEROS at R = 48 000 and 3 UVES at R > 50 000, mainly in the CH⁺ (4232 Å) and CH (3886, 4300 Å) transitions. Detection rates are ~45 per cent for CN and ~67 per cent for the other lines in the POP sample, and ~10-15 per cent for CH⁺ and CH lines in the additional sample. CH and CH⁺ are well correlated between log[N(CH) cm^-2]~12-14, implying that these clouds are CH⁺-like CH and not CN-like CH. CH is also very well correlated with Na I D in the range log[N(Na I cm^-2]) ~12.2-14.2. A few sightlines show tentative velocity shifts of ~2 km s^-1 between CH and CH⁺, which appear to be caused by differences in component strength in blends, and hence do not provide firm evidence for shocks. Finally, we describe a search for ¹³CH⁺ in a sightline towards HD 76341. No ¹³CH⁺ is detected, placing a limit on the ¹³CH⁺ to ¹²CH⁺ ratio of ~0.01. If a formal fit is attempted, the equivalent width ratio in the two isotopes is a factor ~90 but with large errors. 

Energy levels, radiative rates and electron impact excitation rates for transitions in Si II

Energies for the lowest 56 levels, belonging to the 3s² 3p, 3s 3p², 3p³, 3s² 3d, 3s 3p 3d, 3s² 4ℓ and 3s² 5ℓ configurations of Si II, are calculated using the General-purpose Relativistic Atomic Structure Package (GRASP) code. Analogous calculations have also been performed (for up to 175 levels) using the FlexibleAtomicCode (FAC). Furthermore, radiative rates are calculated for all E1, E2, M1 and M2 transitions. Extensive comparisons are made with available theoretical and experimental energy levels, and the accuracy of the present results is assessed to be better than 0.1Ryd. Similarly, the accuracy for radiative rates (and subsequently lifetimes) is estimated to be better than 20 per cent for most of the (strong) transitions. Electron impact excitation collision strengths are also calculated, with the Dirac Atomic R-matrix Code (DARC), over a wide energy range up to 13 Ryd. Finally, to determine effective collision strengths, resonances are resolved in a fine energy mesh in the thresholds region. These collision strengths are averaged over a Maxwellian velocity distribution and results listed over a wide range of temperatures, up to 10^5.5 K. Our data are compared with earlier R-matrix calculations and differences noted, up to a factor of 2, for several transitions. Although scope remains for improvement, the accuracy for our results of collision strengths and effective collision strengths is assessed to be about 20 per cent for a majority of transitions. 

Transit timing variations in WASP-10b induced by stellar activity

The hot-JupiterWASP-10bwas reported by Maciejewski et al. to showtransit timing variations (TTVs) with an amplitude of ~3.5 min. These authors proposed that the observed TTVs were caused by a 0.1M_Jup perturbing companion with an orbital period of ~5.23 d, and hence, close to the outer 5:3 mean-motion resonance with WASP-10b. To test this scenario, we present eight new transit light curves of WASP-10b obtained with the Faulkes Telescope North and the Liverpool Telescope. The new light curves, together with 22 previously published ones, were modelled with a Markov Chain Monte Carlo transit fitting code. Transit depth differences reported forWASP-10b are thought to be due to starspot-induced brightness modulation of the host star. Assuming the star is brighter at the activity minimum, we favour a small planetary radius. We find Rp = 1.039^+0.043 _-0.049R_Jup in agreement with Johnson et al. and Maciejewski et al. Recent studies find no evidence for a significant eccentricity in this system. We present consistent system parameters for a circular orbit and refine the orbital ephemeris ofWASP-10b. Our homogeneously derived transit times do not support the previous claimed TTV signal, which was strongly dependent on two previously published transits that have been incorrectly normalized. Nevertheless, a linear ephemeris is not a statistically good fit to the transit times of WASP-10b. We show that the observed transit time variations are due to spot occultation features or systematics. We discuss and exemplify the effects of occultation spot features in the measured transit times and show that despite spot occultation during egress and ingress being difficult to distinguish in the transit light curves, they have a significant effect in the measured transit times. We conclude that if we account for spot features, the transit times of WASP-10b are consistent with a linear ephemeris with the exception of one transit (epoch 143) which is a partial transit. Therefore, there is currently no evidence for the existence of a companion to WASP-10b. Our results support the lack of TTVs of hot-Jupiters reported for the Kepler sample. 

Expanded iron uta spectra - Probing the thermal stability limits in AGN clouds

The Fe unresolved transition arrays (UTAs) produce prominent features in the 15-17 Å wavelength range in the spectra of active galactic nuclei (AGNs). Here, we present new calculations of the energies and oscillator strengths of inner-shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in AGNs. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration. © 2013. The American Astronomical Society. All rights reserved.

Mechanisms for mhd poynting flux generation in simulations of solar photospheric magnetoconvection

We investigate the generation mechanisms of MHD Poynting flux in the magnetized solar photosphere. Using radiative MHD modeling of the solar photosphere with initial magnetic configurations that differ in their field strength and geometry, we show the presence of two different mechanisms for MHD Poynting flux generation in simulations of solar photospheric magnetoconvection. The weaker mechanism is connected to vertical transport of weak horizontal magnetic fields in the convectively stable layers of the upper photosphere, while the stronger is the production of Poynting flux in strongly magnetized intergranular lanes experiencing horizontal vortex motions. These mechanisms may be responsible for the energy transport from the solar convection zone to the higher layers of the solar atmosphere.

SN 2009ib: a Type II-P supernova with an unusually long plateau

We present optical and near-infrared photometry and spectroscopy of SN 2009ib, a Type II-P supernova in NGC 1559. This object has moderate brightness, similar to those of the intermediate-luminosity SNe 2008in and 2009N. Its plateau phase is unusually long, lasting for about 130 d after explosion. The spectra are similar to those of the subluminous SN 2002gd, with moderate expansion velocities. We estimate the Ni-56 mass produced as 0.046 +/- A 0.015 M-aS (TM). We determine the distance to SN 2009ib using both the expanding photosphere method (EPM) and the standard candle method. We also apply EPM to SN 1986L, a Type II-P SN that exploded in the same galaxy. Combining the results of different methods, we conclude the distance to NGC 1559 as D = 19.8 +/- A 3.0 Mpc. We examine archival, pre-explosion images of the field taken with the Hubble Space Telescope, and find a faint source at the position of the SN, which has a yellow colour [(V - I)(0) = 0.85 mag]. Assuming it is a single star, we estimate its initial mass as M-ZAMS = 20 M-aS (TM). We also examine the possibility, that instead of the yellow source the progenitor of SN 2009ib is a red supergiant star too faint to be detected. In this case, we estimate the upper limit for the initial zero-age main sequence (ZAMS) mass of the progenitor to be similar to 14-17 M-aS (TM). In addition, we infer the physical properties of the progenitor at the explosion via hydrodynamical modelling of the observables, and estimate the total energy as similar to 0.55 x 10(51) erg, the pre-explosion radius as similar to 400 R-aS (TM), and the ejected envelope mass as similar to 15 M-aS (TM), which implies that the mass of the progenitor before explosion was similar to 16.5-17 M-aS (TM).

Optical and near-infrared observations of SN 2011dh - The first 100 days

We present optical and near-infrared (NIR) photometry and spectroscopy of the Type IIb supernova (SN) 2011dh for the first 100 days. We complement our extensive dataset with Swift ultra-violet (UV) and Spitzer mid-infrared (MIR) data to build a UV to MIR bolometric lightcurve using both photometric and spectroscopic data. Hydrodynamical modelling of the SN based on this bolometric lightcurve have been presented in Bersten et al. (2012, ApJ, 757, 31). We find that the absorption minimum for the hydrogen lines is never seen below ~11 000 km s-1 but approaches this value as the lines get weaker. This suggests that the interface between the helium core and hydrogen rich envelope is located near this velocity in agreement with the Bersten et al. (2012) He4R270 ejecta model. Spectral modelling of the hydrogen lines using this ejecta model supports the conclusion and we find a hydrogen mass of 0.01-0.04 M⊙ to be consistent with the observed spectral evolution. We estimate that the photosphere reaches the helium core at 5-7 days whereas the helium lines appear between ~10 and ~15 days, close to the photosphere and then move outward in velocity until ~40 days. This suggests that increasing non-thermal excitation due to decreasing optical depth for the γ-rays is driving the early evolution of these lines. The Spitzer 4.5 μm band shows a significant flux excess, which we attribute to CO fundamental band emission or a thermal dust echo although further work using late time data is needed. Thedistance and in particular the extinction, where we use spectral modelling to put further constraints, is discussed in some detail as well as the sensitivity of the hydrodynamical modelling to errors in these quantities. We also provide and discuss pre- and post-explosion observations of the SN site which shows a reduction by ~75 percent in flux at the position of the yellow supergiant coincident with SN 2011dh. The B, V and r band decline rates of 0.0073, 0.0090 and 0.0053 mag day-1 respectively are consistent with the remaining flux being emitted by the SN. Hence we find that the star was indeed the progenitor of SN 2011dh as previously suggested by Maund et al. (2011, ApJ, 739, L37) and which is also consistent with the results from the hydrodynamical modelling. Figures 2, 3, Tables 3-10, and Appendices are available in electronic form at http://www.aanda.orgThe photometric tables are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/562/A17

Selecting superluminous supernovae in faint galaxies from the first year of the Pan-STARRS1 Medium Deep Survey

The Pan-STARRS1 (PS1) survey has obtained imaging in five bands (grizyP1) over 10 Medium Deep Survey (MDS) fields covering a total of 70 square degrees. This paper describes the search for apparently hostless supernovae (SNe) within the first year of PS1 MDS data with an aim of discovering superluminous supernovae (SLSNe). A total of 249 hostless transients were discovered down to a limiting magnitude of MAB ∼ 23.5, of which 76 were classified as Type Ia supernovae (SNe Ia). There were 57 SNe with complete light curves that are likely core-collapse SNe (CCSNe) or type Ic SLSNe and 12 of these have had spectra taken. Of these 12 hostless, non-Type Ia SNe, 7 were SLSNe of type Ic at redshifts between 0.5 and 1.4. This illustrates that the discovery rate of type Ic SLSNe can be maximized by concentrating on hostless transients and removing normal SNe Ia. We present data for two possible SLSNe; PS1-10pm (z = 1.206) and PS1-10ahf (z = 1.1), and estimate the rate of type Ic SLSNe to be between 3+3−2×10−53+3−2×10−53+3−2×10−5 and 8+2−1×10−58+2−1×10−58+2−1×10−5 that of the CCSN rate within 0.3 ≤ z ≤ 1.4 by applying a Monte Carlo technique. The rate of slowly evolving, type Ic SLSNe (such as SN2007bi) is estimated as a factor of 10 lower than this range.

Late-time spectral line formation in Type IIb supernovae, with application to SN 1993J, SN 2008ax, and SN 2011dh

We investigate line formation processes in Type IIb supernovae (SNe) from 100 to 500 days post-explosion using spectral synthesis calculations. The modelling identifies the nuclear burning layers and physical mechanisms that produce the major emission lines, and the diagnostic potential of these. We compare the model calculations with data on the three best observed Type IIb SNe to-date - SN 1993J, SN 2008ax, and SN 2011dh. Oxygen nucleosynthesis depends sensitively on the main-sequence mass of the star and modelling of the [O I] lambda lambda 6300, 6364 lines constrains the progenitors of these three SNe to the M-ZAMS = 12-16 M-circle dot range (ejected oxygen masses 0.3-0.9 M-circle dot), with SN 2011dh towards the lower end and SN 1993J towards the upper end of the range. The high ejecta masses from M-ZAMS greater than or similar to 17 M-circle dot progenitors give rise to brighter nebular phase emission lines than observed. Nucleosynthesis analysis thus supports a scenario of low-to-moderate mass progenitors for Type IIb SNe, and by implication an origin in binary systems. We demonstrate how oxygen and magnesium recombination lines may be combined to diagnose the magnesium mass in the SN ejecta. For SN 2011dh, a magnesium mass of 0.02-0.14 M-circle dot is derived, which gives a Mg/O production ratio consistent with the solar value. Nitrogen left in the He envelope from CNO burning gives strong [N II] lambda lambda 6548, 6583 emission lines that dominate over Ha emission in our models. The hydrogen envelopes of Type IIb SNe are too small and dilute to produce any noticeable H alpha emission or absorption after similar to 150 days, and nebular phase emission seen around 6550 angstrom is in many cases likely caused by [N II] lambda lambda 6548, 6583. Finally, the influence of radiative transport on the emergent line profiles is investigated. Significant line blocking in the metal core remains for several hundred days, which affects the emergent spectrum. These radiative transfer effects lead to early-time blueshifts of the emission line peaks, which gradually disappear as the optical depths decrease with time. The modelled evolution of this effect matches the observed evolution in SN 2011dh.

PESSTO: Survey description and products from the first data release by the Public ESO Spectroscopic Survey of Transient Objects

Context. The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. PESSTO classifies transients from publicly available sources and wide-field surveys, and selects science targets for detailed spectroscopic and photometric follow-up. PESSTO runs for nine months of the year, January - April and August - December inclusive, and typically has allocations of 10 nights per month. 

Aims. We describe the data reduction strategy and data products that are publicly available through the ESO archive as the Spectroscopic Survey data release 1 (SSDR1). 

Methods. PESSTO uses the New Technology Telescope with the instruments EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5^m for classification. Science targets are selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. We use standard EFOSC2 set-ups providing spectra with resolutions of 13-18 Å between 3345-9995 Å. A subset of the brighter science targets are selected for SOFI spectroscopy with the blue and red grisms (0.935-2.53 μm and resolutions 23-33 Å) and imaging with broadband JHK<inf>s</inf> filters. 

Results. This first data release (SSDR1) contains flux calibrated spectra from the first year (April 2012-2013). A total of 221 confirmed supernovae were classified, and we released calibrated optical spectra and classifications publicly within 24 h of the data being taken (via WISeREP). The data in SSDR1 replace those released spectra. They have more reliable and quantifiable flux calibrations, correction for telluric absorption, and are made available in standard ESO Phase 3 formats. We estimate the absolute accuracy of the flux calibrations for EFOSC2 across the whole survey in SSDR1 to be typically ∼15%, although a number of spectra will have less reliable absolute flux calibration because of weather and slit losses. Acquisition images for each spectrum are available which, in principle, can allow the user to refine the absolute flux calibration. The standard NIR reduction process does not produce high accuracy absolute spectrophotometry but synthetic photometry with accompanying JHK<inf>s</inf> imaging can improve this. Whenever possible, reduced SOFI images are provided to allow this. 

Conclusions. Future data releases will focus on improving the automated flux calibration of the data products. The rapid turnaround between discovery and classification and access to reliable pipeline processed data products has allowed early science papers in the first few months of the survey.

LSQ14bdq: a type Ic Super-luminous Supernova with a double-peaked light curve

We present data for LSQ14bdq, a hydrogen-poor super-luminous supernova (SLSN) discovered by the La Silla QUEST survey and classified by the Public ESO Spectroscopic Survey of Transient Objects. The spectrum and light curve are very similar to slow-declining SLSNe such as PTF12dam. However, detections within ∼1 day after explosion show a bright and relatively fast initial peak, lasting for ∼15 days, prior to the usual slow rise to maximum light. The broader, main peak can be fit with either central engine or circumstellar interaction models. We discuss the implications of the precursor peak in the context of these models. It is too bright and narrow to be explained as a normal ⁵⁶Ni-powered SN, and we suggest that interaction models may struggle to fit the two peaks simultaneously. We propose that the initial peak may arise from the post-shock cooling of extended stellar material, and reheating by a central engine drives the second peak. In this picture, we show that an explosion energy of ∼2 × 10⁵² erg and a progenitor radius of a few hundred solar radii would be required to power the early emission. The competing engine models involve rapidly spinning magnetars (neutron stars) or fallback onto a central black hole. The prompt energy required may favor the black hole scenario. The bright initial peak may be difficult to reconcile with a compact Wolf-Rayet star as a progenitor since the inferred energies and ejected masses become unphysical.

Constraints on explosive silicon burning in core-collapse supernovae from measured Ni/Fe ratios

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of Ni-58/(Fe-54 + Ni-56), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova (SN) that produced a Ni/Fe ratio of 3.4 +/- 1.2 times solar, we find that burning of a fuel with neutron excess eta approximate to 6 x 10(-3) is required. Unless the progenitor metallicity is over five times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. SNe producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of M-ZAMS less than or similar to 13 M-circle dot stars exploding with a delay time of less than one second (M-cut < 1.5 M-circle dot) are able to achieve such silicon-shell ejection. SNe that produce solar or subsolar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic entropy burning outside the iron core, and neutrino-neutronization obtained in electron capture models remains the only viable explanation.

Superluminous Supernovae as Standardizable Candles and High-Redshift Distance Probes

We investigate the use of type Ic superluminous supernovae (SLSN Ic) as standardizable candles and distance indicators. Their appeal as cosmological probes stems from their remarkable peak luminosities, hot blackbody temperatures, and bright rest-frame ultraviolet emission. We present a sample of 16 published SLSN, from redshifts 0.1 to 1.2, and calculate accurate K corrections to determine uniform magnitudes in 2 synthetic rest-frame filter bandpasses with central wavelengths at 400 nm and 520 nm. At 400 nm, we find an encouragingly low scatter in their uncorrected, raw mean magnitudes with M(400) = -21.86 ± 0.35 mag for the full sample of 16 objects. We investigate the correlation between their decline rates and peak magnitude and find that the brighter events appear to decline more slowly. In a manner similar to the Phillips relation for type Ia SNe (SNe Ia), we define a ΔM 20 decline relation. This correlates peak magnitude and decline over 20 days and can reduce the scatter in standardized peak magnitudes to ±0.22 mag. We further show that M(400) appears to have a strong color dependence. Redder objects are fainter and also become redder faster. Using this peak magnitudecolor evolution relation, a surprisingly low scatter of between ±0.08 mag and ±0.13 mag can be found in peak magnitudes, depending on sample selection. However, we caution that only 8 to 10 objects currently have enough data to test this peak magnitudecolor evolution relation. We conclude that SLSN Ic are promising distance indicators in the high-redshift universe in regimes beyond those possible with SNe Ia. Although the empirical relationships are encouraging, the unknown progenitor systems, how they may evolve with redshift, and the uncertain explosion physics are of some concern. The two major measurement uncertainties are the limited numbers of low-redshift, well-studied objects available to test these relationships and internal dust extinction in the host galaxies.