Spectroscopy of SN 1987A at 0.9–2.4 μm: days 1348–3158

We present near-infrared spectroscopic observations of SN 1987A covering the period 1358 to 3158 d post explosion. This is the first time that IR spectra of a supernova have been obtained to such late epochs. The spectra comprise emission from both the ejecta and the bright, ring-shaped circumstellar medium (CSM). The most prominent CSM emission lines are recombination lines of H  i and He  i , and forbidden lines of [S  iii ] and [Fe  ii ]. The ejecta spectra include allowed lines of H  i , He  i and Na  i and forbidden lines of [Si  i ], [Fe  i ], [Fe  ii ] and possibly [S  i ]. The intensity ratios and widths of the H  i ejecta lines are consistent with a low-temperature Case B recombination spectrum arising from non-thermal ionization/excitation in an extended, adiabatically-cooled H envelope, as predicted by several authors. The slow decline of the ejecta forbidden lines, especially those of [Si  i ], indicates that pure non-thermal excitation was taking place, driven increasingly by the decay of ⁴⁴Ti. The ejecta iron exhibits particularly high velocities  (4000–4500 km s^-1)  , supporting scenarios where fast radioactive nickel is created and ejected just after the core bounce. In addition, the ejecta lines continue to exhibit blueshifts with values ∼−200 to −800 km s⁻¹ to at least day 2000. These blueshifts, which first appeared around day 600, probably indicate that very dense concentrations of dust persist in the ejecta, although an alternative explanation of asymmetry in the excitation conditions is not ruled out.     

Solar coronal heating by magnetosonic waves

Solar coronal heating by magnetohydrodynamic (MHD) waves is investigated. ultraviolet (UV) and X-ray emission lines of the corona show non-thermal broadenings. The wave rms velocities inferred from these observations are of the order of 25–60 km s⁻¹ . Assuming that these values are not negligible, we solved MHD equations in a quasi-linear approximation, by retaining the lowest order non-linear term in rms velocity. Plasma density distribution in the solar corona is assumed to be inhomogeneous. This plasma is also assumed to be permeated by dipole-like magnetic loops. Wave propagation is considered along the magnetic field lines. As dissipative processes, only the viscosity and parallel (to the local magnetic field lines) heat conduction are assumed to be important. Two wave modes emerged from the solution of the dispersion relation. The fast mode magneto-acoustic wave, if originated from the coronal base can propagate upwards into the corona and dissipate its mechanical energy as heat. The damping length-scale of the fast mode is of the order of 500 km. The wave energy flux associated with these waves turned out to be of the order of 2.5×10⁵ ergs cm⁻² s⁻¹ which is high enough to replace the energy lost by thermal conduction to the transition region and by optically thin coronal emission. The fast magneto-acoustic waves prove to be a likely candidate to heat the solar corona. The slow mode is absent, in other words cannot propagate in the solar corona.     

Hyperfine splitting of [Al vi] 3.66 μm and the Al isotopic ratio in NGC 6302

The core of planetary nebula NGC 6302 is filled with high-excitation photoionized gas at low expansion velocities. It represents a unique astrophysical situation in which to search for hyperfine structure (HFS) in coronal emission lines from highly ionized species. HFS is otherwise blended by thermal or velocity broadening. Spectra containing  [Al  vi ] 3.66 μm ³P₂←³P₁  , obtained with Phoenix on Gemini South at resolving powers of up to 75 000, resolve the line into five hyperfine components separated by 20–60 km s⁻¹ as a result of the coupling of the   I = 5/2  nuclear spin of ²⁷Al with the total electronic angular momentum J . The isotope ²⁶Al has a different nuclear spin of   I = 5  , and a different HFS, which allows us to place a 3σ upper limit on the ²⁶Al/²⁷Al abundance ratio of 1/33. We measure the HFS magnetic dipole coupling constants for [Al  vi ], and provide the first estimates of the electric quadrupole HFS coupling constants obtained through astronomical observations of an atomic transition.     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

The type IIn supernova 1995G: interaction with the circumstellar medium

We present the photometric and spectroscopic evolution of the type IIn SN 1995G in NGC 1643, on the basis of 4 years of optical and infrared observations. This supernova shows very flat optical light curves similar to SN 1988Z, with a slow decline rate at all times. The spectra are characterized by strong Balmer lines with multiple components in emission and with a P Cygni absorption component blueshifted by only 700 km s⁻¹. This feature indicates the presence of a slowly expanding shell above the SN ejecta as in the case of SNe 1994aj and 1996L. As in other SNe IIn, the slow luminosity decline cannot be explained only with a radioactive energy input, and an additional source of energy is required, most likely that produced by the interaction between supernova ejecta and a pre-existent circumstellar medium (CSM). It was estimated that the shell material has a density   n _H≫10⁸ cm^-3  , consistent with the absence of forbidden lines in the spectra. About 2 years after the burst the low-velocity shell is largely overtaken by the SN ejecta and the luminosity drops at a faster rate.     

Emission lines of Fe xv in spectra obtained with the Solar Extreme-Ultraviolet Research Telescope and Spectrograph

Recent R-matrix calculations of electron impact excitation rates in Mg-like Fe  xv are used to derive theoretical emission-line ratios involving transitions in the 243–418 Å  wavelength range. A comparison of these with a data set of solar active region, subflare and off-limb spectra, obtained during rocket flights by the Solar Extreme-Ultraviolet Research Telescope and Spectrograph (SERTS), reveals generally very good agreement between theory and observation, indicating that most of the Fe  xv emission lines may be employed with confidence as electron density diagnostics. In particular, the 312.55-Å  line of Fe  xv is not significantly blended with a Co  xvii transition in active region spectra, as suggested previously, although the latter does make a major contribution in the subflare observations. Most of the Fe  xv transitions which are blended have had the species responsible clearly identified, although there remain a few instances where this has not been possible. We briefly address the long-standing discrepancy between theory and experiment for the intensity ratio of the  3s^2 1S–3s3p ³P₁  intercombination line at 417.25 Å  to the  3s^2 1S–3s3p ¹P  resonance transition at 284.16 Å.     

Nucleosynthesis of 56Ni in wind-driven supernova explosions and constraints on the central engine of gamma-ray bursts

Theoretically expected natures of a supernova (SN) driven by a wind/jet are discussed. Approximate analytical formulations are derived to clarify basic physical processes involved in the wind/jet-driven explosions, and it is shown that the explosion properties are characterized by the energy injection rate     and the mass injection rate     . To explain observations of Supernova 1998bw associated with gamma-ray burst (GRB) 980425, the following conditions are required:     and     (if the wind Lorentz factor  Γ_w∼ 1  ) or     (if  Γ_w≫ 1  ). In Supernova 1998bw, ⁵⁶Ni  (∼0.4 M_⊙)  is probably produced in the shocked stellar mantle, not in the wind. The expected natures of SNe, e.g. ejected ⁵⁶Ni mass and ejecta mass, vary depending on     and     . The sequence of the SN properties from high     and     to low     and     is as follows: Supernova 1998bw-like – intermediate case – low mass ejecta  (≲1 M_⊙)  where ⁵⁶Ni is from the wind – whole collapse. This diversity may explain the diversity of SNe associated with GRBs. Our result can be used to constrain natures of the wind/jet, which are linked to the central engine of GRBs, by studying properties of the associated SNe.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.     

The circumstellar material around SN IIn 1997eg: another detection of a very narrow P Cygni profile

We report the detection of a very narrow P Cygni profile on top of the broad emission H α and H β lines of the Type IIn Supernova 1997eg. A similar feature has been detected in SN 1997ab, SN 1998S and SN 1995G . The detection of the narrow P Cygni profile indicates the existence of a dense circumstellar material (CSM), into which the ejecta of the supernova is expanding. From the analysis of the spectra of SN 1997eg we deduce (i) that such CSM is very dense  ( n ≳5×10⁷ cm^-3)  , (ii) that it has a low expanding velocity of about 160 km s⁻¹. The origin of such dense CSM can be either a very dense progenitor wind  ( M˙ ∼10^-2 M_⊙ yr^-1)  or a circumstellar shell product of the progenitor wind expanding into a high-pressure environment.     

RR Pic (1925): a Chandra X-ray view

We present the Chandra ACIS-S3 data of the old classical nova RR Pic (1925). The source has a count rate of 0.067 ± 0.002 count s⁻¹ in the 0.3–5.0 keV energy range. We detect the orbital period of the underlying binary system in the X-ray wavelengths. We also find that the neutral hydrogen column density differs for orbital minimum and orbital maximum spectra with values  0.25^+0.23_−0.18× 10²²  and  0.64^+0.13_−0.14× 10²² cm⁻²  at 3σ confidence level. The X-ray spectrum of RR Pic can be represented by a composite model of bremsstrahlung with a photoelectric absorption, two absorption lines centered around 1.1–1.4 keV and five Gaussian lines centered at emission lines around 0.3–1.1 keV corresponding to various transitions of S, N, O, C, Ne and Fe. The bremsstrahlung temperature derived from the fits ranges from 0.99 to 1.60 keV and the unabsorbed X-ray flux is found to be  2.5^+0.4_−1.2× 10⁻¹³ erg  cm⁻² s⁻¹  in the 0.3–5.0 keV range with a luminosity of 1.1 ± 0.2  10³¹ erg s⁻¹  at 600 pc. We also detect excess emission in the spectrum possibly originating from the reverse shock in the ejecta. A fit with a cooling flow plasma emission model shows enhanced abundances of He, C, N, O and Ne in the X-ray emitting region indicating existence of diffusive mixing.     

Fe xiii emission lines in active region spectra obtained with the Solar Extreme-Ultraviolet Research Telescope and Spectrograph

Recent fully relativistic calculations of radiative rates and electron impact excitation cross-sections for Fe  xiii are used to generate emission-line ratios involving 3s²3p²–3s3p³ and 3s²3p²–3s²3p3d transitions in the 170–225 and 235–450 Å wavelength ranges covered by the Solar Extreme-Ultraviolet Research Telescope and Spectrograph (SERTS). A comparison of these line ratios with SERTS active region observations from rocket flights in 1989 and 1995 reveals generally very good agreement between theory and experiment. Several new Fe  xiii emission features are identified, at wavelengths of 203.79, 259.94, 288.56 and 290.81 Å. However, major discrepancies between theory and observation remain for several Fe  xiii transitions, as previously found by Landi and others, which cannot be explained by blending. Errors in the adopted atomic data appear to be the most likely explanation, in particular for transitions which have 3s²3p3d ¹D₂ as their upper level. The most useful Fe  xiii electron-density diagnostics in the SERTS spectral regions are assessed, in terms of the line pairs involved being (i) apparently free of atomic physics problems and blends, (ii) close in wavelength to reduce the effects of possible errors in the instrumental intensity calibration, and (iii) very sensitive to changes in N _e over the range  10⁸–10¹¹ cm⁻³  . It is concluded that the ratios which best satisfy these conditions are 200.03/202.04 and 203.17/202.04 for the 170–225 Å wavelength region, and 348.18/320.80, 348.18/368.16, 359.64/348.18 and 359.83/368.16 for 235–450 Å.     

Neutral hydrogen absorption observations of the central region of NGC 3628

Using MERLIN with 0.2-arcsec resolution we have observed neutral hydrogen absorption against the central region of the starburst galaxy NGC 3628. The central region resolves into ∼16 continuum components at 1.4 GHz. From comparison with published 15-GHz data, we infer that these components are supernova remnants, although three components may be consistent with a weak active galactic nucleus. Neutral hydrogen absorption is seen against the continuum emission with column densities ∼10²² cm⁻². The absorption appears to be from two separate absorbing structures. Assuming a simple morphology, the main velocity structure can be attributed to a ring of neutral gas with a radius 130 pc rotating around a central starburst with a velocity gradient of 1270 km s⁻¹ kpc⁻¹. From simple assumptions, the mass interior to this ring is 0.9 × 10⁹ M_⊙. The second absorption structure may represent outflow from the starburst region or a large-scale galactic structure. Alternatively the absorption structure may be non-axisymmetric, such as a bar.     

The high-velocity molecular stream in W51B: a ring structure with compact H ii regions

We present ¹³ CO J  = 1 − 0 line observations of the H  ii region complex W51B located in the high-velocity (HV) stream. These observations reveal a filamentary and clumpy structure in the molecular gas. The mean local standard of rest (LSR) velocity ∼ + 65 km s⁻¹ of the molecular gas in this region is greater than the maximum velocities allowed by kinematic Galactic rotation curves. The size and mass of the molecular clouds are ∼ 48 × 17 pc² and ∼ 2.4 × 10⁵ M⊙ respectively. In a position–velocity diagram, molecular gas in the southern part comprises a redshifted ring structure with v _LSR=+ 60 to +73 km s⁻¹. The velocity gradient of this ring is ∼ 0.5 km s⁻¹ pc⁻¹, and the mass is ∼ 6.2 × 10⁴ M⊙. If we assume that the ring is expanding with a uniform velocity, the expansion velocity, radius and kinetic energy are ∼ 7 km s⁻¹, ∼ 13 pc and ∼ 3.0 × 10 ⁴⁹ erg respectively. The kinetic energy and mass spectrum of the ring could be explained by an expanding cylindrical cloud with a centrally condensed mass distribution. The locations of two compact H  ii regions, G49.0−0.3 and G48.9−0.3, coincide with the two molecular clumps in this ring. We discuss star formation, and the mechanism that produced the ring structure.     

On invisible plasma content in radio-loud AGNs: the case of TeV blazar Markarian 421

Invisible plasma content in blazar jets such as protons and/or thermal electron–positron ( e ^±) pairs is explored through combined arguments of dynamical and radiative processes. By comparing physical quantities required by the internal shock model with those obtained through the observed broad-band spectra for Mrk 421, we obtain that the ratio of the Lorentz factors of a pair of cold shells resides in about 2 ∼ 20, which implies that the shocks are at most mildly relativistic. Using the obtained Lorentz factors, the total mass density ρ in the shocked shells is investigated. The upper limit of ρ is obtained from the condition that thermal bremsstrahlung emission should not exceed the observed γ-ray luminosity, whilst the lower limit is constrained from the condition that the energy density of non-thermal electrons is smaller than that of the total plasma. Then, we find ρ is 10²–10³ times heavier than that of non-thermal electrons for pure   e ^±  pairs, while 10²–10⁶ times heavier for pure electron–proton ( e / p ) content, implying the existence of a large amount of invisible plasma. The origin of the continuous blazar sequence is briefly discussed and we speculate that the total mass density and/or the blending ratio of e ^± pairs and e / p plasma could be new key quantities for the origin of the sequence.     

Relativistic ionized accretion disc models of MCG–6-30-15

We present BeppoSAX observations of Nova Velorum 1999 (V382 Vel), carried out in a broad X-ray band covering 0.1–300 keV only 15 d after the discovery and again after 6 months. The nova was detected at day 15 with the BeppoSAX instruments which measured a flux F _x≃1.8×10⁻¹¹ erg cm⁻² s⁻¹ in the 0.1–10 keV range and a 2 σ upper limit F _x<6.7×10⁻¹² erg cm⁻² s⁻¹ in the 15–60 keV range. We attribute the emission to shocked nebular ejecta at a plasma temperature kT ≃6 keV . At six months no bright component emerged in the 15–60 keV range, but a bright central supersoft X-ray source appeared. The hot nebular component previously detected had cooled to a plasma temperature kT <1 keV . There was strong intrinsic absorption of the ejecta in the first observation and not in the second, because the column density of neutral hydrogen decreased from N (H)≃1.7×10²³ to N (H)≃10²¹ cm⁻² (close to the interstellar value). The unabsorbed X-ray flux also decreased from F _x=4.3×10⁻¹¹ to F _x≃10⁻¹² erg cm⁻² s⁻¹ .     

A Bayesian periodogram finds evidence for three planets in HD 11964

A Bayesian multiplanet Kepler periodogram has been developed for the analysis of precision radial velocity data. The periodogram employs a parallel tempering Markov chain Monte Carlo algorithm. The HD 11964 data have been re-analysed using 1, 2, 3 and 4 planet models. Assuming that all the models are equally probable a priori, the three planet model is found to be ≥600 times more probable than the next most probable model which is a two planet model. The most probable model exhibits three periods of  38.02^+0.06_−0.05, 360⁺⁴₋₄ and 1924⁺⁴⁴₋₄₃ d  , and eccentricities of  0.22^+0.11_−0.22, 0.63^+0.34_−0.17 and 0.05^+0.03_−0.05  , respectively. Assuming the three signals (each one consistent with a Keplerian orbit) are caused by planets, the corresponding limits on planetary mass ( M sin  i ) and semimajor axis are     respectively. The small difference (1.3σ) between the 360-d period and one year suggests that it might be worth investigating the barycentric correction for the HD 11964 data.     

Predictions of the extent of self-enrichment in oxygen of giant metal-poor H ii regions

In general, H  ii regions do not show clear signs of self-enrichment in products from massive stars  ( M ≥ 8 M_⊙)  . In order to explore why I modelled the contamination with Wolf–Rayet star ejecta of metal-poor  ( Z = 0.001)  H  ii regions, ionized either by a  10⁶ M_⊙  cluster of coeval stars (cluster 1) or by a cluster resulting from continuous star formation at a rate of  1 M_⊙ yr⁻¹  (cluster 2). The clusters have   Z = 0.001  and a Salpeter initial mass function from 0.1 to  120 M_⊙  . Independent one-dimensional constant density simulations of the emission-line spectra of unenriched H  ii regions were computed at the discrete ages 1, 2, 3, 4 and 5 Myr, with the photoionization code cloudy , using as input, radiative and mechanical stellar feedbacks predicted by the evolutionary synthesis code starburst99 . Each H  ii region was placed at the outer radius of the adiabatically expanding superbubble of Mac Low & McCray. For models with thermal and ionization balance time-scales of less than 1 Myr, and with oxygen emission-line ratios in agreement with observations, the volume of the superbubble and the H  ii region was uniformly and instantaneously polluted with stellar ejecta predicted by starburst99 . I obtained a maximum oxygen abundance enhancement of 0.025 dex, with cluster 1, at 4 Myr. It would be unobservable.     

The fading of supernova 1997D

We present a new set of spectroscopic and photometric data extending the observations of SN 1997D to over 400 d after the explosion. These observations confirm the peculiar properties of SN 1997D, such as the very low abundance of ⁵⁶Co (0.002 M_⊙) and the low expansion velocity of the ejecta (∼1000 km s⁻¹). We discuss the implications of these observations for the character of the progenitor and the nature of the remnant, showing that a Crab-like pulsar or an accreting neutron star formed in the explosion of a low-mass progenitor should already have produced a detectable luminosity at this epoch, in contrast with photometric data. On the other hand, the explosion of a high-mass progenitor with the formation of a black hole is consistent with the available observations. The consequences of this conclusion regarding the nature of the explosion and the prospects of directly identifying the black hole are also addressed.     

Calculated spectra for HeH+ and its effect on the opacity of cool metal-poor stars

The wavelength and Einstein A coefficient are calculated for all rotation–vibration transitions of  ⁴He¹H⁺, ³He¹H⁺, ⁴He²H⁺  and  ³ He²H⁺  , giving a complete line list and the partition function for  ⁴HeH⁺  and its isotopologues. This opacity is included in the calculation of the total opacity of low-metallicity stars and its effect is analysed for different conditions of temperature, density and hydrogen number fraction. For a low helium number fraction (as in the Sun), it is found that HeH⁺ has a visible but small effect for very low densities  (ρ≤ 10⁻¹⁰ g cm⁻³)  , at temperatures around 3500 K. However, for high helium number fraction, the effect of HeH⁺ becomes important for higher densities  (ρ≤ 10⁻⁶ g cm⁻³)  , its effect being most important for a temperature around 3500 K. Synthetic spectra for a variety of different conditions are presented.