The Type IIb SN 2011dh: Two years of observations and modelling of the lightcurves

We present optical and near-infrared (NIR) photometry and spectroscopy as well as modelling of the lightcurves of the Type IIb supernova (SN) 2011dh. Our extensive dataset, for which we present the observations obtained after day 100, spans two years, and complemented with Spitzer mid-infrared (MIR) data, we use it to build an optical-to-MIR bolometric lightcurve between days 3 and 732. To model the bolometric lightcurve before day 400 we use a grid of hydrodynamical SN models, which allows us to determine the errors in the derived quantities, and a bolometric correction determined with steady-state non-local thermodynamic equilibrium (NLTE) modelling. Using this method we find a helium core mass of 3.1^+0.7<inf>-0.4</inf> M<inf>⊙</inf> for SN 2011dh, consistent within error bars with previous results obtained using the bolometric lightcurve before day 80. We compute bolometric and broad-band lightcurves between days 100 and 500 from spectral steady-state NLTE models, presented and discussed in a companion paper. The preferred 12 M<inf>⊙</inf> (initial mass) model, previously found to agree well with the observed spectra, shows a good overall agreement with the observed lightcurves, although some discrepancies exist. Time-dependent NLTE modelling shows that after day ∼600 a steady-state assumption is no longer valid. The radioactive energy deposition in this phase is likely dominated by the positrons emitted in the decay of ⁵⁶Co, but seems insufficient to reproduce the lightcurves, and what energy source is dominating the emitted flux is unclear. We find an excess in the K and the MIR bands developing between days 100 and 250, during which an increase in the optical decline rate is also observed. A local origin of the excess is suggested by the depth of the He I 20 581 Å absorption. Steady-state NLTE models with a modest dust opacity in the core (τ = 0.44), turned on during this period, reproduce the observed behaviour, but an additional excess in the Spitzer 4.5 μm band remains. Carbon-monoxide (CO) first-overtone band emission is detected at day 206, and possibly at day 89, and assuming the additional excess to bedominated by CO fundamental band emission, we find fundamental to first-overtone band ratios considerably higher than observed in SN 1987A. The profiles of the [O i] 6300 Å and Mg i] 4571 Å lines show a remarkable similarity, suggesting that these lines originate from a common nuclear burning zone (O/Ne/Mg), and using small scale fluctuations in the line profiles we estimate a filling factor of ≲ 0.07 for the emitting material. This paper concludes our extensive observational and modelling work on SN 2011dh. The results from hydrodynamical modelling, steady-state NLTE modelling, and stellar evolutionary progenitor analysis are all consistent, and suggest an initial mass of ∼12 M<inf>⊙</inf> for the progenitor.

A comparative study of Type II-P and II-L supernova rise times as exemplified by the case of LSQ13cuw

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P.

Mineralogy, chemistry and biological contingents of an early-middle Miocene Antarctic paleosol and its relevance as a Martian analogue

Fossil mesofauna and bacteria recovered from a paleosol in a moraine situated adjacent to the inland ice, Antarctica, and dating to the earliest glacial event in the Antarctic Dry Valleys opens several questions. The most important relates to understanding of the mineralogy and chemistry of the weathered substrate habitat in which Coleoptera apparently thrived at some point in the Early/Middle Miocene and perhaps earlier. Here, Coleoptera remains are only located in one of six horizons in a paleosol formed in moraine deposited during the alpine glacial event (> 15 Ma). A tendency for quartz to decrease upward in the section may be a detrital effect or a product of dissolution in the early stage of profile morphogenesis when climate was presumably milder and the depositing glacier of temperate type. Discontinuous distributions of smectite, laumontite, and hexahydrite may have provided nutrients and water to mesofauna and bacteria during the early stage of biotic colonization of the profile. Because the mesofauna were members of burrowing Coleoptera species, future work should assess the degree to which the organisms occupied other sites in the Dry Valleys in the past. Whereas there is no reasonable expectations of finding Coleoptera/insect remains on Mars, the chemistry and mineralogy of the paleosol is within a life expectancy window for the presence of microorganisms, principally bacteria and fungi. Thus, parameters discussed here within this Antarctic paleosol could provide an analogue to identifying similar fossil or life-bearing weathered regolith on Mars.

The Rossiter-McLaughlin effect reloaded: Probing the 3D spin-orbit geometry, differential stellar rotation, and the spatially-resolved stellar spectrum of star-planet systems

When a planet transits its host star, it blocks regions of the stellar surface from view; this causes a distortion of the spectral lines and a change in the line-of-sight (LOS) velocities, known as the Rossiter-McLaughlin (RM) effect. Since the LOS velocities depend, in part, on the stellar rotation, the RM waveform is sensitive to the star-planet alignment (which provides information on the system’s dynamical history). We present a new RM modelling technique that directly measures the spatially-resolved stellar spectrum behind the planet. This is done by scaling the continuum flux of the (HARPS) spectra by the transit light curve, and then subtracting the infrom the out-of-transit spectra to isolate the starlight behind the planet. This technique does not assume any shape for the intrinsic local profiles. In it, we also allow for differential stellar rotation and centre-to-limb variations in the convective blueshift. We apply this technique to HD 189733 and compare to 3D magnetohydrodynamic (MHD) simulations. We reject rigid body rotation with high confidence (>99% probability), which allows us to determine the occulted stellar latitudes and measure the stellar inclination. In turn, we determine both the sky-projected (λ ≈ −0.4 ± 0.2◦) and true 3D obliquity (ψ ≈ 7+12 −4 ◦ ). We also find good agreement with the MHD simulations, with no significant centre-to-limb variations detectable in the local profiles. Hence, this technique provides a new powerful tool that can probe stellar photospheres, differential rotation, determine 3D obliquities, and remove sky-projection biases in planet migration theories. This technique can be implemented with existing instrumentation, but will become even more powerful with the next generation of high-precision radial velocity spectrographs.

Rates of E1, E2, M1, and M2 transitions in Ni II

Aims. We present rates for all E1, E2, M1, and M2 transitions among the 295 fine-structure levels of the configurations 3d⁹, 3d⁸4s, 3d⁷4s², 3d⁸4p, and 3d⁷4s4p, determined through an extensive configuration interaction calculation. 

Methods. The CIV3 code developed by Hibbert and coworkers is used to determine for these levels configuration interaction wave functions with relativistic effects introduced through the Breit-Pauli approximation. 

Results. Two different sets of calculations have been undertaken with different 3d and 4d functions to ascertain the effect of such variation. The main body of the text includes a representative selection of data, chosen so that key points can be discussed. Some analysis to assess the accuracy of the present data has been undertaken, including comparison with earlier calculations and the more limited range of experimental determinations. The full set of transition data is given in the supplementary material as it is very extensive. 

Conclusions. We believe that the present transition data are the best currently available.

Valence and L-shell photoionization of Cl-like argon using R-matrix techniques

Photoionization cross-sections are obtained using the relativistic DiracAtomic R-matrix Codes (DARC) for all valence and L-shell energy ranges between 27 and 270 eV. A total of 557 levels arising from the dominant configurations 3s²3p⁴, 3s³p5, 3p⁶, 3s²3p³[3d, 4s, 4p], 3p⁵3d, 3s²3p²3d², 3s³p4³d, 3s³p3³d2 and 2s²2p⁵3s²3p⁵ have been included in the targetwavefunction representation of the Ar III ion, including up to 4p in the orbital basis. We also performed a smaller Breit-Pauli (BP) calculation containing the lowest 124 levels. Direct comparisons are made with previous theoretical and experimental work for both valence shell and L-shell photoionization. Excellent agreement was found for transitions involving the ²P^o initial state to all allowed final states for both calculations across a range of photon energies. A number of resonant states have been identified to help analyse and explain the nature of the spectra at photon energies between 250 and 270 eV.

Ultraviolet emission lines of Si II in cool star and solar spectra

Recent atomic physics calculations for Si II are employed within the CLOUDY modelling code to analyse Hubble Space Telescope (HST) STIS ultraviolet spectra of three cool stars, β Geminorum, α Centauri A and B, as well as previously published HST/GHRS observations of α Tau, plus solar quiet Sun data from the High Resolution Telescope and Spectrograph. Discrepancies found previously between theory and observation for line intensity ratios involving the 3s²3p ²PJ-3s3p^{2 4}P^J' intercombination multiplet of Si II at ~ 2335 Å are significantly reduced, as are those for ratios containing the 3s²3p ²PJ-3s3p^{2 2}DJ ~ transitions at ~1816 Å. This is primarily due to the effect of the new Si II transition probabilities. However, these atomic data are not only very different from previous calculations, but also show large disagreements with measurements, specifically those of Calamai et al. for the intercombination lines. New measurements of transition probabilities for Si II are hence urgently required to confirm (or otherwise) the accuracy of the recently calculated values. If the new calculations are confirmed, then a long-standing discrepancy between theory and observation will have finally been resolved. However, if the older measurements are found to be correct, then the agreement between theory and observation is simply a coincidence and the existing discrepancies remain.

ALMA-resolved salt emission traces the chemical footprint and inner wind morphology of VY Canis Majoris

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

An evaluation of the correlation between open solar flux and total solar irradiance

The correlation between the coronal source flux F_{S} and the total solar irradiance I_{TS} is re-evaluated in the light of an additional 5 years' data from the rising phase of solar cycle 23 and also by using cosmic ray fluxes detected at Earth. Tests on monthly averages show that the correlation with F_{S} deduced from the interplanetary magnetic field (correlation coefficient, r = 0.62) is highly significant (99.999%), but that there is insufficient data for the higher correlation with annual means (r = 0.80) to be considered significant. Anti-correlations between I_{TS} and cosmic ray fluxes are found in monthly data for all stations and geomagnetic rigidity cut-offs (r ranging from −0.63 to −0.74) and these have significance levels between 85% and 98%. In all cases, the t is poorest for the earliest data (i.e., prior to 1982). Excluding these data improves the anticorrelation with cosmic rays to r = −0:93 for one-year running means. Both the interplanetary magnetic field data and the cosmic ray fluxes indicate that the total solar irradiance lags behind the open solar flux with a delay that is estimated to have an optimum value of 2.8 months (and is within the uncertainty range 0.8-8.0 months at the 90% level).

The supernova remnant CTB 37B and its associated magnetar CXOU J171405.7-381031: evidence for a magnetar-driven remnant

We discuss the association between the candidate magnetar CXOU J171405.7-381031 and the supernova remnant CTB 37B. The recent detection of the period derivative of the object allowed an estimation of a young characteristic age of only similar to 1000 yr. This value is too small to be compatible even with the minimum radius of the remnant being >= 10 pc, the value corresponding to the lower limit of the estimated distance of 10.2 +/- 3.5 kpc, unless the true distance happens to be even smaller than the lower limit. We argue that a consistent scenario for the remnant`s origin, in which the latter is powered by the energy injected by a young magnetar, is indeed more accurate to explain the young age, and demonstrates its non-standard (i.e. magnetar-driven) nature.

Resonance and chaos: I. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted threebody problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincaré surfaces of section with a mass ratio of 10-3 (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. On theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

Resonance and chaos: II. Exterior resonances and asymmetric libration

The motion of a test particle in the vicinity of exterior resonances is examined in the context of the planar, circular, restricted three-body problem. The existence of asymmetric periodic orbits associated with the 1 : n resonances (where n = 2, 3, 4, 5) is confirmed; there is also evidence of asymmetric resonances associated with larger values of n. A detailed examination of the evolution of the family of orbits associated with the 1:2 resonance shows the sequence that leads to asymmetric libration. On the basis of numerical studies of the phase space it is concluded that the existence of asymmetric libration means that the region exterior to the perturbing mass is more chaotic than the interior region. The apparent absence of 'particles' in 1 : n resonances in the solar system may reflect this inherent bias.

Chaotic diffusion caused by close encounters with several massive asteroids: II. the regions of (10) Hygiea, (2) Pallas, and (31) Euphrosyne

Context. Close encounters with (1) Ceres and (4) Vesta, the two most massive bodies in the main belt, are known to be a mechanism of dynamical mobility able to significantly alter proper elements of minor bodies, and they are the main source of dynamical mobility for medium-sized and large asteroids (D > 20 km, approximately). Recently, it has been shown that drift rates caused by close encounters with massive asteroids may change significantly on timescales of 30 Myr when different models (i.e., different numbers of massive asteroids) are considered. Aims. So far, not much attention has been given to the case of diffusion caused by the other most massive bodies in the main belt: (2) Pallas, (10) Hygiea, and (31) Euphrosyne, the third, fourth, and one of the most massive highly inclined asteroids in the main belt, respectively. Since (2) Pallas is a highly inclined object, relative velocities at encounter with other asteroids tend to be high and changes in proper elements are therefore relatively small. It was thus believed that the scattering effect caused by highly inclined objects in general should be small. Can diffusion by close encounters with these asteroids be a significant mechanism of long-term dynamical mobility? Methods. By performing simulations with symplectic integrators, we studied the problem of scattering caused by close encounters with (2) Pallas, (10) Hygiea, and (31) Euphrosyne when only the massive asteroids (and the eight planets) are considered, and the other massive main belt asteroids and non-gravitational forces are also accounted for. Results. By finding relatively small values of drift rates for (2) Pallas, we confirm that orbital scattering by this highly inclined object is indeed a minor effect. Unexpectedly, however, we obtained values of drift rates for changes in proper semi-major axis a caused by (10) Hygiea and (31) Euphrosyne larger than what was previously found for scattering by (4) Vesta. These high rates may have repercussions on the orbital evolution and age estimate of their respective families. © 2013 ESO.